CN116538996A

CN116538996A - 基于激光雷达的地形测绘系统及方法

Info

Publication number: CN116538996A
Application number: CN202310811599.6A
Authority: CN
Inventors: 李兴辉; 吕则航; 陈锐; 王秀芝
Original assignee: Yunnan Chaotu Geographic Information Co ltd
Current assignee: Yunnan Chaotu Geographic Information Co ltd
Priority date: 2023-07-04
Filing date: 2023-07-04
Publication date: 2023-08-04
Anticipated expiration: 2043-07-04
Also published as: CN116538996B

Abstract

本发明涉及激光雷达技术领域，尤其为基于激光雷达的地形测绘系统及方法，包括：数据采集模块：用于通过激光雷达采集地形三维数据；数据处理模块：用于对采集的地形点云数据进行降噪处理；数据匹配模块：用于对降噪处理后的地形点云数据进行特征提取和地形相似度计算；地形测绘模块：用于根据地形相似度计算结果选取相似度最高的结果进行地形测绘。本发明基于激光雷达采集获取地形点云数据，通过引导滤波对地形点云数据进行降噪处理，不仅能去除小尺度噪声，而且能保持地形点云边缘；并利用激光雷达结合采集的地形点云数据进行数据实时更新匹配，其计算更新的速度快、精度高，能够获得最佳的地形测绘图像。

Description

基于激光雷达的地形测绘系统及方法

技术领域

本发明涉及激光雷达技术领域，尤其是基于激光雷达的地形测绘系统及方法。

背景技术

地形测绘图内包含地面水平情况、地质构成、水文信息等。现今在测量地形时，通常会用到大型水平设备、勘探设备等多个设备，测量装备组成相对复杂，且不同设备间信息格式不同造成信息格式转换过程中数据丢失，导致地形测绘图不够详尽。故现有技术采用激光雷达技术进行地形的测绘。但现有技术基于获取地区的点云数据进行地形测绘，但由于一些人为或环境因素，以及激光雷达扫描装置本身的缺陷，通常会导致获取的点云数据带有一定的噪声，而这些噪声数据会给后续建模和测量中的相关处理带来严重的问题。

故本发明提出基于激光雷达的地形测绘系统及方法，基于激光雷达采集获取地形点云数据，通过引导滤波对地形点云数据进行降噪处理，不仅能去除小尺度噪声，而且能保持地形点云边缘。并利用激光雷达结合采集的地形点云数据进行数据实时更新匹配，获得最佳的地形测绘图像。

发明内容

本发明的目的是通过提出基于激光雷达的地形测绘系统及方法，以解决上述背景技术中提出的缺陷。

本发明采用的技术方案如下：

提供基于激光雷达的地形测绘系统，包括：

数据采集模块：用于通过激光雷达采集地形三维数据；

数据处理模块：用于对采集的地形点云数据进行降噪处理；

数据匹配模块：用于对降噪处理后的地形点云数据进行特征提取和地形相似度计算；

地形测绘模块：用于根据地形相似度计算结果选取相似度最高的结果进行地形测绘。

作为本发明的一种优选技术方案：所述数据采集模块基于三维激光扫描获取地形点云数据与地物信息数据。

作为本发明的一种优选技术方案：所述数据处理模块基于数据滤波算法对采集的地形点云数据进行降噪处理。

作为本发明的一种优选技术方案：所述数据滤波算法中，建立获取的地形点云数据的K-D树的数据结构。

作为本发明的一种优选技术方案：所述数据滤波算法具体如下：

基于建立的K-D树结构选取测试点，基于测试点在半径r内的临近点进行搜索：

计算的质心：

其中，表示在半径r内点的临近点；

假设地形点云数据的局部线性模型为：

其中，、分别为线性模型的系数，为点的第个临近点，为过滤点；

对点在KNN邻域的极小值进行求解：

其中，为平滑参数；

对取极小值时，对和的偏导为零；得到：

得到点的临近点的滤波点云：

遍历所有地形点云数据，得到滤波后的地形点云数据并保存。

作为本发明的一种优选技术方案：所述数据匹配模块中，选取三维坐标点及半径r，遍历滤波点云数据，对于每个滤波点云，判断其与点的距离是否在半径r内，如果点与点的距离在半径r内，将点加入邻域点集，计算邻域点集中点的坐标的平均值，得到邻域内点的均值，基于过滤后生成的地形点云数据，获取三维坐标点选择与其距离在半径r内的点构成其邻域点集，计算三维坐标点与其邻域点集的构成的协方差张量：

其中，表示邻域内三维点的数量，表示邻域内的点坐标，表示邻域内所有坐标的均值，表示数据的转置。

作为本发明的一种优选技术方案：块中，基于协方差张量进行几何特征值提取，并进行特征值提取获得表达因子和，具体如下：计算协方差张量：对于每个点，计算其邻域内点的协方差张量；选择阈值调整特征值，并将特征值从小到大排序；对于每个特征值，计算对应的权重：

其中，为数学常量；

使用特征值和权重来计算表达因子和，具体如下：

其中，和分别是第2个和第3个特征值的权重，和分别是对应的特征向量，，和分别表示第1、2和3个特征值。

作为本发明的一种优选技术方案：所述数据匹配模块的地形相似度计算中，基于局部几何特征提取获取的地形点云数据的表达因子和，令和分别表示地形三维点云图像，其中，为采集的地形三维点云图像，为数据库中存储的地形三维点云图像；则和两个图像分别对于和的表达因子、、、表达如下：

其中，表示采集的地形三维点云图像对于的第个表达因子，表示数据库中存储的地形三维点云图像对于的第个表达因子，表示采集的地形三维点云图像对于的第个表达因子，表示数据库中存储的地形三维点云图像对于的第个表达因子；

表达因子与图像不相似程度如下：

其中，和表示欧式距离，和分别表示和两个图像的空间相似度，且，表示图像的表达因子归一化参数，对表达的因子进行归一化处理后，地形相似度表达如下：

计算地形面积相似度时，以面积重叠率作为衡量其相似度指标，令表示面积相似度；

其中，表示面积重叠率，、分别表示图像重叠位置面积和图像的总面积，设定约束条件如下：

其中，表示图像重叠位置面积和图像的总面积计算过程中的预设误差；

图像面积相似度如下：

令图像和的方向角分别为、，表示和两个图像的方向相似度：

其中，为数学常量。

作为本发明的一种优选技术方案：所述地形测绘模块基于所述数据匹配模块计算的地形相似度、面积相似度和方向相似度分别进行排序，选取相似度最大的结果作为地形测绘匹配结果。

提供基于激光雷达的地形测绘方法，包括如下步骤：

S1：通过激光雷达采集地形三维数据；

S2：对采集的地形点云数据进行降噪处理；

S3：对降噪处理后的地形点云数据进行特征提取和地形相似度计算；

S4：根据地形相似度计算结果选取相似度最高的结果进行地形测绘。

本发明提供的基于激光雷达的地形测绘系统及方法，与现有技术相比，其有益效果有：

本发明基于激光雷达采集获取地形点云数据，通过引导滤波对地形点云数据进行滤波处理，不仅能去除小尺度噪声，而且能保持地形点云边缘。并利用激光雷达结合采集的地形点云数据进行数据实时更新匹配，其计算更新的速度快、精度高，能够获得最佳的地形测绘图像。

附图说明

图1为本发明优选实施例的系统框图；

图2为本发明优选实施例中方法流程图。

图中各个标记的意义为：100、数据采集模块；200、数据处理模块；300、数据匹配模块；400、地形测绘模块。

具体实施方式

需要说明的是，在不冲突的情况下，本实施例中的实施例及实施例中的特征可以相互组合，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

参照图1，本发明优选实施例提供了基于激光雷达的地形测绘系统，包括：

数据采集模块100：用于通过激光雷达采集地形三维数据；

数据处理模块200：用于对采集的地形点云数据进行降噪处理；

数据匹配模块300：用于对降噪处理后的地形点云数据进行特征提取和地形相似度计算；

地形测绘模块400：用于根据地形相似度计算结果选取相似度最高的结果进行地形测绘。

数据采集模块100基于三维激光扫描获取地形点云数据与地物信息数据。

数据处理模块200基于数据滤波算法对采集的地形点云数据进行降噪处理。

数据滤波算法中，建立获取的地形点云数据的K-D树的数据结构。

数据滤波算法具体如下：

计算的质心：

其中，表示在半径r内点的临近点；

假设地形点云数据的局部线性模型为：

对点在KNN邻域的极小值进行求解：

其中，为平滑参数，平滑参数的值通常小于1；

对取极小值时，对和的偏导为零；得到：

得到点的临近点的滤波点云：

数据匹配模块300中，选取三维坐标点和半径r，遍历滤波点云数据，对于每个滤波点云，判断其与点的距离是否在半径r内。如果点与点的距离在半径r内，将点加入邻域点集。计算邻域点集中点的坐标的平均值，得到邻域内点的均值，基于过滤后生成的地形点云数据，获取三维坐标点选择与其距离在半径r内的点构成其邻域点集，计算三维坐标点与其邻域点集构成的协方差张量：

数据匹配模块300中，基于协方差张量进行几何特征值提取，并进行特征值提取获得表达因子和。具体地，本实施例中采用如下步骤：

1.计算协方差张量：对于每个点，使用之前描述的方法计算其邻域内点的协方差张量。

2.设置阈值：选择适当的阈值来调整特征值的影响，将特征值从小到大排序，记作。

3.计算权重：对于每个特征值，计算对应的权重：。

4.计算表达因子：使用特征值和权重来计算表达因子和。使用以下公式：

其中，和分别是第2个和第3个特征值的权重，和分别是对应的特征向量，，和分别表示第1、2和3个特征值。采用阈值和Sigmoid函数来计算权重，实现了自适应调整的能力。通过调整阈值，可以灵活控制特征值对表达因子的影响。这使得表达因子的计算可以根据具体数据的特点进行调整，适应不同地形的几何特征。采用Sigmoid函数作为映射函数，引入了非线性映射的特性。Sigmoid函数的非线性特点使得权重与特征值之间的关系更加灵活，可以处理不同特征值区间的变化。这样可以更好地适应地形数据的局部几何特征，提高表达因子的表征能力。通过权重的计算，可以控制特征值对表达因子的贡献程度。较大的权重意味着对应的特征值在表达因子中的影响更大，而较小的权重则降低了特征值的贡献。这种方式可以灵活地调整不同特征值对表达因子的重要性，使得表达因子更加准确地反映地形的局部几何特征。

数据匹配模块300的地形相似度计算中，基于局部几何特征提取获取的地形点云数据的表达因子和，令和分别表示地形三维点云图像，其中，为采集的地形三维点云图像，为数据库中存储的地形三维点云图像；则和两个图像分别对于和的表达因子、、、表达如下：

表达因子与图像不相似程度如下：

图像面积相似度如下：

其中，为数学常量。

其中，本实施例为了保证相似度计算更符合地形相似度的需求，采用迭代的方法来搜索最佳的归一化参数：

1、首先使用当前的参数值计算相似度，得到相似度值，并设置一个目标相似度值，接下来判断相似度：比较当前的相似度值和目标相似度值：

如果当前相似度值接近目标相似度值，满足设定条件，停止迭代并接受当前的参数值。

如果当前相似度值小于目标相似度值，表明参数需要增加，进入下一步。

如果当前相似度值大于目标相似度值，表明参数需要减小，进入下一步。

2、根据相似度的比较结果来调整参数值。

如果当前相似度值小于目标相似度值，增加的值，，其中是增加的步长。

如果当前相似度值大于目标相似度值，减小的值，，其中是减小的步长。

3、使用更新后的参数值计算新的相似度，并继续迭代，直到满足设定的相似度条件。

其中，目标相似度值通过绘制相似度核密度估计图来分析数据的相似度分布情况，确定合理的相似度阈值作为目标相似度值。

地形测绘模块400基于数据匹配模块300计算的地形相似度、面积相似度和方向相似度分别进行排序，选取相似度最大的结果作为地形测绘匹配结果。

参照图2，提供基于激光雷达的地形测绘方法，包括如下步骤：

S1：通过激光雷达采集地形三维数据；

S2：对采集的地形点云数据进行降噪处理；

本实施例中，数据采集模块100基于安装在无人机上的激光雷达采集某地区地形三维数据，包括该地区的三维点云数据和多种地物信息。数据处理模块200基于获取的地区的三维点云数据搭建K-D树的数据结构，设共有1000个测试点，选取测试点，即第536个测试点，基于测试点在半径r内的临近点进行搜索。

计算的质心：

其中，表示在半径r内点的临近点；

假设地形点云数据的局部线性模型为：

对点在KNN邻域的极小值进行求解：

其中，为平滑参数，平滑参数的值通常小于1；

对取极小值时，对和的偏导为零；得到：

得到点的临近点的滤波点云：

本滤波算法通过利用导引点云的邻域，将滤波后的输出点推导为导引点云对应点的线性模型，通过此类滤波方法不仅能去除小尺度噪声，而且能保持地形点云边缘。

基于滤波后的地形点云数据，计算测试点点与其临近点集的构成的协方差张量，并以此提取局部几何特征；

对于给定的三维坐标点选择与其距离在半径r内的点构成其邻域点集，该邻域点集用于计算协方差张量：

其中，表示邻域内三维点的数量，表示邻域内的点坐标，表示邻域内所有坐标的均值，表示数据的转置。此时，设。

基于协方差张量进行几何特征值提取，并进行几何特征值的提取获得表达因子和。

基于几何特征值的提取获取的地形点云数据的表达因子和，引入地形相似度计算过程中，通过表达因子和衡量地形相似程度，令和分别表示地形三维点云图像，其中，为采集的地形三维点云图像，为数据库中存储的地形三维点云图像；则和两个图像分别对于和的表达因子、、、表达如下：

表达因子与图像不相似程度如下：

其中，和表示欧式距离，和分别表示和两个图像的空间相似度，表达因子与图像不相似程度取值区间为0-1之间，且，对表达的因子进行归一化处理后，地形相似度表达如下：

计算地形面积相似度时，以面积重叠率作为衡量其相似度指标，令表示面积相似度，面积重叠率由表示：

其中，、分别表示图像重叠位置面积和图像的总面积，当面积重叠率数值接近0时，表示两个图像的相似度较大，令表示图像重叠位置面积和图像的总面积计算过程中的预设误差，设定约束条件如下：

图像面积相似度如下：

令图像和的方向角分别为30°、40°，表示两个图像的方向相似度：

其中，为数学常量。

根据地形、面积和方向上的相似度比较和匹配，能够根据点云数据进行实时更新匹配，并且其计算更新的速度快、精度高，能够最大限度地还原地形图像。

地形测绘模块400根据数据匹配模块300计算的地形相似度、面积相似度和方向相似度分别进行排序，选取相似度最大的结果作为地形测绘匹配结果进行地形测绘。

对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。

此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims

1.基于激光雷达的地形测绘系统，其特征在于，包括：

数据采集模块（100）：用于通过激光雷达采集地形三维数据；

数据处理模块（200）：用于对采集的地形点云数据进行降噪处理；

数据匹配模块（300）：用于对降噪处理后的地形点云数据进行特征提取和地形相似度计算；

地形测绘模块（400）：用于根据地形相似度计算结果选取相似度最高的结果进行地形测绘。

2.根据权利要求1所述的基于激光雷达的地形测绘系统，其特征在于：所述数据采集模块（100）基于三维激光扫描获取地形点云数据与地物信息数据。

3.根据权利要求1所述的基于激光雷达的地形测绘系统，其特征在于：所述数据处理模块（200）基于数据滤波算法对采集的地形点云数据进行降噪处理。

4.根据权利要求3所述的基于激光雷达的地形测绘系统，其特征在于：所述数据滤波算法中，建立获取的地形点云数据的K-D树的数据结构。

5.根据权利要求4所述的基于激光雷达的地形测绘系统，其特征在于：所述数据滤波算法具体如下：

计算的质心：

其中，表示在半径r内点的临近点；

假设地形点云数据的局部线性模型为：

对点在KNN邻域的极小值进行求解：

其中，为平滑参数；

对取极小值时，对和的偏导为零；得到：

得到点的临近点的滤波点云：

6.根据权利要求5所述的基于激光雷达的地形测绘系统，其特征在于：所述数据匹配模块（300）中，选取三维坐标点及半径r，遍历滤波点云数据，对于每个滤波点云，判断其与点的距离是否在半径r内，如果点与点的距离在半径r内，将点加入邻域点集，计算邻域点集中点的坐标的平均值，得到邻域内点的均值，基于过滤后生成的地形点云数据，获取三维坐标点选择与其距离在半径r内的点构成其邻域点集，计算三维坐标点与其邻域点集构成的协方差张量：

7.根据权利要求6所述的基于激光雷达的地形测绘系统，其特征在于：所述数据匹配模块（300）中，基于协方差张量进行几何特征值提取，并进行特征值提取获得表达因子和，具体如下：计算协方差张量：对于每个点，计算其邻域内点的协方差张量；选择阈值调整特征值，并将特征值从小到大排序；对于每个特征值，计算对应的权重：

其中，为数学常量；

使用特征值和权重来计算表达因子和，具体如下：

8.根据权利要求7所述的基于激光雷达的地形测绘系统，其特征在于：所述数据匹配模块（300）的地形相似度计算中，基于局部几何特征提取获取的地形点云数据的表达因子和，令和分别表示地形三维点云图像，其中，为采集的地形三维点云图像，为数据库中存储的地形三维点云图像；则和两个图像分别对于和的表达因子、、、表达如下：

表达因子与图像不相似程度如下：

图像面积相似度如下：

其中，为数学常量。

9.根据权利要求8所述的基于激光雷达的地形测绘系统，其特征在于：所述地形测绘模块（400）基于所述数据匹配模块（300）计算的地形相似度、面积相似度和方向相似度分别进行排序，选取相似度最大的结果作为地形测绘匹配结果。

10.基于激光雷达的地形测绘方法，根据权利要求1-9任意一项所述的基于激光雷达的地形测绘系统，其特征在于，包括如下步骤：

S1：通过激光雷达采集地形三维数据；

S2：对采集的地形点云数据进行降噪处理；