提高工业视觉测量系统精度的途径

HEATS

一、引　言

　　工件三维曲面或三维轮廓测量技术广泛应用于工业、科研、国防等领域。汽车车身、飞机机身、轮船船体、汽轮机叶片等加工制造中的在线检测，特别是大型工件的曲面检测一直是生产中的关键技术难题。该类工件在车间条件下一般采用靠模法测量，但可测截面少，测量精度低；在计量室条件下采用三坐标测量机测量虽然精度较高，但数据采集速度慢，测量成本高，且难于实现在线测量。鉴于接触式测量方法的局限性，激光三角法、莫尔投影法、工业视觉测量法等多种非接触测量方法日益受到重视，其应用也渐趋广泛。
　　工业视觉测量技术（或称数字近场摄影测量技术）是一种立体视觉测量技术^［1］，其测量系统结构简单，便于移动，数据采集快速、便捷，操作方便，测量成本较低，且具有在线、实时三维测量的潜力，尤其适合于三维空间点位、尺寸或大型工件轮廓的检测。

/ O2 k" n6 E3 P5 q2 O( N# E

二、测量原理

9 ^ A! M2 t/ v @# t: A" c3 D7 h% a

　　利用CCD摄像机可以获得三维物体的二维图像，即可以实现实际空间坐标系与摄像机平面坐标系之间的透视变换。通过由多个摄像机从不同方向拍摄的两帧（或两帧以上）的二维图像，即可综合测出物体的三维曲面轮廓或三维空间点位、尺寸。
　　为便于说明，设物空间坐标系为O-XYZ，CCD像面的像平面坐标系为o-xy。
　　现以双摄像机为例说明系统的透视变换关系。如图1所示，P为任一空间三维物点，设该点的物空间坐标为P（X，Y，Z），其在摄像机Ⅰ和摄像机ⅡCCD像面上的像点坐标分别为P₁（x₁，y₁）和P₂（x₂，y₂）。

图1　物空间坐标系和双摄像机的像平面坐标系

　　对于摄像机Ⅰ，像点坐标与物点坐标的变换关系为^［2］

       （1）

　　其中w₁为非零参数，a₁，a₂，…，a₁₁为系统变换矩阵的元素，与摄像机Ⅰ的安放位置及成像系统Ⅰ的参数有关，可通过系统定标来确定。
　　对于摄像机Ⅱ，像点坐标与物点坐标的变换关系为

       （2）

4 c# y; y! W' q0 G) i

　　其中w₂为非零参数，b₁，b₂，…，b₁₁为系统变换矩阵的元素，与摄像机Ⅱ的安放位置及成像系统Ⅱ的参数有关，也可通过系统定标来确定。
　　式（1）和式（2）可分别化为

, i1 v" X I t7 U! X. w/ x2 ?

       （3）        （4）

　　其中A＝［a₁，a₂，…，a₁₁］^T；B＝［b₁，b₂，…，b₁₁］^T。a_i和b_i加起来共22个未知参数，利用一个已知靶点和它在两个CCD像面上的像点可建立4个线性方程，欲求出22个未知参数，则至少需要6个已知靶点。利用已知的6个或6个以上靶点，根据上述方程即可求得这些参数。这是测量工作的第一步，称为系统定标，即求出双摄像机组成的测量系统的变换矩阵A和B。
　　测量的第二步是根据被测点在两个CCD像面上的像点坐标P₁（x₁，y₁）和P₂（x₂，y₂）求得未知点P的物空间三维坐标（X，Y，Z）。
　　由式（1）和式（2）还可得到

- d ~& r% J7 z2 H) g

       （5）       （6）

　　由上式可求得3个未知数X，Y，Z，即P点的三维空间坐标值。这样就可以利用双摄像机拍摄的两个二维图像逐点测量物体的三维轮廓或尺寸。
　　工业视觉测量实验系统如图2所示。系统由两个CCD摄像机（分辨率为510（H）×492（V）像素）、一个带有6个以上参考点（已知坐标）的靶标、一块图像采集卡（3路，8位，256灰度阶，分辨率为512×512像素）、被测物体、一台PC机及相应软件组成。

+ @" Q: Q" P" E* M4 c2 ]

图2　工业视觉测量实验系统

) ?$ V& ^. B8 H4 e3 y( v

　　靶标用于对由两个摄像机构成的系统进行定标，其上的各靶点空间位置坐标都是预先经精确测定而已知的。图象采集卡把CCD接收到的图像信号输入计算机，由计算机进行处理。
　　将靶标上特征点的坐标和两个CCD像面上对应像点的坐标代入式（3）和式（4），便可确定两个成像系统的变换矩阵，即完成系统定标。工件上的待测点可以通过粘附高反射标记或用激光投射光点产生。把每个待测点在两个摄像机像面上的像点坐标代入式（5）和式（6），由系统的变换矩阵即可确定该待测点在实际空间坐标系中的坐标值（X，Y，Z）。若采样率足够大，即可逐点测量并重建工件的复杂曲面轮廓。

三、影响精度的因素及提高精度的途径

V, m/ I* f! w* a0 g

　　对于三维空间点位或距离的测量，影响工业视觉测量系统精度的测量误差主要产生于定标和测量两个步骤中。现对影响系统精度的主要因素及提高精度的途径作一分析。

/ @; a, v6 P: s, ^! T5 z

　　1．CCD摄像机和图象采集卡的分辨率
　　由于工业视觉测量系统成像倍率较大，CCD摄像机及图象采集卡的分辨率对点位在空间的测量精度影响极大，因此应用于测量目的的立体视觉系统宜采用尽可能高的分辨率，因硬件系统分辨率有限造成的像面目标点的定位误差可通过软件补偿进一步减小，使其达到亚像素级。

( X/ D; I- ?7 N4 }6 i

　　2．成像系统的畸变
　　在工业视觉测量系统的测量模型中，认为透镜成像处于理想状态，而客观存在的成像系统畸变会导致系统产生定标误差。透镜畸变误差可表示为

f9 \6 h8 X; ]7 q8 E

       （7）

# X4 v5 b5 Z) R# J- S1 y2 s, w5 X

式中　k_i，p_i——与透镜畸变有关的参数
　　　r——像面的点离开光轴的距离
通过对式（1）和式（2）中的x，y进行误差补偿（补偿值为Δx，Δy），然后建立新的非线性方程组，即可解出准确的系统定标参数。对于双摄像机系统，可补偿畸变误差的测量模型^［3］共需求出32个参数，因此定标时至少需要8个已知特征点。

　　3．靶标的设计
　　对于一个摄像机，式（3）和式（4）代表物空间点（X，Y，Z）和其像点（x，y）所确定的一条直线，6个特征点表示有6条直线通过透镜的投影中心，但由于畸变等非线性的影响，6条直线不会完全交汇于同一投影中心，而是在投影中心附近形成一松散的直线束。通过合理安排定标点的位置，使这一直线束越紧密，则表示非线性的影响越小。由此可得出结论：每3个定标点和透镜投影中心非共面时，系统的定标精度可大为提高。一般使靶标上不存在3点共线即可^［4］。
　　有条件时，靶点的分布及其范围应尽可能考虑测量对象的范围和形状，因为被测物体越靠近靶点，测量精度越高。靶标上通常至少需要6个靶点，考虑补偿透镜畸变时则需8个以上靶点。适当增加靶点，在系统定标时可达到较高精度。
　　通常选用圆形靶点，白底黑点或黑底白点均可。靶体应稳定，不易变形，并预先用三坐标测量机精确测量靶点的空间位置，其测量精度对系统定标精度将产生直接影响。

　　4．两摄像机光轴交角
　　成像系统采用透视投影原理，位于一条空间直线上的点都能成像在同一个像点位置。像面上目标点中心的定位误差如在一个圆形区域内，它在物空间则对应一个圆锥形区域，两个摄像机因像面上目标点定位误差产生的纵深方向的物点测量误差将由两个圆锥的交叠部分决定，见图3。

8 n: ^9 O# I" k/ K

图3　纵深误差分析示意图

" ?+ Y2 V7 M& W1 _" ?6 B

　　图中，O为空间物点，AS，A′S′为两摄像机成像面，P和P′为两投影中心，α为两光轴交角，摄像机焦距为f，目标点中心定位误差为Δl，物点到摄像机投影中心的距离为L。根据几何关系，纵深方向测量误差近似为

, \+ C. I. `. S3 k

ΔZ≈LΔl／fsin（α／2）       （8）

' v& J8 o R* v) ~ M

　　例如：Δｓ＝0．001mm，f＝16mm，L＝2000mm，α＝60°，则ΔZ＝0.25mm。
　　随着两摄像机光轴交角的增大，测量时因像面目标点中心定位误差引起的纵深方向测量误差减小。根据放大倍率公式，X，Y方向误差为

ΔX＝ΔY≈LΔl／f≈0.13mm

　　从式（8）还可看出，纵深方向测量误差还与物点距离L成正比，物体距离越远，测量误差越大。

　　5．像面上目标点的定位
　　在系统定标和点位测量两个步骤中，像面目标点像斑中心位置的定位精度是影响最终测量精度的主要原因。影响定位精度的因素有目标点像斑边缘平滑、像斑灰度不均匀、形状不规则等。此外，由于环境光线和电路噪声等因素的影响，还存在图像噪声。
　　通常需要采用图像处理技术对像斑进行边缘锐化等处理。确定中心坐标的算法主要有对称法、几何中心法、能量中心法、二次曲面拟合法^［5］等。其中，对称法受像斑灰度不均匀的影响较大；采用几何中心法时，如像斑边缘确定不准确，则定心误差较大；能量中心法和二次曲面拟合法可达到亚像素级定心精度，但计算量偏大。
　　测量时，靶点或测量点形状应规则，环境照明应尽量均匀，尽量减少杂散光的影响。采用激光束照明产生待测点时，光点形状有时会不规则。

四、空间点位测量实例

　　采用图2所示的工业视觉测量系统进行空间点位测量。被测物体到两摄像机连线的垂直距离为2000mm，两摄像机光轴交角为53°。分别采用6靶点（不考虑透镜畸变）和8靶点（透镜畸变补偿）对系统定标，然后对部分点的三维坐标进行测量。每点测量10次，在X，Y，Z方向均给出测量平均值、标准差、偏差以及测得的点位到实际点位的距离d。测量数据见表。

表　　空间点位测量结果

7 v1 ~8 ]# O, k- u9 x* L

' a7 T4 J) S+ a& D' `% y

|3 H4 N( k& S" N4 m% X+ h" G' Z; H; G# A+ [0 i/ F# A- V, [9 T1 V" S# T* j* P1 c7 e3 A8 D! d! S( L$ I" Z u% o7 u9 R9 X1 U1 V' X' |6 J$ j+ Y$ ^$ o1 }5 _; q; T- I4 ~* L9 o! F6 V1 L% a7 P7 H' {( h; X: B5 Q/ w1 A# _" i+ m8 \' @0 B0 Z$ B0 }) a: W: Y! Z% v, Y# e' ^6 g8 n" ]3 v2 }% S$ O7 l. I. }# G, F! J) [" ^) y' E# o- t3 A: N% Y6 }. C5 _9 s- _6 N# t. k/ H; x# \6 U k+ `+ b2 H0 m0 n1 @) E% N F* N3 Q d4 N# H) U0 I: ^# r( Z5 a5 k+ ^# z$ B0 n K/ n3 R: P6 d8 x4 c/ K6 l# w" q+ D; U. [3 [" n3 X. ?% p; o2 o0 l+ `1 @5 a" s& J. J! c$ F) F8 v# P% d2 F3 j' S0 F0 c/ k% j& o/ O# v% z4 @. e [9 K9 Q2 B5 B( c. D& Y! I! }6 T; H5 t+ b) n( _$ |9 r7 `9 ]3 W% N2 S! G; j9 \. j/ b( L3 U( K2 m2 x9 y7 D5 L# n% ~7 u. @/ I/ K2 y. O' a% D% `9 v) `- l! L) W9 p: v- v6 [# q# K" Z1 G1 q; _. o" i& u9 _, r2 q# T( E. ^$ m5 X7 J- F# @7 i

点号	真值l （mm） （X，Y，Z）	未考虑透镜畸变				考虑透镜畸变
		测量平 均值 （mm）	标准差 σ1	偏差 Δl1＝－l （mm）	平均距 离d1 （mm）	测量平均值 （消畸变） （mm）	标准差 σ2	偏差 Δl2＝－l （mm）	平均距 离d2 （mm）
26	239.622 199.766 75.139	239.721 199.820 75.064	0.0219 0.0525 0.0638	0.099 0.054 －0.075	0.136	239.638 200.779 　75.082	0.014 0.0453 0.0420	－0.016 　0.013 －0.06	0.060
27	199.320 200.645 40.256	199.256 200.568 41.360	0.0435 0.0678 0.0406	－0.064 －0.077 　0.104	0.145	199.286 200.682 40.174	0.0729 0.1117 0.6263	－0.034 　0.037 －0.082	0.096
44	159.313 319.501 40.124	159.347 319.360 40.157	0.0392 0.0299 0.0685	0.034 　－0.141　 　0.033	0.149	159.320 319.566 40.088	0.0435 0.0291 0.0592	0.007 　0.065 　－0.037	0.074
1	0.000 0.000 130.671	－0.098 　0.088 130.820	0.0149 0.0756 0.0444	－0.098　 0.088 0.149	0.184	－0.069　 －0.073　 130.542	0.0043 0.0126 0.0700	－0.069 －0.070 －0.017	0.099
30	39.397 200.814 40.073	39.467 200.861 39.986	0.0151 0.0352 0.0585	0.070 0.047 －0.087	0.121	39.474 200.861 40.111	0.0120 0.0191 0.0372	0.077 0.047 0.038	0.098

" r& \0 ]! l" a

　　实测结果表明，在240×320×90（mm³）空间范围内，测得的点位与实际点位的最大偏差值为0．2mm，补偿透镜畸变后的最大偏差值为0．1mm，测量精度明显提高。当测量范围较大，且位于视场边缘时，提高测量精度的效果更为明显。从标准差可以看出测量重复性较好。

五、结　语

　　工业视觉测量系统可用于空间点位置、空间尺寸或三维型面的测量。为了提高测量精度，应采用高分辨率的CCD摄像机和图像卡；合理安排摄像机距离和光轴交角；靶标设计及靶点预先测量的精度将直接影响系统定标精度；测量范围较大时，必须采用补偿透镜畸变的测量模型；像面上目标点的定位精度也是影响最终测量精度的关键因素，应采用图像处理技术和适当的计算方法提高定位精度。【MechNet】

- U1 Q3 U0 s0 A- z+ T& C- w& Y