特点，在抢险救灾，地质灾害监测等方面发挥了重要作用。无人机航空摄影生成的地形图中的高程点提取是地 形图制作的重要步骤，本文提出了一种自适应的地形图高程点提取算法，相对于传统的算法，具有精确度高，计 算简单，实用性强等特点。该算法的基础是数学形态学中的图像膨胀算法。本文将此算法用于陕西省丹凤县和 紫阳县的无人机航摄生成的地形图中，提取了符合要求的高程点。结果表明，提取的高程点完全满足了地形图
的要求，说明了该算法的准确性和实用性。关键词：无人机；地形图；高程点；算法A fast extraction algorithm for surface peaks of DEM derived from
UAV aerial photographyzhang suAbstract： UAV aerial photography technology has been rapidly developed in recent years. UAV has the advantages of convenient take—off and landing, high resolution and fast response, and plays an important role in disaster relief
and geological disaster mon让oring. The surface peaks extraction in the topographic map generated by aerial photography of UAV is an important step in topographic map production. This paper proposes an adaptive topographic map sur
face peaks extraction 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 paper, the algorithm is applied to the topographic map generated by the UAV aerial 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 fiilly meet the requirements of topographic maps, indicating the accuracy
and practicability of the algorithm.key word: UAV ,DEM ,surface peaks, algorithm1•无人机航摄的应用与发展无人机航空摄影除了可以取得航摄区域的高精度影像, 还可以根据拍摄的航片，通过数据处理软件例如PixGrid,
无人机航空遥感是近年来发展起来的一种遥感技术的 新手段和遥感技术研究的一个热点领域，具有低空云下飞
Pix4d Mapper等软件，得到航摄区域的DEM数据（数字高程 模型）。如果有地面测量的控制点资料,无人机生成的DEM
行、起降方便、分辨率高的特点，是卫星遥感和有人机航空 遥感的重要补充，现在无人机摄影已经广泛用于地质灾害 抢险救灾等领域叭无人机航空数据相对卫星数据来说，受
数据精度可以达到米级。相对于传统的地形图生产制作模 等特点，非常适合小规模区域内的地形图快速制作。2.地形图高程点的提取的基本原理式,无人机航摄数据制作的地形图具有精度高，制作时间短
天气的影响要小,数据快速获取的效率更高。例如，在2018 年10月发生的西藏雅鲁藏布江沿岸山体滑坡及四川与西 藏交界金沙江山体滑坡形成堰塞湖时，无人机都发挥了重
地形图制作中很重要的一个步骤是高程点的提取。传
要作用，在第一时间取得了灾害区域的影像,为抢险救灾提 统的地形图中.高程点即标有高程数值的信息点，通常与等
供了重要的数据影像支持.为救灾的决策提供了重要的依 据冋。而传统卫星由于重访期和云层遮蔽的关系，有可能不
高线配合表达地貌特征的高程信息。高程点分布表达上具
有丰富的影像相关特征,除个体高程属性上的差异外，一方
能获得实时的卫星影像。面与重要地形特征单元（如山峰、鞍部及山谷等）、河流、道路
作者简介：张肃（1982~）,男，四川省仪陇县人，学历：硕士研究生（清华大学土木工程系大地测量专业工学硕士），职称：工程师，研究方向:
遥感影像处理、土地利用分析、环境遥感。146/RESOURCES2()19年第二期 遥感测绘WESTERN RESOURCES等自然或人文要素目标间的分布关系,使得个体高程点在地 理性质上具有不同的重要性意义;另一方面区域地形的起伏
空间分析及运算，步骤较为复杂。本文提岀了一种改进的 自适应的高程点提取算法，可以调整提取的高程点群的空
变化,使得高程点群在空间分布上呈现不同的疏密关系，在 地形起伏变化较大，即地形的变化高频地区,高程点的密度
间分布疏密程度，从而满足不同使用者对于地形图的要求, 并且运算的复杂度和传统的算法相比要简单，有更高的计
也较大。在地形起伏变化较小，即地形的变化低频地区，高 程点的密度也较小。高程点信息和无人机的影像一起给地
算效率和实用性。3.研究方法及技术路线形图使用者提供了重要的地形信息及决策判断依据E。地形图上的高程点一般分为三类：(1) 地形特征高程点。例如山峰、山体鞍部、洼地等地 形特征，山峰是局部区域制高点，并且山峰的视线最为广
阔、鞍部处于两座山峰之间，是连接两座山峰的路径选择的 必经区域、洼地所处的地势是局部最低的，并且是水源存在
的地方，因此地形特征上的高程点提供的高程信息明确地
指岀了这些特殊的地形。(2) 重要地理目标临近的高程点。重要地理目标如道
路、居民地、特殊建筑、河流等。这些地理目标是地图表达 的重要对象.邻近区域内的高程信息更为地形图使用者所
重视，如道路两侧的高程信息是公路建设选线需要考虑的 一个重要因素。河流的高程点信息可以表现河流的平缓程
度，对于水路运输是否可行提供了重要的参考信息。这一 类高程点可以根据影像上的地理目标加以判断，并在关注 的地理目标旁进行高程点的标注。(3) 分布范围上的关键高程点。这一类高程点对于高
程点空间分布的特征有重要的影像，属于特征的重要组成 部分。一旦关键的高程点被删除，整个高程点的空间分布
图1高程点提取流程图本文算法的思路是通过对无人机生成的DEM进行动
态的、全局的分析，并通过搜索局部极值的方式来寻找第一 类高程点。DEM数据可以认为是一组2维的矩阵，山峰在
会发生重大改变。这一类高程点可以利用空间分布的权重
进行分析并加以判断。本文的主要研究目的是研究第一类高程点的提取，传 统的提取方法包括由等高线生成栅格DEM模型，针对每个
影像上的形态可以等同于DEM矩阵上局部最大值。算法
的核心思想是通过局部极大值间距的确定，对2维矩阵进 行图像膨胀运算，通过局部极大值的算法寻找局部极大值
高程点设置以该点为中心的RxR栅格窗口，然后分析中心
并进行标注。该算法的优点在于，可以自定义极大值之间 的距离，从而人工的调节极大值疏密程度，满足不同的使用
栅格点与周围其他栅格点的高程关系进行判断。Woo z的二阶偏导数的正负性来判断山顶点算法，从而提取岀高 程点。刘淑琼等人提出了通过栅格数据的空间分析，窗口 计算功能以及叠置分析提取DEM的山顶点的算法讥陈盼 盼等人提出了 GIS栅格数据邻域分析功能快速提取山顶点 的算法代邹宝裕等人提出了基于封闭等高线套合法提取 大一些。如果地形图使用者对区域内的地形图的地形信息 要求非常的详细准确时，可以极大值的间距调稍小一些，生 成密度较大的高程点群，满足使用者的要求。当寻找完 DEM区域内的高程点后，可以将高程点群的平面坐标和高 程值进行输岀，并将高程点所在矩阵的行列数转换为大地 山顶的算法讥 刘宏建等人以基于流域剖分的山顶点自动 提取方法为基础•详细阐述了流域剖分提取高程点模型的 平面直角坐标系，以便用于高程点在地形图上的标注。对 于寻找洼地等最低值，可以将DEM的高程值取其绝对值的 负数，然后运行此算法，即可将寻找洼地的最低值转换为寻 RESOURCES/147建立过程冏。以上方法均需要在专业软件如Aregis中进行 WESTERN RESOURCES遥感测绘2019年第二期找局部最大值的问题。整个程序在MATLAB软件中进行编 板或掩膜；点最大值；译和运行。寻找高程点的算法流程图如图1所示。② 将核B与图像A进行卷积,计算核B覆盖区域的像素 3.1 DEM教据从地理坐标转换为矩阵坐标无人机处理软件得到的DEM数据一般为.tif格式或者 为.img格式，通常带有DEM的坐标信息，而坐标信息详细的 ③ 将这个最大值赋值给参考点指定的像素。极大值，其原理示意图如图2所示。通过局部极大值的算法，可以寻找到二维矩阵的局部 描述了每一个点所在的平面直角坐标系所处的位置，将 DEM数据转换为二维矩阵，导入到MATLAB软件之后，每 一个DEM点对应矩阵的一个元素,需要将DEM点的地理坐 标转换为矩阵的行列数，地理坐标和矩阵坐标之间的关系 可以由公式（1）进行描述：x_image=y_matlab;y_image=a⑴一x_matlab;x_arcgis=R*x_image+xl;y_arcgis=R*y_image+y 1; （1）其中，x_matlab, y_matlah，表示矩阵的行列数，而x_arc- gis,y_arcgis表Z5了 DEM点的地理坐标，R为DEM的空间分 辨率,xl,yl分别为DEM在地理坐标系中x轴最小值和y轴 最小值。x_image, y_image表亦了 DEM点在一幅普通的影 像上的坐标值。通过公式（1）,即可将DEM的点映射到相 图2局部极大值算法结果示意图应的二维矩阵上。从图2中可以看出，标注为黑色“*”的点即为局部极大 可以寻找山峰和洼地两类高程点。值点，标注为青色“*”极为局部极小值的点，通过该算法就 3.2确定局部极大值之间的间距对于一副地形图，不同的使用者对高程点的疏密程度 有不同的要求。如果为了得到更加详尽的高程点信息，可 以先将局部极大值的间距调整为一个较小的值，本文使用 3.4调整局部圾大值的间距并且得到需要的高程点群根据使用者的需求，调整mindist变量的值，根据不同的 局部极大值的间距的值，得到不同疏密程度的高程点群的 的地形图为1：5000比例尺的地形图，采用了 50m的极大值 间距进行处理。如果发现生成的高程点分布过密,可以加 空间分布，直到该分布满足使用者的预期要求为止.图3显 示了不同的间距对于高程点群的空间分布的影像，可以看 大极大值的间距,重新进行运算，得到相对稀疏的高程群分 布。在程序中，通过设定mindist变量的值进行极大值间距 出，间距为50m时，得到的高程点密度明显要大于间距为 定义,具体语句为：minDist = [25 25];100m时的高程点密度。因为DEM的分辨率为2m,所以实际的极大值间距为 50m o3.3通过局部极大值捜素算法寻找高程点群本文采用的局部极大值的搜索算法为数学形态学的膨 胀算法，膨胀算法主要用途为消除噪声;分割出独立的图像 元素,在图像中连接相邻的元素;寻找图像中明显的极大值 区域;其基本算法步骤如下：问86为50m飾高程点群分布 问距为100m的高程点群令布图3不同间距的高程点群分布情况3.5转换嘉程奴群的坐标并输出高程点丈件高程点在矩阵中的位置为行列数，应该将行列数转换 为平面直角坐标，才能与地形图的坐标系相匹配。本步骤 定义一个卷积核①B核可以是任何的形状和大小，且拥有一个单独定义出 也可以通过公式（1）进行坐标系的转换。得到平面直角坐 标之后,和高程点的高程值一起，用xlswrite命令输出到Ex- 来的参考点-锚点（anchorpoint）;通常和为带参考点的正方形或者圆盘,可将核称为模 148/RESOURCEScel文件中，将.xls文件导入到Arcgis软件之（下转150页）WESTERN RESOURCES遥感测绘2()19年第二期参考文献：[1] 李瑜，王路.浅谈利用DGSinfo绘制剖面图的方法及技巧[J].西部 (a)资源，2015(06):60-61.[2] 中华人民共和国国家质量监督检验检疫总局.GB/T 13908 — 2002固体矿产地质勘查规范总则[S].北京：中国标准出版社, (b)2002.[3] 马延景，李戲，靳浩，朱进守，马文.基于Section软件的矿区实际 图2地质剖面图填充花纹前后的对比图回a.地质剖面图原始形态示意;b.地质剖面图填充花纹后材料图绘制H 现代矿业，2016, 32(09):199-202.[4] 彭桥梁，贺文华，刘瑞，郑时干.基于Section的地球化学剖面光谱 样晶折线图绘制叮现代矿业，2016, 32(03):148-151.4.结论近年基于Mapgis二次开发出来的软件如DGSS.Section 等，为总结实践野外工作提供了诸多便利快捷的功能，以数 据为核心使用这类软件可以轻松绘制各种图件，很大程度 上提高了效率，节约了时间叭[5] 李晓英，唐玉虎.Section在煤炭资源储量估算中的应用[J].科技 经济市场，2015(08):206.[6] 靳晓新.Section在地球化学找矿制图上的应用H现代工业经济 和信息化，2015, 5(01):67-69.[7] 刘永新.利用section软件制作钻孔柱状图[J].内蒙古煤炭经济, 固体矿产地质勘查规范总则中提到鼓励使用计算机辅 助野外采集系统.凡能用计算机成图、成表的资料，都应按 2013(05):45.[8] 徐鹏，方子樊，裴晓东.用Section软件快速绘制区块范围图口资 标准化表格内容的要求填写对需要矢量化图件叫本文通过 多种方式合理使用新的花纹，快速、便捷的使用自制花纹, 源环境与工程，2013, 27(05):706-707.〔9]廖欢，吴辉明.Section软件在绘制煤层底板等高线图中的应用 J],江西煤炭科技，2012(04):44-45.问丁作节省时间成本,大大地提高了工作效率。(上接148页)中，便可以得到高程点文件。4.研究结果步骤和计算过程都相对简单，减小了使用者的操作及学习负 担。该算法可以自由地调整高程点之间的间距.从而得到最 适合使用者需求的高程点群分布。本文通过对不同区域的 无人机航摄生成的地形图进行试验，提取了符合要求的高 本文对陕西省丹凤县和紫阳县无人机航测生成的地形 图进行了高程点的提取，结果如图4所示，地形图使用的比 例尺为1: 10000,该地形图的主要用途的地质灾害遥感解 译。图4分别展示了丹凤县和紫阳县的局部地形.从图中 程点群，说明了此算法的稳定性和实用性。表明了此算法 是日益发展的无人机航摄地形图制作过程的有效工具。 参考文献：[1] 支卫斌.无人机航空摄影测量在地形测绘中的应用出.江西建材, 2015(08):224-225.[2] 汤国安，刘学军，闾国年.数字高程模型及地学分析的原理与方 可以看出，该算法很好的提取了所需要的高程点，较好地完 成了高程点的标注任务。法[M].科学出版社，2005.[3] 金沙江西藏段河道水情基本恢复正常,中国政府网.[4] Wood J,1996,The geomorphological characterization of Digital Ele— vation Models, PhD Thesis Department of Geography, University of Lancaster, UK[5] 刘淑琼，邹时林，邹胜武，et al.基于规则格网DEM提取山顶点的 研究[J].东华理工大学学报(自然科学版)，2013(sl):93-95.[6] 陈盼盼，张友顺，王春，et al.基于DEM的山顶点快速提取技术 丹凤县扃都航摄地昭图 紫阳县扇都就摄地形图[J].现代测绘，2006, 29(02):11-13.图4丹凤县和紫阳县无人机航摄生成地形图高程点提取[7] 邹宝裕，董丞妍，苟娇娇，et al.数字高程模型提取山顶点及空间 5.结论格局差异[J].遥感信息，2016, 31(02):124-128.[8] 刘宏建，刘建忠，蔡中祥，et al. 一种基于流域剖分的山顶点提取 本文提出了一种自适应动态的高程点提取算法，对传 统的高程点提取算法进行了改进，该算法相对于传统算法, 150/RESOURCES方法山.测绘科学技术学报,2014(02):12-15. 因篇幅问题不能全部显示，请点此查看更多更全内容