明提出的移动格网算法能够高效地提取大量初始地面点，能够为后续算法提供 相对精确的初始地表模型； 2. 加密优先级的引入可以保证地面点提取过程中距离已提取到的地面点较近的候选 地面点被优先判定，从而减少了因为地形起伏而造成植被等被错误判定为地面点的情况的 发生； 3. 利用邻近距离判断法进行低矮噪声点判断可以显著提高算法对低矮噪声点的鲁棒 性，进而提_算法的提取精度。
[0022] 综上所述，仅为本发明之较佳实施例，不以此限定本发明的保护范围，凡依本发明 专利范围及说明书内容所作的等效变化与修饰，皆为本发明专利涵盖的范围之内。
【主权项】
1. 一种面向复杂地形机载激光雷达数据的地面点提取方法，其特征在于：包括如下步 骤， 步骤1，获取初始地面点集，利用移动格网法从机载激光雷达数据中提取部分的地面 点，并将该地面点称为初始地面点，而机载激光雷达数据中的其他激光脚点称为候选地面 占 . 步骤2,对步骤1中提取到的初始地面点构建不规则三角网模型，该不规则三角网模型 是当前已经提取到的地面点构成的地表模型； 步骤3,计算每个候选地面点到已经提取到的地面点的最近距离，获取每个候选地面 点的加密优先级； 步骤4,根据加密优先级逐一判断每个候选地面点，如果当前判断的候选地面点满足地 面点判断准则，则将该候选地面点加入到步骤2中的不规则三角网模型中，同时更新该候 选地面点邻域内的候选地面点的加密优先级； 步骤5,确定已经提取到的地面点中的低矮噪声点； 步骤6,如果步骤5中发现有低矮噪声点，则转步骤7,否则转步骤8 ; 步骤7,将步骤5中的低矮噪声点从机载激光雷达数据中剔除，不再参与后续算法步 骤，重新对已经提取到的地面点构建不规则三角网模型，转步骤3 ; 步骤8,将已经提取到的地面点作为算法最终得到的地面点，地面点提取过程结束。2. 如权利要求1所述的一种面向复杂地形机载激光雷达数据的地面点提取方法，其特 征在于：步骤1中所述的移动格网法提取机载激光雷达数据中的地面点进一步包括以下步 骤， ① 获取移动网格步长式中，step为格网移动步长；λ取值为〇. 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. 如权利要求1所述的一种面向复杂地形机载激光雷达数据的地面点提取方法，其特 征在于：所述的步骤3中获取候选地面点的加密优先级进一步包括以下步骤， ① 令已经提取到的地面点的集合为A，候选地面点的集合为B，搜索半径为R，且搜索半 径R设置为5倍的机载激光雷达数据的平均脚点间距； ② 将所有候选地面点的加密优先级设置为〇,对集合A进行八叉树组织； ③ 针对集合B中的每个候选地面点，进行如下操作：假设当前候选地面点为P，查找与 P之间的距离小于R的集合A中的所有地面点，计算这些地面点与P之间的距离，将其中的 最小距离min_distance，按照下列计算出P的加密优先级；其中，式中的Priority (p)为P的加密优先级，min_distance为集合A中所有地面点 距离P之间的最小距离。4. 如权利要求1所述的一种面向复杂地形机载激光雷达数据的地面点提取方法，其特 征在于：所述的步骤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的垂直距离， Θ为连接q和T的三个顶点所构成的边与T所形成的三个夹角中的最小角，％.设置为0.15 米，?设置为15° ;式中，L为候选地面点q到T的三个顶点的最近距离，令Pl为该最近距离对应的T的 顶点，则根据地表模型可以计算出Pl处的地形方向V，δ为q和Pl的连线与V所构成的夹 角，Tli设置为3倍的机载激光雷达数据的平均脚点间距，τ,设置为10° ; ⑥ 将q加入到集合A中，并更新地表模型，然后重新计算q在半径R内的候选地面点的 加密优先级，如果计算出来的新加密优先级大于原加密优先级，则更新其在加密队列PQ中 的位置； ⑦ 如果加密队列PQ为空，说明当前搜索半径R内的候选地面点已判断完毕，转步骤⑧， 否则转步骤④； ⑧ 扩大搜索半径R，R=kXR，其中k=l. 5,转步骤⑦。5. 如权利要求2所述的一种面向复杂地形机载激光雷达数据的地面点提取方法，其特 征在于：所述的步骤5中查找机载激光雷达数据中低矮噪声点的邻近距离判断步骤进一步 还包括以下步骤， ① 令Pg为当前待判断的地面点，计算内的k个邻近地面点丨巧…為），将/?和 !+IV"Aj构成一个新的集合，其中A (〇幻^ )包含的基本属性有三维空间坐标，对集合 丨根据其三维空间坐标的Z分量从小到大进行排序，得到有序数组V = %, 而k值大小可以根据以下公式计算得到其中，area为搜索面积，area设置为5平方米，average_interval为机载激光雷达数 据的平均脚点间距； ② 计算有序数组V中相邻地面点间的高差D ,其中③ 从j=〇开始到k-Ι，根据以下公式判断是否为真，如果为真，则丹将从地面点中剔 除。6. 其中，f为最大允许的高程，设置为2米，Wma5i为噪声点集的最大点数， 设置为5。
【专利摘要】一种面向复杂地形机载激光雷达数据的地面点提取方法，包括如下步骤，获取初始地面点集，利用移动格网法提取机载激光雷达数据中的地面点；对初始地面点构建不规则三角网模型，该不规则三角网模型是已经提取到的地面点构成的地表模型；计算每个候选地面点到已经提取到的地面点的最近距离，获取每个候选地面点的加密优先级；根据加密优先级逐一判断每个候选地面点，更新该脚点领域内候选地面点的加密优先级；将已经提取到的地面点作为最终得到的地面点。本发明针对现有地面点提取算法在处理复杂地形机载激光雷达数据存在较大误差的问题，提供了一种面向复杂地形机载激光雷达数据的地面点提取方法，以达到提高复杂地形区域的地面点提取精度的目的。
【IPC分类】G01S17/89, G01S7/48
【公开号】CN104898106
【申请号】CN201410725773
【发明人】高广, 张良, 凡孝举, 潘振华, 胡高文
【申请人】武汉天擎空间信息技术有限公司
【公开日】2015年9月9日
【申请日】2014年12月4日