烟台市30m DEM与shp数据处理实战：从解压到地形分析

简介：本资源为山东省烟台市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分辨率又是什么概念

DEM（Digital Elevation Model）就是数字高程模型，说白了是一张用数值记录地面高度的栅格图。每个像素存一个海拔值，单位一般是米。你把它想象成一个巨大的Excel表，每一格对应地面上一个小方块，格子里填的是这个小方块的平均海拔高度，投影后放到地图上，就成了一张能反映地形起伏的灰度图——高亮的地方高，暗的地方低。

30m分辨率的意思就是：地面上一个像素点对应30米×30米的实际面积。这个概念对后续所有分析都至关重要，因为它直接决定了你能做什么、不能做什么。举个例子，烟台市辖区面积大约是1.39万平方公里，30米分辨率的话，每个像元是900平方米，一景数据大约有1500多万个像元。这个精度对区域性地形分析完全够用，但你要是想分析一条10米宽的冲沟，那就不行了——它根本显示不出来。

那为什么30m这么流行？原因很实际：它是SRTM和ASTER GDEM这两大全球公开数据源的常见分辨率，后来Copernicus DEM（GLO-30）也是30米。很多省市级基础地理信息平台对外发布的DEM也是30m。再往上是12.5m（ALOS PALSAR）、10m（Copernicus 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都极快，做简单的裁剪、格式转换非常方便。
• 纯Python（geopandas + 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 51N（EPSG: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值：不同来源的DEM，NoData值可能不一样，有的是-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是地理坐标系（经纬度），直接算坡度得到的结果是错的（单位混乱），必须先投影成米制坐标系再算。

等高线（Contour）
AutoCAD用户比较需要这个。从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 Attributes；QGIS的Split vector layer），把市级合并图斑按区县字段拆分出多个shp。

这里额外提一个很多新手问的问题：“shp文件导出的时候没有cpg文件是怎么回事？”cpg文件记录的是dbf属性表的字符编码，如果缺了，最常见的结果是打开shp后，中文属性值显示成乱码。解决办法很简单：缺失cpg时，手动指定编码方式（UTF-8或GBK）重新读取。如果手头有多个shp都缺cpg，可以批量写一个空的.cpg文件，内容写UTF-8或GBK保存即可。

4. 常见问题与排查技巧实录

我挑了这些年关于这类数据出现频率最高的几个问题，整理成速查表，每条都是实操验证过的。

这个表里的问题，我几乎每一条都真实踩过。最惨的一次是接了一个外地项目，对方发来一个分卷zip，命名还改过，我对着文件大小比对了一个多小时才拼凑出完整数据。所以在这里再提醒一次：数据包传输时最好保留原始文件名，分卷不要改名，下载完先校验大小和完整性再解压。

还有一个高频痛点是“file is not a zip file”。很多人第一反应是下载坏了，但很多时候其实是文件被第三方软件自动改了后缀（比如下载工具把.vnd格式修改了），或者微信/QQ传输时，文件被加上了类似于.zip(1)的后缀。你只要把后缀改回.zip再解压，往往就能解开。这个坑大概有三分之一的人不知道。

5. 从烟台DEM数据包延伸到更大的数据生态

5.1 开源DEM数据源怎么选

既然聊到30m DEM，顺便把主流开源DEM数据源做个对比，这对以后自己找数据会很有用：

我的建议是：新项目优先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较劲的朋友们一些实际的帮助。

本文还有配套的精品资源，点击获取