一种三维点云的快速配准方法技术

一种三维点云的快速配准方法，包括：获得三维参考点云A和待配准点云B；建立点云A的投影坐标系，并对其进行透视投影，得到二维点云A’；基于三角剖分算法，将点云A’剖分成若干个三角形；在所述投影坐标系中，对点云B进行透视投影，得到二维点云B’；判断点云B’中的点是否落入所述三角形中，以此判别得到若干个落入点o；计算出落入点o与其所在三角形三个顶点之间的插值系数；根据插值系数和三角插值算法，计算出落入点o在点云A中的对应点点q的坐标值，以得到点q与落入点o在点云B中的对应点点p所构成的点对<qn，pn>，n≥1；计算出点对<qn，pn>的最小刚体变换值T；基于变换值T，更新待配准点云B中各点的坐标值。

【技术实现步骤摘要】
一种三维点云的快速配准方法
本专利技术涉及计算机视觉与模式识别
，具体来说，涉及一种三维点云的快速配准方法。
技术介绍
三维图像采集设备一次只能得到物体一个侧面的三维数据，为了得到物体的整体三维数据，需要从多个角度获取三维数据，并对这些三维数据进行配准。常规的点云配准过程为：先通过采样，提取不同视角的三维数据，然后，基于一些算法，如迭代最近点（ICP）算法等，通过计算，建立起点云间点与点之间的联系，最后，基于这些联系，通过对某一初始位置的三维数据，进行逐层迭代，最终建立起完整的三维图像模型。传统的ICP算法是通过将待配准点云中的点逐个与参考点云中的每一个点进行距离长短的比较，以寻找最近点，从而建立起点云间点与点的联系，这种计算方式计算速度太慢，效率太低，而且也存在较大的误差匹配的问题，因此，之后，许多研究者在其基础之上，提出了很多改进方式，例如，《ImageVisionCompute》1992年第10期第3卷发表的"Objectmodellingbyregistrationofmultiplerangeimages"一文中，提出运用法矢方向的距离来代替点到点距离作为匹配的评价函数，这样两个视图中的点无需一一对应，但这种方法要求解非线性最小二乘问题，计算效率仍然较低；又如，《InternationalJournalofComputerVision》1994年第12期第13卷发表的"Iterativepointmatchingforregistrationoffree-formcurvesandsurfaces"一文中，提出一种通过修改KD树来加快搜索最近点的方法，然而，建立类似KD树之类的空间检索结构需要消耗额外的空间和时间，故只能在事先的预处理中进行，因此，仅适用于静态三维点云，而不适合动态三维点云采集过程；再如，CN103150747A中公开了一种基于拓扑特征的点云配准算法，其通过提取点云的拓扑特征点(边界特征点和突出特征点)作为配准元素选择，建立点云间的联系，然而，同样存在配准效率低，无法满足对连续获得的三维点云进行快速、即时配准的要求。
技术实现思路
为了解决上述问题，本专利技术提出一种三维点云的快速配准方法，其可以实现三维点云的动态实时、快速配准，配准效率高。实现本专利技术所提供的配准方法的第一种具体技术方案如下：一种三维点云的快速配准方法，包括以下步骤：（1）获得三维参考点云A和待配准点云B；（2）建立参考点云A的投影坐标系，将参考点云A投影到二维平面上，得到二维点云A’；（3）基于Delaunay三角剖分算法，将点云A’剖分成若干个三角形，形成Delaunay三角网；（4）在参考点云A的投影坐标系中，将待配准点云B投影到所述二维平面上，得到二维点云B’；（5）判断点云B’中的点是否落入所述三角形中，以此判别得到若干个落入点o；（6）基于三角插值算法，计算出落入点o与其所在三角形三个顶点之间的插值系数；（7）利用所述插值系数，对所述三个顶点所对应的参考点云A中的三个点进行坐标值的插值计算，得到所述落入点o在参考点云A中的对应点点q的坐标值；所述点q与所述落入点o在待配准点云B中的对应点点p构成点对<qn，pn>，n≥1；（8）计算出点对<qn，pn>的最小刚体变换值T；（9）基于所述最小刚体变换值T，更新待配准点云B中各点的坐标值，以使待配准点云B与参考点云A相配准。其中，所述步骤（6）中，计算插值系数时，所采用的计算公式如下：u*x1+v*x2+w*x3=x0u*y1+v*y2+w*y3=y0u+v+w=1式中，(x0,y0)为所述落入点o的坐标值；(x1,y1)、(x2,y2)和(x3,y3)分别为所述落入点o所在三角形三个顶点的坐标值；u、v、w为所述落入点o依次对应(x1,y1)、(x2,y2)、(x3,y3)三个顶点的插值系数。所述步骤（7）中，插值计算点q的坐标值时，所采用的计算公式如下：u*X1+v*X2+w*X3=Xu*Y1+v*Y2+w*Y3=Yu*Z1+v*Z2+w*Z3=Z式中，(X1,Y1,Z1)、(X2,Y2,Z2)和(X3,Y3,Z3)依次为所述(x1,y1)、(x2,y2)、(x3,y3)三个顶点所对应的参考点云A中的三个点的坐标值；(X,Y,Z)为所述落入点o在参考点云A中的对应点点q的坐标值。此外，所述参考点云A和待配准点云B优选为相互间有重叠部分的三维点云数据，并且，在通过诸如三维扫描仪等采集设备，获取所述参考点云A和待配准点云B的过程中，可以优选采用动态实时扫描采集的方式进行。所述步骤（2）和步骤（4）中，相应投影生成点云A’和点云B’的方式，优选采用透视投影方式。实现本专利技术所提供的配准方法的第二种具体技术方案如下：一种三维点云的快速配准方法，其特征在于，包括以下步骤：（1）获得三维参考点云A和待配准点云B；所述参考点云A中的点均带有同侧法向方向值。（2）建立参考点云A的投影坐标系，将参考点云A投影到二维平面上，得到二维点云A’；（3）基于Delaunay三角剖分算法，将点云A’剖分成若干个三角形，形成Delaunay三角网；（4）在参考点云A的投影坐标系中，将待配准点云B投影到所述二维平面上，得到二维点云B’；（5）判断点云B’中的点是否落入所述三角形中，以此判别得到若干个落入点o；（6）基于三角插值算法，计算出落入点o与其所在三角形三个顶点之间的插值系数；（7）利用所述插值系数，分别对所述三个顶点所对应的参考点云A中的三个点进行坐标值和法向方向值的插值计算，得到所述落入点o在参考点云A中的对应点点q的坐标值和法向方向值；（8）基于所述点q的坐标值和法向方向值，确定出所述落入点o在参考点云A中对应点点q所处的平面S，以最终得到所述落入点o在待配准点云B中的对应点点p与所述平面S所构成点-面对<pn，Sn>，n≥1；（9）计算出点-面对<pn，Sn>的最小刚体变换值T；（10）基于所述最小刚体变换值T，更新待配准点云B中各点的坐标值，以使待配准点云B与参考点云A相配准。其中，所述步骤（6）中，计算插值系数时，所采用的计算公式如下：u*x1+v*x2+w*x3=x0u*y1+v*y2+w*y3=y0u+v+w=1式中，(x0,y0)为所述落入点o的坐标值；(x1,y1)、(x2,y2)和(x3,y3)分别为所述落入点o所在三角形三个顶点的坐标值；u、v、w为所述落入点o依次对应(x1,y1)、(x2,y2)、(x3,y3)三个顶点的插值系数，均为常数值。所述步骤（7）中，插值计算点q的坐标值时，所采用的计算公式如下：u*X1+v*X2+w*X3=Xu*Y1+v*Y2+w*Y3=Yu*Z1+v*Z2+w*Z3=Z式中，(X1,Y1,Z1)、(X2,Y2,Z2)和(X3,Y3,Z3)依次为所述(x1,y1)、(x2,y2)、(x3,y3)三个顶点所对应的参考点云A中的三个点的坐标值；(X,Y,Z)为所述落入点o在参考点云A中的对应点点q的坐标值。以及，插值计算点q的法向方向值时，所采用的计算公式如下：式中，、和依次为所述(x1,y1)、(x2,y2)、(x3,y3)三个顶点所对应的参考点云A中的三个点的法向方本文档来自技高网...

【技术保护点】
一种三维点云的快速配准方法，其特征在于，包括以下步骤：（1）获得三维参考点云A和待配准点云B；（2）建立参考点云A的投影坐标系，将参考点云A投影到二维平面上，得到二维点云A’；（3）基于Delaunay三角剖分算法，将点云A’剖分成若干个三角形，形成Delaunay三角网；（4）在参考点云A的投影坐标系中，将待配准点云B投影到所述二维平面上，得到二维点云B’；（5）判断点云B’中的点是否落入所述三角形中，以此判别得到若干个落入点o；（6）基于三角插值算法，计算出落入点o与其所在三角形三个顶点之间的插值系数；（7）利用所述插值系数，对所述三个顶点所对应的参考点云A中的三个点进行坐标值的插值计算，得到所述落入点o在参考点云A中的对应点点q的坐标值；所述点q与所述落入点o在待配准点云B中的对应点点p构成点对<qn，pn>，n≥1；（8）计算出点对<qn，pn>的最小刚体变换值T；（9）基于所述最小刚体变换值T，更新待配准点云B中各点的坐标值，以使待配准点云B与参考点云A相配准。

【技术特征摘要】
1.一种三维点云的快速配准方法，其特征在于，包括以下步骤：（1）获得三维参考点云A和待配准点云B；（2）建立参考点云A的投影坐标系，将参考点云A投影到二维平面上，得到二维点云A’；（3）基于Delaunay三角剖分算法，将点云A’剖分成若干个三角形，形成Delaunay三角网；（4）在参考点云A的投影坐标系中，将待配准点云B投影到所述二维平面上，得到二维点云B’；（5）判断点云B’中的点是否落入所述三角形中，以此判别得到若干个落入点o；（6）基于三角插值算法，计算出落入点o与其所在三角形三个顶点之间的插值系数；（7）利用所述插值系数，对所述三个顶点所对应的参考点云A中的三个点进行坐标值的插值计算，得到所述落入点o在参考点云A中的对应点点q的坐标值；所述点q与所述落入点o在待配准点云B中的对应点点p构成点对<qn，pn>，n≥1；（8）计算出点对<qn，pn>的最小刚体变换值T；（9）基于所述最小刚体变换值T，更新待配准点云B中各点的坐标值，以使待配准点云B与参考点云A相配准。2.如权利要求1所述的一种三维点云的快速配准方法，其特征在于，所述步骤（6）中，计算插值系数时，所采用的计算公式如下：u*x1+v*x2+w*x3=x0u*y1+v*y2+w*y3=y0u+v+w=1式中，(x0,y0)为所述落入点o的坐标值；(x1,y1)、(x2,y2)和(x3,y3)分别为所述落入点o所在三角形三个顶点的坐标值；u、v、w为所述落入点o依次对应(x1,y1)、(x2,y2)、(x3,y3)三个顶点的插值系数。3.如权利要求2所述的一种三维点云的快速配准方法，其特征在于，所述步骤（7）中，插值计算点q的坐标值时，所采用的计算公式如下：u*X1+v*X2+w*X3=Xu*Y1+v*Y2+w*Y3=Yu*Z1+v*Z2+w*Z3=Z式中，(X1,Y1,Z1)、(X2,Y2,Z2)和(X3,Y3,Z3)依次为所述(x1,y1)、(x2,y2)、(x3,y3)三个顶点所对应的参考点云A中的三个点的坐标值；(X,Y,Z)为所述落入点o在参考点云A中的对应点点q的坐标值。4.如权利要求1所述的一种三维点云的快速配准方法，其特征在于，所述参考点云A和待配准点云B为相互间有重叠部分的三维点云数据。5.一种三维点云的快速配准方法，其特征在于，包括以下步骤：（1）获得三维参考点云A和待配准点云B；所述参考点云A中的点均带有同侧法向方向值；（2）建立参考点云A的投影坐标系，将参考点云A投影到二维平面上，得到二维点云A’；（3）基于Delauna...

【专利技术属性】
技术研发人员：刘晖，
申请(专利权)人：长沙中科院文化创意与科技产业研究院，刘晖，
类型：发明
国别省市：湖南;43