一种基于多尺度噪声的激光雷达点云去噪方法技术

本发明专利技术提供一种基于多尺度噪声的激光雷达点云去噪方法，涉及三维点云去噪技术领域，包括根据不同尺度的噪声点进行分类处理；获取DBSCAN算法参数，将点云数据集分为不同形状的簇；使用Kd

【技术实现步骤摘要】
一种基于多尺度噪声的激光雷达点云去噪方法


[0001]本专利技术属于3D点云数据处理
，具体涉及一种基于多尺度噪声的激光雷达点云去噪方法。

技术介绍

[0002]随着激光雷达遥感技术的快速发展，点云数据由于其高精度、高效率、不受时间限制等优点，在地理测绘、灾害评估、人工智能等诸多研究领域得到了广泛应用。但是三维点云数据在获取时，由于设备的硬件缺陷、环境的不稳定性、人为操作等因素影响，导致点云数据不可避免地会存在大量的噪声点，直接影响到三维模型表面的精度，对后续的数据处理造成了困难，因此需要对采集的三维点云模型执行去噪的处理。
[0003]由于点云数据点散乱、庞大，常常伴随有奇异点、干扰点等，如何保持点云数据高频特征信息并提供去噪方法效率，是点云去噪的一大难度。传统点云滤波方法主要通过局部表面拟合、局部或非局部的平均或依赖于底层噪声模型的统计假设，主要有不规则三角网加密滤波、基于坡度的滤波、基于形态学的滤波和基于移动曲面的滤波等，但是每种算法有相应的应用场景，目前还没有一种算法能够适应所有场景。

技术实现思路

[0004]本专利技术的目的在于提供一种基于多尺度噪声的激光雷达点云去噪方法，以解决现有技术的不足，对于点云数据中存在的大尺度噪声点和小尺度噪声点，该方法能够在去除不同尺度噪声的同时保留模型的细节特征，提高去噪效率。
[0005]为实现上述目的，本专利技术提供了一种基于多尺度噪声的激光雷达点云去噪方法，包括以下步骤：
[0006]S1、输入点云数据；
[0007]S2、获取DBSCAN算法中Eps和MinPts参数；
[0008]S21、采用欧氏距离公式对未知数据的特征与数据集D中的每一个数据相对应的特征进行计算，得到数据集D的距离分布矩阵。然后对距离分布矩阵进行升序排列得到KNN分布矩阵，并按列取平均值为K值的平均距离，得到Eps参数列表；
[0009]以上所述欧式距离公式定义为：
[0010][0011]其中，Dist
n
×
n
为n
×
n的实对称矩阵，n为数据集D中的对象个数；dist
(i,j)
为数据集D中第i个对象到第j个对象的距离；
[0012]以上所述数据集D的距离分布矩阵定义为：
[0013]Dist
n
×
n
＝{dist(i,j)|1≤i≤n,1≤j≤n}
[0014]其中，Dist
n
×
n
为n
×
n的实对称矩阵，n为数据集D中的对象个数；dist
(i,j)
为数据集D中第i个对象到第j个对象的距离；
[0015]S22、给定数据集D，以Eps为邻域求出核心点邻域内数据点数量的期望值，采用数学期望法生成MinPts参数；
[0016]以上所述数学期望法定义为：
[0017][0018]S3、建立Kd
‑
tree加速邻域搜索，对大尺度噪声进行DBSCAN聚类去噪；
[0019]S31、将阈值参数代入DBACAN聚类算法，对点云数据集D进行读取，以Eps为邻域遍历核心点，寻找与核心点密度可达的数据点，进行聚类成簇，当簇数稳定，密度阈值最小时K值所对应的Eps、MinPts参数为最优参数；
[0020]S32、创建Kd
‑
tree加速搜索邻域，保留最优参数下DBSCAN算法生成的聚类簇，删除边界点和噪声点，得到去除大尺度噪声后的点云数据；
[0021]由上，步骤S32中数据点在K维空间中进行构造，树中的每一层对应一个维度，Kd
‑
tree通过数据点的属性个数循环的构造K维二叉搜索树.Kd
‑
tree中左子树在给定维度上的值小于父节点，而右子树则大于父节点。
[0022]S4、设定移动最小二乘法阶数，进行法线估计，确定拟合半径；
[0023]S41、建立K近邻邻域及局部区域的拟合函数，设定拟合函数中基函数的阶数为2；
[0024]在拟合区域的一个局部子域上，拟合函数f(x)定义为：
[0025][0026]其中，p
T
(x)＝[p1(x)p2(x)
…
p
m
(x)]为基函数，α
i
(x)(i＝0,1,...,n)为待求的系数矢量，系数对应的空间坐标为x＝[x,y,z]的函数；
[0027]S42、然后确定权函数为三次样条权函数，对每个数据点的邻域拟合移动最小二乘曲面；
[0028]所选取三次样条权函数定义为：
[0029][0030]其中，r＝x
‑
x
i
，h
i
为第i个节点权函数影响区域的范围，β为影响系数，取值范围通常为1.5～2.5，一般为2.0；
[0031]S43、计算局部K邻域范围内的曲面，根据曲面进行法线估计；
[0032]S44、创建Kd
‑
tree搜索邻域，确定拟合的K近邻半径；
[0033]S45、进行曲面重建，去除小尺度噪声；
[0034]S5、移动最小二乘法平滑去噪；
[0035]S6、输出点云数据。
[0036]上述点云去噪方法中，根据点云数据噪声点分布的特点，与单一的滤波方法相比，
DBSCAN算法具有很强的适应性，能够对点云数据进行聚类成簇，将点云数据中的核心点、边界点、噪声点进行分类处理，去除明显的离群点，针对在模型周围仍然存在一些小尺度的噪声，使用移动最小二乘法进行二次去噪。该方法对点云数据中存在的多尺度噪声能够有效去除，并保持了点云模型的细节部分，解决了两种算法单独使用的局限性。
附图说明
[0037]图1是本专利技术实施例提供的方法的流程示意图；
[0038]图2是本专利技术实施例中基于改进DBSCAN算法的具体去噪流程图；
[0039]图3是本专利技术实施例中点云数据原始加噪图；
[0040]图4是本专利技术实施例中基于改进DBSCAN算法的具体去噪效果图；
[0041]图5是本专利技术实施例中基于移动最小二乘法的去噪流程图；
[0042]图6是本专利技术实施例中基于移动最小二乘法的去噪效果图。
具体实施方式
[0043]下面结合附图和具体实施方式对本专利技术进行详细说明；
[0044]本专利技术基于多尺度噪声的激光雷达点云去噪方法，如图1所示，具体按照
[0045]以下步骤实施：
[0046]S1、输入点云数据；
[0047]S2、获取DBSCAN算法中Eps和MinPts参数；
[0048]S21、采用欧氏距离公式对未知数据的特征与数据集D中的每一个数据相对应的特征进行计算，得到数据集D的距离分布矩阵。然后对距离分布矩阵进行升序排列得到KNN分布矩阵，并按列取平均值为K值的平均距本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于多尺度噪声的激光雷达点云去噪方法，其特征在于，包括以下步骤：S1、输入点云数据；S2、获取DBSCAN算法中Eps和MinPts参数；S3、建立Kd
‑
tree加速邻域搜索，对大尺度噪声进行DBSCAN聚类去噪；S4、设定移动最小二乘法阶数，进行法线估计，确定拟合半径；S5、移动最小二乘法平滑去噪；S6、输出点云数据。2.根据权利要求1所述的基于多尺度噪声的激光雷达点云去噪方法，其特征在于，步骤S2中，DBSCAN算法是一种基于密度的空间的数据聚类方法；该方法将簇定义为密度相连的点的最大集合，把具有足够密度的区域划分为簇，并在有噪声的空间数据集中发现任意形状的簇。3.根据权利要求1所述的基于多尺度噪声的激光雷达点云去噪方法，其特征在于，步骤S2中，获取DBSCAN算法参数的具体方式为：S21、如果数据集在特征空间中的K个最相似的数据集中的大多数属于某一个类别，则该样本也属于这个类别；采用欧氏距离公式对未知数据的特征与数据集D中的每一个数据相对应的特征进行计算，得到数据集D的距离分布矩阵；然后对距离分布矩阵进行升序排列得到KNN分布矩阵，并按列取平均值为K值的平均距离，得到Eps参数列表；所述欧式距离公式定义为：其中，d表示距离，xi表示未知数的特征，yi表示训练集中的每一个数据相对应的特征；所述数据集D的距离分布矩阵定义为：Dist
n
×
n
＝{dist(i,j)|1≤i≤n,1≤j≤n}其中，Dist
n
×
n
为n
×
n的实对称矩阵，n为数据集D中的对象个数；dist
(i,j)
为数据集D中第i个对象到第j个对象的距离；S22、给定数据集D，以Eps为邻域求出核心点邻域内数据点数量的期望值，采用数学期望法生成MinPts参数；所述数学期望法定义为：其中，n为数据集D中核心对象的总数量，Pi为第i个对象的Eps邻域的核心对象的数量。4.根据权利要求1所述的基于多尺度噪声的激光雷达点云去噪方法，其特征在于，步骤S3中，DBSCAN聚类去噪具体实现方式为：S31、将阈值参数代入DBACAN聚类算法，对点云数据集D进行读取，以Eps为邻域遍历核心点，寻找与核心点密度...

【专利技术属性】
技术研发人员：杨家志，沈洁，陈梦强，于广旺，周国清，
申请(专利权)人：桂林理工大学，
类型：发明
国别省市：