无人机航摄生成DEM的高程点快速提取算法

1、无人机航摄生成DEM的高程点快速提取算法A fast extraction algorithm for surface peaks of DEM derived from UAV aerial photography zhang suAbstract： UAV aerial photography technology has been rapidly developed in recent years. UAV has the advantag- es of convenient take-off and landing, high resolution and fast response

2、and plays an important role in disaster relief and geological disaster monitoring. The surface peaks extraction in the topographic map generated by aerial photogra- phy of UAV is an important step in topographic map production. This paper proposes an adaptive topographic map surface peaks extraction

3、 algorithm, which has higher accuracy and simple calculation than traditional algorithms. And this method is more practical. The basis of this algorithm is the image expansion algorithm in mathematical morphology. In this pape, the algorithm is applied to the topographic map generated by the UAV aer

4、ial photography in Danfeng County and Ziyang County of Shaanxi Province, and the surface peaks that meet the requirements are extracted. The results show that the extracted surface peaks fully meet the requirements of topographic maps, indicating the accuracy and practicability of the algorithm.key

5、word：UAV DEM surface peaks algorithm无人机航摄的应用与发展无人机航空遥感是近年来发展起来的一种遥感技术的新手段和 遥感技术研究的一个热点领域，具有低空云下飞行、起降方便、分辨 率高的特点，是卫星遥感和有人机航空遥感的重要补充，现在无人机 摄影已经广泛用于地质灾害抢险救灾等领域。无人机航空数据相对卫 星数据来说，受天气的影响要小，数据快速获取的效率更高。例如， 在2021年10月发生的西藏雅鲁藏布江沿岸山体滑坡及四川与西藏 交界金沙江山体滑坡形成堰塞湖时，无人机都发挥了重要作用，在第 一时间取得了灾害区域的影像，为抢险救灾提供了重要的数据影像支 持，为救灾的决

6、策提供了重要的依据。而传统卫星由于重访期和云层 遮蔽的关系，有可能不能获得实时的卫星影像。无人机航空摄影除了可以取得航摄区域的高精度影像，还可以根 据拍摄的航片，通过数据处理软件例如PixGrid，Pix4d Mapper等软 件，得到航摄区域的DEM数据（数字高程模型）。如果有地面测量的 控制点资料，无人机生成的DEM数据精度可以达到米级。相对于传 统的地形图生产制作模式，无人机航摄数据制作的地形图具有精度高， 制作时间短等特点，非常适合小规模区域内的地形图快速制作。地形图高程点的提取的基本原理地形图制作中很重要的一个步骤是高程点的提取。传统的地形图 中，高程点即标有高程数值的信息点，通常与

7、等高线配合表达地貌特 征的高程信息。高程点分布表达上具有丰富的影像相特征，除个体 高程属性上的差异外，一方面与重要地形特征单元（如山峰、鞍部及 山谷等）、河流、道路等自然或人文要素目标间的分布关系，使得个 体高程点在地理性质上具有不同的重要性意义;另一方面区域地形的 起伏变化，使得高程点群在空间分布上呈现不同的疏密关系，在地形 起伏变化较大，即地形的变化高频地区，高程点的密度也较大。在地 形起伏变化较小，即地形的变化低频地区，高程点的密度也较小。高 程点信息和无人机的影像一起给地形图使用者提供了重要的地形信 息及决策判断依据。地形图上的高程点一般分为三类：（1）地形特征高程点。例如山 峰、山体

8、鞍部、洼地等地形特征，山峰是局部区域制高点，并且山峰 的视线最为广阔、鞍部处于两座山峰之间，是连接两座山峰的路径选 择的必经区域、洼地所处的地势是局部最低的，并且是水源存在的地 方，因此地形特征上的高程点提供的高程信息明确地指出了这些特殊 的地形。重要地理目标临近的高程点。重要地理目标如道路、居民 地、特殊建筑、河流等。这些地理目标是地图表达的重要对象，邻近 区域内的高程信息更为地形图使用者所重视，如道路两侧的高程信息 是公路建设选线需要考虑的一个重要因素。河流的高程点信息可以表 现河流的平缓程度，对于水路运输是否可行提供了重要的参考信息。 这一类高程点可以根据影像上的地理目标加以判断，并在关

9、注的地理 目标旁进行高程点的标注。分布范围上的关键高程点。这一类高程点对于高程点空间 分布的特征有重要的影像，属于特征的重要组成部分。一旦关键的高 程点被删除，整个高程点的空间分布会发生重大改变。这一类高程点 可以利用空间分布的权重进行分析并加以判断。本文的主要研究目的是研究第一类高程点的提取，传统的提取方 法包括由等高线生成栅格DEM模型，针对每个高程点设置以该点为 中心的RxR栅格窗口，然后分析中心栅格点与周围其他栅格点的高 程关系进行判断。Wood J.提出了从地形特征点的几何形态特征研究 了 DEM提取山顶点的算法，提出在局部区域内用x方向和y方向关 于高程z的二阶偏导数的正负性来判断

10、山顶点算法，从而提取出高程 点。刘淑琼等人提出了通过栅格数据的空间分析，窗口计算功能以及 叠置分析提取DEM的山顶点的算法。陈盼盼等人提出了 GIS栅格数 据邻域分析功能快速提取山顶点的算法。邹宝裕等人提出了基于封闭 等高线套合法提取山顶的算法。刘宏建等人以基于流域剖分的山顶点 自动提取方法为基础，详细阐述了流域剖分提取高程点模型的建立过 程。以上方法均需要在专业软件如Arcgis中进行空间分析及运算，步 骤较为复杂。本文提出了一种改进的自适应的高程点提取算法，可以 调整提取的高程点群的空间分布疏密程度，从而满足不同使用者对于 地形图的要求，并且运算的复杂度和传统的算法相比要简单，有更高 的计

11、算效率和实用性。研究方法及技术路线本文算法的思路是通过对无人机生成的DEM进行动态的、全局 的分析，并通过搜索局部极值的方式来寻找第一类高程点。DEM数 据可以认为是一组2维的矩阵，山峰在影像上的形态可以等同于DEM 矩阵上局部最大值。算法的核心思想是通过局部极大值间距的确定， 对2维矩阵进行图像膨胀运算，通过局部极大值的算法寻找局部极大 值并进行标注。该算法的优点在于，可以自定义极大值之间的距离， 从而人工的调节极大值疏密程度，满足不同的使用者对于地形图信息 的需求。当地形图使用者仅仅是需要整个区域内的地形大致分信息时， 可以将极大值的间距调稍大一些。如果地形图使用者对区域内的地形 图的地形

12、信息要求非常的详细准确时，可以极大值的间距调稍小一些， 生成密度较大的高程点群，满足使用者的要求。当寻找完DEM区域 内的高程点后，可以将高程点群的平面坐标和高程值进行输出，并将 高程点所在矩阵的行列数转换为大地平面直角坐标系，以便用于高程 点在地形图上的标注。对于寻找洼地等最低值，可以将DEM的高程 值取其名倒对值的负数，然后运行此算法，即可将寻找洼地的最低值转 换为寻找局部最大值的问题。整个程序在MATLAB软件中进行编译 和运行。寻找高程点的算法流程图如图1所示。3.1DEM数据从地理坐标转换为矩阵坐标无人机处理软件得到的DEM数据一般为.tif格式或者为.img格 式，通常带有DEM的

13、坐标信息，而坐标信息详细的描述了每一个点 所在的平面直角坐标系所处的位置，将DEM数据转换为二维矩阵， 导入到MATLAB软件之后，每一个DEM点对应矩阵的一个元素，需 要将DEM点的地理坐标转换为矩阵的行列数，地理坐标和矩阵坐标 之间的关系可以由公式(1)进行描述：x_image=y_matlab;y_image=a (1) -x_matlab;x_arcgis二R*x_image+x1;y_arcgis二R*y_image+y1; (1)其中，x_matlab，y_matlab，表示矩阵的行列数，而x_arc-gis， y_arcgis表示了 DEM点的地理坐标，R为DEM的空间分辨率，x

14、l, yl分别为DEM在地理坐标系中x轴最小值和y轴最小值x_image, y_image表示了 DEM点在一幅普通的影像上的坐标值。通过公式(1)， 即可将DEM的点映射到相应的二维矩阵上。3.2确定局部极大值之间的间距对于一副地形图，不同的使用者对高程点的疏密程度有不同的要 求。如果为了得到更加详尽的高程点信息，可以先将局部极大值的间 距调整为一个较小的值，本文使用的地形图为1：5000比例尺的地 形图，采用了 50m的极大值间距进行处理。如果发现生成的高程点 分布过密，可以加大极大值的间距，重新进行运算，得到相对稀疏的 高程群分布。在程序中，通过设定mindist变量的值进行极大值间距

15、定义，具体语句为：minDist=；因为DEM的分辨率为2m，所以实际的极大值间距为50m。3.3通过局部极大值搜索算法寻找高程点群本文采用的局部极大值的搜索算法为数学形态学的膨胀算法，膨 胀算法主要用途为消除噪声;分割出独立的图像元素，在图像中连接 相邻的元素;寻找图像中明显的极大值区域;其基本算法步骤如下：1定义一个卷积核B核可以是任何的形状和大小，且拥有一个单独定义出来的参考点 -锚点(anchorpoint);通常和为带参考点的正方形或者圆盘，可将核称为模板或掩膜；2将核B与图像A进行卷积，计算核B覆盖区域的像素点最大 值；3将这个最大值赋值给参考点指定的像素。通过局部极大值的算 法，

16、可以寻找到二维矩阵的局部极大值，其原理示意图如图2所示。从图2中可以看出，标注为黑色“*”的点即为局部极大值点，标注 为青色“* ”极为局部极小值的点，通过该算法就可以寻找山峰和洼地两 类高程点。3.4调整局部极大值的间距并且得到需要的高程点群根据使用者的需求，调整mindist变量的值，根据不同的局部极 大值的间距的值，得到不同疏密程度的高程点群的空间分布，直到该 分布满足使用者的预期要求为止，图3显示了不同的间距对于高程点 群的空间分布的影像，可以看出，间距为50m时，得到的高程点密 度明显要大于间距为100m时的高程点密度。3.5转换高程点群的坐标并输出高程点文件高程点在矩阵中的位置为行列数，应该将行列数转换为平面直角 坐标，才能与地形图的坐标系相匹配。本步骤也可以通过公式（1） 进行坐标系的转换。得到平面直角坐标之后，和高程点的高程值一起， 用xlswrite命令输出到Ex-cel文件中，将.xls文件导入到Arcgis软件 之中，便可以得到高程点文件。研究结果本文对陕西省丹凤县和紫阳县无人机航测生成的地形图进行了 高程点的提取，结果如图4所示，地形图使用的比例尺为1：*， 该地形图的主要用途的地质灾害遥感解译。图4分别展示了