简介本资源为河南省焦作市30米分辨率数字高程模型DEM地理信息数据包面向GIS初学者、地理信息专业学生、城乡规划及环境评估从业者解决区域地形分析、坡度坡向计算、流域划分与三维可视化等基础空间分析需求。压缩包共12个文件包含核心高程数据文件焦作市dem.tif、配套市级行政边界Shapefile.shp/.shx/.dbf/.prj等及辅助文件.ovr缩略图、.tfw地理配准、.xml元数据全面支持ArcGIS、QGIS等主流平台直接加载与分析。资源大小仅11.79MB轻量易用结构规范符合GIS数据标准组织方式。目前已有647人学习下载用户可直接调用dem.tif开展地形分析结合shp文件实现空间叠加与属性关联无需额外数据预处理是开展焦作市尺度地貌研究与教学实践的即用型基础数据集。 当初在找河南省内地级市的数字高程数据时我第一个想到的就是焦作。这座城市北依太行、南临黄河地形从山区到平原的过渡非常典型做水文分析、坡度分级、城市选址模拟都特别合适。而一份“河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip”名字看上去就是一个标准的数据包——30米分辨率DEM栅格加上市级边界的矢量范围解压之后既能直接出图也能用来裁剪更大范围的原始数据。这篇文章我就围绕这种数据包从内部结构、读取方式、坐标系处理、实际分析场景到常见报错完整过一遍。无论你是刚接触GIS的学生还是经常跟栅格数据打交道的从业者都可以照着操作。1. 拆开压缩包之前先弄清楚DEM和shp各自是什么角色1.1 30m分辨率到底意味着什么DEMDigital Elevation Model数字高程模型是一张用栅格像元记录地面高度的数据表。每个像元里存一个数值代表该位置的海拔单位一般是米。30m分辨率的意思是一个像元对应地面上30米乘30米的真实区域也就是每个像元覆盖900平方米。放到焦作这种地形起伏明显的区域30m足以把太行山前的丘陵、山间谷地、黄河冲积平原的轮廓刻画得比较清楚。相比之下90m分辨率的SRTM数据在平缓地区会显得“平”小沟壑和低矮岗地经常被平滑掉而12.5m或更高精度的数据虽然细节更好但数据量成倍增加处理起来对电脑配置要求也高。所以30m是一个性价比很高的平衡点——既能看清地形骨架又不至于让一张全市范围的影像把内存占满。1.2 市级范围shp文件是干什么用的shp是ESRI Shapefile矢量格式它存的不是连续的高度场而是离散的边界、点、线、面。这份压缩包里带上“市级范围shp文件”意思就是你不用自己去网上找焦作市的行政边界直接拿它当分析范围使用。在GIS软件里你可以把这个shp当作一个“蒙版”用来裁剪DEM只保留焦作市辖区内的地形用来统计每个区县的平均海拔、最高点、坡度分布用来叠加道路、河流、居民点等其他矢量数据限定分析只在市界内进行。许多人在处理DEM数据时最容易忽略的一件事就是原始DEM往往覆盖了比目标区域大得多的范围。比如从公开数据源下载的全球或全国拼接数据一打开可能就是整块几千乘几千像素的大影像。直接用会非常卡而且分析结果也包含范围之外的区域。用shp去裁剪是标准工作流里的第一步。1.3 焦作的地形特点对这份数据的使用价值焦作市地处河南省西北部地势总趋势是北高南低。北部是太行山余脉海拔可以达到一千米以上山势陡峻沟谷发育中部是山前丘陵和洪积扇群南部则是黄河北岸的冲积平原海拔只有几十米。这种“一市之内落差过千米”的地形格局让焦作的DEM数据在各类实验中显得特别有代表性。如果你做坡度分析能同时看到陡峭的山区和平缓的平原如果做水文分析从北部山区发源的河流向南流经平原能清楚展现汇水路径和河网的发育过程。很多教程里用沿海城市或者纯平原城市的数据效果往往不如焦作这种“半山半原”的结构鲜明。2. 打开数据包的正确步骤解压、编码、坐标系统统别踩坑2.1 稳妥的解压方式与常见解压报错拿到“河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip”这个文件后第一步自然是解压。Windows用户可以右键选择“全部解压缩”也可以安装7-Zip或Bandizip。Linux服务器上一行命令就能搞定unzip 河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip -d jiaozuo_dem/如果是在命令行环境下很多朋友会遇到“file is not a zip file”或者“End-of-central-directory signature not found”这类报错。这通常不是解压软件的问题而是文件本身没下载完整。尤其是用浏览器直接下载较大文件时网络中断或者代理缓存可能导致文件缺失尾部索引块。判断方法是先看文件大小——如果压缩包显示的大小和下载页面标注的源文件大小明显不一致八成就是下载中断了需要重新下载。另外还有一种特殊情况文件后缀虽然叫zip但实际可能只是把其他格式直接改了名这时候用file命令看一下真实类型file 河南省焦作市DEM数字高程数据30m含市级范围shp文件.zip如果输出显示是“Zip archive data”说明是正常的压缩包如果显示HTML文档或纯文本那多半是下载到了一个错误页面。2.2 shp文件不是“一个文件”而是“一组文件”很多第一次接触shp的人会困惑解压后怎么冒出来一堆同名但后缀不同的文件从严格意义上讲一个完整的Shapefile应由这几个文件共同组成后缀作用是否可缺失.shp要素几何信息真正存储点线面坐标不可缺失.shx几何索引加速读取不可缺失.dbf属性表存储要素的字段值不可缺失.prj坐标系描述文件建议保留.cpg属性表字符编码声明建议保留.sbn / .sbx空间索引非必需可缺失这个机制里最常被忽略的是.prj和.cpg。没有.prjArcGIS和QGIS会不知道这个面边界是经纬度还是投影坐标导致加载后位置错乱没有.cpgdbf属性表里的中文或特殊字符容易显示成乱码。所以拿到数据包后先检查一下这几个文件是不是都在。如果发现缺失可以用其他途径补充坐标信息——比如对比数据包说明文档里写的坐标系然后在软件里手动指定。2.3 DEM栅格读取时最容易忽略的NoData值DEM栅格在存储时并不是每一个像元都有效。比如焦作市范围之外、或原始数据拼接时出现的无值区域软件会用NoData标记。很多人在ArcGIS里加载tif后看到一大片黑色或者异常值第一反应是数据损坏其实只是没有正确设置NoData。在ArcMap或ArcGIS Pro里可以在图层属性的“栅格信息”里查看NoData值是多少在QGIS里则可以在图层属性“透明度”选项卡里把额外的NoData值设置为透明。这个细节处理不好会直接影响后续分析。例如做坡度分析时如果NoData没有被识别这些“伪像元”会被当作极高的异常数值参与计算导致坡度结果里出现莫名其妙的红色区域。做填洼、流向分析时更是如此一个NoData点就可能切断整个汇水网络。所以拿到数据包的第一件事不是急着出图而是先确认有效值范围和NoData设置。3. 在GIS软件里加载这份焦作数据两种主流工作流3.1 ArcGIS系列从目录窗口拖入图层ArcGIS是老牌GIS平台处理这种“一个DEM加一个shp”的组合非常方便。打开ArcMap或者ArcGIS Pro后在Catalog面板里找到解压后的文件夹直接把tif文件拖进地图视图再把shp拖进来。这里我需要提醒一点务必将shp和DEM都加到同一张地图里再对比位置关系。如果shp能和DEM的山形轮廓贴合那说明坐标系正确如果shp跑到了很远的地方或者形状被拉伸得不成样子那就是坐标系不匹配。对于焦作市的数据如果DEM是地理坐标系GCS_WGS_1984而shp是投影坐标系比如CGCS2000 3-degree Gauss-Kruger zone 38两者直接叠加很可能出现位置偏移。解决方案是在后期处理时统一坐标系在ArcToolbox中找到“Data Management Tools - Projections and Transformations - Raster - Project Raster”将DEM投影到和shp一致的坐标系或者反过来用“Define Projection”先给缺失坐标系的文件定义正确的参考系统再做投影转换。处理完坐标系后可以做一件很有用的事——用shp裁剪DEM。找到“Spatial Analyst Tools - Extraction - Extract by Mask”输入DEM和焦作范围shp输出就得到了一张刚好覆盖市界、范围外全是NoData的裁剪后栅格。后续所有分析都用这个裁剪后的数据运算速度会明显提升。3.2 QGIS方案免费开源也完全够用如果你不想折腾授权QGIS是一个能打的开源选择。直接“图层 - 添加栅格图层”选择tif“添加矢量图层”选择shp。QGIS对坐标系的判断通常比ArcGIS更“透明”——如果发现图层范围、单位不对可以从图层属性里看到当前坐标系详细信息再使用“导出 - 另存为”对话框中重新指定CRS方式完成坐标转换。在QGIS里裁剪DEM也很快找到工具箱里的“裁剪 - 按掩膜图层裁剪栅格”输入层选DEM掩膜层选焦作市shp点击运行即可。QGIS的批处理功能比ArcGIS更直观假如你想把市界的多个分区shp分别裁剪DEM可以右键工具选择“批量处理”一次配置所有参数自动循环运行。3.3 出图前先做一版漂亮的渲染默认加载进来的DEM往往是一张灰度图很多没接触过地形渲染的人会觉得“这有什么好看的”。其实只要做两步地形立刻立体起来第一步生成山体阴影Hillshade。在ArcGIS的空间分析工具里找“Hillshade”或者在QGIS里用“地形分析 - 山体阴影”。山体阴影本质是模拟一个光源照射地表面根据坡向和坡度计算明暗程度。第二步把原始DEM和山体阴影叠加显示。在图层属性里把DEM的符号系统设置为“拉伸”或“分级着色”再把透明度调到70%-80%然后叠加在Hillshade之上。这样既有颜色表达海拔高低又有阴影表达地形起伏视觉效果接近测绘部门的标准晕渲图。这套操作我几乎每次处理DEM都会做因为不管是做PPT汇报还是写技术报告一张高清晕渲图都比一堆数字有说服力得多。4. 焦作DEM数据的典型应用从坡度分析到水文提取4.1 坡度与坡向分析山地开发的基础指标焦作北部山区面积不小坡度是判断土地是否适合建设、农业开发潜力如何、地质灾害风险高低的核心因子。在ArcGIS里找“Spatial Analyst Tools - Surface - Slope”输入裁剪后的DEM输出量纲选择“DEGREE”就能得到每一个像元的坡度角。产生的坡度栅格与行政区shp叠加后可以用“Zonal Statistics as Table”分区统计每个乡镇的平均坡度、最大坡度、坡度大于25度的面积占比。这种数据对城乡规划、林业部门、自然资源调查都很有参考价值。坡向分析同样简单用“Aspect”工具即可生成从0到360度的坡向栅格。结合焦作的地形特征北坡和南坡的日照、蒸发差异很大山区农业生产中经常要用到坡向来选择种植区。导出坡向栅格后还可以重分级为平坦、北、东、南、西等九个方向方便统计各方向的面积。对于30m的数据坡度和坡向结果在小范围内可能不如高精度LiDAR数据准确但做区域尺度上的宏观分析完全够用。需要特别留意的是在平原地区坡向会显得比较“碎”因为高程差太小噪音信号都被放大了这时可以对坡度结果做一次低通滤波或者先对DEM做轻度平滑。4.2 等高线提取老测绘人熟悉的操作DEM的一个经典应用是生成等高线。用“Spatial Analyst Tools - Surface - Contour”工具设置等间距比如20m或50m就能从焦作DEM中提取出一整套等高线矢量数据。提取后的等高线可以直接转成CAD格式在AutoCAD里做进一步的规划制图这个流程在工程领域应用极广。有一个实操经验值得分享如果发现提取出的等高线在平原地区过于密集、互相交叉或者呈锯齿状通常是因为原始DEM在平缓地区存在微小的台阶效应。遇到这种情况可以在提取前先做一步“Focal Statistics”选择一个3x3或5x5的矩形邻域统计类型为“MEAN”对DEM做一个平滑处理。平滑之后的等高线会顺眼很多当然代价是局部微小地形被弱化所以到底做不做平滑取决于出图用途。4.3 水文分析从填洼到河网提取的完整链条焦作市南部紧邻黄河北部山区的降水会通过山间沟谷汇入河流。利用DEM做水文分析可以自动提取河网、划分汇水区这对山洪灾害风险评估、水资源管理都非常有价值。ArcGIS的Hydrology工具集是一个相当经典的工作流填洼Fill把DEM中的凹陷区域填平保证水流能连续向下流动流向Flow Direction用D8算法计算每一个像元的水流方向汇流累积量Flow Accumulation统计有多少上游像元的水会流经当前像元设置阈值提取河网比如把汇流累积量大于1000的像元提取为河道栅格河网转矢量Stream to Feature得到河流线文件。这里最容易出错的是“填洼”这一步。很多人在ArcGIS里做填洼后发现生成的河网在平原地区变成一条直线或者乱窜这是因为平原地区高程差异太小水流方向的判定几乎取决于DEM里的微小误差。对此可以使用“Cost”填充法或者结合现有河网数据进行“burn-in”——把真实河道的高程人为降低几个像元逼迫流向算法沿着真实河道走。对于30m数据我个人建议在提取河网时将“流域分析”和“河网提取”的阈值调高一些比如取5000或10000否则提取出来的河道太多了和真实水系对比会很乱。4.4 通视分析与视域分析选址与景观评价的利器如果再往深走一点30m DEM还可以用于通视分析。比如在焦作的山区选择通信基站站址、观景平台位置或者评估风力发电机的可视影响范围。“Viewshed”工具就是为此设计的。输入观察点位置和DEM软件会计算该点能看到的区域范围。实际操作时观察点高度参数很重要——如果是人眼观察一般设置1.6米到2米如果是通信塔则要设置几十米甚至上百米。对于30m分辨率的数据通视分析的结果比更精细的数据要粗糙一些但在宏观选址阶段已经足够。如果你发现算出来的可视范围大得离谱通常是因为没有考虑地球曲率和大气折射在软件里开启“使用地球曲率”选项可以改善。4.5 将DEM与shp结合做面积-高程积分如果做地貌学研究面积-高程积分是一个非常有用的指标用来判断一个区域的地貌发育阶段。操作思路是先按shp范围裁剪DEM然后在GIS里把DEM重分类为若干等间距的高程带统计每一带占区域总面积的比例再画出高程-面积曲线。焦作这种山区到平原的过渡地带通常能算出明显的“少年期-壮年期”地貌特征。虽然这个分析不复杂但在论文和报告中用起来很加分。5. 进阶玩法Python批处理、格式转换、以及几个绕不开的坑5.1 用Python读DEM和shprasterio与geopandas组合当数据量变大、或者需要重复处理多个城市的DEM时光靠鼠标点GIS界面效率太低了。用Python可以做一个稳定的批处理流水线。读取这份焦作DEM可以用rasterioimport rasterio from rasterio.mask import mask import geopandas as gpd # 读取DEM with rasterio.open(jiaozuo_dem_30m.tif) as src: dem src.read(1) profile src.profile bounds src.bounds print(像素大小:, profile[width], x, profile[height]) print(地理范围:, bounds) # 读取焦作市shp jiaozuo gpd.read_file(jiaozuo_city.shp) print(jiaozuo.crs) # 按shp裁剪DEM with rasterio.open(jiaozuo_dem_30m.tif) as src: out_image, out_transform mask(src, jiaozuo.geometry, cropTrue) out_profile src.profile.copy() out_profile.update({ height: out_image.shape[1], width: out_image.shape[2], transform: out_transform }) with rasterio.open(jiaozuo_dem_clipped.tif, w, **out_profile) as dst: dst.write(out_image)这段代码的逻辑很直接先读出DEM的基本信息再读取shp矢量最后用shp的几何范围去裁剪DEM。如果你有十几个地级市的数据需要批量处理写一个循环就能一次跑完剩下的时间大可用来调试出图样式。5.2 导出为其他格式从tif到3DTilesGIS领域里把DEM和shp做成三维场景越来越常见。比如用“shp转3dtiles”的思路把焦作市的建筑边界或行政区边界转换成3DTiles瓦片然后在Cesium或本地三维地球里加载。30m的DEM可以转化成带地形高度的Tile模型配合倾斜摄影或白模一起展示。转换工具有很多常见的是用CesiumLab或基于cesium的定制脚本如果是纯代码路线可以考虑用py3dtiles处理点云和模型不过学习成本稍高。需要提醒的是DEM转3DTiles后文件体积往往比原始tif大好几倍而且瓦片金字塔的层级选择非常关键。层级过多会让加载变慢层级过少则细节丢失。以焦作全市范围为例30m数据通常生成到14-16级就足够再高就是浪费。5.3 DSM和DEM的区别什么时候不要混用热搜词里有“dsm生成dem”这个说法这里顺便理一下。DSM数字表面模型记录的是地物表面高度包含建筑物、树冠等DEM记录的是裸地高程。很多开源数据源默认给的是DSM比如部分卫星立体像对提取的高程数据。如果你要用DSM做水文分析得先做滤波把地表覆盖物从高度场中移除。这个过程没有统一标准常用的手段包括形态学开运算、局部最小值滤波、或者结合土地覆盖数据分类剔除。所以如果你下载的数据源没有明确标注“bare earth”或“terrain”最好先判断一下它到底是DEM还是DSM。对于本主题里的这份焦作30m数据如果是从经典DEM数据源如SRTM、Copernicus GLO-30派生而来那么基本可以当作裸地高程使用否则要谨慎做验证。5.4 中文路径与编码问题一份带中文标题的zip包解压后的文件夹名和文件名也往往包含中文。在ArcGIS里打开中文路径下的tif或shp通常问题不大但有些第三方工具、Python库在读取含中文路径的文件时会直接报错比如常见的“No such file or directory”。解决办法也很直接在使用命令行或脚本前先把文件重命名为纯英文路径。比如解压后把文件夹改成mv 河南省焦作市DEM数字高程数据30m含市级范围shp文件 jiaozuo_dem mv jiaozuo_dem/*.tif jiaozuo_dem/dem.tif mv jiaozuo_dem/*.shp jiaozuo_dem/jiaozuo.shp这样不仅避免Python跨平台时的编码问题也方便在各类开源工具之间自由切换。5.5 如何判断一份DEM数据包可不可信最后聊一个很多人忽略的问题从网上下载的DEM质量怎么把关我拿到任何DEM数据之后会做三件验证加载后读取几个已知地物点的高程值。比如焦作北部山区的一些山峰在在线地形工具或纸质地形图上有明确海拔对比一下是否接近生成一版山体阴影肉眼观察山脊线和沟谷是否自然检查数据范围内有没有明显的接边痕迹或条带噪声。如果在影像里看到成片的斜纹、方块状边界可能是原始数据经过了较差的重采样使用时要谨慎。对于30m数据来说焦作这种地形的参考验证比较容易——因为太行山前的地形特征很明显随便找一个山谷或者山峰和测绘部门公开的高程做对比就能判断出数据是否合格。选择DEM数据和shp边界数据我现在的习惯是数据下载后先做一次完整性检查再做一次坐标系确认最后才进入正式分析流程。这个过程不复杂却能省掉后面大量的返工时间。如果你手头也有一份类似命名的数据压缩包可以参考我的流程逐步操作从解压到应用一次走通。本文还有配套的精品资源点击获取