首页 理论教育免费教程:用2步法轻松实现摄像机标定

免费教程:用2步法轻松实现摄像机标定

【摘要】:先利用线性变换方法或透视变换方法求解摄像机参数,再以求得的参数作为初始值,考虑畸变因素,利用非线性优化方法进一步提高标定精度,就形成了所谓的两步法。两步法计算量适中、精度较高。摄像机坐标系定义为原点在摄像机光学中心,z轴与光轴重合。

摄像机线性模型除包括线性模型中的全部参数外,还包括非线性径向畸变参数k1k2和切向畸变参数p1p2。线性模型的参数和非线性模型的参数一起构成了非线性模型的摄像机内部参数。

非线性模型求解是由非线性优化方法来估计,如广泛采用的Levenberg Marquardt(LM)优化算法和Newton Raphson(NR)优化算法。这些方法都涉及非线性方程求解,或需假设摄像机部分内参数可由其他方法测出;或者用线性模型首先计算出线性模型的参数,作为近似初值,再用迭代方法计算精确解。

针对三维立体靶标上的特征点,用线性模型首先计算出线性模型的参数,作为近似值,再用迭代的方法计算精确解,是非线性模型摄像机标定较为有效的方法。

先利用线性变换方法或透视变换方法求解摄像机参数,再以求得的参数作为初始值,考虑畸变因素,利用非线性优化方法进一步提高标定精度,就形成了所谓的两步法。

典型两步法是基于径向排列约束RAC的两步法:第一步是利用最小二乘法解超定线性方程,求得外部参数;第二步求解内部参数,先假定摄像机无镜头畸变,求解超定线性方程组给出初值。再结合非线性优化方法获得全部参数。两步法计算量适中、精度较高。

1.两步法采用的摄像机模型

如图5.16所示,设(XYZ)是三维世界坐标系中某目标点P的三维坐标。摄像机坐标系定义为原点在摄像机光学中心,z轴与光轴重合。(xyz)是同一点P在摄像机坐标系中的三维坐标。x′o′y′是原点在光轴z与图像平面的交点o′点平行于xy轴的像平面坐标系。有效焦距f是图像平面和光学中心的距离。(x′y′)是小孔成像模型下世界坐标点P的像平面坐标。(x*y*)是由于摄像机镜头畸变引起偏移的实际像平面坐标。(uv)是以像素为单位的图像坐标。(u0v0)是像主点坐标。与前文类似,从点P世界坐标(XYZ)通过投影变换映射为图像坐标由以下步骤完成。

(1)由世界坐标系到摄像机坐标转换

978-7-111-34687-6-Chapter05-157.jpg

式中,R为3×3旋转矩阵T为3×1平移矩阵。

(2)成像变换

978-7-111-34687-6-Chapter05-158.jpg

(3)畸变模型

造成成像坐标偏差畸变有径向畸变、偏心畸变和薄棱镜畸变。只考虑径向畸变并用二阶多项式近似描述径向畸变,有

978-7-111-34687-6-Chapter05-159.jpg

理想像平面坐标和实际像平面坐标的关系如下:

978-7-111-34687-6-Chapter05-160.jpg

式中,k为径向畸变系数。

(4)实际像平面坐标到图像坐标的转换

978-7-111-34687-6-Chapter05-161.jpg

式中,μ为尺度因子;Sx=SxNcx/Nfx为像素行方向的等效间距。所谓等效间距是指相邻两个像素对应的CCD阵列中的感光单元间距。假定CCD感光阵列行相邻单元间距为Sx,CCD每一行的感光单元个数为Ncx。从感光阵列到像素阵列的处理过程如下:图像获取硬件首先对每一行感光单元利用零阶采样保持器获得模拟波,然后再对该波进行采样,得到一行Nfx个像素。这样理想情况下像素行方向上的等效距离为Sx=SxNcx/Nfx,由于图像扫描是逐行进行的,像素列方向上的等效间距与感光阵列列方向的实际间距相同,均为Sy

通常情况下,CCD生产商提供的SxSy可以达到亚微米精度,NcxNfx对于特定CCD和图像采集卡为已知值,进而SxSy均为已知值。但在实际情况下,图像获取硬件和摄像机扫描硬件间存在的同步差或摄像机硬件本身时钟的不精确会造成一些不确定因素,直接影响像素行方向的等效间距值,因此引入尺度引子μ来表示这种不确定性。(www.chuimin.cn)

2.两步法原理

摄像机标定中,需标定外部参数6个,内部参数6个,即u0v0SxSyfk。其中,u0v0SxSy通过预标定确定。采用两步标定法,可以先求解外部6个参数,再计算2个内部参数。

978-7-111-34687-6-Chapter05-162.jpg

图5.16 摄像机成像示意图

根据摄像机成像模型,理想的透视投影关系如图5.16所示。空间点PXYZ)在摄像机像平面上的理想像素点为p′x′y′),由于摄像机镜头的径向畸变,实际像素点为p*x*y*)。它与PXYZ)之间不符合投影透视关系。

由图5.16可以看出,978-7-111-34687-6-Chapter05-163.jpg978-7-111-34687-6-Chapter05-164.jpg的方向一致,978-7-111-34687-6-Chapter05-165.jpg是摄像机光轴上的一点,其z坐标与物点在摄像机坐标系下的坐标值相同。

1)假设摄像机镜头畸变只是径向的,那么无论畸变如何变化,从图像中心点o′到图像点p*x*y*)的矢量978-7-111-34687-6-Chapter05-166.jpg的方向将保持不变,并且978-7-111-34687-6-Chapter05-167.jpg978-7-111-34687-6-Chapter05-168.jpg平行。

2)摄像机的等效焦距fx*y*产生同样的影响,f的大小不影响矢量978-7-111-34687-6-Chapter05-169.jpg的方向。

3)世界坐标系沿着x*轴和y*轴放置和平移,使得每一点有978-7-111-34687-6-Chapter05-170.jpg平行于978-7-111-34687-6-Chapter05-171.jpg,对x*y*产生的影响相同,从而o978-7-111-34687-6-Chapter05-172.jpg的方向保持不变。

4)在每一点处矢量978-7-111-34687-6-Chapter05-173.jpg978-7-111-34687-6-Chapter05-174.jpg的平行约束条件,与径向畸变表达式的系数、等效焦距f、平移分量Tz无关。而RTxTy受该约束制约。

因此,运用径向排列约束,存在如下关系:

978-7-111-34687-6-Chapter05-175.jpg

978-7-111-34687-6-Chapter05-176.jpg根据坐标转换公式,则

978-7-111-34687-6-Chapter05-177.jpg

根据RAC含义,有

978-7-111-34687-6-Chapter05-178.jpg

整理得

978-7-111-34687-6-Chapter05-179.jpg

行矢量978-7-111-34687-6-Chapter05-180.jpg是已知的,由标定模板和图像处理确定,列矢量978-7-111-34687-6-Chapter05-181.jpg是待求参数。

对每一标定模板目标点Pi,已知其XYZx*y*就可以列出如上的一个方程。R有9个参数,但是其正交性规定了R仅有三个自由度,再加上TxTy正好五个独立变量,故上式可惟一确定旋转矩阵RTxTy

求出RTxTy后,再利用最优化方法求解有效焦距fTTz分量和透镜畸变系数k