一种三维点云匹配方法组成比例

本发明专利技术公布了一种三维点云匹配方法。在逆向工程、人脸识别和现代制造业中，解决通过光栅投影到物体表面由摄像机获得物体两组三维点云的对应匹配问题。将两组三维点云输入到目标函数中，将目标函数中的对应矩阵进行归一化处理，再将目标函数中的旋转和平移用对偶数表示形成新的目标函数，然后通过对该新目标函数进行最小化目标，找到点云之间的对应矩阵，即完成三维点云的匹配，使用该方法能够提高三维点云匹配的准确率。

【技术实现步骤摘要】

本专利技术涉及计算机视觉和模式识别领域，特指数字图像处理中三维点云匹配。
技术介绍
三维点云匹配是研制视觉三维测量系统的第一关键技术步骤，点云匹配配准可以为后续的三维重构提供精确的匹配，提高视觉三维测量的精度，用于开发实用的三维测量系统。一方面是用于逆向工程等工业生产的漫反射表面的准静态测量三维测量系统，另一方面是用于人脸等生活中的准静态三维测量系统，还可以开发多种复杂零部件测量专用设备，达到降低成本，提升制造业水平的目标。因此，找到一种三维点云匹配准确率高的方法，成为当前该研究领域急待解决的问题。
技术实现思路
本专利技术的目的在于：提供一种仿射变换下三维点云匹配的新方法，在逆向工程、人脸识别和现代制造业中，解决通过光栅投影到物体表面由摄像机获得物体两组三维点云的对应匹配问题。该专利技术的技术方案如下：通过光栅投影到物体表面，由摄像机拍摄获得物体的两组三维点云{Xj

【技术保护点】
一种三维点云匹配方法，其特征在于：该方法具体为：通过光栅投影到物体表面，由摄像机拍摄获得物体的两组三维点云{Xj}和{Yk}，第一组三维点云{Xj}在坐标系X′Y′Z′下，具体表示为(x′(j),y′(j),z′(j))，第二组三维点云{Yk}在坐标系X″Y″Z″下，具体表示为(x″(k),y″(k),z″(k))，其中j,k均为正整数，由于光栅投影和摄像机拍摄时受到噪声的影响，使得三维点云{Xj}与三维点云{Yk}中对应点的位置发生了变化，即三维点云{Xj}到三维点云{Yk}中点的位置对应关系发生了变化，这个变化是由平移、尺度、旋转和变形产生的；将三维点云{Xj}与三维点云{Yk}输入到下面的目标函数式(1)中：E3D(m,T,R)=Σj=1JΣk=1Kmjk||RXj+T-Yk||2-αΣj=1JΣk=1Kmjk---(1)]]>在式(1)中，E3D(m,T,A)表示关于变量m,T,A的期望值；T＝[tx,ty,tz]i是平移矩阵，tx、ty、tz分别表示在x轴、y轴、z轴上的平移，i表示矩阵的转置，tx、ty、tz的值来源于均匀分布U[2.5,7.5]，R=rx000ry000rz]]>表示旋转，rx、ry、rz分别表示在x轴、y轴、z轴上的旋转，其值来源于均匀分布U[20,70]；α与匹配目标有关的参数，其取值为0.2≤α≤0.5；是一个从三维点云{Xj}到三维点云{Yk}的对应矩阵，矩阵m的行不等式约束为矩阵m的列不等式约束为式(1)中的旋转和平移写成如下形式：RT0001]]>利用对偶数的性质将式(1)中的旋转和平移可以分别记作：R=(r42-(rx,ry,rz)i(rx,ry,rz))I+2(rx,ry,rz)(rx,ry,rz)i+2r4K(rx,ry,rz)]]>其中，I表示单位矩阵，K(R)=0-rzryrz0-rx-ryrx0]]>表示变形；进一步，将旋转和平移写成如下形式：R＝W(r)iQ(r)，T＝W(r)is其中，s=(d/2)cos(θ/2)n+sin(θ/2)(Xj×n)-(d/2)sin(θ/2),]]>r=sin(θ/2)ncos(θ/2),]]>W(r)=r4I-K(R)R-Rir4,]]>Q(r)=r4I+K(R)R-Rir4,]]>d为坐标系X′Y′Z′和X″Y″Z″的原点在世界坐标系下的距离，n为旋转轴的单位方向向量，θ为点{Xj}在沿着向量n的旋转角度；因此，将目标函数式(1)写成下式：E3D(m,T,R)=Σj=1JΣk=1Kmjk||W(r)iQ(r)xj+W(r)is-yk||2-αΣj=1JΣk=1Kmjk---(2)]]>式(2)中，xj＝(Xj,0)i，yk＝(Yk,0)i分别表示式(1)中的点Xj和Yk；利用矩阵变换的性质，目标函数式(2)写成下式：E3D＝riC1r+siC2s+siC3r+λ1(rir‑1)+λ2(sir)  (3)其中：C1=-Σj=1JΣk=1KmjkQ(yk)iW(xj),]]>C2=12Σj=1JΣk=1KmjkI,]]>C3=Σj=1JΣk=1Kmjk(W(xj)-Q(yk)),]]>λ1，λ2为拉格朗日乘数；式(3)中的对应矩阵mjk有三种取值：第一种，mjk＝1，表示三维点云{Xj}中的一个点与三维点云{Yk}中最多只有1个点是完全匹配；第二种，0＜mjk＜1，表示三维点云{Xj}中的一个点与三维点云{Yk}中的点是部分匹配；第三种，mjk＝0，表示三维点云{Xj}中的一个点与三维点云{Yk}中的点是空匹配；目标函数式(3)的对应矩阵m中每个mjk值的过程如下：步骤一：通过对矩阵m进行初始化；步骤二：对目标函数式(3)中的矩阵m进行归一化处理，得到新的矩阵m；步骤三：将步骤二中得到的新的矩阵m代入到目标函数式(3)中，对目标函数式(3)中的参数r、s、λ1和λ2求偏导数，得到r、s、λ1和λ2的新值，并将r、s、λ1和λ2的新值代入到目标函数式(3)中；步骤四：重复进行步骤二～步骤三，直到对应矩阵m中的元素mjk＝1或mjk＝0，即此时对应矩阵m中只有完全匹配和空匹配，已经没有了部分匹配；因此，对应矩阵m完成了两组三维点云的匹配。...

【技术特征摘要】
1.一种三维点云匹配方法，其特征在于：该方法具体为：
通...

【专利技术属性】
技术研发人员：乔付，周波，刘忠艳，车向前，边莉，
申请(专利权)人：乔付，
类型：发明
国别省市：黑龙江;23