基于tsai方法多项畸变模型单目摄像机标定.docVIP

基于tsai方法多项畸变模型单目摄像机标定 　　【摘要】在基于tsai一阶径向畸变标定方式的基础上，提出了一种改进算法。算法结合传统标定原理首先线性求解摄像机模型内外参数，作为初值，应用于带有二阶径向畸变，偏心畸变和薄棱镜畸变的图像矫正模型，通过最小二乘法求解6项畸变尺度因子。充分利用OpenCV视觉函数库对图像分步校正，避免了重复计算并判断透镜非线性映射。实验结果表明，改进算法提高了系统标定精度，降低了计算时间复杂度，可以实现较高精度的单目视觉标定及立体视觉匹配的应用需求。 　　【关键词】单目标定；畸变校正；立体匹配；tsai方法 　　计算机立体视觉技术的实现主要包括图像获取，摄像机标定，特征提取，立体匹配和三维重建五部分[1]。摄像机标定作为空间三维图景到成像二维平面的首要条件，是立体视觉系统图像深度测量的基础部分。而夹角，都不可避免的作用在摄像机透视成像过程上。对于以上镜头径向，切向畸变，目前标定方法主要分制造中所引起的光学镜头径向曲率变化，装配中导致 的多镜头光轴不同轴以及镜头与摄像机像平面存在的为传统标定方法[2]，自标定方法[3]以及基于主动视觉的标定方法。典型的传统标定方法主要有abdel-aziz和karara提出的直接线性变换法（DLT），hallert提出的非线性优化法，tsai的两步法以及张正友的平面模板法[4]。这其中tsai[5]两步标定法的主要贡献是在RAC条件下，用DLT方法求解透镜矩阵系统参数，并以此作为初始迭代值，非线性优化估计畸变参数，减少了求解参数空间维数，提高了线性求解速度和非线性求解精度，并广泛应用于实时性较高的机器视觉场合。tsai的方法中只引入了一级径向畸变，对于一些精度要求高的场合，往往需要引入切向畸变。本文在tsai方法基础上，引入二级径向畸变，偏心畸变，薄棱镜畸变，全面分析摄像机透镜模型，结合OpenCV视觉函数库标定摄像机内外参数，并给出畸变校正后的棋盘格图像。 　　1.透镜成像及畸变校正 　　图1为广泛应用于工业非接触尺寸测量的基于径向畸变的摄像机针孔模型[6]，O-XWYWZW为全局坐标系，摄像机透镜光心与坐标系O-XCYCZC的原点重合，xc轴和yc轴分别与xw轴yw轴平行，zc轴与zw轴平行，O-XY为像平面坐标系，以mm为单位，O-UV为像素坐标系，以pixels为单位，点P在像素坐标系O-UV下的离散坐标（u，v）恒为整数值，?表示摄像机焦距，标定点P的成像过程由以下几步坐标系转换完成：（1）全局坐标系到摄像机坐标系的刚体变换；（2）像平面坐标系到像素坐标系的变换；（3）摄像机坐标系到像平面坐标系的变换；经过上述坐标系齐次变换，得到点P全局坐标与像素坐标关系如下[7]： 　　式（1）中，M1为摄像机内参数矩阵，ax，ay分别为x轴y轴上尺度焦距，单位为pixels，u0，v0为图像中心坐标，M2为摄像机外参数矩阵，为3×3正交旋转矩阵，平移矩阵，M为摄像机投影矩阵。摄像机标定的目的就是精确求解内参数?x，?y，u0，v0和外参数矩阵R，T。由于无法避免透镜畸变，（1）式中的待求解参数并不能代表实际的透镜参数，因此需要引入非线性畸变校正摄像机投影矩阵各项参数。成像点 的坐标畸变转换如下[8]： 　　式中，（xu，yu）和（xd，yd）表示点p在像平面坐标系O-XY下的无畸变坐标和畸变坐标，ζx，ζy为投影畸变在像平面x轴和y轴上的映射，tsai的标定方法中只考虑了一阶径向畸变，限制了复杂背景下对摄像机标定的高精度要求，扩展带有二阶径向畸变，偏心畸变和薄棱镜畸变的校正模型如下： 　　式中r2=xd2+yd2，k1，k2为径向畸变尺度，p1，p2为切向畸变尺度，ε1，ε2为薄棱镜畸变尺度在成像平面上的径向和切向分量。通过求解（1），（3）式摄像机参数，空间点的三维坐标得以重建。 　　2.改进的标定方法 　　首先线性求解摄像机投影矩阵M，然后以M各项投影参数为初始迭代值，结合三项畸变模型，求解畸变系数，最后以中心点坐标区域局部最小值为数据拟合函数，全局优化所求摄像机内外参数初值和畸变因子。具体求解方法如下：（1）求解摄像机内外参数。对于空间i=1，2，3...n个标定坐标点，由式（1）可以得到： 　　其中，（Xwi，Ywi，Zwi，1）表示标定靶第i个点的三维坐标的齐次坐标， 表示第（ui，vi，1）个点的图像坐标的齐次坐标，为投影元素。为了消除摄像机坐标Zci的影响，可以将（4）式用线性方程表示如下： 　　根据线性模型的成像原理，每个已知的空间点全局坐标和像素坐标都不同的对应一个式（5）所表示的两个线性方程，由（5）式可知，线性方程一共有12个待求投影参数，通常认为空间两坐标系的变换并不受相乘系数的影响，因此对于（5）式，可令m12=1，当有