DEM地形操作(geotools方式与nga方式)

Posted 2023-05-02

参考技术A 由于项目中需要用到大范围tiff的图像（全中国30米分辨率的dem影像），并且需要单点获取高程，以及实现部分范围的dem裁切与获取趋于范围极值，当时在网上查找的部分，很多都不满足预期，或者计算结果与实际并不够契合，因此单独开一篇专门讲这块内容。

引入基础的依赖包，使用中间组件完成所需功能。

由于组件之间存在交互关系，因此将其全部放到一起。

由于接下去的操作都需要依赖同一份数据源，因此在最开始的时候，需要先封装关于数据源的相关操作。

通过获取gridCoverage2D 的的投影信息，加载投影并使用经纬度进行计算位置点高程

总体来说，整个的实现流程基本上可以概括为[引入依赖]->[加载数据源]->[业务数据读取]。读取DEM的高程流程相对简单，获取范围内的高程极值操作也不复杂，主要是需要算出范围对应的格网，以及获取每个格网像元的高程值；较为复杂的是DEM的图像导出，由于使用geotools导出的dem图像存在一定的偏移与投影问题，导致其他平台软件无法正常读取dem的一些属性信息，因此在操作上，先使用geotools自定义的导出，先输出一份DEM，再使用正常平台导出的dem文件，通过读取它的属性信息，赋值给空的dem影像，再将像元的高程值写入，最后输出成一份符合规范的dem。此操作较为繁琐，但基础满足功能所需，此为抛砖引玉，如有更好的办法实现，望诸位不吝赐教，将感激不尽！

QGIS入门实战精品教程10.1：QGIS基于DEM数据的地形分析案例教程

本文讲解QGIS中基于DEM数据的地形分析方法，包括：坡度分析、坡向分析、山体阴影、地貌分析、强度指数(地形复杂性)。

文章目录

• 一、加载DEM
• 二、坡度分析
• 三、坡向分析
• 四、山体阴影
• 五、地貌分析
• 六、强度指数(地形复杂性)

---

QGIS提供了地形分析的工具，位于工具箱→Raster terrain analysis，如图所示：

一、加载DEM

二、坡度分析

坡度指过地面人任意一点的切平面与水平地面之间的夹角。

计算结果：