大数据量、少特征的产品点云数据快速全局优化配准方法与流程

本发明属于机器视觉技术领域，特别涉及一种大数据量、少特征的产品点云数据快速全局优化配准方法。

背景技术：

三维点云配准是计算机视觉中的一个热点问题，所谓的点云配准，即给定不同坐标系中的两组点集，或者是相同坐标系下的两个不同姿态的点集，将这两组点集通过旋转、平移转换到同一个坐标系下或是将两组点集重叠在一起。三维点云配准在机器视觉领域应用非常广泛，包括3d建模、对象识别、姿态估计、机器人导航等。点云配准主要分为两类：一是刚体变换的配准，二是非刚体变换的配准，本发明针对的是刚体变换的配准。由于数据采集设备的视场限制、物体自身的遮挡现象等原因，单次测量只能采集到物体的部分点云，需要通过不断配准才能获得所需的完整模型，在配准过程中，存在单次测量点云与大部分点云或完整点云(已知cad模型情况下)的配准问题，这种部分对整体的配准存在数据量大、初值位置偏差大、全局优化配准速度慢的问题。

点云配准主要流程一般可以分成两个阶段，点云粗配准和点云精配准。

点云粗配准阶段，通常情况下都是通过求解点云或者关键点的局部形状几何信息，比如点与点之间的法向量夹角、点云的曲率、点之间的欧几里得距离、点云深度值等等，然后利用这些信息来构成点云的局部特征描述子，利用局部特征描述子进行配准，因此局部特征描述子在配准中起着关键的作用。一个好的特征描述子应该是具有高度描述性的，以便提供全面、明确的局部几何形状。为了确保精确和高效的点云配准，特征描述子还应该具有计算效率，紧凑性和鲁棒性。目前很多学者提出了许多局部特征描述子，包括快速点特征直方图描述子(fpfh)、点签名方向直方图描述子(shot)，旋转投影统计描述子(rops)，局部特征统计图(lfsh)等。这些都提供了相当详细的局部形状几何描述。因此，它们具有高维特征向量，准确性高但计算量大。相对而言曲率、法向等信息描述特征具有低描述性但易于计算。

点云精配准阶段，目前最著名的方法是迭代最近点(icp)算法，这种算法速度较慢且对点云的初始位置要求较高，所以一般都是用于点云的精配准。很多学者对icp算法进行了改进，jost和hugli利用多分辨率来对icp算法进行改进，提高了icp算法的速度；dai等人采用两个新的结构约束对icp算法进行改进。这样，不可靠的匹配点对就会被移除，高质量的点对将会留下来用来计算转换关系；等人通过嵌套bnb算法和icp算法实现转换矩阵的全局收敛，但是该方法速度慢，而且仅对旋转矩阵进行全局优化。

技术实现要素：

本发明的目的在于克服现有技术的云配准过程存在数据量大、初始位置不精确、全局配准速度慢的问题，提供一种大数据量、少特征的产品点云数据快速全局优化配准方法。

本发明解决其技术问题所采用的技术方案是：

一种大数据量、少特征的产品点云数据快速全局优化配准方法，包括：

基于超体素聚类算法对滤波后的cad模型点云进行重叠区域的提取；对每个点云块的质心点提取相对应的fpfh特征描述子，利用这些点云特征描述子进行粗配准；

基于iss算法得到的关键点并结合全局优化算法实现精配准。

优选的，基于超体素聚类算法对滤波后的cad模型点云进行重叠区域的提取；对每个点云块的质心点提取相对应的fpfh特征描述子，利用这些点云特征描述子进行粗配准，具体包括：

a1、基于超体素聚类算法实现点云的重叠区域提取，即对cad模型点云提取与测量点云大致重叠的区域；包括：基于超体素聚类算法将点云分块，并提取出各分块点云的质心点；求出cad模型点云质心点及测量点云质心点的主曲率、查询质心点的法向量与查询质心点到其领域点重心连线的夹角，并利用这两个元素作为重叠区域相似性提取判别准则，得到重叠区域对应的质心点，利用质心点索引将重叠区域点云提取出来；将cad模型重叠区域质心点集和测量点云质心点集加以保存，作为点云粗配准的输入点云；

a2、利用质心点集进行粗配准；包括：计算cad模型重叠区域质心点集和测量点云质心点集的快速点特征fpfh描述子，将其作为点云的点特征；基于fpfh特征描述子建立匹配点对，应用随机采样一致性ransac算法去除错误匹配点对；计算对应点对的协方差矩阵并进行奇异值分解svd，获得初始转换矩阵，完成点云的粗配准。

优选的，所述基于超体素聚类算法实现点云的重叠区域提取，具体包括：

a1-1、首先利用超体素聚类算法将测量点云与cad模型点云进行分块，将种子点均匀分布在点云空间clouddiagonal中，然后对体素的距离即voxel的分辨率rvoxel和种子点的分辨率rseed进行定义，采用aabb包围盒求得点云的最小包围盒，再利用最小包围盒的对角线长的一半作为体素分辨率和种子点分辨率的因子，实现自适应分辨率提取：

rvoxel＝0.01×clouddiagonal+0.8(1)

rseed＝0.037×clouddiagonal+3.0(2)

式(1)为体素分辨率的取值公式，式(2)为种子点分辨率的取值公式；

a1-2、对超体素聚类分块后的各点云块提取其质心点，获得cad模型的各点云块的质心点集p{m1，m2...mn}，以及测量点云的质心点集q{n1，n2...nn}；

a1-3、对a1-2中提取到的质心点进行判定筛选，实现重叠区域的提取；评判准则为质心点集p、q中每个质心点的主曲率和质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角；

a1-4、对a1-3中的质心点集p、q中每个质心点的主曲率进行计算，采用离散曲率估算法计算主曲率，包括最小曲率和最大曲率，采用基于局部最小二乘的法估计向量；

a1-5、对a1-3中的质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角进行求解；首先计算k邻域的重心点centerpσ(i)和centerqσ(j)，求得查询点质心点到其重心的连线||centerpσ(i)-pi||和||centerqσ(j)-qj||；令np(i)和nq(j)分别为点集p、q各查询点的法向量，则它们之间的夹角分别为α＝accos＜np(i)，||centerpσ(i)-pσ(i)||＞、β＝accos＜nq(j)，||centerqσ(j)-pqσ(j)||＞；

a1-6、质心点集p、q中每个质心点的最大主曲率分别记为pik1和qjk1，最小主曲率分别记为pik2和qjk2；质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角记为α和β；判别准则不等式为：||pik1-qjk1||＜ε1、||pik2-qjk2||＜ε2、||α-β||＜ε3，其中ε1、ε2、ε3取值分别为0.03、0.03、0.9。

优选的，步骤a2中利用质心点而非全部点云集进行粗配准。

优选的，所述基于iss算法得到的关键点并结合全局优化算法实现精配准，具体包括：

a3、采用iss关键点提取算法提取滤波后的cad模型点云及测量点云的关键点，对提取到的关键点采用分支定界法bnb嵌套icp算法进行点云的精配准，即可得到全局优化的转换矩阵。

优选的，步骤a3中利用iss提取的关键点作为匹配对象，分支定界法bnb嵌套icp算法进行点云的精配准。

本发明与现有相关技术相比，具有以下优点：

(1)本发明对部分对整体的点云配准，提出了一种基于超体素聚类分块并提取重叠区域的方法，并直接利用重叠区域中提取到的每块点云的质心点作为点云粗配准数据源，大大减少了计算量，有效提高了粗配准的效率。

(2)本发明对重叠区域提取过程中的参数进行了自适应操作，提高了自动化程度，避免了由于人工输入参数而消耗的时间；

(3)本发明利用iss关键点提取算法提取关键点，提出基于关键点进行全局优化配准(bnb嵌套icp算法)，在避免局部优化的基础上，提高了配准的速度和精度。

以下结合附图及实施例对本发明作进一步详细说明；但本发明的一种大数据量、少特征的产品点云数据快速全局优化配准方法不局限于实施例。

附图说明

图1是本发明的流程图；

图2是原始cad模型与测量点云图；其中图2(a)为原始cad模型，图2(b)为测量点云图；

图3是本发明实施例中重叠区域提取的效果图；

图4是本发明实施例中粗配准的效果图；

图5是本发明实施例中精配准的效果图。

具体实施方式

参见图1所示，本发明提出了一种大数据量、少特征的产品点云快速全局优化配准方法，包括如下步骤：

a1、基于超体素聚类算法实现点云的重叠区域提取，即对cad模型点云提取与测量点云大致重叠的区域。首先基于超体素聚类算法将点云分块，并提取出各分块点云的质心点；求出cad模型点云质心点及测量点云质心点的主曲率、查询质心点的法向量与查询质心点到其领域点重心连线的夹角，并利用这两个元素作为重叠区域相似性提取判别准则，得到重叠区域对应的质心点，利用质心点索引将重叠区域点云提取出来；将cad模型重叠区域质心点集和测量点云质心点集加以保存，作为点云粗配准的输入点云；

a2、利用质心点集进行粗配准。计算cad模型重叠区域质心点集和测量点云质心点集的快速点特征(fpfh)描述子，将其作为点云的点特征；基于fpfh特征描述子建立匹配点对，应用随机采样一致性(ransac)算法去除错误匹配点对；计算对应点对的协方差矩阵并进行奇异值分解(svd)，获得初始转换矩阵，完成点云的粗配准；

a3、采用iss关键点提取算法提取滤波后的cad模型点云及测量点云的关键点，对提取到的关键点采用分支定界法(bnb)嵌套icp算法进行点云的精配准，即可得到全局优化的转换矩阵。

本实施例中，步骤a1之前还包括：对cad模型点云及测量点云进行统计滤波，可以去除异常点及噪声点，消除噪声点对点云配准的影响。

本实施例中，所述的a1中基于超体素聚类算法实现点云的重叠区域提取的步骤为：

a1-1、首先利用超体素聚类算法将测量点云与cad模型点云进行分块。将种子点均匀分布在点云空间(clouddiagonal)中，然后对体素的距离即voxel的分辨率(rvoxel)和种子点的分辨率(rseed)进行定义，采用aabb包围盒求得点云的最小包围盒，再利用最小包围盒的对角线长的一半作为体素分辨率和种子点分辨率的因子，实现自适应分辨率提取：

rvoxel＝0.01×clouddiagonal+0.8(1)

rseed＝0.037×clouddiagonal+3.0(2)

式(1)为体素分辨率的取值公式，式(2)为种子点分辨率的取值公式。

a1-2、对超体素聚类分块后的各点云块提取其质心点，获得cad模型的各点云块的质心点集p{m1，m2...mn}以及测量点云的质心点集q{n1，n2...nn}；

a1-3、对a1-2中提取到的质心点进行判定筛选，实现重叠区域的提取。评判准则为质心点集p、q中每个质心点的主曲率和质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角；

a1-4、对a1-3中的质心点集p、q中每个质心点的主曲率进行计算，采用离散曲率估算法计算主曲率，包括最小曲率和最大曲率，采用基于局部最小二乘的法估计向量。

a1-5、对a1-3中的质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角进行求解。首先计算k邻域的重心点centerpσ(i)，centerσ(j)，求得查询点质心点到其重心的连线||centerpσ(i)-pi||，||centerqσ(j)-qj||。令np(i)nq(j)为点集p、q各查询点的法向量。则它们之间的夹角为α＝accos＜ηp(i)，||centerpσ(i)-pσ(i)||＞、β＝accos＜nq(j)，||centerqσ(j)-pqσ(j)||>。

a1-6、质心点集p、q中每个质心点的最大主曲率分别记为pik1和qjk1，最小主曲率分别记为pik2和qjk2；质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心连线之间的夹角记为α和β；判别准则不等式为：||pik1-qjk1||＜ε1、||pik2-qjk2||＜ε2、||α-β||＜ε3，其中ε1、ε2、ε3取值分别为0.03、0.03、0.9。

所述的a2中利用质心点而非全部点云集进行粗配准。

所述的a3中利用iss提取的关键点作为匹配对象，分支定界法(bnb)嵌套icp算法进行点云的精配准。

如下将以具体的cad模型点云与测量点云为输入进行配准，具体包括：

步骤一：首先输入cad模型点云与测量点云。参见图2所示，cad模型点云包含约298599个点，测量点云包含约145813个点；再分别对cad模型点云与测量点云进行统计滤波。

步骤二：对滤波后的点云基于超体素聚类算法实现点云的重叠区域提取。首先基于超体素聚类算法将点云分块，接着提取出各分块点云的质心点，然后通过求得cad模型点云质心点及测量点云质心点的主曲率以及查询质心点的法向量与查询质心点到其领域点重心连线的夹角。利用这两个元素作为重叠区域相似性提取判别准则，得到重叠区域对应的质心点，然后利用质心点索引将重叠区域点云提取出来，效果如图3所示。最后将cad模型重叠区域质心点集和测量点云质心点集加以保存，并作为点云粗配准的输入点云；

所述的点云质心点的主曲率以及查询质心点的法向量与查询质心点到其领域点重心连线的夹角这两个因素的计算方法为：

(1)主曲率的提取采用的是离散曲率估算法中的最常用的解析法，主曲率包括最小曲率和最大曲率，它们代表着法曲率的极值。要求得主曲率值，就得先求得质心点的法向量。通常情况下点云数据都伴随着噪声，采用基于局部最小二乘的法向量估计。

(2)计算质心点集p、q中每个质心点的法向量与查询质心点到k领域的重心的连线之间的夹角。首先计算k邻域的重心点centerpσ(i)，centerqσ(j)，求得查询点质心点到其重心的连线||centerpσ(i)-pi||，||centerqσ(j)-qj||。令np(i)nq(j)为点集p、q各查询点的法向量。则它们之间的夹角为α＝accos＜np(i)，||centerpσ(i)-pσ(i)||＞、β＝accos＜nq(j)，||centerqσ(j)-pqσ(j)||>。

(3)通过以上(1)、(2)步骤求得的点云质心点的主曲率以及查询质心点的法向量与查询质心点到其领域点重心的连线的夹角来建立提取相似重叠区域的判别不等式：

||pik1-qjk1||＜ε1

||pik2-qjk2||＜ε2

||α-β||＜ε3

其中，pik1和qjk1为两点云的最大主曲率，pik2和qjk2为两点云的最小主曲率。其中ε1、ε2、ε3取值分别为0.03、0.03、0.9。

步骤三：利用质心点集进行粗配准。计算cad模型重叠区域质心点集和测量点云质心点集的快速点特征(fpfh)描述子，利用fpfh描述子作为点云的点特征。

步骤四：基于fpfh特征描述子建立匹配点对，然后再应用随机采样一致性(ransac)算法去除错误匹配点对，最终得到精确的匹配点对；

步骤五：计算对应点对的协方差矩阵并进行奇异值分解(svd)从而得到初始转换矩阵。至此，完成了点云的粗配准。效果如图4所示；

步骤六：采用iss关键点提取算法提取滤波后的cad模型点云及测量点云的关键点，对提取到的关键点采用分支定界法(bnb)嵌套icp算法进行点云的精配准，可以得到全局优化的转换矩阵。

对转换后的测量点云与cad模型点云的位置进行分析，可以得到均方差为0.652mm。整个的配准时间为53.258s。

按照文献的方法得到的均方差为3.457mm。整体时间为80.233s。

由上述实施例1的实验数据可知，本发明提出的大数据量少特征的点云快速全局优化配准方法从实质上减少了数据量，从而提高了配准的速度。

以上仅为本发明实例中一个较佳的实施方案。但是，本发明并不限于上述实施方案，凡按本发明所做的任何均等变化和修饰，所产生的功能作用未超出本方案的范围时，均属于本发明的保护范围。