QGIS如何制作等高线

网优随笔

于 2025-03-13 16:40:56 发布

阅读量1.1k

点赞数 31

文章标签： QGIS 开源软件

本文链接：https://blog.youkuaiyun.com/qq_27987847/article/details/146235376

版权

一、等高线基础理论与应用价值

1.1 等高线测绘原理

等高线作为三维地形二维表达的核心要素，本质是海拔相同的相邻点投影连线。其数学基础可表示为：

z=f(x,y)=C(C为高程常数)z=f(x,y)=C(C为高程常数)

通过间距Δh控制地形表现精度，Δh越小地形细节越丰富。现代测绘中，计算机辅助生成已取代传统手工绘制，DEM数据成为主要来源。

1.2 应用场景扩展

工程勘测：道路选线时，通过等高距分析确定最大纵坡
灾害预警：滑坡体识别需结合等高线曲率变化
军事应用：三维地形重构支持战场环境模拟
生态研究：微地形分析揭示物种分布规律

二、DEM数据获取与预处理

2.1 多源数据获取方案

数据源	分辨率	获取方式	适用场景
OpenTopography插件	30m	QGIS插件在线下载	大区域概略分析
ALOS World 3D	30m	JAXA官方申请	全球地形建模
国家基础地理信息数据库	5m/10m	测绘部门授权使用	工程级精度要求
无人机航测	0.1-1m	Pix4D/Photoscan生成	小范围精细建模

通过QGIS插件OpenTopography DEM Downloader获取数据时需注意：

注册获取API密钥，建议创建多账户应对下载限额
坐标系强制转换工具推荐使用GDAL Warp
大数据量处理时启用PyQGIS批处理脚本

2.2 数据质量控制流程

高程突变

空洞修补

原始DEM

异常值检测

形态学滤波

IDW插值

坡度一致性检查

重采样优化

最终DEM

三、QGIS等高线生成技术体系

3.1 GRASS工具链工作流

模块调用：通过Processing Toolbox启动r.contour
参数设置：
- 基准高程(Offset)：设定起始等高线值
- 等高距(Step)：根据地貌复杂度动态调整
- 简化阈值(Simplify)：控制顶点密度
并行计算：启用--overwrite参数避免内存溢出

3.2 GDAL等高线生成对比

特性	GRASS r.contour	GDAL Contour
算法核心	改进型Marching Squares	标准Marching Cubes
拓扑检查	内置拓扑验证	需后期处理
大数据支持	分块处理	单线程处理
属性字段	自动生成坡度属性	仅高程值
执行速度(1GB DEM)	2分18秒	3分45秒

实验表明，复杂地形宜选用GRASS工具，简单地形可优先GDAL。

四、参数优化与质量控制

4.1 等高距动态设置模型

Δh={0.5m当 S≤5∘1.0m当 5∘<S≤15∘2.0m当 S>15∘Δh=⎩⎨⎧0.5m1.0m2.0m当 S≤5∘当 5∘<S≤15∘当 S>15∘

其中S为区域平均坡度，通过r.slope计算获得。

4.2 断线修复技术

形态学修复：使用v.clean工具消除微小间隙
流线追踪：基于D8算法重建山谷线
人工编辑：QGIS拓扑检查器定位拓扑错误

五、三维建模与高级分析

5.1 三维可视化实现

# QGIS Python控制台脚本示例
layer = iface.activeLayer()
extent = layer.extent()
rlayer = QgsRasterLayer("dem.tif", "DEM")
QgsProject.instance().addMapLayer(rlayer)

params = {
    'MAP_CONFIG': {
        'coord': 'EPSG:4326',
        'layers': [rlayer],
        'resolution': 1,
    },
    'CAMERA_CONFIG': {'altitude': 5000}
}
controller = Qgs3DMapCanvasWidget(params)
controller.show()

运行

该脚本实现DEM数据的三维即时渲染。

5.2 地形因子关联分析

地形指标	计算方法	应用领域
地表曲率	`r.param.scale`	水土流失评估
地形湿度指数	`r.topidx`	生态分区规划
可视域分析	`r.viewshed`	通信基站选址
侵蚀强度	`r.uslek`	地质灾害预警

六、常见问题诊断矩阵

故障现象	可能原因	解决方案
等高线交叉	DEM数据噪声	使用`r.neighbors`中值滤波
高程值跳跃	坐标系不匹配	检查CRS一致性并重投影
边缘锯齿明显	重采样方法不当	改用双三次卷积插值
属性字段丢失	输出格式兼容性问题	指定ESRI Shapefile格式输出
运行内存不足	未启用分块处理	设置GRASS_REGION环境变量