(10)申请公布号 CN 102915560 A
(43)申请公布日 2013.02.06C N 102915560 A
*CN102915560A*
(21)申请号 201210353952.2
(22)申请日 2012.09.21
G06T 17/00(2006.01)
G01S 7/48(2006.01)
(71)申请人中国石油大学（华东）
地址266580 山东省青岛市经济技术开发区
长江西路66号地球科学与技术学院
(72)发明人曾喆 万剑华 黄荣刚 刘善伟
(74)专利代理机构东营双桥专利代理有限责任
公司 37107
代理人
侯华颂
(54)发明名称
及装置
(57)摘要
本发明涉及一种数字化制图技术领域，特别
涉及一种阈值无关的机载激光雷达点云滤波方法
及装置。其技术方案是由以下步骤实现：步骤1）
构建；步骤2）点云-地面高差数据生成；步骤3）
基于点云-地面高差的统计滤波处理。有益效果
是：利用这种整体统计特征，本发明的方法及装
置可以不采用与点云数据的高差、坡度等几何特
征直接相关的阈值来实现机载Lidar 点云数据的
滤波过程，这样就降低了整个滤波数据处理过程
人工干预量，提高了整个机载Lidar 点云数据处
理的工作效率。
(51)Int.Cl.
权利要求书2页 说明书6页 附图2页
(19)中华人民共和国国家知识产权局(12)发明专利申请
权利要求书 2 页 说明书 6 页 附图 2 页
1.一种阈值无关的机载激光雷达点云滤波方法，其特征是由以下步骤实现：
步骤1）点云数据分割及初始地面不规则三角网的构建：主要是以矩形格网规则划分Lidar点云数据，将每个规则格网中的最低点作为地面点，并根据这些点构建初始地面不规则三角网，也称TIN网；控制单元首先将点云数据加载到点云数据存储器中，并在分割器模块中实现步骤1），控制单元则协助分割器完成点云数据存储器与分割器模块之间的数据传输；
步骤2）点云-地面高差数据生成：主要通过步骤1得到的初始地面TIN网，计算出点云数据的高差数据；其主要在统计数据生成器模块中完成，控制单元协助完成统计生成器与点云数据存储器之间的数据传输控制；
步骤3）基于点云-地面高差的统计滤波处理：实现Lidar点云数据的高差统计滤波处理，其主要是将步骤2得到的点云高差值通过计算其统计特征Tsallis熵，并利用熵最大化原则将其在概率分布上分离为两部分，即：一部分为地面点云A，一部分为非地面点云B，步骤3主要在统计滤波器模块中完成主要计算操作，控制单元协助完成统计滤波器与点云数据存储器之间的数据传输控制。
2.根据权利要求1所述的阈值无关的机载激光雷达点云滤波方法，其特征是：步骤1）的详细步骤如下：
1.1）点云数据平面格网分割：点云数据坐标信息主要有（x,y,h）三个方向，其中x，y 为水平方向，h为高程方向，分割器利用点云数据存储器中所有点云数据的x，y方向坐标，计算整个点云数据在x，y平面上的最小外接矩形，对最小外接矩形以矩形规则格网划分，分割器比较Lidar点云数据中的所有点的x，y坐标与划分好的矩形格网范围，将所有的Lidar点划分到其x，y坐标所属的矩形格网内，并将所有的点云数据根据其所属矩形格网按Morton码的顺序在点云数据存储器中重新存储；
1.2）搜索每个格网中点云的最小高程点：分割器对步骤1.1中排好序的点云分块数据遍历每个矩形格网，对在每个矩形格网中的点云数据实施下面两个步骤：
1.2.1）对每个矩形区域中的所有点云数据按其高程值h排序；
1.2.2）排好序后，取高程值最小为该矩形格网内初始地面点；
1.3）初始地面点的不规则三角网构建：对于在步骤1.
2.2)中得到的所有初始地面点构建TIN网，采用计算几何中的凸包法来构建初始地面点的TIN网，将TIN结果存放入点云数据存储器中。
3.根据权利要求1所述的阈值无关的机载激光雷达点云滤波方法，其特征是：步骤2）的详细步骤如下：
2.1) 内插地面高程值：遍历Lidar点云数据中的每个点云数据点实施如下操作：
2.1.1) 取点云数据点P0（x0,y0,h0）以及TIN网中所有顶点的x与y两个水平坐标，这样得到P0在水平面上的对应点P1（x0,y0）及TIN网在水平面上的对应三角网N1，在水平面上计算N1的梯形图，根据梯形图计算出水平面上点P1在N1的三角形A1内，并给出A1对应的TIN网中三角形A0；
2.1.2）计算点云数据点水平坐标（x0，y0）在A0所在平面内其所对应点（x0，y0，h1）；
2.1.3）将h1设置为P0点对应的地面高程值；
2.2) 计算点云数据高差：遍历Lidar点云数据中的每个点云数据点，将其高程值与步
骤2.1中得到的地面高程值求差，得到每个点云数据点的高差值；
2.3) 边界情况的高差计算：在边界矩形格网内的一些点云数据点无法得到步骤2.1.1的计算结果A0时，计算这些点云数据点的高程与其所在的边界矩形格网内的最低点高程之差，将其设为该点云数据点的高差数据。
4.根据权利要求1所述的阈值无关的机载激光雷达点云滤波方法，其特征是：步骤3）的详细步骤如下：
3.1）对步骤2中得到的点云-地面高差数据做概率分布统计：将步骤2中得到的所有Lidar点云高差值等分为N个区段，并计算点云数据点在每个高差值区段的概率为（取值为）；
3.2）以t（t为整数且取1~N）为分离点计算Tsallis熵，取t的初始值为，在高差值区段中t对应的高差值将地面点云和非地面点云分离，即小于等于t区段的高差值为地面点云，大于的则为非地面点云，分别按公式计算地面点云高差分布的Tsallis熵，非地面点云高差分布的Tsallis熵；
3.3）最大熵的优化计算：
3.3.1) 对t从初始值以1为步长进行增加，分别计算地面点云高差的Tsallis熵
S A 、非地面点云高差的Tsallis熵S
B
以及整体高差的Tsallis熵S
t
，
经过N-t
B
次迭代，可以
得到N-t
B 个整体高差的Tsallis熵值S
t
以及其对应的N-t
B
个分离点t值；
3.3.2)对步骤3.3.1)中计算得到的所有整体高差的Tsallis熵值S
t
排序，得到其最
大值，并获取最大值所对应的t值t
max
；
3.4）滤波分离地面与非地面点云：根据步骤3.3)中得到的最大Tsallis熵所对应的
t max 对点云进行滤波处理，将高差概率分布于t
max
右侧的点云滤除为非地面点云，即其高差
值大于t
max
区段的点云，剩余点云作为地面点云，并通过控制单元将地面点云放入点云数据存储器。
5.一种阈值无关的机载激光雷达点云滤波装置，其特征是：主要由控制单元、点云数据存储器、分割器、统计数据生成器和统计滤波器组成，其中，控制单元主要控制整个自动化点云滤波的过程以及各个模块之间的数据传输；点云数据存储器主要用于存放与点云数据相关的初始化数据，滤波过程中生成的中间数据以及最后得到的地面与非地面点云结果数据；分割器主要实现对点云数据的一些初始化及分割操作；统计数据生成器主要完成统计滤波中一些统计数据的生成，统计滤波器采用熵优化的原理实现点云数据的自动滤波处理。
6.根据权利要求5所述的阈值无关的机载激光雷达点云滤波装置，其特征是：所述的控制单元可以采用ARM处理器芯片来实现；分割器、统计数据生成器以及统计滤波器可以采用FPGA器件来实现。
7.根据权利要求6所述的阈值无关的机载激光雷达点云滤波装置，其特征是：所述的ARM处理器芯片采用SAMSUNG的S3C4510B；所述的FPGA器件采用Altera公司的EP1S80F1020C5。
一种阈值无关的机载激光雷达点云滤波方法及装置
技术领域
[0001] 本发明涉及一种数字化制图技术领域，特别涉及一种阈值无关的机载激光雷达点云滤波方法及装置，适合于针对机载Lidar点云实施自动化处理而将地面与非地面点云分离的一种统计滤波方法及装置。
背景技术
[0002] 近年来，机载激光雷达技术已成为获取高精度数字高程模型（DEM）的一种重要手段。通过机载激光雷达获取到的点云数据往往包含地面点云数据和非地面点云数据。为了构建数字高程模型(DEM)，通常要将Lidar中地面点云数据与非地面点云数据分离，也就是要把非地面点云数据通过滤波方法从原始Lidar点云数据中滤除，现有的Lidar点云数据滤波方法在整个滤波过程中需要人为设定一些经验阈值才能做滤波处理。人为经验阈值的调整往往导致在Lidar点云数据生产DEM过程中的工作量剧增。
[0003] LIU X. Airborne LiDAR for DEM generation: some critical issues [J]. Progress in Physical Geography, 2008, 32(1): 31-49.
以上这篇文章，给出了机载Lidar点云建立DEM的过程，其中点云数据的滤波处理是其中最难且最关键的一步，该文给出了目前广泛采用的滤波方法，这些方法中大多数需要一些经验知识给出参数或者阈值，因此，如何自动地通过滤波方法将地面和非地面点云分离是一个难点。
[0004] 专利(中国专利申请号201010262366.8 )在步骤3.2中采用预设的高程差阈值和面积比阈值实施滤波操作，因而这种经验阈值会导致整个方法的一种人工干预和调优。采用这种阈值相关的方法，往往存在一个潜在的问题：对于不同地形的Lidar点云数据滤波过程，人工经验阈值的设定调整将导致工作量剧增。
[0005] 从上面所述的关于Lidar点云数据滤波的现有技术方法来看，这些方法在实际的滤波过程往往会需要人工根据经验设定一些阈值，而设定这些阈值往往需要大量实践积累出经验，并根据这种经验来挑选阈值，因而实际上带来Lidar点云数据滤波工作量剧增以及处理过程的效率低下。
发明内容
[0006] 本发明的目的就是针对现有技术存在的上述缺陷，提供一种阈值无关的机载激光雷达点云滤波方法及装置。
[0007] 一种阈值无关的机载激光雷达点云滤波方法，其技术方案是由以下步骤实现：步骤1）点云数据分割及初始地面不规则三角网（TIN）的构建：主要是以矩形格网规则划分Lidar点云数据，将每个规则格网中的最低点作为地面点，并根据这些点构建初始地面TIN网；控制单元首先将点云数据加载到点云数据存储器中，并在分割器模块中实现步骤1，控制单元则协助分割器完成点云数据存储器与分割器模块之间的数据传输；
步骤2）点云-地面高差数据生成：主要通过步骤1得到的初始地面TIN网，计算出点
云数据的高差数据；其主要在统计数据生成器模块中完成，控制单元协助完成统计生成器与点云数据存储器之间的数据传输控制；
步骤3）基于点云-地面高差的统计滤波处理：实现Lidar点云数据的高差统计滤波处理，其主要是将步骤2得到的点云高差值通过计算其统计特征Tsallis熵，并利用熵最大化原则将其在概率分布上分离为两部分，即：一部分为地面点云A，一部分为非地面点云B，步骤3主要在统计滤波器模块中完成主要计算操作，控制单元协助完成统计滤波器与点云数据存储器之间的数据传输控制。
[0008] 步骤1）的详细步骤如下：
1.1）点云数据平面格网分割：点云数据坐标信息主要有（x,y,h）三个方向，其中x，y 为水平方向，h为高程方向，分割器利用点云数据存储器中所有点云数据的x，y方向坐标，计算整个点云数据在x，y平面上的最小外接矩形，对最小外接矩形以矩形规则格网划分，分割器比较Lidar点云数据中的所有点的x，y坐标与划分好的矩形格网范围，将所有的Lidar点划分到其x，y坐标所属的矩形格网内，并将所有的点云数据根据其所属矩形格网按Morton码的顺序在点云数据存储器中重新存储；
1.2）搜索每个格网中点云的最小高程点：分割器对步骤1.1中排好序的点云分块数据遍历每个矩形格网，对在每个矩形格网中的点云数据实施下面两个步骤：
1.2.1）对每个矩形区域中的所有点云数据按其高程值h排序；
1.2.2）排好序后，取高程值最小为该矩形格网内初始地面点；
1.3）初始地面点的不规则三角网构建：对于在步骤1.
2.2)中得到的所有初始地面点构建TIN网，采用计算几何中的凸包法来构建初始地面点的TIN网，将TIN结果存放入点云数据存储器中。
[0009] 步骤2）的详细步骤如下：
2.1) 内插地面高程值：遍历Lidar点云数据中的每个点云数据点实施如下操作：
2.1.1) 取点云数据点P0（x0,y0,h0）以及TIN网中所有顶点的x与y两个水平坐标，这样得到P0在水平面上的对应点P1及TIN网在水平面上的对应三角网N1，在水平面上计算N1的梯形图，根据梯形图计算出水平面上点P1在N1的三角形A1内，并给出A1对应的TIN网中三角形A0；
2.1.2）计算点云数据点水平坐标（x0，y0）在A0所在平面内的对应点（x0，y0，h1）；
2.1.3）将h1设置为P0点对应的地面高程值；
2.2) 计算点云数据高差：遍历Lidar点云数据中的每个点云数据点，将其高程值与步骤2.1中得到的地面高程值求差，得到每个点云数据点的高差值；
2.3) 边界情况的高差计算：在边界，矩形格网内的一些点云数据点无法得到步骤2.1.1的计算结果A0时，计算这些点云数据点的高程与其所在的边界矩形格网内的最低点高程之差，将其设为该点云数据点的高差数据。
[0010] 步骤3）的详细步骤如下：
3.1）对步骤2中得到的点云-地面高差数据做概率分布统计：将步骤2中得到的所有Lidar点云高差值等分为N个区段，并计算点云数据点在每个高差值区段的概率为（取值为）；
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议