【gis系列】从等高线到地形分析：dem生成与高程、坡度、坡向解析
1. 从等高线到DEM：数据准备与处理
搞地形分析的朋友们都知道，DEM（数字高程模型）是我们的基础粮草。但很多人卡在第一步——怎么把原始等高线数据变成可用的DEM？我当年第一次做这个的时候，踩了不少坑，现在把经验都分享给你们。
先说数据来源。常见的有几种方式：专业测绘数据、公开DEM数据（比如ASTER GDEM）、或者从地图软件导出。我常用的是QGIS搭配OpenStreetMap数据，免费又方便。拿到数据后，第一步是检查数据质量：
- 等高线完整性：有没有断线？等高距是否一致？
- 高程值准确性：检查Z值是否合理，避免出现平地突然冒出一座山的情况
- 拓扑错误：比如等高线交叉、重叠等
处理CAD数据时有个小技巧：先把所有高程点和高程注记都放到同一个图层。用AutoCAD的"LAYMRG"命令就能快速合并图层。记得删除无关元素，只保留等高线和高程点，这样导入GIS时不会带一堆垃圾数据。
# 用PyQGIS检查等高线数据的示例代码 layer = iface.activeLayer() for feature in layer.getFeatures(): geom = feature.geometry() if geom.isMultipart(): print("发现多段线！可能需要分解") z_values = [vertex.z() for vertex in geom.vertices()] if min(z_values) == max(z_values): print(f"ID {feature.id()} 的高程值异常")2. DEM生成实战：TIN与栅格的转换艺术
有了干净的等高线数据，接下来就是重头戏——生成DEM。这里推荐先用TIN（不规则三角网）作为中间产物，再转成栅格DEM。为什么？因为TIN能更好地保留地形特征，特别是在陡峭区域。
在ArcGIS Pro里操作很简单：
- 打开"3D Analyst工具"
- 选择"TIN管理→创建TIN"
- 设置好输出路径和坐标系统
- 关键参数：高程字段选"Elevation"，类型选"硬替换"
等处理完成后，你会得到一个看起来像三角网拼成的曲面。这时候右键图层选"符号系统"，改用高程渐变色显示，就能直观看到地形起伏了。
# GDAL创建DEM的命令行方案（适合批量处理） gdal_contour -a elevation -i 10.0 input_dem.tif output_contour.shp gdal_grid -zfield "elevation" -outsize 1000 1000 -ot Float32 contour.shp dem.tif转换栅格时要注意三个参数：
- 像元大小：一般用10-30米，精度要求高可以到5米
- 插值方法：自然邻域法适合复杂地形，反距离权重法运算更快
- 输出格式：GeoTIFF兼容性最好，IMG格式支持更大文件
3. 高程分析：不只是看高度那么简单
拿到DEM后，第一件事当然是看高程分布。但别急着出图，先做数据检查：
- 用栅格计算器排除异常值：
Con(("dem" > -50) & ("dem" < 9000), "dem") - 统计基本参数：最小值、最大值、平均值、标准差
- 生成高程分级图：建议用自然断点法(Jenks)分5-7类
进阶技巧来了——高程剖面分析。想规划一条徒步路线？用"剖面图"工具沿着预设路径提取高程变化。我做过一个项目，用这个方法发现了地图上没标注的悬崖，帮客户避免了重大安全隐患。
# 使用rasterio进行高程分析示例 import rasterio import numpy as np with rasterio.open('dem.tif') as src: elevation = src.read(1) print(f"最高点：{np.nanmax(elevation)}米") print(f"最低点：{np.nanmin(elevation)}米") # 生成高程分级 bins = [0, 200, 500, 1000, 2000, np.inf] classified = np.digitize(elevation, bins)4. 坡度分析：安全评估的关键指标
坡度图是工程规划中最重要的衍生品之一。在ArcGIS里生成坡度图很简单，但有几个细节要注意：
- Z因子：如果水平单位和垂直单位不同（比如经纬度vs米），需要设置Z因子转换
- 坡度表示法：百分比坡度(0-100%)还是度数(0-90°)? 工程常用百分比，科研常用度数
- 重分类：根据应用场景划分坡度等级。比如：
- 0-8%：适合建设
- 8-15%：需简单处理
- 15-25%：需工程措施
25%：不建议开发
有个实际案例：某开发商想在山区建度假村，我们做的坡度分析显示75%区域坡度超过25%，最后建议他们调整了选址，省下了大笔地基处理费用。
5. 坡向分析：太阳光照与生态研究利器
坡向图看起来就是个彩色圆盘，但信息量很大。正北是0°，顺时针增加到360°。在QGIS里生成时要注意：
- 使用"坡向"工具时勾选"三角函数法"，结果更精确
- 输出建议用16位整型，节省存储空间
- 可视化时用"圆形色带"样式，直观显示各方向
坡向分析最实用的三个应用：
- 太阳能板布置：南坡（北半球）接收更多日照
- 农业规划：不同坡向的作物选择
- 生态研究：分析植被分布与坡向的关系
有个小技巧：用栅格计算器结合坡向和高程，可以生成"日照时长模拟图"。公式类似：("aspect" > 135 & "aspect" < 225) & ("elevation" > 500)
6. 常见问题排查与性能优化
做了上百个项目后，我整理出DEM生成的七大典型问题：
- 锯齿状等高线：通常是像元大小设置过大，试试缩小2-3倍
- 平坦区域出现噪点：检查原始数据是否有零散高程点
- 边缘缺失：处理前先做缓冲区分析
- 处理速度慢：分块处理，或者改用TIN后再转换
- 坐标系统错误：始终检查垂直基准面是否一致
- 内存不足：改用64位软件，或调整处理范围
- 结果异常：先用小范围测试参数
性能优化方面，我的经验是：
- 超过1GB的数据建议用PostGIS+GDAL方案
- 定期清理临时文件
- 关闭不必要的图层和符号渲染
- 对于超大数据，考虑使用LAS点云数据替代
最后提醒：所有分析结果都要实地验证。有次我们发现某区域坡度分析结果异常，到现场才发现是新建的高架桥影响了DEM数据。地形分析从来都不是纯数字游戏，实地经验同样重要。