三、单目标定
从之前的
从世界坐标到像素坐标可以得到:
Zc⎣⎡uv1⎦⎤=⎣⎡dx1000dy10u0v01⎦⎤⎣⎡f000f0001000⎦⎤[R0T1]⎣⎢⎢⎡XwYwZw1⎦⎥⎥⎤
相机标定也是利用了上面这个公式,等式右边前两个矩阵乘积为相机内参,第三为相机外参(也成为Toc)。
各个参数的含义:dx代表一个像素的宽度(x方向),与x同单位,x/dx表示x轴上有多少个像素,同理y/dy表示y轴上的像素个数,(u0,v0)是图像平面中心。相机模型类似小孔成像模型,自行百度,也可参考前面的坐标转换关系相机坐标系与图像坐标系与世界坐标系与相机坐标系。
举个实际的例子,相机内参为:
K=⎣⎡dxf000dyf0u0v01⎦⎤
1、如何进行标定?(理想情况下无畸变)
单目标定的目的就是为了得到dx,dy,f,R,T这几个参数(暂时不考虑畸变参数)。而最开始接触标定只是按照不同姿势地摆放标定板,前提当然要固定相机位置,为什么要这么做呢。
设有N个角点,K个棋盘,则可以列出2NK个方程,未知参数有4+6K个,所以从理论上来说,方程组数大于未知数个数就能解出未知数,即2NK>=6K+4,即(N-3)K>=2,但对于同一个平面来说,4个角点确定了一个与投影平面共面的四边形,在四个方向同时延伸就可以变成任意四边形,所以N的约束即为4,那么就只有K>=2这个条件,理论上至少需要两个视场列得到的2×4×2=16个方程才能求得4+6×2=16个参数。但是考虑噪声和数值稳定性,图片上的角点可以取多后采用最小二乘法、ransac或其他方法取最优参数。一般标定取10幅7×7或7×8角度不同的图片。到此为止,求出了dx,dy,f,R,T这几个参数。
当然,标定还有另外一种方式,对于eye in hand 这种case,可以固定标定板,相机在机械手上,随着机械手变化位置,去不同的点位拍照,(当硬件结构确定情况下的大批量生产情况采用这种自动标定方式可以节约人力成本,)提前预设好几个位置,让机械手随机在附近进行拍照,完成内外参的标定。
2、实际情况下(有畸变)
理想的针孔模型透过小孔的光线少,导致了相机的曝光时间慢,实际中为了加快图像的生成使用透镜,但同时也引入了畸变。所以畸变发生在相机内部,所需要增加的相应计算多了几个畸变参数(会在下面定义),计算方法参照的方程式仍然是从世界坐标系到像素坐标的转换Tpw,常见的畸变有(1)径向畸变,(2)切向畸变。
2.1、径向畸变
径向畸变的现象有桶形失真与枕形失真。在成像仪光轴中心的畸变为0,沿着镜头半径方向向边缘移动,径向位移增大,畸变越来越严重。畸变的数学模型用参数k1,k2,k3(或更多)描述,用主点的泰勒级数展开式前三项(或更多)表示:
{xcorrect=x(1+k1r2+k2r4+k3r6)ycorrect=y(1+k1r2+k2r4+k3r6)
2.2、切向畸变
切向畸变是由于透镜的安装偏差产生,导致透镜本身与相机成像平面不平行。切向畸变模型用参数p1,p2描述:
{xcorrect=x+[2p1xy+p2(r2+2x2)]ycorrect=y+[2p2xy+p1(r2+2y2)]
(r为该点距离成像中心距离)
2.3、总结
综合以上两种畸变,得到畸变矫正前后的坐标对应(理想图像坐标到实际图像坐标)关系:
[xcyc]=(1+k1r2+k2r4+k3r6)[xpyp]+[2p1y+p2(r2+2x2)2p2y+p1(r2+2y2)]
此时,世界坐标系到像素坐标系的关系就多了一个畸变转换:
Zc⎣⎡uv1⎦⎤=⎣⎡dx1000dy10u0v01⎦⎤[distortransform]⎣⎡f000f0001000⎦⎤[R0T1]⎣⎢⎢⎡XwYwZw1⎦⎥⎥⎤
畸变转换至少包含五个参数k1,k2,k3,p1,p2,标定需要求得的未知数就有4+5+6K个,所以,为了标定结果的准确性,前面提到的标定取10幅7×7或7×8角度不同的图片是很有必要的。
3、疑问
3.1、标定的世界坐标从何而来
给予标定板的建立世界坐标的xy轴,z轴垂直标定板射出,原点位置定义为标定板的圆孔的第一个点。
3.2、一张图像也能得到标定结果
N的约束4还要保留吗?
一张图片确实可以得到标定结果,会有很多解。
4、自动标定
由于标定过程比较繁杂,无非就是尽可能多的在工作平面附近拍摄标定板的图片进行计算。使用机械手自动走点位进行拍照的方法较为节省人力成本。
本次分析参考:自动标定分析
eye in hand
原始位置:
rTo=rTt⋅tTc⋅cTo
新位置:
rTo=rTt′⋅tTc⋅cTo′
cTo′=cTo⋅T
结合以上两个公式,移项化简得:
rTt′=rTo⋅T−1⋅oTt
eye to hand
与eye in hand同理
原始位置:
rTc=rTt⋅tTo⋅oTc
新位置:
rTc=rTt′⋅tTo⋅oTc′
oTc′=T−1⋅oTc
结合以上两个公式,移项化简得:
rTt′=rTo⋅T⋅oTt
以上的T均假设移动的object的偏移量。