一种面向复杂地形机载激光雷达数据的地面点提取方法
【技术领域】
[0001] 本发明属于空间信息技术领域，具体涉及一种面向复杂地形机载激光雷达数据的 地面点提取方法。
【背景技术】
[0002] 机载激光雷达（Airborne Light Detection And Ranging, LiDAR)是一种快速发展 的地形测绘技术。机载激光雷达的工作原理是根据脉冲从发射到接收地表反射回的脉冲信 号的传播时间来计算激光出射点到地表的距离，然后利用GPS (Global Position System, 全球定位系统）提供的发射点的三维空间信息以及姿态测量系统获得的飞行平台的姿态信 息来确定地表反射点的坐标。该技术在一些使用传统摄影测量方法进行地形测绘比较困难 的地区，如植被覆盖区、海岸带、岛礁地区等等，发挥着不可替代的作用。
[0003] 提取机载激光雷达数据中的地面点是指从机载激光雷达数据中将地面点与非地 面点分离开，它是提取数字高程模型（Digital Elevation Model，DEM)的核心技术。目 前常见的地面点提取方法，如基于形态学的算法、基于线性预测的算法、基于坡度的算法 等，主要是用于解决较平坦城区机载激光雷达数据的地面点提取问题，而专门针对复杂地 形，如地形起伏较大地区的地面点提取问题的算法相对较少。国际摄影测量与遥感协会 (International Society for Photogrammetry and Remote Sensing, ISPRS)特别指出现 有的地面点提取算法在地形起伏、陡坡陡坎多，同时分布有植被的区域，存在较大的分类误 差，因此本发明提出的面向复杂地形机载激光雷达数据的地面点提取方法具有较强的应用 价值。
【发明内容】
[0004] 为解决上述技术问题，本发明针对现有地面点提取算法在处理复杂地形机载激光 雷达数据时存在较大误差的问题，提供了一种面向复杂地形机载激光雷达数据的地面点提 取方法，以达到提高复杂地形区域的地面点提取精度的目的。
[0005] 实现上述目的，本发明之一种面向复杂地形机载激光雷达数据的地面点提取方 法，其特征在于：包括如下步骤， 步骤1，获取初始地面点集，利用移动格网法从机载激光雷达数据中提取部分的地面 点，并将该地面点称为初始地面点，而机载激光雷达数据中的其他激光脚点称为候选地面 占. 步骤2,对步骤1中提取到的初始地面点构建不规则三角网模型，该不规则三角网模型 是当前已经提取到的地面点构成的地表模型； 步骤3,计算每个候选地面点到已经提取到的地面点的最近距离，利用该距离和优先级 计算公式得到每个候选地面点的加密优先级； 步骤4,根据加密优先级逐一判断每个候选地面点，如果当前判断的候选地面点满足地 面点判断准则，则将该候选地面点加入到步骤2中的不规则三角网模型中，同时更新该候 选地面点邻域内的候选地面点的加密优先级； 步骤5,利用邻近距离判断法确定已经提取到的地面点中的低矮噪声点； 步骤6,如果步骤5中发现有低矮噪声点，则转步骤7,否则转步骤8 ; 步骤7,将步骤5中的低矮噪声点从机载激光雷达数据中剔除，不再参与后续算法步 骤，重新对已经提取到的地面点构建不规则三角网模型，转步骤3 ; 步骤8,将已经提取到的地面点作为算法最终得到的地面点，地面点提取过程结束。
[0006] 在上述方案基础上优选，步骤1中所述的移动格网法提取机载激光雷达数据中的 地面点进一步包括以下步骤， ① 获取移动网格步长
式中，step为格网移动步长；入取值为0. 8 ;average_interval为机载激光雷达数据 的平均脚点间距，也称为机载激光雷达数据的空间分辨率； ② 三维空间坐标是机载激光雷达数据所包含的最基本的属性之一，三维空间坐标的X、 Y、Z分量分别对应于WGS-84坐标系中的X、Y、Z轴的坐标值，统计机载激光雷达数据在X轴 和Y轴上的最小值和最大值，分别为minx、miny、maxx和maxy ; ③ 确定格网的坐标原点（orgx, orgy);其中，
④ 令maX_length为机载激光雷达数据中所包含的建筑物的最大边长，利用格网的坐 标原点（〇rgx，orgy)以及算法参数max_length对机载激光雷达数据进行二维格网划分，将 二维格网中每个网格内的最低点判定为地面点； ⑤ 如果orgx大于maxx或者orgy大于maxy，则步骤1结束，否则转步骤⑥； ⑥ 更新格网的坐标原点，计算公式如下，转步骤④；
在上述方案基础上优选，所述的步骤3中获取候选地面点的加密优先级进一步包括以 下步骤， ① 令已经提取到的地面点的集合为A，候选地面点的集合为B，搜索半径为R，且搜索半 径R设置为5倍的机载激光雷达数据的平均脚点间距； ② 将所有候选地面点的加密优先级设置为〇,对集合A进行八叉树组织； ③ 针对集合B中的每个候选地面点，进行如下操作：假设当前候选地面点为P，查找与 P之间的距离小于R的集合A中的所有地面点，计算这些地面点与P之间的距离，将其中的 最小距离min_distance，按照下列计算出P的加密优先级；
其中，式中的Priority (p)为P的加密优先级，min_distance为集合A中所有地面点 距离P之间的最小距离。
[0007] 在上述方案基础上优选，所述的步骤4中基于加密优先级的地面点提取方法进一 步包括以下步骤， ① 经过步骤3处理后得到现阶段已经提取到的地面点的加密优先级，令加密队列为 PQ，PQ具有如下性质：它能够保证加密优先级大的候选地面点被优先判断，PQ的具体实现 形式为最大堆，步骤4需要两个算法参数，分别为搜索半径R和增长比率k，搜索半径R的大 小为5倍的机载激光雷达数据的平均脚点间距，k=l. 5,令已经提取到的地面点的集合为A， 候选地面点的集合为B ; ② 如果集合B为空，则步骤4结束，否则转步骤③； ③ 将加密队列PQ置空，然后依次插入集合B中的候选地面点，每个候选地面点中携带 有对应的加密优先级，PQ根据加密优先级的大小对候选地面点从大到小排序； ④ 从加密队列PQ中取出加密优先级最大的候选地面点q，并从集合B中将q剔除，转步 骤⑤； ⑤ 判断候选地面点q是否满足以下地面点条件，如果不满足地面点条件，转步骤⑦，否 则转步骤⑥；
式中，令T为候选地面点q投影到地表模型时所位于的三角形，d为q到T的垂直距 离，0为连接q和T的三个顶点所构成的边与T所形成的三个夹角中的最小角，设置为 〇. 15米，7；设置为15。；
式中，L为候选地面点q到T的三个顶点的最近距离，令P1为该最近距离对应的T的 顶点，则根据地表模型可以计算出P1处的地形方向V，5?为q和P1的连线与V所构成的夹 角，凡设置为3倍的机载激光雷达数据的平均脚点间距，T5.设置为1〇° ; ⑥ 将q加入到集合A中，并更新地表模型，然后重新计算q在半径R内的候选地面点的 加密优先级，如果计算出来的新加密优先级大于原加密优先级，则更新其在加密队列PQ中 的位置； ⑦ 如果加密队列PQ为空，说明当前搜索半径R内的候选地面点已判断完毕，转步骤⑧， 否则转步骤④； ⑧ 扩大搜索半径R，R=kXR，其中k=l. 5,转步骤⑦。
[0008] 在上述方案基础上优选，所述的步骤5中查找机载激光雷达数据中低矮噪声点的 邻近距离判断步骤进一步还包括以下步骤， ①令丹为当前待判断的地面点，计算的k个邻近地面点将列和 〔A…巧J构成一个新的集合
包含的基本属性有三维空 间坐标，对集合丨~丹:根据其三维空间坐标的2分量从小到大进行排序，得到有序数组
，而k值大小可以根据以下公式计算得到
其中，area为搜索面积，area设置为5平方米，average_interval为机载激光雷达数 据的平均脚点间距； ② 计算有序数组V中相邻地面点间的高差
，其中
③ 从j=〇开始到k-1，根据以下公式判