第31卷第1期 2 0 1 4年2月 沈阳航空航天大学学报 Journal of Shenyang Aerospace University VO1.31 No.1 Feb.2 0 1 4 文章编号:2095—1248(2014)O1—0072—05 基于直线成像特征的视觉测量系统综合标定方法 王静 ,蔡清华 ,任维 ,赵萍 (1.中航(沈阳)高新科技有限公司研发中心,沈阳110141;2.辽宁省交通高等专科学校机械电子工程系,沈阳110122) 摘要:针对视觉测量系统存在复杂系统误差的特点,提出一种基于直线成像特征的系统综合标定 方法,有效地保证了测量精度。该方法通过对1等量块的平行直线边缘提取,建立量块边缘空间 物点位置和像点位置的原始对应关系,利用量块边缘的理想直线特征,通过统计计算对测量系统 进行综合标定,建立描述空间物点位置和像点位置的相互对应关系的二元三次标定函数。对在测 量范围内不同方位的量块进行多次测量,表明使用这种标定方法标定的视觉测量系统的测量精度 能够达到。 关键词:视觉测量;摄像机标定;正交多项式;边缘检测 中图分类号:TP391.4 文献标志码:A doi:10.3969/j.issn.2095—1248.2014.01.015 Vision measurement system comprehensive calibration based on linear imaging characteristics WANG Jing ,CAI Qing—hua ,REN Wei ,ZHAO Ping (1.R&D Center,Avic(Shenyang)High Technology Co.,LTD.,Shenyang 110141; 2.Department of Mechanical and Electronic Engineering,Liaoning Provincial College of Communications,Shenyang 1 10122) Abstract:Aimed at the features of various system errors in vision measuring system,the paper proposed a comprehensive system calibration method based on linear imaging characteristics with measurement accura— cy,According to the extraction of parallel linear edges of 1 grade gauge blocks,the original correspondence between the space point positions and the image point positions of the edges is built.Using the ideal linear features of gauge block edges,comprehensive calibration of measurement system is carried out through staffs— tical calculation,and a two-parameter second—order polynomial which describes the corresponding relation between the space point positions and the image point positions is established.Adopting the general engineer- ing gauge block as the calibration tool,the calibration method is simple to operate and is easy to realize and has good versatility.It can comprehensively correct all kinds of system errors such as the optical distortion, he perspectitve error,the sensor position error and location error of edge detection algorithm.Bythe calibra— tion method,multiple measurement experiments of different azimuth gauge blocks in the field show that the precision of the measurement system can reach. Key words:vision measurement;camera calibration;orthogonal polynomials;edge detection 在机器视觉测量中,为了从图像中定量提取、 二维视觉测量系统,在理想情况下,实际物体与物 体的成像之间只相差一个放大倍数,但是在测量 过程中,存在多种系统误差,如镜头畸变、由于像 测量空间物体的几何信息,必须建立图像中像点 位置和空间物点位置的相互对应关系…。对于 收稿日期:2013—07—16 基金项目:国家863计划课题(项目编号:2012AA041310) 作者简介:王静(1970一),女,吉林长春人,工程师,主要研究方向:复杂曲面测量、制造及数字化化装配等,E-mail:jingwang@ avicgxkj.com。 第1期 王静,等:基于直线成像特征的视觉测量系统综合标定方法 73 机光轴与测量平面不垂直引起图像透视剪切变 形,边缘检测算法的定位误差使得测量点位置偏 移。因而需要对系统进行标定,建立像点位置和 空间物点位置的相互对应关系,消除或尽可能小 地减小系统误差对测量精度的影响。 人们主要通过对摄像机参数的标定来实现对 视觉测量系统的标定 J,这种标定方法大体可 以分为两类:传统标定方法和自标定方法。传统 标定方法是利用一个标准参照物的一些特征点 平行于图像坐标系 轴的直线作为测量坐标系的 x轴,将该点且平行于图像坐标系Y轴的直线作 为测量坐标系的l,轴,坐标方向与图像坐标系相 同。将一组量块依次按照如图2所示的位置定位 装夹到i贝0量平台上,提取量块测量面的边缘,并根 据测量范围选取采样点(图中标志为・的点)。 (或特征线)和它们的像素坐标,然后通过非线性 优化的方法求取摄像机模型的内部参数和外部参 数。自标定方法不需要标定参照物,它是利用摄 像机在运动过程中周围环境的图像及图像问的对 应关系对摄像机进行标定,但是该类方法求解不 稳定,且标定精度受摄像机运动精度的影响¨]。 总的来说,通过对摄像机的标定的精度达不到微 米级,主要原因是:①误差模型只考虑了镜头畸变 误差,而不是整个系统的误差;②目前常用的标定 参照物,如黑白方块阵列和正交网格制作精度很 难达到0.1级;③参照物与被测物的成像条件 (成像物距、光照条件、拍摄参数)不一致。天津 大学的孙双花设计了可用于背光条件的精密网格 靶标,采用分区线性标定方法对二维视觉系统进 行标定,标定的最大误差为5.1l9 J,但受靶标的限 制,该方法很难作为一种通用方法。 为了实现高精度的标定,本文提出了一种基 于直线成像特征的系统综合标定方法,利用通用 的1等量块的平行直线边缘对测量系统进行综合 标定。建立了描述像点位置和空问物点位置的相 互对应关系的标定函数,并利用量块测量面成像 的直线特征确定标定函数的系数。 1标定原理 用于视觉测量系统标定的参照物一般具有以 下三个基本条件:①标定参照物精度高于待标系 统至少一个数量等级;②参照物特征点的坐标容 易求取;③参照物与被测物的成像条件(成像物 距、光照条件、拍摄参数)一致。 按JJG2056—90《长度计量器具(量块部分) 检定系统》规范,1 1/llTI~10 mm的1等量块测量 面的测量精度可以达到0.05,且边缘简单,有利 于高精度提取边缘,使用背光源条件采集量块的 图像如图1所示,图像对比度高,边缘锐利,因此 本文使用1等量块对系统进行标定。设定与图像 中心对应的点作为测量坐标系原点,将过该点且 图1量块图像 在理想情况下,直线成像后仍然是直线,因此 量块成像后两个测量面之间的像素距离所对应的 空间距离等于量块的标称尺寸。这样对于图2 (a),采样点的x坐标是确定的,对于图2(b),采 样点的】,坐标是确定的。如果采用2 rnnl、4 mm、 6 1TI1TI、8 rnnl、10 mm这五个量块标定,则采样点 的x坐标从左到右依次为{一10,一8,一6,一4, 一2,2,4,6,8,10},l,坐标从上到下依次为{一10, 一8,一6,一4,一2,2,4,6,8,10},单位为mm。根 据标定函数,可以分别建立x、y坐标与图像坐标 ( ,Y)的关系。 1.1标定函数的建立 根据中心透视投影的成像原理,在理想情况 下,测量平面上的点(X,Y)与其对应理想像点 (X ,Y )的关系为: fX= ( 一c )一Mx(X。一c )=Mx(X 一 。)… 【y= (Y 一c )一My(X。一c )=My(Y 一Y。) 式中,(C ,C )为光心的坐标,(X。,Y。)为测量 坐标原点在图像坐标系中的坐标。 、 分别 为测量系统x、y向的尺寸当量,单位为mm/像 素。 但实际情况是,由于受各种系统误差影响,物 体的图像不能反映其实际的形态。图3是对 7 mm量块和直径为19.951 l/iar的圆样板(精密 磨削)的边缘提取结果,可见直线边缘和圆边缘 都有较大的变形,经计算量块的直线度为2.935 像素,圆样板的圆度误差为1.732像素,在当前测 量条件下,尺寸当量为0.01375 mm/像素,可见都 超过了实际值。 为此,本文综合考虑系统中存在的各类误差, 建立理想像点(X ,Y )与实际像点( ,Y)之间的关 系为: 74 沈阳航空航天大学学报 第31卷 ~ ———————+X I , 一一● l一● I—I 一 l一● ●一I l—I I一● l—I I—I I一● }一I 一 l—I l—I I—I ●一● l一● r '1 r—— 一 X l—I I—I ●一I I—I ●一I I一1 I—I I一一 l一● l一一 I—I l一1 l一1 l--q I-- V y (a) 矬标的标定 (b)y坐标的标定 图2标定原理 x/pixel x/pixel (a)量块的亚像素边缘 图3亚像素边缘提取结果 f c=ao+al +a2Y+a3x +a4xy+asY + (b)圆样板的亚像素边缘 定位块的A面上,如图4所示。 / /{I Y bc:=a。 b06x+3 + b1 a7+x +2y2Y+ +a +bsx3Yx:2 b+ +a94yx3y+ +bsY +c 2 6 + 7 Y+b8xY +b9y 基准线 /一 量块 fX=A0+Al +A2Y+A3 +A4 +A5Y + l L , 、 y l A6 +A7 Y+A8xy +A9Y l y=曰0+B1 +B2Y+B3 +B4xy+B5Y + L B6X +B7 y+日8j +B9Y 式中,Ao=Mx(a0一 ),A =Mxa。,B0= 图4量块与定位块的关系 (bo—Yo),B =M b ,其中i=1,2,…,9。 式(3)描述了像点位置与已消除了误差影响 的物点位置的相互对应关系,利用该式就可根据 物体的成像求得其尺寸和位置。 =但在实际操作中,定位块A面的成像边缘很 难正好调整到定位基准线(x向标定时,基准线为 1 224,Y向标定时,基准线为Y=1 025)上,此 外量块直接靠在A面上也存在位置偏差,可能会 出现如图5所示的情况。为了解决这个问题,本 文利用坐标平移和旋转来校正量块定位时的偏移 和倾斜。 2.2量块的定位 对于本文标定方法,量块的定位精度很大程 度上影响了系统标定的精度。为此,本文先用定 位块对量块进行粗定位,把定位块的A面作为定 位面,反复调整定位块,使定位面A的成像边缘 位于基准线的位置,然后锁紧,标定时将量块靠在 提取作为定位基准的测量面边缘,将边缘点 进行线性拟合,获得直线Y=ax+b,根据a、b确 定量块偏离定位基准线的截距b和夹角0。这样 第1期 王静,等:基于直线成像特征的视觉测量系统综合标定方法 fX=CO+ClX+c2),+c3(・ 一2)+c4 + 75 量块上采样点的坐标( ,Y)修正为 =1223 作为定位基准 的测量面边缘 fy _2)+c6 一 )+c ( _2)y+ { yI。一2)x+c9 一 ) y: + ),+ +以(), 一2)+以 + ‘5 / v:1025 / df5(x2_2)+d6(Y。一 ) ),。-2) + l d8(x2—2)y+a9( ) 对于上式,分别用最小二乘法求取拟合系数C 、 d ,即得 图5量块偏离定位基准线 f =口cos +(),一易)sin 【Y =一asin0+(Y—b)cos0 负号。 (4) cf ∑Pi( ,y)X f=1 — ————~ 式中,0=±arctana,a>0时,取正号,a<0时,取 ∑P ( ,y) =1 (6) 由于通过定位块的粗定位,量块偏离定位基 l【 ∑Pi( ,),) —~ 准线的程度较小,因此采样点上携带的误差信息 不会改变太大,不会对标定精度产生较大的影响。 I d = ∑p ( ,y) 2.3标定函数的拟合 式(3)所描述的标定函数是两个关于 、Y的 二元三次多项式,由于采样点的个数m远大于待 定系数的个数,因此用最小二乘法来求这两个多 项式的系数。为了避免解方程组时出现病态系数 矩阵,先用正交多项式组作基底来重构多项式。 式(3)可以写成 3实验及分析 本文选用德国basler公司的piA2400—17gc 全帧型面阵CCD黑白像机,日本Computar公司 的M5018一MP2手动光圈镜头以及basler piA2400—17gc像机自带的网络接口(Gigabit Jf x=∑ciP ,),) 。 N(4) 【l,=∑diPj( ,y) 式中,C 、d 是拟合系数,N=9,P 是基底函数。 采用2个一维Chebyshev多项式正交多项式 Ethernet)数字图像采集卡搭建视觉测量系统,利 用本文提出的标定方法对该系统进行标定。 为了验证标定后系统的测量精度,使用该系 统测量如图6所示的5 mm和8 rain量块的尺寸 和直线度,截取图像中量块测量面边缘的中段部 集的张量积来生成P ,则式(4)可写为 分进行测量,其在测量平面内的位置如图6所示。 fa)5 mm量块 (b)5 mm量块 (c)8 mtil量块 (d)8 mm量块 图6测试量块图像 将图7中各量块基准边缘上的测量点进行最 小二乘线性拟合,得到评定基线L ,计算测量边缘上各测量点到基准边缘的距离d ,采用统计的 方法计算,量块的尺寸D为 D: .d ‘ (7) 76 沈阳航空航天大学学报 第31卷 其中n为测点的个数。尺寸误差△为计算 尺寸与量块的标称尺寸之差。 将测量边缘上的测量点进行最小二乘线性拟 合,得到评定基线 ,计算该边缘上各测量点到 的距离 ,,则量块的直线度误差d 为 d = 一 i (8) 按照上述方法,计算的尺寸和直线度误差如 表1所示 _ 1 / 2 / 0 —5 0 5 】 X|mm (a)5 mm量块 图7量块中段部分的位置 表1测试量块的尺寸与直线度(单位nlin) 尺寸 直线度 图像—— D △ d 图7(a)5.00043 0.00043 n 0o3139一n0O329Q00I蚴 图7(b)5.IXf2462 0.009-462 0.0o1778一 0( 6n0o6038 图7(C)8.004339 n 0043339 n哪!8—0.0o149 0.便I2208 图7(d)8.IXr2586 n( ̄Y2586 n0o4818—0.0038 0.08618 由表1可见,对于按照统计方法计算的尺寸 和直线度可以达到±5 m。 4结论 基于机器视觉的测量系统是一个复杂的系 统,存在多种系统性误差:如光学畸变误差、透视 误差、边缘检测算法的定位误差等。因而需要对 系统进行标定来消除或尽量减小这些误差对测量 的影响。为了实现高精度的标定,本文提出并实 现了一种基于直线成像特征的系统综合标定方 法。该方法通过对1等量块的平行直线边缘提 取,建立量块边缘空间物点位置和像点位置的原 始对应关系,利用量块边缘的理想直线特征,通过 统计计算对测量系统进行综合标定,建立描述空 间物点位置和像点位置的相互对应关系的二元三 次标定函数。这种标定方法操作简便,容易实现, 通用性好,可以综合修正各项系统性误差。通过 对量块的多次测试证明,采用这种标定方法,视觉 测量系统的测量精度能够达到±5 m。 参考文献(Referenees): [1]于起峰,尚洋.摄像测量学原理与应用研究[M].北 京:科学出版社,2007. [2]郭羽,杨红,杨照金,等.CCD摄像系统镜头的畸变测 量[J].应用光学,2008,29(2):279—282. [3]郁春潮,陈韶华.广角镜头的摄像机非线性标定技术 [J].湖北大学学报(自然科学版),2005,27(4): 351—354. [4]王子亨,穆森,邱桂苹.摄像机非线性标定方法[J]. 计算机工程与设计,2010,31(15):3526—3529. [5]张玉发,孙晓泉.一种基于同心圆环的图像畸变校正 方法[J].光电技术应用,2007,22(2):63—65. [6]朱日红,李建欣.光学成像系统中非线性畸变的数字 校正方法[J].南京理工大学学报,2004,28(4): 414—416. [7]朱铮涛,黎绍发.镜头畸变及其校正技术[J].光学技 术,2005,31(1):136—138. [8]周富强,胡坤,张广军.基于共线特征点的摄像机镜 头畸变校正[J].机械工程学报,2006,42(9):174— 177. [9]孙双花.视觉测量关键技术及在自动检测中的应用 [D].天津:天津大学,2007. (责任编辑:刘划英文审校:刘敬钰)