激光与光电子学进展, 2020, 57 (4): 041514, 网络出版: 2020-02-20   

基于改进粒子群算法的相机内参优化方法 下载: 1097次

Optimization Method for Camera Intrinsic Parameters Based on Improved Particle Swarm Algorithm
作者单位
1 南京林业大学机械电子工程学院, 江苏 南京 210037
2 南通职业大学机械工程学院, 江苏 南通 226007
摘要
相机标定是机器人视觉系统中实现精确定位的重要前提,针对传统相机标定精度不高的问题,提出基于改进粒子群算法的相机标定优化方法。该方法以张正友标定方法获得相机内参初始值,在不同迭代阶段实现对惯性参数非线性自适应调整,以平衡局部和全局搜索能力;对社会和自身学习率采用不同迭代阶段正余弦变化的动态自调整策略,进一步提高全局搜索能力与后期搜索精度;在粒子群快要陷入局部最优时,采用驱散机制扩大粒子群所在空间范围,避免算法过早收敛。实验结果表明,所提相机标定方法与传统标定方法相比具有较高的精度和较好的可重复性。
Abstract
Camera calibration is an important premise for accurate positioning in robot machine vision systems. To solve the problem of low accuracy of traditional camera calibration, this paper proposes a camera calibration optimization method based on an improved particle swarm optimization algorithm. This method uses Zhang Zhengyou calibration method to obtain the initial value of camera intrinsic parameters and realizes the nonlinear self-adaptive adjustment of inertial weight parameters in different iteration stages, balancing the local and global search capabilities. Dynamic self-adjusting strategies of sines and cosines changes in different iteration stages are adopted for global and local learning factors to further improve the global search ability further and late search accuracy. When a particle swarm is about to fall into the local optimum, the dispersing mechanism is used to enlarge the spatial range of the particle swarm to avoid premature convergence of the algorithm. Experimental results show that the proposed method has better precision and repeatability as compared with the traditional methods.

1 引言

基于视觉伺服的机器人操作在很多行业中的应用越来越广泛,其中机器人末端定位作业是目前研究较多的功能[1-3],需要对目标物体进行识别和精确定位,这与相机内参密切相关。另外,机器人与相机组成的视觉系统在应用前首先需要确定相机坐标系与机器人基坐标系或末端坐标系之间的关系,即完成手眼标定,而相机内参精度是手眼标定精度的重要影响因素。综上所述,提高相机标定精度是提高机器人视觉作业精度的重要基础手段,而相机标定内参优化是提高相机标定精度的一种重要方法。

目前应用广泛的相机标定方法主要有张正友标定法[4]和Tisa两步标定法[5],Tisa两步标定法是先求出外部参数再求出内部参数,过程简单,但对设备要求高,成本高。张正友标定法所需的标定模板简单,只需获取多幅不同角度的标定模板图像,便可求解参数,但标定精度偏低。近年来很多学者提出各种改进算法,主要包括:利用线结构光完成相机标定的方法[6],此类方法操作简单,但设备成本高;将模拟退火算法与遗传算法相结合的方法[7],此类方法精度有所保证,但算法计算量大;相机标定线性法[8],此类方法计算过程简单,但标定精度有待提高。还有方法将粒子群算法融入传统相机标定中[9-11],得到的标定精度有所提高,但总体来说这些算法还是不能有效避免陷入局部最优,标定精度仍有待提高。

本文提出一种基于改进粒子群算法(IPSO)的相机标定内参优化方法,先确立粒子群优化算法的初始值,再对基本粒子群算法中的惯性参数和学习率进行优化调整,采用驱散机制扩大粒子群所在空间范围,以避免陷入局部最优,提高全局搜索能力,从而更好地提高算法的标定精度。

2 单目相机成像模型

相机模型可表征场景目标空间点与对应成像平面点的投影关系。本研究所使用的相机成像模型为针孔成像几何模型,成像关系呈线性变化。如图1所示,场景空间中任一目标点Pw在成像平面xoy上的投影为该点和光心oc的连线与该成像平面的交点Ppix图1还包含了世界坐标系(ow-xwywzw)、相机坐标系(oc-xcyczc)、图像物理坐标系(o-xy)和图像像素坐标系(opix-uv)。

世界坐标系下的目标点Pw与其在像素坐标系下的投影点Ppix之间的关系可表示为

zcuv1=fx0u00fyv0001r11r12r13t1r21r22r23t2r31r32r33t3xwywzw1=MRTxwywzw1,(1)

式中:fxfy分别为相机在xy方向的焦距,单位为pixel;M为内参矩阵,主要与相机内部结构有关; RT为相机坐标系相对于世界坐标系位姿的外参矩阵。

图 1. 单目相机模型

Fig. 1. Monocular camera model

下载图片 查看所有图片

在实际应用中,所有相机都会存在一定的畸变,而单目相机模型是理想的线性模型,并没有涉及相机的畸变。因此,必须将相机畸变考虑进去,则有[12]

xuyu=1+l1(x2+y2)+l2(x2+y2)2xy+2p1xy+p2(y2+3x2)p1(x2+3y2)+2p2xy,(2)

式中:(xu,yu)为图像坐标系畸变点的位置;(x,y)为图像坐标系理想点的位置;l1l2p1p2均为畸变系数。(2)式等号右边第一项描述径向畸变,第二项描述切向畸变。

3 相机内参优化算法初始值的确立

相机内参优化算法初始值确立的好坏直接影响优化算法的执行结果。相机内参的标定需要通过采集图像,根据物体在世界坐标系下的已知坐标矩阵[ xwywzw1]T和物体在图像坐标系下的像素坐标矩阵[ uv1]T,计算相机的内参矩阵M。本研究采用标准的标定板来定义世界坐标系,在标定板定义的世界坐标系中,每个坐标点均分布在标定板平面内,即zw=0,则有

zcuv1=MRTxwywzw1=Mr1r2txwyw1=Hxwyw1,(3)

式中:H= h11h12h13h21h22h23h31h32h33;t为位置向量;r1r2分别为旋转矩阵R的第一列向量和第二列向量。

将(3)式改写为

xwyw1000-uxw-uyw-u000xwyw1-vxw-vyw-vh=0,(4)

式中:hH矩阵拉直的列向量。给定一个标定板的图像,得到

h1h2h3=λMr1r2t,(5)

式中:λ为尺度系数;h1h2h3分别为矩阵H的第一列向量、第二列向量和第三列向量。

因为r1r2具有正交关系: r1Tr1= r2Tr2=1且 r1Tr2=0,有

hT1M-TM-1h2=0hT1M-TM-1h1=h2-1M-TM-1h2(6)

通过上述两个基本约束关系可求解相机内参,即

B=M-TM-1=1fx20-u0fx201fy2-v0fy2-u0fx2-v0fy2u02fx2+v02fy2+1(7)

矩阵B是对称矩阵,由6个参数确定:

b=b11b12b22b13b23b33T(8)

H的第i列为hi= hi1hi2hi3T,有

hTiBhj=qTijb=hi1hj1hi1hj2+hi2hj1hi2hj2hi3hj1+hi1hj3hi3hj2+hi2hj3hi3hj3Tb11b12b22b13b23b33,(9)

式中:qij为矩阵H的元素运算后组成的列向量。由正交约束关系(6)式可得:

qT12q11-q22Tb=0(10)

对于(10)式,可以采用最小二乘来求解b,从而获得相机模型的内参:

fx=λb11fy=λb11b11b22-b122u0=-b13fx2λv0=b12b13-b11b23b11b22-b122(11)

4 基于改进的粒子群算法的相机内参优化

4.1 基本粒子群算法

基本粒子群算法(SPO)是一种智能进化算法,已在问题优化过程中广泛应用。其基本思路是先随机初始化一组粒子,之后在每次迭代周期中更新粒子[13]。如在一个n维的搜索空间中,第i个粒子当前位置为Xi=(Xi1,Xi2,…,Xin)T,它的广义速度为Vi=(Vi1,Vi2,…,Vin)T,则粒子更新方程可表示为

Vi(k+1)=wVi(k)+c1e1pbest-Xid(k+1)+c2e2gbest-XikXi(k+1)=Xi(k)+Vi(k+1),(12)

式中:pbest为该粒子的历史最优位置;gbest为整个粒子群的最优位置;k为当前迭代次数;w为惯性参数,用以调节收敛方向;c1c2分别为自身和社会学习率,用于调节迭代的步距;e1e2为(0,1)区间内互相独立的随机数。

相机标定的目标函数可表示为

fobj=mini=1nj=1mp^ij-p(M,k1,k2,p1,p2,Ri,Ti),(13)

式中: p^ij为第j个标定点在第i次采样时候的像素坐标。

4.2 惯性参数调整

惯性参数w是SPO算法的重要参数之一,它的大小决定了当前迭代粒子的飞行速度对上一次迭代粒子飞行影响的程度[14-15]w取较大值时,全局搜索能力增强,但不易收敛,尤其在后期不易于寻找最优精度;w取较小值时,局部搜索能力增强,但容易陷入局部最优,算法过早收敛。国内外学者对此进行了大量研究,经典的线性递减权值 (LDW)方法可控制线性变化、计算过程简单,但对复杂或非线性问题的求解性能一般。经过分析研究,发现在不同迭代阶段利用e指数可实现对惯性参数的非线性自适应调整,以达到预期效果。具体关系式为

w=[exp(-k)(wmax-wmin)(kmax-k)/kmax]+wmin,(14)

式中:wmax为惯性参数最大值;wmin为惯性参数最小值;kmax为最大迭代次数;k为当前迭代次数。

由(14)式可知,w的值与收敛程度有关,迭代初期,w的值相对较大,可充分扩大搜索空间范围,提高全局寻优能力。迭代后期,w的值相对较小,以精细的搜索达到快速收敛的目的。

4.3 学习率调整

由粒子群算法特性可知:迭代初期的个体粒子最优位置对粒子当前位置影响程度较大,大的c1和小的c2更有利于提高全局搜索能力;在迭代后期,粒子群的最优位置主要取决于粒子当前最优位置,小的c1和大的c2更有利于精确收敛于全局最优值。经过分析研究发现,对全局和局部学习率采用不同迭代阶段正余弦变化的动态自调整策略更符合算法进化规律,符合要求的关系式可表示为

c1=(c1max-c1min)cosπk2kmax+c1min,(15)c2=(c2max-c2min)sinπk2kmax+c2min,(16)

式中:c1maxc1min分别表示自身学习率c1的最大值和最小值;c2maxc2min分别表示社会学习率c2的最大值和最小值。

4.4 驱散机制

根据前述对粒子群算法的惯性参数和学习率调整,可以在一定程度上避免陷入局部最优,但并不能完全避免算法过早收敛。因此,如果在粒子群算法趋于局部最优的过程中,引入驱散机制触发驱散动作来扩大当前粒子群搜索空间的范围,即跳出局部最优重新在一定范围内搜索,以保障全局最优搜索能力。

传统的驱散机制在整个迭代过程中一直存在,而这会导致收敛变慢[16]。只有在陷入局部最优时引入驱散机制,才能在收敛速度不受影响的情况下避免陷入局部最优。在迭代初期,相邻两次的全局最优位置相差比较大,随着迭代次数的增加,相邻几次的全局最优位置之差减小,这就很容易陷入局部最优值,可用最近的相邻两次迭代差值的平均值衡量随迭代次数而变化的收敛程度,具体关系式可表示为

s(k)=1za=0z[gbest(k-a)-gbest(k-a-1)],(17)

式中:s(k)表示迭代k次的收敛程度;gbest(k-a)表示迭代(k-a)次的全局最优位置;z一般取3~4。由(17)式可知,通过判断s(k)是否小于某一值ε,就可判断是否需要触发驱散动作。触发驱散动作后粒子群被分配到新的空间范围,随着迭代次数的增加,每次驱散后的全局搜索空间范围越来越小,因此驱散动作程度也相应减小。驱散动作程度可表示为

d=402+exp(k/30),(18)

式中:d表示驱散动作程度;k表示触发驱散动作时的迭代次数。粒子的新位置关系可表示为

Xi(k+1)=dXi(k),(19)

式中:Xi(k)为第i个粒子当前位置;Xi(k+1)为粒子被驱散后的新位置。

4.5 相机内参优化流程

基于改进粒子群算法的内参优化标定流程如下:

1) 粒子初始化。随机生成N个包含位置与速度信息的粒子,初始化粒子群优化算法中使用的参数wc1c2等,确定最大迭代次数kmax,并初始化当前迭代次数,即k=0。

2) 相机内参初始化。根据前述相机内参初始值的确立,计算相机标定内参优化初始值fxfyu0v0l1l2p1p2

3) 选择适应度函数。优化目标函数是利用角点提取获得的实际像素坐标点(u'i,v'i)与世界坐标系中通过成像数学模型计算得到的像素坐标点(ui,vi)之间的距离,适应度函数的定义为

f=mini=1m(u'i-ui)2+(v'i-vi)2(20)

4) 更新惯性参数和学习率。对惯性参数采用非线性自适应调整策略,对全局和局部学习率采用正余弦变化的动态自调整策略。

5) 选择粒子。采用适应度函数计算每个粒子的得分,并与其之前历史最好状态pbest进行对比,如果更好,则更新pbest为当前状态;同时对每个粒子与所有全局粒子的最好状态gbest进行对比,如果更好,则更新gbest为当前状态。

6) 更新粒子。根据前述粒子更新公式[(12)式],在每次迭代中更新粒子。

7) 判断是否局部最优。判断当前位置是否可能陷入局部最优,若是则转入步骤8),若不是转入步骤9)。

8) 驱散动作。触发驱散动作,扩大当前粒子群的搜索空间范围,即跳出局部最优在一定范围内重新搜索。

9) 重复步骤4)~8),直到最大迭代次数或者达到预设的适应度值则结束。

5 实验结果与分析

基于前述相机标定原理与流程,用标准的棋盘格作为标定板。将单目相机固定在机械臂末端,控制机械臂末端的位姿使其带着待标定的相机从不同视角拍摄标定板,得到20幅采样图像。每幅采样图像中每个棋盘格角点在标定板平面上的坐标固定不变且已知,而棋盘格角点在图像平面上的像素坐标可以直接采用Harris角点提取计算。根据前述相机内参初始值的确定方法求得内参初始值,并将其附近区间设定为初始值优化范围,具体为 fxfyu0v0±80(80为初始值上下浮动的区间,下同), l1l2p1p2±10。设定粒子群的粒子数量为80,惯性系数wmax=0.8,wmin=0.3,自身学习率c1max=1.5,c1min=0.7,社会学习率c2max=2,c2min=1.5,最大迭代次数kmax取400次。

表1中分别列出了张正友标定法(Zhang)、基本粒子群算法(PSO)和改进的粒子群算法(IPSO)的相机内参标定结果。为了验证这3种方法的标定精度,计算了标定数据的重投影误差(reprojection error),图2所示为张正友标定法的重投影误差,其平均重投影误差为0.250178 pixel。图3所示为基本粒子群算法的重投影误差,其平均重投影误差为0.223765 pixel,比张正友方法的标定精度有所提高。图4所示为改进的粒子群算法的重投影误差,其平均重投影误差为0.085126 pixel,重投影误差与前两种方法相比有所减小,标定精度得到一定提高。

图 2. 张正友方法的重投影误差

Fig. 2. Reprojection error of Zhang Zhengyou calibration method

下载图片 查看所有图片

表 2. 可重复性实验的重投影误差

Table 2. Reprojection error of repeatable experiments

ExperimentalgroupAmountof imagesZhang methodPSO methodIPSO method
Meanvalue /pixelStandarddeviation /pixelMeanvalue /pixelStandarddeviation /pixelMeanvalue /pixelStandarddeviation /pixel
1100.2453230.2145360.2231670.1677360.0842330.054263
2250.2531780.2875520.2365860.1899210.1123580.082621
3400.2884360.3142530.2633520.2543730.1213260.153323
4550.2988320.3021770.2681210.2477360.1355370.151852
5700.2766390.2989270.2308550.2066250.1052330.127735
6850.2566320.2799820.2279630.1867290.0983120.063298
71000.2543860.2656740.2216580.1756730.0979450.043652

查看所有表

图 3. 基本粒子群算法重投影误差

Fig. 3. Reprojection error of PSO method

下载图片 查看所有图片

表 1. 相机内参标定结果

Table 1. Calibration results of camera intrinsic parameters

MethodZhangPSOIPSO
fx297.45560297.65350298.14325
fy297.32785297.61227298.13513
u0252.65332253.01923251.35418
v0210.35674210.74318211.23546
l10.132600.123720.12540
l2-0.43980-0.42615-0.42913
p1-0.00556-0.00672
p2--0.006870.00543

查看所有表

图 4. 改进粒子群算法的重投影误差

Fig. 4. Reprojection error of IPSO method

下载图片 查看所有图片

为了验证改进粒子群算法的可重复性,重新进行了多次采样实验。利用上述三种方法进行标定,记录重投影误差并计算平均值与标准差,各组数据结果如表2所示,重投影误差对比如图5所示,其中改进粒子群优化算法得到的平均重投影误差比张正友方法减小了0.15 pixel左右,比基本粒子群算法减小了0.1 pixel左右,基本维持在0.1 pixel左右,优于前两种方法。

图 5. 重复性实验平均投影误差

Fig. 5. Mean reprojection errors of repeatable experiments

下载图片 查看所有图片

6 结论

提出一种用于提高相机标定精度的优化方法,首先建立了存在畸变影响的相机模型,再利用张正友标定方法确定相机内参优化算法的初始值;然后对粒子群算法中的参数作了调整并增加了全局寻优的保障机制,即对惯性参数作了非线性动态调整,对学习率采用正余弦变化的自适应调整策略,这些都提高了全局寻优能力与后期搜索精度,并采用驱散机制根据驱散动作程度重新分配粒子群空间范围,使其跳出局部最优;最后将几种传统方法与所提改进方法进行对比实验,结果表明改进后的相机标定方法重投影误差减小,基本维持在0.1 pixel左右,标定精度得到一定的提高,并且具有较好的可重复性。在机器视觉应用系统中,所提方法具有较高的有效性和可靠性。

参考文献

[1] 马庭田. 基于机器人视觉的多类型工件识别与定位问题研究[D]. 南京: 南京航空航天大学, 2018: 8- 18.

    Ma TT. Research on identification and positioning problems for a variety of work pieces on binocular stereo vision[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2018: 8- 18.

[2] 李双. 基于机器视觉的叶片夹取机器人研究[D]. 哈尔滨: 东北林业大学, 2018: 42- 47.

    LiS. Research on blade recognition robot based on machine vision[D]. Harbin: Northeast Forestry University, 2018: 42- 47.

[3] 桑孟雷. 基于机器视觉的圆形工件位姿估计[D]. 杭州: 浙江大学, 2018: 11- 24.

    Sang ML. Machine vision based location and attitude estimation for circular workpiece[D]. Hangzhou: Zhejiang University, 2018: 11- 24.

[4] Zhang Z. A flexible new technique for camera calibration[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2000, 22(11): 1330-1334.

[5] Tsai R. A versatile camera calibration technique for high-accuracy 3D machine vision metrology using off-the-shelf TV cameras and lenses[J]. IEEE Journal on Robotics and Automation, 1987, 3(4): 323-344.

[6] 解则晓, 张安祺. 超大尺度线结构光传感器内外参数同时标定[J]. 光学学报, 2018, 38(3): 0315001.

    Xie Z X, Zhang A Q. Simultaneous calibration of the intrinsic and extrinsic parameters of ultra-large-scale line structured-light sensor[J]. Acta Optica Sinica, 2018, 38(3): 0315001.

[7] 游江, 唐力伟, 邓士杰. 基于改进遗传模拟退火的相机标定方法[J]. 计算机工程与设计, 2017, 38(3): 819-824.

    You J, Tang L W, Deng S J. Camera calibration method based on genetic algorithm and simulated annealing[J]. Computer Engineering and Design, 2017, 38(3): 819-824.

[8] Zhang Z, Zhao R J, Liu E H, et al. A single-image linear calibration method for camera[J]. Measurement, 2018, 130: 298-305.

[9] 黄伟光, 董安国. 基于粒子群算法的摄像机自标定[J]. 计算机应用与软件, 2015, 32(5): 216-219, 233.

    Huang W G, Dong A G. Camera self-calibration based on particle swarm optimisation[J]. Computer Applications and Software, 2015, 32(5): 216-219, 233.

[10] 郭彤颖, 李宁宁, 刘雍. 基于粒子群算法的摄像机内参数优化方法[J]. 激光与光电子学进展, 2017, 54(11): 111504.

    Guo T Y, Li N N, Liu Y. Optimization of camera internal parameters based on particle swarm algorithm[J]. Laser & Optoelectronics Progress, 2017, 54(11): 111504.

[11] 王道累, 胡松. 基于量子粒子群优化算法的摄像机标定优化方法[J]. 激光与光电子学进展, 2018, 55(12): 121502.

    Wang D L, Hu S. Optimization method of camera calibration based on quantum-behaved particle swarm optimization algorithm[J]. Laser & Optoelectronics Progress, 2018, 55(12): 121502.

[12] 秦瑞康, 杨月全, 李福东, 等. 基于全参数自适应变异粒子群算法的单目相机标定[J]. 东南大学学报(自然科学版), 2017, 47(z1): 193-198.

    Qin R K, Yang Y Q, Li F D, et al. Monocular camera calibration based on particle swarm algorithm with all parameter adaptive mutation mechanism[J]. Journal of Southeast University(Natural Science Edition), 2017, 47(z1): 193-198.

[13] 李少波, 张成龙, 郑凯. 基于改进粒子群算法的制造车间WSN部署优化[J]. 仪表技术与传感器, 2017( 10): 101- 104.

    Li SB, Zhang CL, ZhengK. Node deployments optimization in wireless sensor networks based on improved particle swarm algorithm[J]. Instrument Technique and Sensor, 2017( 10): 101- 104.

[14] 赵丰明, 樊艳芳, 钱福如, 等. 关于新能源并网供电无功性能优化设计研究[J]. 计算机仿真, 2018, 35(11): 118-122.

    Zhao F M, Fan Y F, Qian F R, et al. Research on optimized design of reactive power performance of new energy grid connected power supply simplified particle swarm optimization[J]. Computer Simulation, 2018, 35(11): 118-122.

[15] 杜江, 袁中华, 王景芹. 动态改变惯性权重的新模式粒子群算法[J]. 安徽大学学报(自然科学版), 2018, 42(2): 60-66.

    Du J, Yuan Z H, Wang J Q. New model of particle swarm optimization algorithm with dynamically changing inertia weight[J]. Journal of Anhui University(Natural Science Edition), 2018, 42(2): 60-66.

[16] 游佳丽, 周志勇, 章程, 等. 基于自适应驱散机制的粒子群优化算法[J]. 计算机工程与应用, 2017, 53(7): 41-48, 103.

    You J L, Zhou Z Y, Zhang C, et al. Adaptive dispersion mechanism based particle swarm optimization algorithm[J]. Computer Engineering and Applications, 2017, 53(7): 41-48, 103.

徐呈艺, 刘英, 肖轶, 曹健. 基于改进粒子群算法的相机内参优化方法[J]. 激光与光电子学进展, 2020, 57(4): 041514. Chengyi Xu, Ying Liu, Yi Xiao, Jian Cao. Optimization Method for Camera Intrinsic Parameters Based on Improved Particle Swarm Algorithm[J]. Laser & Optoelectronics Progress, 2020, 57(4): 041514.

本文已被 5 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!