平面靶标测头中心的两步法标定 下载: 846次
1 引言
基于手持靶标的视觉坐标测量系统因具有体积小、重量轻、便携等优点而得到广泛的研究与应用[1-4]。该系统使用具有一系列特征靶点且各靶点位置关系已知的手持靶标,通过其前端的可拆卸测头接触被测目标表面,摄像机对靶标上点亮的靶点进行成像,分析靶点与测头的位置关系求解出接触点的三维坐标,实现接触和非接触测量的结合。为实现高精度的三维坐标测量,必须对手持靶标视觉坐标测量系统的关键参数进行标定[5-7]。对于已有机械加工精度保证的手持靶标,其各靶点的相对位置关系在测量过程中不会改变,可通过离线方式,并借助特殊的测量设备来获取这些参数[8]。但是为了满足现场不同的测量环境以及保护测头的要求,对于测头的拆卸和更换是不可避免的。即使是使用同一测头,由于安装力度的不同,也无法保证其位置不变。因此,如何快速精确地对测头中心进行现场标定是亟待解决的问题。
文献[ 2]利用靶标标记点与测球球心之间的距离信息提出旋转轴求交法,借助两组图像信息通过奇异值分解(SVD)和最小二乘法求得球心坐标,整个算法过程分为三步,且测头中心坐标在x、y、z方向的稳定性为(0.178 mm,0.188 mm,0.221 mm),并不能很好地满足工业要求。文献[ 8]根据10点3列的立体靶标测量模型和位置不变原理,通过非线性方程组最小二乘解对目标函数进行优化求解,求取测头中心在靶标坐标系下的坐标,该方法可取得较高的坐标精度,可达(0.033 mm,0.030 mm,0.043 mm),但难以标定近似平面的靶标。当靶标相对于摄像机成像面以不同倾角进行测量时,平面靶标测头中心的z向误差对测量精度有影响,可以将文献[ 8]的方法改进为平面靶标标定方法,但无法直接区分测头中心轴向坐标变化及测量距离变化,即无法精确实现z向标定。
考虑到平面靶标视觉坐标测量的广泛研究与应用[9-15],本文提出一种适用于平面靶标测头中心两步标定的方法。该方法基于平面点阵测头测量数学模型和位置不变原理,建立以测头中心为坐标原点的靶标坐标系。标定过程分为两步,其基本思想是通过两组不同姿态的靶标图像建立靶标坐标系和摄像机坐标系之间的转换关系,先标定出靶标坐标系下测头中心x、y方向最优化的修正值,再利用靶标相对于成像面倾斜不同角度的图像信息求出z方向的误差值,进而利用单摄像机和平面靶标现场准确地标定测头中心。
2 测头中心标定方法
2.1 平面靶标视觉坐标测量原理
手持平面靶标视觉坐标测量系统主要由手持平面靶标、摄像机、计算机处理单元以及被测目标组成,平面靶标上至少有两个发光二极管(LED)靶点不在同一直线上,示意图如
2.2 标定模型及测头中心参数求解
如
对于每一位置γ、β对应的图像k,设各靶点在靶标坐标系下的坐标为(xki,yki,zki)及其对应的经过畸变修正的图像像面坐标为(Uki,Vki)时,其中i=1,2,…,j,j为靶标上靶点的个数,则依据透视变换成像原理有:
式中f为摄像机有效焦距,Tk为两坐标系间的平移矩阵,其物理意义是靶标坐标系原点P与摄像机坐标系原点OC之间的平移位置关系,则Tk即为测头中心在摄像机坐标系下的坐标。在平面靶标各靶点共面情况下,其z向坐标zki为0,设
变量代换得到未知量alk的线性方程组为:
当2j≥8时,则可解出未知量alk的最小二乘解,再由Rk矩阵的正交约束性,可得:
即对于图像k可解出唯一的测头中心在摄像机坐标系下的坐标Tk。
2.3 标定的两步法
为了能够准确标定测头中心,先获取一组靶标绕z轴旋转不同角度的图像,再获取一组绕y轴旋转不同角度的靶标图像,通过相关图像处理和标定方法实现测头中心的标定。标定的两步法具体过程如下:
第一步:测头中心x、y方向的标定。
1) 标定过程中,摄像机固定,摄像机与靶标的相对距离一定,靶标靶点可以清晰成像。当倾角β近似等于0时,通过多次改变倾角γ的大小,即获取一组靶标平行于摄像机成像面的图像,并进行图像处理得到图像靶点中心坐标(Uki,Vki)。
2) 先假定测头中心不存在误差,根据(1)~(4)式可得该组第k幅图像对应的测头中心坐标的理论值(Txk,Tyk,Tzk),则取初始值X0=(0,0,Tx1,Ty1,Tz1)。
3) 实际测头中心可能存在误差,设测头中心x、y方向在靶标坐标系下的修正值为(Δx,Δy),待求测头中心在摄像机坐标系下的坐标为(x0,y0,z0),则各靶点在靶标坐标系下的坐标应修正为(xki+Δx,yki+Δy,0),记X=(Δx,Δy,x0,y0,z0),第k幅图像求解的测头中心在摄像机坐标系下实际与理论位置的偏差为:
建立以该组所有靶标图像对应的测头中心在摄像机坐标系下的坐标误差平方和为最小的目标函数:
以X0作为X的优化初值,采用Levenberg-Marquardt优化算法解出X,即靶标坐标系下测头中心最优化的修正值(Δx,Δy)和在摄像机坐标系下的测头中心坐标(x0,y0,z0)。
4) 依据所提方法原理,得到修正值(Δx,Δy)后,则各靶点在靶标坐标系下坐标修正为(xki+Δx,yki+Δy,0),可按照(2)~(4)式求出该组每幅图像对应的摄像机坐标系下实际测头中心坐标(T'x,T'y,T'z)。根据位置不变原理,坐标数据应具有较好的一致性,因此以该组测头中心坐标数据的稳定性来衡量所提方法的精度。
第二步:测头中心z方向的标定。
1) 在第一步的标定系统环境基础上,令倾角γ近似等于0,通过多次改变倾角β的大小,得到一组靶标相对成像面的图像,并进行图像处理得到图像靶点中心坐标(Uki,Vki)。
2) 由第一步可知此时各靶点在靶标坐标系下的坐标对应修正为(xki+Δx,yki+Δy,0),代入(2)~(4)式求出该组每幅图像对应的摄像机坐标系下实际测头中心坐标(T″x,T″y,T″z)。当测头中心z方向不存在误差时,则该组数据应具有较好的一致性。
3) 若测头中心z方向存在误差Δz,依据所建立的坐标系关系,随着倾角β的改变,Δz将会影响测头中心在摄像机坐标系下的坐标T″y,即:
利用不同倾角下坐标T″y的偏差便可得到z方向的误差值Δz,即:
最终可以得到标定后靶标坐标系下的测头中心坐标。
3 标定方法的仿真验证
为了验证所提方法的可行性,在2.3节标定两步法的基础上利用Matlab对其进行仿真分析,仿真验证采用的摄像机参数和靶点参数与标定实际实验的参数基本一致。设接触点P在摄像机坐标系下的坐标为平移矩阵T=(0,150,1500),通过设置不同β、γ(靶标竖直时γ=β=0),得到摄像机相对于靶标的空间姿态,即旋转矩阵R,利用透视变换成像原理求出像面坐标系下靶点的理想图像坐标,详细推导可参考文献[ 5]。考虑图像处理误差,给理想图像坐标增加符合正态分布的随机噪声,噪声标准差为σ1=0.05 pixel,得到含噪声的图像坐标(Uki,Vki)。假设靶标坐标系下测头中心x、y、z方向存在误差且(Δx,Δy,Δz)=(1 mm,1 mm,-1 mm)。
靶标以不同方位测量同一点,设置倾角γ在0°~40°范围内随机变化,由仿真可得在摄像机坐标系下测头中心坐标多次测量的结果。
图 3. 测头中心标定前与修正后坐标T'x和T'y的仿真结果。(a) T'x;(b) T'y
Fig. 3. Simulation results of coordinate T'x and T'y of probe tip center before and after calibration. (a) T'x; (b) T'y
使用文献[
8]扩展到平面靶标情形的方法对测头中心进行20次仿真标定,可得每组图像对应的测头中心x、y、z方向的修正值变化结果如
图 4. 测头中心修正值Δx、Δy和Δz的仿真结果
Fig. 4. Simulation results of corrected value Δx, Δy and Δz of probe tip center
设倾角β为0,给倾角β增加符合正态分布的随机噪声,噪声标准差为2°,按所提方法第一步对测头中心进行20次仿真标定,可得第一步中每组图像对应的测头中心x、y方向的修正值变化结果如
图 5. 测头中心修正值Δx和Δy的仿真结果
Fig. 5. Simulation results of corrected values Δx and Δy of probe tip center
再进行第二步的标定仿真。当倾角β从0°到30°变化时,给靶标各倾角增加符合正态分布的随机噪声,噪声标准差为2°,运行仿真程序可得测头中心在摄像机坐标系下的坐标T″y与倾角β的关系如
图 6. 测头中心坐标T″y与靶标倾角β的关系
Fig. 6. Relationship between coordinate T″y of probe tip center and angle tilt β of target
4 标定方法的实验验证
实验采用JAI GO-5000M-USB工业摄像机,摄像机分辨率为2560 pixel×2048 pixel,像元大小为5 μm,摄像机标定后焦距为52.5 mm。平面靶标安有8个亮度自适应的红外LED,其中有两个LED靶点不在同一直线上,示意图参考
装在靶标上,采用的测头的中心在摄像机坐标系下的坐标参考值为(-20 mm,0,-1 mm)。靶标坐标系下各靶点的坐标如
表 1. 靶标坐标系下各靶点坐标
Table 1. Coordinates of targets under target coordinate systemmm
|
实验在满足靶标靶点全部在视场内情况下,取γ∈(-30°,30°),β∈(0°,30°)。先使靶标在xPy平面内旋转,β近似为0,倾角γ从图像最左侧旋转11次到最右侧变化,得到第一组的11幅图像;再令γ近似为0,靶标在xPz平面内从β近似为0旋转11次到β近似为30°的位置,得到第二组的11幅图像,对两组图像进行图像处理得到靶点中心图像坐标。借助第一组图像信息,分别使用文献[
8]扩展到平面靶标情形的改进方法和本研究所提方法对测头中心进行标定,测头标定后坐标计算结果如
表 2. 第一组图像测头中心标定的结果
Table 2. Results of probe tip center calibration with the first set of imagesmm
|
借助第一组图像,通过所提方法第一步可以求得测头中心x、y方向的修正值Δx和Δy分别为19.271,0.234 mm,进而求得摄像机坐标系下测头中心实际坐标如
借助第二组的图像信息,求出每张图像对应摄像机坐标系下实际的测头中心坐标和倾角β,
式中σβ表示角度β的标准差。假设角度β的标准差为2°,代入已知数据可得σz≈0.140 mm,即测头中心z方向的重复性精度为0.140 mm。
综上所述,标定后靶标坐标系下的测头中心坐标为(-19.271 mm,-0.234 mm,-0.967 mm),在1900 mm成像距离和靶标相对成像面倾斜最大30°的情况下,测头中心在x、y、z方向的重复性精度可达0.077,0.035,0.140 mm,验证了所提标定方法的有效性和测头中心位置的准确性。理论上x、y方向精度一致性较好,这里是因为实验系统采用的平面靶标宽度分布范围较小所以影响了x方向的精度。
5 结论
提出了测头中心两步标定算法,该方法较为简便,适用于不同类型的平面靶标测头,不仅可以高精度地标定测头中心x、y方向,而且可以利用靶标相对于摄像机成像面倾斜来标定测头中心z方向,现场测量仅需借助两组不同位姿的靶标图像信息即可实现测头中心的有效标定。实验结果表明,利用本方法,测头中心在x、y、z方向的重复性精度分别可达0.077,0.035,0.140 mm,有效地解决了手持平面靶标测头中心现场标定的难题。
[3] 陈洋, 林嘉睿, 高扬, 等. 视觉与倾角传感器组合相对位姿测量方法[J]. 光学学报, 2015, 35(12): 1212003.
[4] 马国鹭, 曾国英. 基于无衍射探针的空间坐标测量及其不确定度分析[J]. 光学学报, 2015, 35(5): 0512005.
[5] 赵敏, 郑兴纯, 黄秋红, 等. 深度空间小视角摄像机标定[J]. 激光与光电子学进展, 2017, 54(9): 091502.
[6] 苏建东, 齐晓慧, 段修生. 基于单目视觉和棋盘靶标的平面姿态测量方法[J]. 光学学报, 2017, 37(8): 0815002.
[7] 黄风山, 刘书桂, 彭凯, 等. 光笔视觉测量系统结构参数的自标定与仿真[J]. 光电子·激光, 2006, 17(6): 705-708.
[9] 徐巧玉, 车仁生. 基于光学测棒的立体视觉坐标测量系统的研究[J]. 光学学报, 2008, 28(11): 2181-2186.
[10] 尹英杰, 徐德, 张正涛, 等. 基于单目视觉的平面测量[J]. 电子测量与仪器学报, 2013, 27(4): 347-352.
[11] 冯萍, 魏振忠. 光笔式大视场三维视觉测量系统[J]. 光学精密工程, 2013, 21(9): 2217-2224.
[12] 解则晓, 孙洪磊, 王晓敏. 光笔式双摄像机三维坐标视觉测量系统[J]. 光学技术, 2012, 38(4): 459-464.
[13] 张元元, 张丽艳, 杨博文. 基于双目立体视觉的无线柔性坐标测量系统[J]. 仪器仪表学报, 2010, 31(7): 1613-1619.
[14] 于之靖, 孙大钊, 陶洪伟. 基于编码测棒的单目视觉测量系统[J]. 半导体光电, 2014, 35(6): 1080-1083.
郑兴纯, 赵敏, 冯少华. 平面靶标测头中心的两步法标定[J]. 激光与光电子学进展, 2018, 55(1): 011201. Zheng Xingchun, Zhao Min, Feng Shaohua. Two-Step Calibration of Probe Tip Center of Planar Target[J]. Laser & Optoelectronics Progress, 2018, 55(1): 011201.