简介本资源为山东省烟台市30米分辨率数字高程模型DEM地理信息数据集面向GIS初学者、地理信息专业学生及城乡规划、环境分析等领域的实践者用于开展地形可视化、坡度坡向计算、水文建模与三维地形渲染等基础空间分析任务。压缩包共12个文件包含核心栅格数据烟台市dem.tif及其辅助文件.tfw、.ovr、.aux.xml等以及完整市级行政边界的Shapefile矢量数据.shp、.shx、.dbf、.prj等确保坐标系统准确、属性结构完整、GIS软件开箱即用。资源大小31.8MB格式标准兼容ArcGIS、QGIS等主流平台。目前已有884人学习下载提供从数据导入、投影定义到基础地形分析的全流程支撑特别适合课堂实训、课程设计及中小尺度区域地理分析项目使用。 拿到这个数据包的时候我第一反应是这玩意儿太典型了。很多做GIS、做规划、做地质、做农业区划的朋友日常工作里打交道最多的就是这类“XX省XX市DEM数据30m含shp范围”但真正能把这套数据用明白的人其实不多。有人解压完就不知道下一步该干嘛有人在ArcGIS里打开发现一片黑有人裁剪完发现坐标对不上还有人连shp文件缺了cpg导致中文乱码都搞不清楚。这篇我拿烟台这个包当例子把从zip解压、shp读取、DEM裁剪到坡度坡向、等高线、山体阴影这些实际产品制作链路完整捋一遍顺便把这几年踩过的坑一并吐出来。先说清楚这套数据能干什么烟台市30米分辨率的数字高程模型配上本市级行政边界shp是典型的“一个范围一张高程底图”组合。DEM解决的是“地表起伏长什么样”的问题shp解决的是“在哪里做分析”的问题两者一配地形分析、可视域分析、淹没模拟、工程选址、坡度分级、等高线生成、三维地形可视化全都能做。适合的人群很明确GIS初学者拿它练手规划/测绘从业者拿它做前期底图搞环境建模的人拿它当地形输入还有做三维可视化的人拿它生成地形网格。下面我按实际操作的顺序来聊从打开zip到最终出成果每一步的“为什么”都说清楚。1. 先搞清楚你手里的数据到底是什么1.1 DEM是什么30m分辨率又是什么概念DEMDigital Elevation Model就是数字高程模型说白了是一张用数值记录地面高度的栅格图。每个像素存一个海拔值单位一般是米。你把它想象成一个巨大的Excel表每一格对应地面上一个小方块格子里填的是这个小方块的平均海拔高度投影后放到地图上就成了一张能反映地形起伏的灰度图——高亮的地方高暗的地方低。30m分辨率的意思就是地面上一个像素点对应30米×30米的实际面积。这个概念对后续所有分析都至关重要因为它直接决定了你能做什么、不能做什么。举个例子烟台市辖区面积大约是1.39万平方公里30米分辨率的话每个像元是900平方米一景数据大约有1500多万个像元。这个精度对区域性地形分析完全够用但你要是想分析一条10米宽的冲沟那就不行了——它根本显示不出来。那为什么30m这么流行原因很实际它是SRTM和ASTER GDEM这两大全球公开数据源的常见分辨率后来Copernicus DEMGLO-30也是30米。很多省市级基础地理信息平台对外发布的DEM也是30m。再往上是12.5mALOS PALSAR、10mCopernicus GLO-10再往上是1m、0.5m这种商业级LiDAR数据但获取成本和对机器性能的要求完全是另一个量级。所以30m是“免费开源数据里最能打的分辨率收费数据里最基础的分辨率”属于性价比最均衡的选择。1.2 烟台这个数据包的结构拆解从标题来看这个zip里至少包含两组东西30m分辨率的DEM栅格文件大概率是tif或img格式烟台市本市级范围的shp文件包含面要素用于限定行政边界“本市级范围”这个说法很关键。烟台市在行政区划上是地级市下辖多个区县芝罘、福山、牟平、莱山、蓬莱、龙口、莱阳、莱州、招远、栖霞、海阳以及长岛等。如果你拿到的是“本市级范围”而不是“县域范围”意味着shp里可能是多个面要素每个区县一个图斑也可能全部合并成一个整体轮廓。这个细节在后续裁剪DEM时会直接影响操作方式——你到底是按整个市域一刀切还是按每个区县分别提取得看具体需求。另外还有一层隐含信息标题里只写了shp但shp文件其实不是一个文件。一个完整的shapefile至少包括.shp几何、.shx索引、.dbf属性表还可能包括.prj投影、.cpg字符编码、.sbn/.sbx空间索引、.xml元数据等。很多新手只把.shp拷走结果打开报错这个我在后面问题排查部分细聊。1.3 为什么需要shp范围文件有人可能会问DEM本身就是一张覆盖整个区域的栅格图为什么还要额外的shp范围答案在于分析效率和数据组织方式。第一DEM往往覆盖的范围比你要研究的范围大。比如你拿的是覆盖烟台市整市的30m DEM但你要做的项目可能只涉及福山区一个工业园区的选址这时候你直接对整个DEM做坡度分析、填洼计算吃的是不必要的算力。用shp裁剪出研究区再后续处理计算量可能直接缩小到原来的十分之一甚至更少。第二shp范围是做掩膜提取的基础。在ArcGIS里Extract by Mask、Clip栅格版都需要一个范围要素。没有shp你只能按经纬度手动框选范围又慢又不精确。第三shp本身还携带属性信息比如区县名、行政区代码、面积等。这些字段可以作为分析结果的标签也可以用来做分级统计、分区统计Zonal Statistics比如“按区县统计平均海拔”“按乡镇统计最大坡度”这种分析没有shp根本做不了。所以你别小看这个附带shp文件它是整个数据包的“分析框架”DEM只是在这个框架里填地形属性。2. 从zip到可用数据的解锁过程2.1 解压Zip时最容易翻车的几个细节拿到一个zip压缩包第一反应肯定是右键解压。但GIS数据包的解压有一个不少人踩过坑的地方——文件夹路径。尽量把解压路径放在纯英文、无空格的目录下比如D:\data\yantai_dem\不要放在C:\Users\张三\桌面\新建文件夹 (2)\这种路径里。原因是ArcGIS、QGIS这类软件在读取数据时如果路径包含中文、空格、特殊字符偶尔会抽风。这不是玄学是底层文件读取库对非ASCII路径的支持不完善。尤其是老版本的ArcMap,对中文路径的容忍度非常差经常报The supported source data formats are...之类让人摸不着头脑的错误。另外一个常被忽略的点解压前检查zip包是否完整。资料在网盘间转存、QQ微信传输时经常损坏最典型的表现就是解压到某个进度时突然报错“文件末端缺少结束符”或者直接提示file is not a zip file。如果你的zip是分卷压缩.z01、.z02加.zip需要把所有分卷放在同一目录下再解压缺一个分卷就解不出来。我自己习惯的做法是解压之前先看一眼zip的大小和压缩包内文件列表心里有数解压后第一时间检查每个文件的大小是不是0字节0字节的基本就是坏了。如果发现损坏优先尝试修复工具比如Linux系的zip -FF damaged.zip --out repaired.zip命令可以在一定程度上重建zip索引但前提是压缩包本体没有大面积损坏。如果是分卷压缩千万别单独解压某一个分卷那只会得到一堆乱码。2.2 shp文件家族为什么它不是一个文件解压完之后你大概率会看到一坨后缀不同的同名文件yantai_city.shp、yantai_city.shx、yantai_city.dbf、yantai_city.prj、yantai_city.cpg可能还有.sbn、.sbx、.qpj。这时候很多人会困惑我要“打开shp文件”到底双击哪个答案是在GIS软件里你只需要选择那个.shp文件即可软件会自动读取同目录下的其他附属文件。但如果你把.shp单独拷走、把其他文件留在原目录那么打开后要么只有图形没有属性要么干脆打不开。所以在拷贝或移动shp文件时忠告就是永远成组拷贝别单独拿.shp。这里分享一个我自己的管理习惯我经常把shp的全部配套文件打成一个zip再传输或者干脆用文件夹整体复制。如果某些软件只让你上传单个文件那我建议把shp压成zip上传因为很多在线GIS平台支持zip自动解包读取shp。还有一个容易忽略的文件是.prj。它记录了坐标系信息。没有它虽然shp能打开但软件会默认它是未知坐标系在叠加其他数据时会错位得一塌糊涂。所以我打开shp后的第一个动作就是看它的坐标系这个习惯帮我至少省掉了十次麻烦。2.3 用什么工具打开这一堆文件打开shp和DEM的工具有很多我按使用频次排一下ArcGIS Pro / ArcMap国内用得最多的商业GIS软件功能全教程多但授权成本高。QGIS免费开源GIS软件完全够用我最近几年日常分析基本在QGIS里完成ArcGIS更多用来出图。Global Mapper轻量级GIS查看转换工具打开DEM、shp都极快做简单的裁剪、格式转换非常方便。纯Pythongeopandas rasterio适合批量处理和数据流水线对写代码的人最友好。如果你是第一次接触这类数据我建议直接从QGIS开始因为它免费、跨平台、插件生态丰富而且对中文编码的支持比ArcGIS还要好一点。3. 核心实操把DEM和shp组合出可用成果3.1 第一步检查坐标系统一投影打开shp和DEM之后第一件事不是急着裁剪而是先检查坐标系。这一步太重要了我单独拎出来说。DEM这类栅格数据如果是公开来源SRTM、Copernicus等原始坐标系一般是WGS84地理坐标系EPSG:4326单位是度。但烟台市的shp文件尤其是国内制作的基础地理数据投影往往是CGCS2000 / Gauss-Kruger投影比如EPSG:4548是CGCS2000 / 3-degree Gauss-Kruger CM 120E烟台在东经120度附近很可能落到这个带或者WGS 84 / UTM zone 51NEPSG:32651。这两个坐标系不统一的话后果很直观你把shp叠加到DEM上边界和地形完全对不上或者裁剪出来的结果是一片空白。所以你先分别查看两个数据的属性看坐标系是什么、单位是什么。如果坐标系不一致处理方案不是“把DEM转换成shp的坐标系”就是反过来。但有个原则分析时尽量让所有数据处于同一个投影坐标系。因为地理坐标系经纬度虽然在定位上没问题但在计算面积、距离、坡度这些量时会导致数值不准因为一个经度的长度在不同纬度下不一样。而投影坐标系是平面直角坐标单位是米计算距离和面积是天生的优势。所以推荐流程是先把DEM投影到和shp一致的投影坐标系如果shp的投影适合该区域或者把shp重投影到DEM的坐标系如果DEM本身已经是投影坐标。对于烟台这种中纬度地区我一般倾向用CGCS2000高斯克吕格投影或UTM具体用哪个取决于shp的已有坐标系尽量少做一次转换就少做一次避免精度损失。3.2 第二步按shp范围裁剪DEM坐标系统一好后正式进入第一个实际任务按烟台市范围把30m DEM裁剪出来。在ArcGIS Pro里可以用“Extract by Mask”工具输入栅格是原始DEM输入掩膜数据是shp输出就是裁剪好的DEM。在QGIS里可以用“裁剪栅格以掩膜图层”工具Raster clip by mask layer。原理上都是利用shp的边界几何将落在边界外的栅格像元设为NoData边界内的保留。裁剪的时候有几个参数要注意输出像元大小保持和输入DEM一致30m不要随意改大改小。改大了损失精度改小了纯属占存储空间。NoData值不同来源的DEMNoData值可能不一样有的是-9999有的是0有的是255。裁剪后如果发现边缘出现异常数值记得检查NoData值设置。是否保留外部区域一般裁剪掉但如果后边要做洪水淹没这种需要外部缓冲区的分析就不要裁太紧留一点缓冲更好。裁剪完成后建议你顺手做个“栅格统计”看看最小值、最大值、平均值是否合理。烟台的地形特征是南高北低南部有昆嵛山、牙山等山地最高海拔大约900多米北部是滨海平原最低接近海平面0米。所以统计结果如果最大值是惊人的3000米那基本可以判断数据源有问题或者投影转换出错。3.3 第三步从DEM提取你想分析的地形因子裁剪好DEM之后接下来能做的事就多了。我列几个最常用的都是ArcGIS/QGIS直接能跑的坡度Slope和坡向Aspect坡度就是地表倾斜程度单位是度或百分比。坡向是坡面朝向哪个方向9个方向平坦、北、东北、东、东南、南、西南、西、西北。这俩是DEM最基础、最常用的派生数据。搞工程建设要算坡度是否适合开发搞农业要知道哪些坡向光照好搞林业要判断坡向对植被分布的影响。计算原理不复杂对每个像元用周围3×3窗口的8个邻域像元计算最大变化率。ArcGIS里直接用Slope工具和Aspect工具QGIS里用“坡度”“坡向”工具。重点提示如果原始DEM是地理坐标系经纬度直接算坡度得到的结果是错的单位混乱必须先投影成米制坐标系再算。等高线ContourAutoCAD用户比较需要这个。从DEM提取等高线ArcGIS里用Contour工具QGIS里用“等高线”工具。生成时注意设置合适的等高距地形平缓区域10米就够密了山地丘陵区域可能得20米或50米。生成后可以导出为DXF格式给CAD使用。这里有个实操经验直接生成的等高线会比较毛糙尤其是30m分辨率的DEM提取出来的等高线经常带锯齿。建议生成前先对DEM做一次低通滤波比如3×3均值滤波或者重采样到较大像元虽然细节少一点但线条干净很多出图效果好。山体阴影Hillshade这个主要是做可视化底图用的。ArcGIS里是Hillshade工具需要设置太阳方位角默认315度西北方向、太阳高度角默认45度。生成后叠加在DEM上能做出立体感很强的地形图。实际上我更常用的做法是把山体阴影设置成半透明叠加在带颜色的DEM渲染层上这样既有地形高度色彩又有立体阴影效果可比单独一个灰度好看多了。水文分析填洼、流向、累积流量、河网提取这个属于进阶用法。DEM在处理时会有一些虚假的洼地spurious sinks如果不填掉水流方向会卡在洼地里导致河网提取断裂。所以用ArcGIS的Fill工具先填洼然后流向Flow Direction、流量累积Flow Accumulation、栅格计算器设阈值提取河网。30m分辨率提取的河网用于区域级的流域划分是可以的但用于精细的河道分析精度不够。我之前用烟台DEM做过一次流域提取的试验——不填洼直接提取河网在中段断了好几处填洼之后主河道完整多了。所以做水文分析填洼那一步千万别跳过。三维可视化DEM最直观的用途是三维显示。ArcGIS Pro的3D场景里直接把DEM作为高程表面叠加影像或shp就能看到烟台市的立体地形。还有两个方向很热门一是把DEM导出为TIN、3D模型再转成3D Tiles发布到Web端让用户在浏览器里看三维地形二是用Blender等三维软件加载DEM高度图生成可渲染的地形mesh这种常用于游戏场景、虚拟仿真。这里要提一下热搜词里的“dsm生成dem”和“copernicus dem glo 30下载”。DSM数字表面模型和DEM的区别在于——DSM包含了地物高度树、房子DEM是纯地表高程。如果只有DSM想做DEM需要做滤波或分类把非地面点去掉操作比较复杂。而Copernicus DEM GLO-30是欧洲航天局发布的全球30米DEM目前在开源数据里精度和时效性都很能打比老一代SRTM可靠做烟台这种30m DEM数据时Copernicus GLO-30本身就是很好的替代数据源。3.4 第四步shp文件的其他玩法说到shp不能只讨论“按此裁剪”。这个烟台市级shp还有其他用途转KML在Google Earth里看边界或者在航拍软件里叠加显示用QGIS右键图层导出即可。转GeoJSON前端Leaflet/Mapbox做Web地图时最常用的格式。QGIS里另存为GeoJSON几秒钟搞定。转3D Tiles把shp的边界挤出高度比如按属性字段拉伸生成三维建筑物或区域模型然后转成3D Tiles在Cesium里展示。这个流程涉及到shp转3dtiles时下三维GIS很热这个方向可以深入研究。渔网分割做方格网统计时把烟台市范围按指定格网大小打散比如5km×5km的渔网再统计每个格网内的平均海拔、最大坡度等。ArcGIS的Create Fishnet工具或QGIS的“创建网格”都可以。按需拆分如果想按区县分别出图可以利用“按属性分割”功能ArcGIS Pro的Split by AttributesQGIS的Split vector layer把市级合并图斑按区县字段拆分出多个shp。这里额外提一个很多新手问的问题“shp文件导出的时候没有cpg文件是怎么回事”cpg文件记录的是dbf属性表的字符编码如果缺了最常见的结果是打开shp后中文属性值显示成乱码。解决办法很简单缺失cpg时手动指定编码方式UTF-8或GBK重新读取。如果手头有多个shp都缺cpg可以批量写一个空的.cpg文件内容写UTF-8或GBK保存即可。4. 常见问题与排查技巧实录我挑了这些年关于这类数据出现频率最高的几个问题整理成速查表每条都是实操验证过的。问题现象可能原因排查手段和解决办法解压时报“file is not a zip file”zip文件损坏、下载不完整、或文件被改后缀重新下载用zip -FF尝试修复检查文件头是否为PK开头解压分卷.zip时失败分卷不全或分卷名被改确认所有.z01/.z02/.zip在同一目录且命名一致按顺序解压第一分卷ArcGIS/QGIS无法识别shp缺少.shx或.dbf或路径含中文把整套shp组件放同一目录迁移到纯英文路径检查.shp文件大小是否异常打开shp后中文属性乱码缺少.cpg或编码与软件默认不一致添加.cpg文件并指定GBK/UTF-8在QGIS中重设编码DEM在ArcGIS里全黑或全白显示拉伸方式不对或数据范围未识别右键图层属性符号系统选择“拉伸”“百分比裁剪”显示检查栅格统计看是否有极端值DEM和shp位置对不上坐标系不一致分别查看两数据的坐标系统一重投影到同一投影坐标裁剪后DEM一片空白掩膜范围与栅格范围相交很少或坐标系不一致先投影统一检查shp边界是否超出DEM有效范围检查NoData设置直接算坡度数值离谱输入是地理坐标系经纬度先用投影工具转成米制投影坐标系再计算坡度、面积、距离导出DXF给CAD后线条太密或太毛糙等高距设置不合理或DEM未滤波调大等高距对DEM做平滑滤波后再提取等高线打开zip里某文件提示“could not find EOCD”分卷损坏或zip结构不完整用zip -FF重建索引或用7-Zip尝试从损坏压缩包中提取文件这个表里的问题我几乎每一条都真实踩过。最惨的一次是接了一个外地项目对方发来一个分卷zip命名还改过我对着文件大小比对了一个多小时才拼凑出完整数据。所以在这里再提醒一次数据包传输时最好保留原始文件名分卷不要改名下载完先校验大小和完整性再解压。还有一个高频痛点是“file is not a zip file”。很多人第一反应是下载坏了但很多时候其实是文件被第三方软件自动改了后缀比如下载工具把.vnd格式修改了或者微信/QQ传输时文件被加上了类似于.zip(1)的后缀。你只要把后缀改回.zip再解压往往就能解开。这个坑大概有三分之一的人不知道。5. 从烟台DEM数据包延伸到更大的数据生态5.1 开源DEM数据源怎么选既然聊到30m DEM顺便把主流开源DEM数据源做个对比这对以后自己找数据会很有用数据源分辨率覆盖范围时间特点SRTM30m1弧秒/90m全球北纬60°~南纬56°2000年采集最经典但部分地区有空洞和噪点ASTER GDEM30m全球1999-2011年覆盖全球山区细节一般Copernicus DEM GLO-3030m全球2011-2015年基于TanDEM-X精度高推荐使用ALOS PALSAR RTC12.5m全球2006-2011年免费里分辨率最高的但数据量大NASADEM30m全球2000年SRTM的修复版做了很多修正我的建议是新项目优先Copernicus GLO-30它填补了很多SRTM的坑比如山区空洞同时高程精度比ASTER更可靠。如果对细节要求高可以考虑ALOS 12.5m不过处理起来机器负担会大不少。国内的DEM有时还可以从地理空间数据云、国家基础地理信息中心下载数据源更贴合国内坐标系。5.2 DEM数据在不同领域能干什么拿到烟台市的DEM后能做的远不止“看地形”这么简单。我给你列几个真实场景城市规划坡度大于25%的区域通常不适合大规模建设利用DEM提取坡度并叠加用地规划图可以快速识别不适合开发的土地为“三区三线”划定提供辅助决策。水利与洪涝用填洼后的DEM做流向分析和汇流累积量可以生成水系为洪水风险评估提供地形基础。30m分辨率在区域级洪水淹没分析中完全够用。交通选线道路纵断面设计需要精确的高程剖面DEM可以直接在任意两点之间生成剖面线辅助选线。通信基站选址利用DEM做视距分析Viewshed判断山顶或高层建筑顶部的信号覆盖范围。农业区划坡度、坡向、海拔是农作物适宜性评价最基础的三要素。比如烟台是中国重要的苹果产区苹果适合种植在海拔100-800米、坡度15度以下的丘陵地用DEM数据做适宜性分析分分钟把合适区域筛出来。生态评价通过DEM提取的地形湿润指数TWI常用来识别潜在湿地和水分汇集区这在生态红线划定中经常用到。5.3 数据落地展示从DEM到三维Web场景最后提一个现在很常见的需求方向把DEM和shp搬到Web端展示。这个流程涉及“shp转3dtiles”和“DEM生成地形瓦片”两类技术。简单说下思路先用GIS工具把DEM重采样为适合Web显示的层级比如Cesium terrain format或者生成quantized-mesh地形瓦片shp则先生成三维要素带拉伸高度再通过工具转成3D Tiles。前端用CesiumJS加载这两个数据就能在一个网页里看到烟台市的三维地形和行政边界。这个方向的详细实操细节比较多这里不展开但可以明确告诉大家从我们这套数据出发是完全可以落地到这个程度的。6. 一些实用的个人建议写到最后分享几个我在实操中养成的习惯不算全面但真的能提升效率和少踩坑第一所有原始数据zip原包务必单独归档不要解压后就把压缩包删了。数据出了问题要重新解压时原包就是你最后的保底。第二建立自己的目录规范。我的习惯是project/ 01_raw/ # 初始数据不做任何修改 02_intermediate/ # 中间结果如裁剪后的DEM、投影后的shp 03_final/ # 最终成果如坡度图、等高线、分析结果 04_scripts/ # 处理脚本Python、ModelBuilder 05_output/ # 图表、报告、配图这样看起来多此一举但在项目周期长、迭代多次时能避免“最终版”“最终版2”“最终版最终版”这种灾难现场。第三每次处理完数据随手记下坐标系、投影参数、分辨率、NoData值、处理时间。这个信息密度不一定很高但在复盘或出问题时极其有用。第四涉及到坐标转换时不要盲目用默认参数。国内CGCS2000与其他坐标系之间的转换往往需要七参数或四参数随便做一个近似转换结果可能会偏差好几十米。如果项目精度要求高最好从正规途径获取转换参数或者用专业的坐标转换工具处理。第五DEM数据的质量验证看似简单其实很容易忽略。拿到数据先跑一下统计检验再叠加影像目视检查几个关键地形位置基本能筛掉90%的明显错误。我个人的体会是像“山东省烟台市DEM数字高程数据30m含本市级范围shp文件.zip”这种数据包几乎每个GIS从业者都接触过但很多人只停留在“打开看一眼”的层面。真正把它的价值释放出来靠的是把“范围文件高程底图”这个组合用熟用透——从一张灰度图开始既能出坡度坡向、等高线、山体阴影这些基础产品也能延展到水文分析、三维可视化、选址评价这些复杂场景。这套能力是通用的拿到任何城市的数据包都能用同样的方法处理。希望这篇内容能给正在跟DEM和shp较劲的朋友们一些实际的帮助。本文还有配套的精品资源点击获取