1.一种基于局部变换一致的非刚性点集配准方法,其特征在于:通过计算机,并按如下步骤进行:步骤1,向计算机输入由三维点云获取设备获取的源点云和目标点云,获取这两个点云的关键点集:源关键点集S和目标关键点集O;具体如下:步骤1.1,点云获取:获取两个不同视角的点云,利用三维几何处理软件将这两个点云保存为只包含三维坐标信息的点云,得到源点云和目标点云;
步骤1.2,关键点提取:根据步骤1.1得到的源点云和目标点云,在源点云数据范围内建立一个网格,所述网格由三维体素组成,通过将源点云包围在三维体素网格中,每个体素中包含源点云的三维点,三维点在一个以上,为简化点的数量,则在每个三维体素中选一个点,这样在该三维体素内所有点就可以用这一个点表示,完成了对点云数据的简化;对网格中所有体素处理后得到的点云即为源点云对应的源关键点集S;同理,目标点云对应一个目标关键点集O步骤2,由步骤1中获取的两个关键点集,转换获取关键点集之间的对应关系;
步骤3,建立源关键点集S的邻域索引矩阵;邻域索引矩阵的获取方法为:采用被广泛应用的kd-tree最近邻搜索算法,将源关键点集S中的每一个点作为查询点,检索在kd-tree树中与查询点距离最近的K个相邻点,计算这K个相邻点在源关键点集S中的位置索引,根据索引建立源关键点集S的邻域索引矩阵Idy=[idy1,···,idyi,···,idyM]T;M表示的是源关键点集S的点的个数;
步骤4,根据源关键点集S及其邻域索引矩阵,定义源关键点集S中每个点的空间变换与其K个相邻点的空间变换的差异程度的函数式,即约束项;该约束项由源关键点集S中每个点的空间变换分别与其K个相邻点的空间变换的差值之和组成,用来约束相邻点之间的变换保持一致;数据项表示的是源关键点集S和目标关键点集O中每一对对应点之间的距离;约束项与数据项组成目标函数式的作用是最小化目标函数式,获取此时对应的空间变换矩阵,并将此空间变换矩阵作用到源关键点集S,获取新的源关键点集S';
步骤5,设定迭代的参数值,迭代执行步骤2到步骤4,当达到参数值时停止迭代,将此时的源关键点集S'和目标关键点集O作为配准结果并输出。
2.根据权利要求1所述的一种基于局部变换一致的非刚性点集配准方法,其特征在于:在步骤1中,获取两个不同视角的点云,对这些点云进行关键点提取,获取各自的关键点集,即源关键点集S和目标关键点集O;
在步骤2中,对步骤1中获取的两个关键点集,计算源关键点集S中的每个点在目标关键点集O中的对应点,获取这两个关键点集之间的对应关系;
在步骤3中,对源关键点集S中的每一个点进行K近邻搜索,即对关键点集S中每一个点,在整个源关键点集S中找到与其欧几里得距离最近的K个点,这K个点便称为它的相邻点,并根据这些相邻点在源关键点集S中的位置信息,建立源关键点集S的邻域索引矩阵;
在步骤4中,利用相邻点的空间变换一致,对非刚性变换进行局部约束,即对源关键点集S中的每一个点的变换都进行局部约束,计算每一个点的空间变换分别与其K个相邻点的空间变换的差异,在保证差异较小的情况下获得较优的空间变换矩阵,以此获得变换后的源关键点集S';
在步骤5中,设定一个最大迭代次数和一个参数阈值,迭代执行步骤2到步骤4,参数值超过阈值或者迭代次数超过最大迭代次数,停止迭代,此时经空间变换矩阵变换后的源关键点集S'及目标关键点集O作为配准结果并输出。
3.根据权利要求1或2所述的一种基于局部变换一致的非刚性点集配准方法,其特征在于:详细步骤如下:在步骤1中,点云获取与提取点云关键点:获取仅包含三维坐标信息的两个点云,即源点云和目标点云,分别对这两个点云采用降采样的方式提取关键点,得到源关键点集S和目标关键点集O;所述的降采样是指对点云进行稀疏化处理,源关键点集S和目标关键点集O分别对应的是源点云和目标点云进行降采样后的数量不大于1万个点的两个点集,这两个点集分别保留了源点云和目标点云的形状特征及空间结构信息;
在步骤2中,计算源关键点集S与目标关键点集O之间的对应关系:在步骤1获取关键点集的基础上,对源关键点集S中的每一个点,计算它与目标关键点集O中的每一个点的欧式距离,依据对应的欧式距离确定两个关键点集中每个点与点之间的对应概率;
在步骤3中,建立源关键点集S的邻域索引矩阵:对步骤1获取的源关键点集S中的每一个点,在整个源关键点集S范围内计算与其欧式距离最近的K个相邻点,并记录这些相邻点在源关键点集S中的位置索引,根据已记录的索引建立源关键点集S的邻域索引矩阵;
在步骤5中,迭代执行步骤2到步骤4:设置一个最大迭代次数和参数σ2的阈值,每循环执行一次步骤2到步骤4,便重新计算一下参数σ2的值,迭代次数加一;σ2表示每个高斯分量的协方差;σ表示标准差;
当参数值大于阈值或者迭代次数达到最大迭代次数,停止迭代,此时将空间变换矩阵作用在源关键点集S上,得到新的源关键点集S',将其和目标关键点集O,作为配准结果并输出;
当参数值不大于阈值、迭代次数未达到最大迭代次数时,说明对应点之间的距离并没有接近于0,没有实现将两个点集统一到一个坐标系下,即没有完成两个关键点集的配准,返回步骤2。
4.根据权利要求1或2所述的一种基于局部变换一致的非刚性点集配准方法,其特征在于:步骤2具体如下:步骤2.1,在步骤1.2得到的源关键点集S和目标关键点集O的基础上,计算源关键点集S中的每个点与目标关键点集O中每个点之间的欧几里得距离:假设ym为源关键点集S中的一个点,xn为目标关键点集O中的一个点,则两点之间的距离公式如下:d(xn,ym)=||xn-ym||2 (2-1)
步骤2.2,计算高斯混合模型概率密度函数:将源关键点集SM×D=(y1,···,yM)T表示为高斯混合模型的质心,将目标关键点集ON×D=(x1,···,xN)T作为由高斯混合模型生成的数据点,高斯混合模型简称GMM,在步骤2.1得到的点与点之间的欧几里得距离基础上,建立高斯混合模型概率密度函数公式如下:
其中,M和N分别表示源关键点集S和目标关键点集O的点的个数,m表示高斯混合模型的第m个高斯分量,n表示目标关键点集O中点的下标,D表示两个关键点集的维数,取值为3,exp表示以自然常数e为底的指数函数;σ2表示每个高斯分量的协方差,初始值为:
假设GMM的所有分量都是独立同分布的,则联合高斯混合模型概率密度函数公式如下:
其中,M表示的是源关键点集S的点的个数,N表示的是目标关键点集O的点的个数,P(m)=1/M表示每个高斯分量的隶属概率;
步骤2.3,计算源关键点集S中的每个点与目标关键点集O中每个点之间的对应关系:基于高斯混合模型的点集配准方法,当源关键点集S和目标关键点集O对齐时,对于一个已知点xn,其与点ym之间的对应关系利用最大化高斯混合模型的后验概率获得;在步骤2.2的基础上,计算高斯混合模型的后验概率,获得两个点集之间的对应概率;高斯混合模型的后验概率公式如下:P(m|xn)=P(m)p(xn|m)/p(xn) (2-5)P(m)表示每个高斯分量的隶属概率;m表示高斯混合模型的第m个高斯分量;xn为已知点;p(xn|m)高斯混合模型概率密度函数公式。