一种高效高精度计算刚性点云曲面之间距离的方法技术

本发明专利技术公开了一种高效高精度计算刚性点云曲面之间距离的方法，包括以下步骤：S1、对点云进行区域划分，得到磨损区域和非磨损区域；对待测点云进行处理，得到纯净的待测点云；S2、将非磨损区域待测点云和参考点云进行基于点云不确定度的ICP迭代配准方法的刚性配准，获取非磨损区域待测点云的不确定系数和不确定度；S3、高精度磨损测量，获取点云曲面之间的距离。本发明专利技术解决在测量点云和目标点云存在形貌差别的情况下如何进行高精度的配准问题，解决在测量点云和目标点云进行配准时如何保证配准不陷入局部解并达到全局最优的问题，解决目前点云测量中测量精度低，无法满足高精度测量的问题。的问题。的问题。

【技术实现步骤摘要】
一种高效高精度计算刚性点云曲面之间距离的方法


[0001]本专利技术属于图像处理
，具体涉及一种高效高精度计算刚性点云曲面之间距离的方法。

技术介绍

[0002]随着计算机视觉技术尤其是三维立体视觉技术的快速发展，对三维点云测量的需求日益提升。三维点云的应用领域包括车辆自主导航、工业自动化、无人机空间定位等各种领域。由于点云的采集可能源于不同的传感器，因此他们的参考坐标系各不相同，如何将他们的坐标系进行快速统一，对于后续的目标识别、姿态估计具有重大的指导意义。另外，在实际的更加精密要求的工业生产中有时候需要测量同一目标的外形形貌变化情况。然而目前的测量方法仍存在测量精度低、计算复杂度高、误配准等问题，无法满足高精度的测量要求，为后续的工业生产带来严重的安全隐患，因此急需一种快速的高精度的点云曲面距离计算方法。
[0003]点云曲面的距离测量设计的主要技术有点云配准和点云计算。在点云配准方面，目前常用方式是将两个刚性点云进行配准，常用的点云配准经常ICP迭代进行点云配准，由于ICP迭代以欧式距离作为优化的目标，未考虑点云中点分布的不确定度，因此可能导致配准的精度不够理想。此外，由于测量点云和目标点云并非完全一致，某些位置存在形貌变化，会导致配准精度下降。测量点云存在大量的测量噪声，这也会对配准造成干扰，导致配准陷入局部解，甚至无法收敛。最后，优化方法对于初始化的角度和位置的具有严格的要求，大范围的角度和位置变化会导致配准失败。
[0004]在点云测量方面，现有的方法多适用于同一点云密度，对于目标点云较为稀疏的情况测量精度会下降。另外，点云测量仅依据目标点云和测量点云的最近对应点的距离，测量的精度低，无法满足更高精度的要求。
[0005]中国专利技术专利，公开号为CN107123134A，专利名称为一种基于特征的危岩体滑坡监测方法，公开了一种局部曲面拟合法的点云数据形变量计算方法，采用通用的局部平面拟合方法，形变量确定是通过计算参考点云数据中的点与圆柱体内待测点云数据的重心点之间的距离。
[0006]中国专利技术专利，公开号为CN113160129A，专利名称为组合式精简点云数据的快速配准方法，公开了通过非边界点高斯曲率判断特征区域和非特征区域，后通过非特征区域的PCA粗配准和ICP精配准。
[0007]中国专利技术专利，公开号为CN113205486A，专利名称为基于三维点云的表面缺陷检测方法、装置及存储介质，公开了共面四点的粗配准的改进方法，但并未对ICP精匹配进行优化。
[0008]但上述专利并不能实现高精度计算刚性点云曲面之间距离。

技术实现思路

[0009]为解决现有技术中存在的技术问题，本专利技术的目的在于提供一种高效高精度计算刚性点云曲面之间距离的方法。
[0010]为实现上述目的，达到上述技术效果，本专利技术采用的技术方案为：
[0011]一种高效高精度计算刚性点云曲面之间距离的方法，包括以下步骤：
[0012]S1、点云区域划分
[0013]对点云进行区域划分，得到磨损区域的点云和非磨损区域的点云；对点云进行处理，得到纯净的待测点云；
[0014]S2、点云与参考点云配准
[0015]通过改进的基于点云不确定度的ICP迭代配准将非磨损区域的点云和参考点云进行高精度配准：利用非磨损区域的点云的不确定系数和不确定度，对ICP迭代的点云欧式距离进行加权求和，从而降低噪声点对配准的干扰，提高配准精度；
[0016]S3、通过高精度点云间曲面距离计算，最终获取点云曲面之间的距离。
[0017]步骤S2中，所述改进的基于点云不确定度的ICP迭代配准中，将点云进行降采样2
L
倍，作为迭代的初始值，并将每次进行的降采样2
i
倍的待测点云和参考点云均进行保存，随着迭代的过程，每10*i迭代采用2
L
‑
i
降采样结果作为配准的点云集合，直到恢复实际的点云集合大小，以加速点云配准时间，同时，在每次迭代过程中，剔除置信度低的噪声，提高配准精度。
[0018]进一步的，在改进的基于点云不确定度的ICP迭代配准中，在第10*j次迭代过程中，待测非磨损区域的点云表示为其中，
[0019]设点云的总数为N，则降采样后的点云总数为N/2
(L
‑
j)
；
[0020]非磨损区域的参照点云表示为Q＝{q
i
}，其中，
[0021]进一步的，在改进的基于点云不确定度的ICP迭代配准中，在每次迭代过程中，确定非磨损区域的待配准点云的不确定度，可靠点的置信度高，噪声因置信度低而被剔除，得到点云欧式距离为：
[0022][0023]T
f
＝[R
f
|t
f
]＝arg min E(R
f
|t
f
)
[0024]其中，q
j
是相对应的最近点；∑
k
表示由不确定度组成的协方差矩阵；R
f
、t
f
分别表示每次ICP迭代配准中获取的最佳旋转矩阵和平移矩阵，均通过求解公式E(R|t)的最小值获得。
[0025]进一步的，非磨损区域的待配准点云的不确定度采用以下步骤得到：
[0026]通过最远点采样法对非磨损区域的点云进行点云采样，获得距离点云表面最近的K个关键点q
K
；
[0027]通过Ransac方式从待配准点云上选择任意两个点，并分别计算所述两个点与关键点q
K
的单元向量，所得两个单元向量的交点为点U
K，i
，重复上述步骤N遍，得到不确定点集合{U
K，i
|i＝1，...，N}，进而得到不确定系数：
[0028][0029]其中，q表示非磨损区域的参照点云的点云集合Q的其中一点，ρ表示信任性阈值，不确定点集的单元向量越接近，不确定度的系数越大，置信度越高；
[0030]进而得到关键点q
K
的三维空间分布，得到关键点q
K
的均值u
k
及不确定度
[0031][0032][0033]步骤S3中，将待测点云和参考点云进行点云投影，计算投影后待测点云邻域对应的参照点云的相关点的索引；通过索引值计算邻域参照点云的局部平面法向量，从而得到待测点云到参照曲面距离测量结果；通过移动最小二乘算法和点云重采样去除剩余点云噪声，获取高精度的三维点云曲面的点云。
[0034]步骤S3中，高精度磨损的测量具体包括以下步骤：
[0035](1)获取邻域相关点的索引集：
[0036]将磨损区域待测点云和参照区域点云沿高度方向投影为二维点云，将投影后的待测磨损点作为中心点，以r为半径查找左右参照磨损区域点云的所有点，并储存对应的索引结合；
[0037](2)确定参照点云平面法向量本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种高效高精度计算刚性点云曲面之间距离的方法，其特征在于，包括以下步骤：S1、点云区域划分对点云进行区域划分，得到磨损区域的点云和非磨损区域的点云；对点云进行处理，得到纯净的待测点云；S2、点云与参考点云配准通过改进的基于点云不确定度的ICP迭代配准将非磨损区域的点云和参考点云进行高精度配准：利用非磨损区域的点云的不确定系数和不确定度，对ICP迭代的点云欧式距离进行加权求和，从而降低噪声点对配准的干扰，提高配准精度；S3、通过高精度点云间曲面距离计算，最终获取点云曲面之间的距离。2.根据权利要求1所述的一种高效高精度计算刚性点云曲面之间距离的方法，其特征在于，步骤S2中，所述改进的基于点云不确定度的ICP迭代配准中，将点云进行降采样2
L
倍，作为迭代的初始值，并将每次进行的降采样2
i
倍的待测点云和参考点云均进行保存，随着迭代的过程，每10*i迭代采用2
L
‑
i
降采样结果作为配准的点云集合，直到恢复实际的点云集合大小，以加速点云配准时间，同时，在每次迭代过程中，剔除置信度低的噪声，提高配准精度。3.根据权利要求2所述的一种高效高精度计算刚性点云曲面之间距离的方法，其特征在于，在改进的基于点云不确定度的ICP迭代配准中，在第10*j次迭代过程中，待测非磨损区域的点云表示为其中，设点云的总数为N，则降采样后的点云总数为N/2
(L
‑
j)
；非磨损区域的参照点云表示为Q＝{q
i
}，其中，q
i
∈R3×
(M/2(L
‑
j))
。4.根据权利要求2所述的一种高效高精度计算刚性点云曲面之间距离的方法，其特征在于，在改进的基于点云不确定度的ICP迭代配准中，在每次迭代过程中，确定非磨损区域的待配准点云的不确定度，可靠点的置信度高，噪声因置信度低而被剔除，得到点云欧式距离为：T
f
＝[R
f
|t
f
]＝arg min E(R
f
|t
f
)其中，q
j
是相对应的最近点；Σ
k
表示由不确定度组成的协方差矩阵；R
f
、t
f
分别表示每次ICP迭代配准中获取的最佳旋转矩阵和平移矩阵，均通过求解公式E(R|t)的最小值获得。5.根据权利要求4所述的一种高效高精度计算刚性点云曲面之间距离的方法，其特征在于，非磨损区域的待配准点云的不确定度采用以下步骤得到：通过最远点采样法对非磨损区域的点云进行点云采样，获得距离点云表面最近的K个关键点q
K
；通过Ransac方式从待配准点云上选择任意两个点，并分别计算所述两个点与关键点q
K
的单元向量，所得两个单元向量的...

【专利技术属性】
技术研发人员：李栋，李苏祺，黄磊，沈国锋，陆晓隽，
申请(专利权)人：江苏集萃智能光电系统研究所有限公司，
类型：发明
国别省市：