简介本资源为河南省焦作市30米分辨率数字高程模型DEM地理信息数据包面向GIS初学者、城乡规划从业者、地理科研人员及遥感分析学习者支撑地形分析、坡度坡向计算、流域划分、灾害风险评估等实际应用。压缩包共12个文件包含核心高程数据文件焦作市dem.tif、配套市级行政边界Shapefile含.shp/.shx/.dbf/.prj等标准组件及辅助文件.ovr缩略图、.tfw地理配准、.xml元数据等完整覆盖GIS空间分析所需的数据结构与坐标系统信息总大小11.79MB。已有647人学习下载资源结构规范、开箱即用无需额外处理即可在ArcGIS、QGIS等平台直接加载并开展三维可视化、地形因子提取与空间叠加分析是理解区域地貌特征与开展基础地理建模的实用型入门级实测数据集。 直接说结论这份“河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip”属于典型的地方级开源地理数据包里面装的东西就两样——焦作市30米分辨率的DEM栅格数据外加一个市级行政区划边界shp文件。对做地形分析、规划选址、淹没模拟、视线分析、土地利用研究的人来说这个包相当实用。尤其是焦作这种北依太行、南临黄河、地势落差大的城市30米DEM能做的事情太多了。但问题在于很多人拿到zip包之后要么不知道怎么校验数据完整性要么不清楚DEM坐标参考和shp边界不匹配怎么统一甚至有人连“zip解压失败”“shp打开是空图层”这类基础坑都踩过。网上一搜全是碎片化问答缺一篇能把数据格式、工具操作、实战流程串起来的完整文章。所以这篇东西不打算讲教科书理论只围绕这份数据包把从解压到出成果的全过程掰开揉碎讲清楚。适合谁看如果你正在做焦作市域或河南境内的地形相关项目或者刚接触DEM、shp这类空间数据这篇能帮你省下大量试错时间。1. 内容整体设计与思路拆解为什么选30米DEM为什么附带shp边界拿到这份数据先别急着解压。先想清楚两个关键问题DEM为什么是30米分辨率shp边界文件又有什么用。这两点决定了你后续处理数据的整个思路。1.1 30米分辨率到底意味着什么先说分辨率。30米DEM意味着每个像元代表地面30×30米的区域。这个尺度放在地形分析里属于“中等精度”的通用配置。比它更粗的SRTM 90米适合大区域宏观分析精度不够看细节比它更细的12米、5米甚至1米数据虽然能看到更精细的地形变化但数据量呈几何级数上涨处理成本也高。焦作市域面积约4071平方公里30米分辨率大概对应450万到500万个像元这个规模在ArcGIS Pro、QGIS里跑起来非常流畅普通办公电脑完全带得动。为什么不是15米或更高从数据源角度看30米这个级别在国内外开源数据里特别常见像ASTER GDEM、ALOS AW3D30、Copernicus GLO-30都提供了30米左右的产品。这份数据大概率脱胎于公开的全球或全国DEM再经过裁剪、重投影处理形成焦作市级范围。换句话说它不是原始测绘成果而是经过加工、适合项目使用的成品数据。那30米能做什么比较典型的有全市尺度的坡度坡向分析、汇水区提取、洪涝淹没模拟、光伏选址的朝向评估、交通选线的纵断面概算。这些场景对精度要求没那么苛刻30米足够用。但要去做精确到村级的耕地平整设计或者一栋楼的日照遮挡分析30米就太糙了需要换更高精度数据。1.2 shp边界文件在数据包里的角色很多人忽略了这个附带shp边界文件的作用。这个市级范围shp可不是让你画地图好看的。它最重要的用途有三个。第一统一分析范围。DEM原始数据往往覆盖一个矩形区域包含了焦作周边的新乡、济源、晋城部分地区。如果你按DEM栅格直接做统计统计结果会把市域外的地形也算进去导致结果失真。正确做法是先加载shp边界用“按掩膜提取”把DEM严格裁剪到焦作市行政边界内再去算面积、算体积、算平均坡度这样结论才有行政口径上的意义。第二作为成果出图的裁剪框。做坡度图、坡向图、山体阴影图时不加边界直接出图图面会显得很散而且审阅者会问“你这是焦作市吗怎么边界都不对”。把shp边界叠加在图面顶层或者用来设置制图范围出图效果专业很多。第三做空间统计时的分组依据。如果你要按县区统计地形指标比如“解放区平均坡度多少修武县最大高程多少”那这份shp还能再细分。虽然这份shp是市级范围不一定带县界字段但你可以用“面转线”再结合其他数据做叠加至少保证市域边界是准确的。1.3 这份数据解决的核心痛点在没有这份数据的情况下你需要自己去下载全国或全球DEM然后做投影转换、范围裁剪、边界配准这一整套流程。自己动手不是不行但步骤繁琐而且很容易出问题——比如下载的DEM是WGS84坐标而焦作的shp是CGCS2000两者叠在一起差了几百米不经过投影转换直接分析结果全是错的。这份数据把最关键的一步“按市级范围裁好”做了拿到手可以直接进分析流程。不过别高兴太早我拿到这种包之后从来不会直接信任它的元数据。压缩包里的文件可能被人重投影过也可能坐标系统信息缺失甚至裁切边界和DEM范围之间存在半像元偏移。这些风险后面都要靠实际检查来排除。2. 数据包解压与基础检查从zip到可用数据每一步都不能省绝大多数人拿到zip的第一反应是双击打开把里面的东西拖出来完事。如果只是普通文档这么干没问题。但空间数据不同解压过程和数据完整性直接挂钩。这里把解压、文件识别、坐标检查三步一次讲透。2.1 linux和windows环境下正确解压zip如果你用的是Windows直接右键“全部解压缩”就行。但我更建议装一个7-Zip尤其是这个数据包里的DEM文件命名可能比较长Windows自带的解压工具偶尔会在长文件名上出错解压出来文件后缀是好的但栅格数据集文件.tif.aux.xml之类丢失导致后续加载出问题。7-Zip对中文文件名和长路径的支持更稳定。如果你在Linux服务器上做处理那解压命令要敲对。核心命令是unzip常规用法unzip 河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip -d jiaozuo_dem-d参数指定解压目录。如果文件较多用unzip -o强制覆盖用-unzip -q静默模式减少输出。还有个大坑有些数据包是分卷压缩的后缀是.z01、.z02加最后的.zip这种要先确保所有分卷放在同一目录然后只对最后那个.zip执行解压工具会自动识别分卷。如果在Linux下碰到“file is not a zip file”或者“could not find eocd”这类报错十有八九是文件没下载完整或者从QQ、网盘传输时被改了后缀。解决办法是用file命令看一眼真实文件类型file 河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip如果输出显示“Zip archive data”说明文件是好的如果显示“HTML document”或者“data”那基本可以断定它根本不是zip——要么下载链接跳到了错误页面要么传输截断了。至于网上流传的zip密码移除、超人zip解密助手之类的工具在这个数据包场景下用不上开源DEM数据不会加密如果解压要求输入密码先怀疑文件来源是否正规。2.2 解压后的文件组成与格式识别解压完成后你会看到两类核心文件。DEM部分大概率是GeoTIFF格式后缀.tif可能还带.tfw世界文件、.prj投影信息、.xml元数据等附属文件。记住GeoTIFF的所有附属文件必须和主文件在同一目录且保持文件名前缀一致缺一个都可能导致坐标系识别失败或加载异常。shp边界部分则是一组同名不同后缀的文件.shp几何数据、.shx索引、.dbf属性表、.prj投影、.cpg字符编码。这组文件里只有.shx和.dbf同时存在QGIS、ArcGIS才能正常打开缺了.shx软件会重新构建索引但会报警告缺了.dbf属性表就没了缺了.prj坐标系统乱了缺了.cpg属性表中文乱码。很多人问“shp文件导出的时候没有cpg文件是怎么回事”——这种情况常见于从某些老系统或AutoCAD导出的数据没有.cpg文件不代表不能用但中文属性字段可能出现乱码。真正要注意的是.prj缺失那才是大问题。解压完之后在QGIS里按顺序加载DEM和shp两个图层先看空间位置对不对。2.3 坐标参考一致性检查最容易翻车的一步把DEM和shp都加载进地图后第一步不是做分析而是检查两个图层的坐标参考。按下图层的“属性”查看“信息”里的“坐标系”一栏。理想状态是两者一致如果都是CGCS2000_3_Degree_GK_Zone_38或者都是WGS 84 / UTM zone 49N那好办直接干活如果DEM是WGS 84地理坐标系经纬度而shp是CGCS2000投影坐标系米两者虽然能叠在一起但DEM会被压扁成一条窄带或者两个图层相距很远这就是典型的坐标系不匹配。处理方式很简单在QGIS里对DEM执行“重投影”工具目标坐标系选shp的坐标系。不要用“设置图层CRS”这种强制改定义的方式那是把错误信息硬掰过来数值没变坐标系统标签变了结果更乱。正确的重投影是真正重采样每个像元输出新的GeoTIFF。还有一个必须要查的像元大小。原始DEM可能是度约0.000277778°或者米30m为单位。如果是度那“30米”这个说法就是近似值在焦作的纬度上一度经度大约90公里一度纬度约111公里要投影后才能得到标准的30米像元。这点直接决定坡度、坡向计算的正确性。如果发现像元单位是度务必先投影再用“重采样”工具把像元大小设成30米。3. 核心细节解析与实操要点DEM分析前必须做的数据预处理很多人拿到DEM就直接扔进“坡度”工具里跑然后跑出一堆奇怪结果——坡度全在85度以上或者山脊线完全变形。这不是数据错了而是少做了几个预处理步骤。这里把栅格计算前必须处理的三件事讲清楚无效值处理、填洼与水文分析的选择性操作、以及像元对齐。3.1 无效值NoData处理为什么你的DEM有洞DEM原始数据在个别区域会存在数据空洞主要是因为云遮挡、水体反射异常或者传感器盲区。在栅格数据里这些区域以NoData值存在。问题在于很多分析工具对NoData的处理策略不同有的直接忽略导致统计结果面积偏小有的把NoData当0计算出来的高程和坡度完全离谱。打开DEM属性查看“统计信息”里有没有出现比如-32768、-9999这类异常最小值。如果有必须先把它们转成NoData。在QGIS里可以用“栅格计算器”公式写((dem1 -9999) * dem1)意思是把所有大于等于-9999的像元保留原值小于这个值的都变成NoData。ArcGIS里对应的脚本是SetNull(dem -9999, dem)有些数据包自带的DEM是用浮点数存储的NoData可能是“nan”用整数存储的NoData常见值是-9999或-32768。做完转换后再统计一次确认最小值在合理范围内。焦作的地势最低处约80米黄河流域最高处约1300米太行山如果最小值出现负数或者超过2000的异常值那要么是NoData没处理干净要么是投影转换时出的问题。3.2 填洼还是填平决定你后面能不能顺利做水文分析在水文分析中填洼是一个必经步骤。真实地表存在很多坑洼这些洼地在DEM里表现为局部低点水流进去就出不来导致提取出来的河网断断续续。如果用这个DEM直接做流向分析你会得到一堆环形流向或者断流河段。所以做河网提取、汇水区划分之前必须“填洼”。QGIS里填洼工具叫“填洼wang–liu”ArcGIS里叫“填”。填洼分两步先填一次对比填洼前后的高程差如果最大填洼深度超过合理范围焦作地区一般不会超过100米除非是人工矿坑就要检查是不是原始DEM里有错误像元。焦作是矿区历史上煤炭开采形成了一些采空区和沉陷区这些地方的DEM可能表现出大量局部洼地填洼后地形起伏被抹平分析出来的汇水区和实际偏差会很大。这种情况下建议把填洼后的DEM和原始DEM做差值生成“填洼深度栅格”在出图时标注哪些区域是人工干预过的。还有一个细节填洼工具只处理NoData以外的区域所以第3.1小节的NoData清理得先做否则填洼结果会出现接缝。顺序错一步结果全废。3.3 像元对齐DEM和shp叠加溢出的隐形坑如果你把DEM重投影到投影坐标系后直接做“按掩膜提取”偶尔会发现边界处出现一圈无数据带或者提取出来的栅格范围比市界大了几十米。这是因为DEM像元网格和shp矢量边界没有精确对齐。解决办法是用“对齐栅格”功能。QGIS里在“栅格-对齐-对齐栅格”中勾选DEM设定参考网格的范围和像元大小强制让DEM的像元网格与shp边界的外包矩形对齐。实际操作时有一个技巧先做一个空的矢量网格实际上是矢量面让它覆盖焦作市界并向外扩展半个像元15米然后以这个网格的范围作为DEM裁剪范围这样裁剪后边界最干净。ArcGIS里则用“裁剪栅格”工具勾选“使用输入要素裁剪几何可选”几何范围输出后就不会有多余外扩。说到底这一步做的就是“确保分析范围的每一条边界线都恰好落在像元边框线上”为后面的坡度、坡向、统计计算打好基础。别嫌麻烦这步省了后面做面积统计时误差能达到几公顷。4. 实操过程与核心环节实现以焦作市为例跑一次完整地形分析现在假设数据已经解压、坐标检查完毕、NoData也清理好了接下来直接进入核心实操环节。这一节围绕一份真实的DEM数据走一遍坡度坡向分析、山体阴影生成、以及等高线提取的完整流程。所有步骤都用QGIS和ArcGIS双版本对照给出方便你在自己熟悉的环境里复现。4.1 坡度坡向分析参数怎么设才合理坡度分析的意义在于从DEM中计算每个像元的最大坡度变化率坡向则给出该像元面向的方向。在QGIS中工具位于“栅格-分析-坡度”。关键参数有三个坡度单位选“度”不是“百分比”Z因子填1这个参数是垂直方向与水平方向的单位比例因子如果DEM是以米为单位的投影坐标系Z因子保持1就对了算法选择“Horn”它比Zevenbergen-Thorn算法更平滑适用于自然地表。以焦作为例北部太行山区坡度可达30到50度中部山前丘陵坡度5到15度南部平原基本小于3度。跑完之后图例色带建议用黄到红渐变这样山区和平原的区分度很高。坡向分析同理工具就在坡度旁边输出结果中0度代表正北90度是正东360度又回到正北。坡向在焦作的应用场景很典型光伏板的朝向选择和农业区划中避风向阳坡的识别都需要坡向数据参与叠加分析。ArcGIS里对应操作是“3D Analyst工具-栅格表面-坡度”输出测量单位可选“DEGREE”Z因子填1。有一个小技巧在ArcGIS中给结果图层启用“山体阴影渲染”来叠加显示坡向分布会一眼看清。跑完坡度坡向后把结果和原始DEM一起按第3.3小节的掩膜边界裁剪确保斜坡面积统计严格落在焦作市域范围内。用“分区统计”工具QGIS里是“分区统计”ArcGIS里是“区域分析-以表格显示分区统计”按shp边界统计平均坡度就能得到全市平均坡度、最大坡度、平均坡向这些数值直接导出成csv后续写报告用得上。4.2 山体阴影与可视化出图让地形图有“立体感”出图是DEM分析里最能体现专业度的一环。山体阴影Hillshade模拟太阳光照射下的地表明暗变化可读性极强。QGIS工具在“栅格-分析-山体阴影”核心参数是太阳方位角默认315度即西北方向和太阳高度角默认45度。为什么默认值要选315度因为西北方向来的光照能最大程度突出正南北走向的山脊线这是地形可视化里经验值不要随意改成0度正北光照那样山脊线阴影会缺失。高度角45度用来保持明暗对比适中。如果你追求更柔和的效果可以生成三个不同方位角270度、315度、360度的山体阴影然后按50%透明度叠加效果接近专业地图产品的“多向山体阴影”。ArcGIS的“山体阴影”工具参数相同输出时建议勾选“模拟阴影”让阴影区域更清晰。山体阴影生成后把它放在DEM彩色渲染图层的下方把DEM的透明度调到60%-70%叠加效果就是一个立体的彩色地貌图。放大到焦作北部太行山区时能看到明显的断层崖和峡谷地貌放大到南部平原则是一片平缓的绿色。这种图直接截图放进PPT或报告里比纯色DEM图专业得多。还有一招适合做汇报展示把DEM图层混合模式设置成“叠加”或“正片叠底”放在山体阴影上方。这样既能保留高程色彩信息又有阴影立体感整个过程不涉及任何额外计算但视觉提升非常大。4.3 等高线提取与shp转3D应用从栅格到矢量的一次跨越如果报告需要等高线图直接从DEM提取即可。QGIS里用“栅格-提取-等高线”间隔参数要看图幅比例。焦作市域1:50万图用50米等高距1:20万图用20米如果只做太行山区的细部图10米等高距比较合适。ArcGIS对应工具是“3D Analyst工具-栅格表面-等值线”。等高线提取完成后是一堆shp线要素。有时候需要结合dxf在CAD里继续制图那就执行“导出-保存要素为”格式选AutoCAD DXF。ArcGIS里则用“转换工具-转为CAD-导出CAD”。网上常说的“批量把多个shp转为CAD”其实也走这条路在ArcGIS里先把多个shp合并工具叫“合并”再统一导出DXF。CAD打开后要注意比例尺和坐标单位如果导入CAD后图离得太远多半是投影坐标系在CAD里没被识别用“MOVE”命令把图移到原点附近或者用“MAPCSLIBRARY”工具重新指定坐标系。如果想更进一步把DEM转成3DTiles用于Web端三维展示推荐用QGIS插件“qgis2threejs”先看效果再到CesiumLab里做切片。这里有一个容易被吐槽的坑直接拿DEM生成的3DTiles因为30米精度有限在大比例尺下会看到明显的“马赛克台阶”所以Web端三高展示时要结合LOD策略远处显示DEM近处叠加倾斜摄影模型效果更真实。4.4 基于DEM的流域提取与淹没模拟扩展做完基础分析之后利用这份DEM做一次小流域提取体验会更好。按水文分析的流程走一遍填洼后的DEM计算流向D8算法再算流量累积设置流量阈值为1000个像元后提取河网再以出水口为点生成汇水区。焦作境内的主要河流像大沙河、山门河它们的源头就在北部太行山区利用DEM自动提取的河网走向和实际水系对比匹配度很高。淹没模拟则借助“积水区”分析或者“GIS-SWAT”插件来做。30米DEM做淹没模拟的精度有限但对大范围洪泛区识别是够用的。模拟结果出来后叠加shp边界内的村庄点位就能判断哪些村庄位于洪水风险区。市面上还有把DEM用来做“太阳辐射分析”的插件如QGIS的SAGA“日辐射”工具这个在光伏选址中非常好用输入DEM和纬度焦作市约北纬35.2度输出全年太阳辐射总量栅格能直接指导光伏板倾角设计。5. 常见问题与排查技巧实录处理这份DEM时最容易踩的坑数据本身再干净工具链上该出的问题一个都跑不掉。把这些年处理地方DEM数据的经验总结成一份速查表再挑几个高频问题展开说。5.1 打开shp却一片空白是什么原因shp加载后地图窗口没有任何图层显示多半是几何类型或坐标系出了问题。先确认shp在QGIS里是否自动识别出WKT投影信息如果没识别到设置里手动指定。还有一种情况文件本身是空几何用QGIS的“修复几何”工具可以尝试重建几何但如果是原始矢量数据就是空的那没办法只能找数据源重下。另外检查shp配套的.dbf文件如果属性表能打开但不显示几何那可能是.shp和.shx索引不匹配用“Shapefile修复工具”跑一遍十有八九能救回来。5.2 加载DEM颜色怪异全是黑白灰块怎么办DEM加载后没有自动拉伸显示成一片灰这是很多新手会遇到的。QGIS里右键图层属性在“符号系统”选项卡选择“单波段假彩色”配色方案选“Spectral”或者“Terrain”然后点“分类”软件会自动计算最小最大值。ArcGIS里右键图层属性-符号系统-拉伸类型选“百分比截断”直方图拉伸选“2%到98%”这样可以避免个别异常像元把整个配色带偏。如果拉伸后仍然异常检查统计信息里有没有出现极大的孤立值比如一座山上出现9999这种说明NoData没清理干净回到第3.1小节重新处理。5.3 “invalid zip archive: could not find eocd”和“file is not a zip file”怎么破这两个报错是zip解压领域的经典问题。EOCDEnd Of Central Directory记录位于zip文件末尾解压器找不到它基本可以认定压缩包不完整。常见场景有两个一是从网盘下载了文件但下载过程中断或者保存时被浏览器当成HTML页面下载下来二是QQ/微信传文件时下载到一半中断文件大小比原始文件小了几KB到几MB。排查方法很简单用file命令检测文件头或者把后缀改成.rar试试看是不是传输工具把压缩格式改了。还有一种特殊情况如果包里有.z01分卷主zip包里缺少EOCD需要把全部分卷放在同一目录再尝试解压。解决不了就找数据源头重新下载不要浪费时间尝试修复。5.4 中文文件名和路径引发的连环坑这份数据文件名里有“河南省焦作市”这样的中文字符很多工具在读取时不会报错但在中间过程中一旦生成临时文件路径就可能出问题。特别是ArcGIS里的“输出位置”如果选择了带中文或空格的目录工具经常弹“创建要素类失败”之类莫名的错。经验是在开始处理前把解压后的目录统一改成纯英文路径比如D:\gisdata\jiaozuo\dem.tif不要有中文目录名。文件名也可以简化把DEM改名为jz_dem_30m.tifshp改名为jz_boundary.shp。这不是强迫症是省时间的预判。涉及后续Python自动化分析时更是要全程保证ASCII路径否则arcpy或gdal大概率在第一步就挂掉。5.5 与坐标和高程基准相关的隐蔽错误最后提醒一个最隐蔽的坑高程基准问题。我国使用的国家高程基准是1985国家高程基准正常高系统而全球尺度的DEM数据如Copernicus GLO-30默认的是EGM2008大地水准面或者WGS84椭球高。两者在部分地区可能相差几十米。焦作这种地处中原、地形起伏相对平缓的地区差距通常在10到20米之间但如果你拿这份DEM和现场RTK测量值对比会发现系统偏差。处理方式是在报告中注明高程基准或者根据实测控制点计算改正数对整个DEM进行批量平移。很多人把精力放在平面坐标上忽略了高程基准最后做土方计算时才发现全错了那时候返工代价就大了。6. 从这份数据还能往哪些方向扩展拿到手的DEM和shp只是起点真正有意思的是基于它们可以延伸出哪些应用方向。最后分享三个我实际做过、并且觉得费效比很高的扩展方向。第一个是地形起伏度分析。在QGIS里用“焦点统计”计算DEM在3×3窗口内的最大最小高差输出地形起伏度栅格然后按shp边界统计各分区平均起伏度。这个指标对城市扩张评估、农业机械化可行性判断特别有用。焦作南部平原起伏度基本在10米以内北部山区可能超过100米一套数据直接支撑了“宜建区”和“不宜建区”的空间划分。第二个是结合多源数据做地质灾害易发性初判。把DEM生成的坡度、坡向、地形湿度指数TWI叠加断层分布、降雨量、植被覆盖度用加权叠加法跑一个易发性分区图。这份30米DEM可以作为基础数据源虽然精度不高但用来做全市尺度的风险识别完全够用。第三个方向比较冷门但实用把DEM转成3D打印用的高度图。QGIS里把DEM导出为灰度PNG再用图像处理软件生成STL模型焦作的地形就能被“打印”出来。这对于规划展示、教育基地的实体沙盘制作非常有用。这些扩展方向不需要新的数据源一份DEM、一个shp边界、一个开源GIS软件就能启动。数据本身就是工具用得好不好全看你愿意在上面花多少心思。本文还有配套的精品资源点击获取