CN109373901A - 位于平面上的孔的中心位置的计算方法 - Google Patents

Abstract

本发明提供一种位于平面上的孔的中心位置的计算方法，包括如下步骤：1)利用单目视觉传感器采集孔所在平面包含孔区域的灰度图，标记此区域为区域I；2)向区域I分别投射横、纵向的三频四相相移光栅，采集图片，畸变校正，再分别沿横向、纵向解算相机平面上的绝对相位；解算投影仪平面上的横、纵绝对相位；3)取一系列相机平面和投影仪平面所对应的同位置点，计算单应性矩阵；4)提取灰度图中圆孔边缘的亚像素坐标，依据单应性矩阵计算出其在投影仪平面中的坐标，求取圆孔边缘各点的三维坐标，计算圆孔中心在平面中的位置。该方法将三频四相法与单目传感器相结合，能快速、准确的获取孔的中心位置；同时设备体积减小，便于测量。

Description

位于平面上的孔的中心位置的计算方法

技术领域

本发明涉及大型交通工具制造、安装领域，具体涉及一种位于平面上的孔的中心位置的计算方法。

背景技术

孔位置检测在汽车、轮船、飞机等制造领域的自动装配工作中具有非常重要的作用，实现孔位置的准确检测有利于提高在线安装或在线检测工作的准确性，继而提高产品的质量。目前，孔类中心位置的解算方式主要有以下三种：1)三坐标测量机、2)线激光器和3)双目立体视觉。三坐标测量机能够获得较高的测量精度，但其测量速度较慢，难以满足白车身在线测量的要求。线激光的测量方式又可分为以下三种：二维影像结合激光测距检测方案、线激光机器人扫描检测方案、线激光单光刀方案。其中，二维影像结合激光测距检测方案需要额外的LED光源进行打光，且对光照强度、光照角度、摄像机采集角度都有严格的要求；线激光机器人扫描检测方案需要机器人在小范围内进行运动并进行检测，机器人运动绝对精度较低，小范围的轨迹难以精确控制，容易在测量过程中引入机器人运动误差；线激光单光刀方案拍摄图像少，测量速度快，但其拍摄时需使摄像机光轴垂直于待测圆孔所在平面，且只提取两个光刀点进行计算，可用数据较少，稳定性较差。双目立体视觉匹配较为困难，难以获得理想的测量精度。综上所述，现有方法均不能准确、快速的获取孔的中心位置。

发明内容

为了解决上述技术问题，本发明提供一种位于平面上的孔的中心位置的计算方法，其将三频四相法与单目传感器相结合，能快速、准确的获取孔的中心位置；同时设备体积减小，便于测量。

为此，本发明的技术方案如下：

一种位于平面上的孔的中心位置的计算方法，包括如下步骤：

1)利用单目视觉传感器采集孔所在平面包含孔区域的灰度图，标记此区域为区域I；

2)向包含区域I的区域分别投射横向的三频四相相移光栅和纵向的三频四相相移光栅，同时采集图片，对得到的图片进行畸变校正，再分别沿横向、纵向解算相机平面上的绝对相位；

解算投影仪投射出的横向的三频四相相移光栅和纵向的三频四相相移光栅的绝对相位信息，得到投影仪平面上的横、纵绝对相位；

3)取一系列相机平面、投影仪平面横向、纵向绝对相位均相等的点，计算相机平面与投影仪平面之间的单应性矩阵；

4)提取步骤1)得到的灰度图中的圆孔边缘的亚像素坐标，依据步骤3)得到的单应性矩阵计算出其在投影仪平面中的坐标，求取圆孔边缘各点的三维坐标，计算圆孔中心在平面中的位置。

进一步，所述平面为孔上端面的四周至少1cm处与孔的上端面处于同一个平面。

进一步，步骤4)提取步骤1)得到的灰度图中的圆孔边缘的亚像素坐标的步骤为：

①先利用canny边缘检测法提取圆孔的整像素边缘；

②利用亚像素曲面拟合法求解圆孔边缘的亚像素坐标，其中曲面拟合法利用二元二次多项式进行拟合：

G(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

将拟合窗口内的数据代入上式，得到六个系数：a₀,a₁,a₂,a₃,a₄,a₅；

再利用解算出圆孔边缘的亚像素坐标。

进一步，步骤4)中计算圆孔中心在平面中的位置时，先将获取的圆孔边缘的三维坐标点利用最小二乘法，拟合出圆柱，获取圆柱的轴线与圆孔所在平面的交点，即为圆孔中心在平面中的位置。

本发明提供的位于平面上的孔的中心位置的计算方法将三频四相法与单目传感器相结合，通过横向、纵向分别投射的三频四相相移光栅而获取24张图片，计算相机平面与投影仪平面之间的单应性矩阵，将图片上孔的位置信息在投影仪平面上得以反馈，继而计算得到圆孔中心的位置。该方法能快速、准确的获取孔的中心位置，且使用的设备尺寸小，扩大了领域。

附图说明

图1为本发明一个实施例利用相机采集的孔所在区域的灰度图；

图2为图1中相同区域被投射横向三频四相相移光栅后得到的图片；

图3为图1中相同区域被投射纵向三频四相相移光栅后得到的图片；

图4为投影仪平面上的横向绝对相位图；

图5为投影仪平面上的纵向绝对相位图；

图6为相机平面上的横向绝对相位图；

图7为相机平面上的纵向绝对相位图；

图8为采用本发明提供的计算方法对图1采集孔的亚像素边缘坐标进行拟合得到的结果；

图9为本发明一个实施例得到的孔的中心位置的结果展示图。

具体实施方式

以下结合附图对本发明的技术方案进行详细描述。

一种位于平面上的孔的中心位置的计算方法，包括如下步骤：

1)利用单目视觉传感器采集孔所在平面包含孔区域的灰度图，如图1所示，标记此区域为区域I；

2)向包含区域I的区域分别投射横向的三频四相相移光栅和纵向的三频四相相移光栅，同时采集图片(如图2、3)，对得到的图片进行畸变校正，再分别沿横向、纵向解算相机平面上的绝对相位(得到图6、7)；

解算投影仪投射出的横向的三频四相相移光栅(得到图4)和纵向的三频四相相移光栅(得到图5)的绝对相位信息，得到投影仪平面上的横、纵绝对相位；

3)取一系列相机平面、投影仪平面横向、纵向绝对相位均相等的点，计算相机平面与投影仪平面之间的单应性矩阵

设

则可得：

其中[u_p v_p 1]^T为投影仪像平面坐标的齐次表示；[u_c v_c 1]^T为相机像平面坐标的齐次表示；

(1)式展开可得：

将(2)式化简并表示为矩阵形式：

设共有N(N≥4)组点，共可列出2N个方程，故可求解出单应性矩阵求出的应归一化，使得h₃₃＝1；

4)提取步骤1)得到的灰度图中的圆孔边缘的亚像素坐标(如图8)，步骤如下：

①先利用canny边缘检测法提取圆孔的整像素边缘；

②取一个3×3的拟合窗口，将该部分视为一个微小的曲面，该曲面的极值点就为待求的亚像素坐标。采用二元二次函数G(x,y)来拟合该曲面，其表达形式如下：G(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

函数G(x,y)对x,y求偏导，可得：

化简可得：

表1中取了6个点，分别列出了canny边缘检测的坐标与曲面拟合法提取的坐标

表1Canny边缘检测坐标与曲面拟合坐标比较

将拟合窗口内的数据代入上式，得到六个系数：a₀,a₁,a₂,a₃,a₄,a₅；

然后，利用解算出圆孔边缘的亚像素坐标；其次，依据步骤3)得到的单应性矩阵计算出其在投影仪平面中的坐标，根据单目传感器标定所获得的内外参数(内参数包括相机与投影仪的焦距、相机与投影仪的像主点坐标、相机与投影仪的畸变系数；外参数包括相机坐标系至投影仪坐标系的旋转矩阵R与平移矩阵T)，以及圆孔边缘点在相机像平面与投影仪像平面的对应点坐标，即可计算出圆孔边缘在相机坐标系下的三维坐标；解算的公式如下：x_c＝u_c·z_c、y_c＝v_c·z_c、

其中，u_p,u_c,v_c坐标都已归一化并进行了畸变校正；

代表相机坐标系至投影仪坐标系的旋转矩阵；

T＝[t₁ t₂ t₃]^T代表相机坐标系至投影仪坐标系的平移矩阵；

利用如下方法计算圆孔中心在平面中的位置：先将获取的圆孔边缘的三维坐标点利用最小二乘法，拟合出圆柱，要想唯一确定一个圆柱，需要7个参数，分别为该圆柱中心轴线得方向向量(a,b,c)和直线上得某一点坐标(x₀,y₀,z₀)，以及圆柱半径r，由此可得圆柱的中心轴线

其中，t为参数变量；

空间圆柱面得方程可表示为：

其中，(x,y,z)为空间圆柱上得点，(a,b,c)为单位方向向量，化简可得：

列出误差方程，令v＝r²-R²

则可得：

v＝[a(z_i-z₀)-c(x_i-x₀)]²+[b(x_i-x₀)-a(y_i-y₀)]²+[c(y_i-y₀)-b(z_i-z₀)]²-R² (9)

将式(9)进行线性化，并将所求出的圆孔边缘点坐标代入，按最小二乘法即可解算出七个未知参数，继而求出圆柱的轴线所在的直线方程，从而可计算出轴线与圆孔所在平面的交点，获得圆孔在该平面中的位置信息，最终的结果如图9所示。

以上实施例以白车身在线检测中螺孔中心点的检测为例，对该本发明的实际应用进行描述，该计算方法将三频四相法与单目传感器相结合，通过横向、纵向分别投射的三频四相相移光栅而获取24张图片，计算相机平面与投影仪平面之间的单应性矩阵，将图片上孔的位置信息在投影仪平面上得以反馈，继而计算得到圆孔中心的位置。该方法能快速、准确的获取孔的中心位置，且使用的设备尺寸小，扩大了领域。

Claims (4)

1.一种位于平面上的孔的中心位置的计算方法，其特征在于包括如下步骤：

1)利用单目视觉传感器采集孔所在平面包含孔区域的灰度图，标记此区域为区域I；

2)向包含区域I的区域分别投射横向的三频四相相移光栅和纵向的三频四相相移光栅，同时采集图片，对得到的图片进行畸变校正，再分别沿横向、纵向解算相机平面上的绝对相位；

解算投影仪投射出的横向的三频四相相移光栅和纵向的三频四相相移光栅的绝对相位信息，得到投影仪平面上的横、纵绝对相位；

3)取一系列相机平面、投影仪平面横向、纵向绝对相位均相等的点，计算相机平面与投影仪平面之间的单应性矩阵；

4)提取步骤1)得到的灰度图中的圆孔边缘的亚像素坐标，依据步骤3)得到的单应性矩阵计算出其在投影仪平面中的坐标，求取圆孔边缘各点的三维坐标，计算圆孔中心在平面中的位置。

2.如权利要求1所述位于平面上的孔的中心位置的计算方法，其特征在于：所述平面为孔上端面的四周至少1cm处与孔的上端面处于同一个平面。

3.如权利要求1所述位于平面上的孔的中心位置的计算方法，其特征在于：步骤4)提取步骤1)得到的灰度图中的圆孔边缘的亚像素坐标的步骤为：

①先利用canny边缘检测法提取圆孔的整像素边缘；

②利用亚像素曲面拟合法求解圆孔边缘的亚像素坐标，其中曲面拟合法利用二元二次多项式进行拟合：

G(x,y)＝a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²

将拟合窗口内的数据代入上式，得到六个系数：a₀,a₁,a₂,a₃,a₄,a₅；

再利用解算出圆孔边缘的亚像素坐标。

4.如权利要求1所述位于平面上的孔的中心位置的计算方法，其特征在于：步骤4)中计算圆孔中心在平面中的位置时，先将获取的圆孔边缘的三维坐标点利用最小二乘法，拟合出圆柱，获取圆柱的轴线与圆孔所在平面的交点，即为圆孔中心在平面中的位置。