一种基于垂直度约束的大物体面形测量的拼接方法与流程

本发明涉及一种光学三维测量方法，属于光电检测技术领域。尤其涉及一种基于垂直度约束的大物体面形测量的拼接方法。

背景技术：

随着科学技术交叉、多融合式的发展，基于光学、图像处理和计算机视觉的三维测量方法因高速、高精度、非接触、全场测量、快速信息获取等优点，已广泛应用于工业制造、文物保护、医疗等诸多领域，在工业生产和日常生活中发挥着重要作用。随着我国生产结构的工业化升级，对大型生产器件三维测量的需求逐渐增加。在航空航天、大型船舶、汽车车身、高速铁路等加工制造环节均离不开快速精密测量来保证每一个零件乃至最终产品的质量。

然而，在对大物体三维面形测量的实际过程中，影响大物体三维测量的速度和精度的主要因素均在于大物体面形数据的拼接技术。由于大型物体几何尺寸的限制，基于面结构光的三维测量技术无法单次完成物体表面整体三维面形的测量，需要进行三维图像拼接融合，在拼接融合过程中会出现迭代计算和累计误差等问题，严重影响了大物体整体面形的三维测量速度和精度。因此，在基于面结构光对大物体三维测量的基础上，研究出一种适用于大物体面形测量数据的拼接方法具有重要的科研意义和现实的经济价值。本发明就是针对大型物体三维测量数据拼接的速度和精度，提出一种基于垂直度约束的大物体面形测量的拼接方法。

技术实现要素：

为解决大型物体三维测量数据拼接的速度和精度的问题，本发明在基于面结构光对大物体三维测量的基础上，提出一种基于垂直度约束的大物体面形测量的拼接方法。其主要分为四个过程：

(1)根据harris角点检测算法检测相机所拍摄的棋盘格特定点的像素坐标值的差值来调整ccd相机光轴与参考面达到垂直。(2)然后利用传统的基于面结构光对大物体面形三维重建的方法得到单幅面形点云数据。(3)需在垂直度标定基础下，运用三标志点所确定平面的法向量与误差校正后的y轴垂直的法向量分量求解出旋转角度进而得到旋转平移矩阵r、t，然后将不同视角的点云数据统一到同一世界坐标系下，完成拼接。(4)运用点距阈值消除法和加权平均融合算法对拼接之后的数据进行无缝融合。

在步骤(1)中，用来调整ccd相机光轴与参考面达到垂直的具体方法为：取视场中央区域某一点，再取离中央点适当距离且近似等距的上下左右四个点构成十字型，当四个点到中央点的距离相等时，则ccd相机光轴与参考面垂直，否则需要继续调整相机方向直到满足条件为止。

在步骤(4)中，增加垂直度标定，以误差校正后的y轴为旋转轴，误差主要是由于参考面不能完全与载物台垂直，标志点所确定平面一般不与旋转轴平行，因此，需要求出与旋转轴垂直的法向量分量。

本发明的优点：(1)本发明运用一种基于垂直度标定的系统使得求解旋转平移矩阵r、t的数据量减少，提高了数据拼接速度。(2)增加了对y轴方向上的旋转轴的误差校正，使得测量大型物体时，拼接精度明显提高。(3)用一种新三标志点法，采用先旋转、后平移的方式，能有效避免基于面结构光对大物体测量时柱面拼接方法的轴心定心不准问题。

附图说明

图1为本发明的三维测量系统示意图。

图2为本发明的标定效果图。

图3为本发明的旋转轴误差校正图。

图4为本发明的旋转角求解原理图。

图中：1.投影仪、2.ccd相机、3.工作站、4.载物台、5.待测物体、6.参考平面。

具体实施方式

以下结合附图说明对本发明的实施例作进一步详细描述，但本实施例并不用于限制本发明，凡是采用本发明的相似结构及其相似变化，均应列入本发明的保护范围。

一种基于垂直度约束的大物体面形测量的拼接方法包括以下步骤：

本发明的三维测量系统俯视图如图1所示，包括投影仪、ccd相机、工作站、测量支架、参考平面、载物台和待测物体；投影仪和ccd相机放在测量支架上；投影仪、ccd相机分别通过数据线与工作站相连接；待测物体放在参考平面前的载物台上；工作站内包含图像采集卡、投影软件、测量软件。投影仪将带有特征信息的条纹聚焦投射到被测物体表面，由ccd相机采集条纹信息，经过工作站处理后提取出特征信息，并按照特定算法进行三维重建。投影仪光轴和ccd相机光轴相交于o点。dlp投影仪和ccd相机为同一高度，它们之间的距离为d，它们到参考平面的距离为l。被测物体某点h的高度计算公式为：

其中f为参考平面上生成的正弦条纹的频率，为物体表面图像和参考平面图像对应点的连续相位差。

(1)ccd相机垂直度标定。标定过程为：在测量之前，紧贴着参考面正放一块棋盘格标定板，将相机采集到的棋盘格图片进行harris检测算法检测棋盘格上每一个角点的像素坐标。为避免ccd相机畸变的影响，取视场中央区域某一点，再取离中央点适当距离且近似等距的上下左右四个点构成十字型，当四个点到中央点的距离相等时，则ccd相机光轴与参考面垂直，否则需要继续调整相机方向直到满足条件为止。标定完成后，需要固定ccd相机的位置。标定效果如图2所示。实验过程中只需稍加调整，标定误差不超过一个像素。

(2)由于对大物体的测量区域比较大，且投影仪采用倾斜投影的方式，因此会出现周期性正弦条纹在参考面上投影不均匀的条纹展宽现象，先采用基于条纹周期校正的条纹四步相移法，再采用基于条纹周期校正的时间相位展开法即可精确地测量出单幅三维大物体的面形数据：xw,yw,zw。

(3)选用标准的圆形标志点，为了使标志点采集清晰，又能保证测量精度，标志点的直径大小尽量为测量区域y轴长度的1/20。三个标志点尽量避免构成狭长三角形且在移动前后的公共视场中，通过霍夫变换算法检测出圆形标志点圆心的二维坐标，再利用被测物体的高度计算公式和相机畸变校正模型确定标志点圆心所对应的世界坐标。物体移动前为第一视角，标志点坐标为：p1(x1,y1,z1)、p2(x2,y2,z2)、p3(x3,y3,z3)；物体移动后为第二视角，标志点坐标为：q1(x4,y4,z4)、q2(x5,y5,z5)、q3(x6,y6,z6)。

(4)令第一视角三标志点所确定平面为p，第二视角三标志点所确定平面为q。则p、q平面的法向量求解如下：

其中则：

同理：

(5)由于p、q平面一般不与载物台垂直，且参考面会在与参考面垂直方向上出现一定偏差角δ，因此要想借助旋转轴y轴来确定旋转角度就必须进行误差校正，如图3所示。由垂直度标定得，z方向的世界坐标全部沿着垂直于参考面的方向，y方向的世界坐标全部沿着平行于参考面方向，而载物台上物体的移动方向却是沿着载物台平面的，因此，旋转轴的误差方向只与y、z轴方向有关，设原旋转轴为则校正之后的旋转轴为因为该偏差角一般较小，所以偏差角δ余弦值可由如下公式求解：

(6)旋转角度求解原理图如图4所示，旋转角度可由p、q平面的法向量在旋转轴方向的垂直分量来确定，则p、q两个法向量与y轴的夹角余弦值分别为：

因此两个法向量与y轴垂直分量分别为：

最后根据这两个与y轴垂直的法向量分量求解出旋转角度的余弦值与正弦值。

由法向量分量得：

为保证视场公共区域的范围，θ一般为锐角，因此：

(7)由上可得，绕误差校正后的y轴旋转的旋转矩阵r为：

令p平面三标志点旋转之后的坐标为：p1′(x1′,y1′,z1′)、p2′(x2′,y2′,z2′)、

p3′(x3′,y3′,z3′)

则平移矩阵为：

其中：pi′＝pi·r(i＝1,2,3)

因此可由公式：

pw′(xw′，yw′，zw′)＝pw(xw，yw，zw)·r+t(10)

完成拼接，其中：pw(xw，yw，zw)为上述第一视角的点云世界坐标。

(8)最后运用点距阈值消除法和加权平均融合算法对拼接之后的数据进行无缝融合。

点距阈值消除法：

其中，(xp,yp,zp),(xq,yq,zq)分别是p视角和q视角的重叠区域某点，重叠区域可由前面拼接可得；ε是设定的某一阈值。当上式成立，则去掉p图(或者q图)上该点。

由于ε阈值的设置始终存在，导致重叠区域的数据融合会出现一定的裂缝，此时需进行裂缝区域的数据加权平均融合：

其中，kz为融合数据，d为裂缝宽度，w宽度中间值，p为裂缝区域内某一点。可达到无缝拼接效果，有效的提高了拼接数据融合的误差。