基于多个标志点自动识别的点云数据拼接方法技术

本发明专利技术公开了一种基于多个标志点自动识别的点云数据拼接方法，在实际点云数据拼接中，由于存在测量误差，两块点云数据中的标志点的相对位置不可能完全一致，因此得到的变换矩阵将导致被变换物体变形。为求得最佳的拼接效果，从实际应用的角度出发，针对点云数据巨大散乱的特点，为了提高灵活性和精度，实现点云拼接的自动化，考虑在输入标志点首先控制在一定精度范围的前提下，应取尽可能多的标志点作为标志点群，在前面自动搜索识别算法的基础上，自动对准标志点，实现特征群目标函数最小二乘拼接的算法，通过三维空间中的多个点求解转换矩阵，计算方法简单，拼接精度也比较高。

【技术实现步骤摘要】

本专利技术属于反求工程领域中的点云数据拼合技术，具体涉及到一种基于多坐标点自动识别的点云数据的拼接方法。
技术介绍
研究与实现的是在反求工程中的数据预处理，包括点云的处理，三维空间识别，多视点云的拼接，点云的优化处理等，这一部分的处理决定了三维点云数据的精度与质量。CAD模型建立的前提就是具备完整的高质量三维点云数据，三维点云数据的完整性和精度很大程度上决定了 CAD模型所能达到的数据质量。常规三维检测与转换过程中得到的只是单个的点云数据，完整数据模型需要多幅点云的拼接和处理，其精度也主要依赖于三维点 云的拼接系统这一环节的处理。三维点云拼接与处理不仅仅是承上启下的作用，而且直接决定了三维重构的成功与质量，因此，这一部分的处理在反求工程中是非常关键和重要的一部分。反求工程三维检测中，常需要把物体表面分成多个局部相互重叠的子区域，从多个角度获取零件不同方位的表面信息，以完成对整个物体模型的测量。从各个视觉分块测量得到多个独立的点云，称为多视点云。由于在测量不同的区域时，都是在测量位置对应的局部坐标系下进行的，多次测量所对应的局部坐标系是不一致的，必须把各次测量对应的局部坐标系统一到同一坐标系，并消除两次测量间的重叠部分，以得到被测物体表面的完整数据，此处理过程即为多视点云的拼接。根据获得的多视角点云数据，可以采取直接点集配准的方法，拼接点云数据。多视标签定位点云拼接方法一般包括点位法、固定球法、平面法等。目前这些方法采用的原理都是一种简单、实用的多视点云拼接算法，即基于三基准点的拼合算法。具体做法是根据每个测量位置至少可以看见三个标志点的原则，在待测物体相对平坦的表面以及标定平板上贴一些大小合适的测量专用标签，这样在每块测量数据中，分别测量指定的三个基准点，记为A、B、C。基准点的布置原则是①基准点必须是在两次测量中都可以测到的点指定的基准点既可以是人工设定的，也可以是测量对象固有的，但必须是能够在测量数据点集中准确识别、匹配和定位的点伽、B、C三点不能共线，应避免Λ ABC成为狭长三角形，其包围面积应尽可能大，以减少测量误差对变换的影响。分别以两次测量所得基准点的坐标值建立坐标系，通过坐标平移和旋转变换，就实现了坐标系统一。理论上说，只要知道两块点云数据中对应的3个标志点就可进行拼接，但在实际点云数据拼接中，由于存在测量误差，两块点云数据中的标志点的相对位置不可能完全一致，因此得到的变换矩阵将导致被变换物体变形。目前还有很多方法是采用三点定位组合变换。ICP方法可以多次迭代求取，进一步满足精度的坐标转换。P. Bsel和D. McKay等人在文献“A method for registration of3-D shapes”中提出了 ICP算法(Iterative Closet Point Algorithm),它是首先找到与对象点集中的每一个点在参考点集中的一个距离最近的点，建立一种点对点映射关系，然后使用最小二乘法计算最优坐标变换。但实际应用中，很多情况下都缺乏明确对应关系，因此寻找点集每个点对点的映射关系是比较困难的，使得计算速度跟不上，实际应用起来很不方便。Besl和Mckay 1992年提出了快速最近点迭代ICP算法，迭代地最小化两个点集之间的平均距离，可以较好地解决6个自由度的数据点和模型匹配的问题，不过速度也比较慢，本身不能保证收敛到全局(甚至局部)最小值，并且需要为算法提供一个合理的刚体变换的初始猜测。国内目前也已经在开展点云数据拼接方法的研究，并初步取得一些成果。例如，东南大学达飞鹏等人，在文献“三维扫描系统中点云拼接用标志点匹配方法”中，提出了一种能够对点云拼接标志点进行快速匹配的方法。采用动态分层处理和预估处理。标志点数据实现方式采用线性链表，同时采用数组存储必要的信息，具有明显的计算效率。浙江大学方慧等人，在文献“三维点云数据获取用的基准体及点云合成方法”中，提供了一种利用基准体来确定转台的中心转轴，根据计算得到的中心转轴拼接局部三维点云，并合成为完整的三维点云图像的方法。为了避免在实际点云数据拼接中，由于两块点云数据中的标志点的相对位置不完全一致，致使得到的变换矩阵引起被变换物体变形，并针对实际测量情况，以提高运算速度，达到实用要求为目的，本专利技术提供了一种。
技术实现思路
本专利技术提供一种运算速度较快，拼接精度较高的多视点云拼接时多个标志点自动识别的运算方法。一种，包括如下步骤 步骤I:首先读入两块点云数据，设置目标图与待拼接图，初始化点云拼接数据； 步骤2 :根据标志点的空间几何关系采取特征搜索与识别算法，由程序自动找出对应标志点对，若标志点对大于3对，则将其分别存入矩阵A、B，否则不能进行拼接，如果搜索过程结束后仍然没有解，那么适宜采用递归的方法进行搜索，从外到内开展搜索，在搜索并识别出标志特征点后结束；权利要求1.一种，其特征在于包括如下步骤 步骤I:首先读入两块点云数据，设置目标图与待拼接图，初始化点云拼接数据； 步骤2 :根据标志点的空间几何关系采取特征搜索与识别算法，由程序自动找出对应标志点对，若标志点对大于3对，则将其分别存入矩阵A、B，否则不能进行拼接，如果搜索过程结束后仍然没有解，那么适宜采用递归的方法进行搜索，从外到内开展搜索，在搜索并识别出标志特征点后结束；aIl a2l aM ,^li hi -hi I ^ Φ ±5aSI]β η α2 a^s LPinhn L h hi hi -,。 令^ + + + U可縛线性方程乱AXs=IalCi=LJ) ^ . hJ' .—' 4= ι —ιΓ 令三维变换矩阵T3d= ，并能够使得最小二乘目标函数最小，即T3d= A+ B ； 步骤3 :根据T3d= A+B，求出T3d，其中，； 步骤4 :将一块点云数据经过T3d坐标变换，加入另一块点云数据中，3D坐标变换矩阵1U iIa hi : iW可甬132 ,2J ilM褒示.从变換功缝上货，Ti!}可分为4个子矩 hi hi %3 *_____________Λ__ M % ' iM.阵-.其φ，Lc in i=广主比5IH- ■切几两变ft:L1JI hi 产生平移变换；产生投影变换；产生整体比例变换；由于采用的坐标变换是刚性变换，故有[ 14 tu 34Γ=r, = ，若点P在第一块点云数据中的坐标为A=C^i,J^1,Z1, D，则在第二块点云数据中的坐标为々2= (Ζ2，72，Α，I)，则两者关系可表不为PlT3D=P2 ; 步骤5 :若拼接未完，转为第一步，拼接完成消除重叠点，拼接结束。2.如权利要求I所述的，其特征在于所述步骤2中的标志点为4个或4个以上,标志点对越多,拼接精度越高。全文摘要本专利技术公开了一种，在实际点云数据拼接中，由于存在测量误差，两块点云数据中的标志点的相对位置不可能完全一致，因此得到的变换矩阵将导致被变换物体变形。为求得最佳的拼接效果，从实际应用的角度出发，针对点云数据巨大散乱的特点，为了提高灵活性和精度，实现点云拼接的自动化，考虑在输入标志点首先控制在一定精度范围的前提下，应取尽可能多的标志点作为标志点群，在前面自动搜索识别算法的基础上，自动对准标志点，实现特征群目标函数最小二乘拼接的本文档来自技高网...

【技术保护点】
一种基于多个标志点自动识别的点云数据拼接方法，其特征在于包括如下步骤：步骤1：首先读入两块点云数据，设置目标图与待拼接图，初始化点云拼接数据；步骤2：根据标志点的空间几何关系采取特征搜索与识别算法，由程序自动找出对应标志点对，若标志点对大于3对，则将其分别存入矩阵A、B，否则不能进行拼接，如果搜索过程结束后仍然没有解，那么适宜采用递归的方法进行搜索，从外到内开展搜索，在搜索并识别出标志特征点后结束；令三维变换矩阵？T3D=？[X1？？X2？？X3？？X4]，则有AT3D=B；假设Xi1满足？||？AXi1_bi？||2？=？min？||？AXi1_bi？||2？(i=1,2,3,4)，则称Xi1为方程组AXi=？bi的最小二乘解，此方程组的最小二乘通解为Xi=A+bi+？(I=A+A)Y,Y？Cn；方程组中长度最小的解即为？Xi=A+bi，是唯一解，可称作极小二乘解：其中A+是A的广义逆，将求出的极小二乘解合并可求得T3D=[X1？？X2？？X3？？X4]，并能够使得最小二乘目标函数最小，即T3D=？A+？B；步骤3：根据T3D=？A+B，求出T3D，其中，A+=AT(AAT)？1；步骤4：将一块点云数据经过T3D坐标变换，加入另一块点云数据中，3D坐标变换矩阵[t41？？t42？？t43]产生平移变换；[t14？？t24？？t34]产生投影变换；[t44]产生整体比例变换；由于采用的坐标变换是刚性变换，故有[t14？？t24？？t34]T=[0？？0？？0]T，[t44]=[1？？0]，若点P在第一块点云数据中的坐标为p1=(x1,y1,z1,1)，则在第二块点云数据中的坐标为p2=(x2,y2,z2,1)，则两者关系可表示为：P1T3D=P2？；步骤5：若拼接未完，转为第一步，拼接完成消除重叠点，拼接结束。459dest_path_image001.jpg,931375dest_path_image002.jpg...

【技术特征摘要】

【专利技术属性】
技术研发人员：余东满，史增芳，李晓静，杨峰，苏静，高志华，张玉华，王笛，户燕会，孙育竹，
申请(专利权)人：河南工业职业技术学院，
类型：发明
国别省市：