发明名称 基于正交补偿的大视场摄像机标定方法
摘要 本发明公开了一种基于正交补偿的大视场摄像机标定方法,属于计算机视觉检测和图像处理领域,涉及大型锻件尺寸测量系统中摄像机内、外参数的标定方法。通常主动视觉标定过程中电控平台两平移导轨并非完全正交,而标定只是盲目选择精度高的电控平台近似认为其正交,因此提出一种基于正交补偿的摄像机标定方法通过分别求取主点和尺度因子、平移向量和旋转矩阵获取摄像机内、外参数。本发明考虑在标定过程中加入导轨正交度因素,即引入正交补偿,以降低导轨正交度对标定结果带来的误差。
申请公布号 CN104167001B 申请公布日期 2017.02.15
申请号 CN201410432768.6 申请日期 2014.08.27
申请人 大连理工大学 发明人 贾振元;刘巍;王亮;刘阳;杨景豪;王灵丽;张弛
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 关慧贞;梅洪玉
主权项 一种基于正交补偿的标定方法,其特征在于如下步骤:步骤1:变焦距法求主点用摄像机在不同焦距下拍摄二维靶标,当焦距改变时,光心沿光轴移动,在摄像机视场中二维靶标将会等比例扩缩,得到不同视场下靶标图像,不同图像中有一点的图像保持不变,叫做扩缩中心,在不同视场的图像中扩缩中心即为主点,因此将靶标各个特征点在不同焦距下的连线的最小二乘意义下的交点即为主点;步骤2:利用四维电控平台的主动视觉标定方法获取靶标图像控制四维电控平台带动摄像机做四组平面内的正交运动,每组正交运动含有两次平移,摄像机在每次正交运动中停留三个位置,进而得到三张图像,而且在确保靶标在视场内的前提下,每组正交运动之间摄像机均有不小于5度的俯仰角度或扫视角度,通过四组的两次纯平移运动得到共12张靶标图像;步骤3:求取尺度因子平行直线与无穷远平面相交于同一个无穷远点,即隐消点,将畸变校正后的12张图像分组进行消隐点估计,每组3张图像中对应点连线得到2个消隐点,4组共8个消隐点;将得到的8个消隐点坐标、主点坐标以及导轨间角度带入下式中:<maths num="0001"><math><![CDATA[<mrow><mi>cos</mi><mi>&theta;</mi><mo>=</mo><mfrac><mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>T</mi></msubsup><msub><mi>wv</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mrow><mrow><msqrt><mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>T</mi></msubsup><msub><mi>wv</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></mrow></msqrt><msqrt><mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>T</mi></msubsup><msub><mi>wv</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mrow></msqrt></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo><mn>4</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001135871660000011.GIF" wi="1126" he="191" /></maths>其中,θ为电控平台两根平移导轨之间的角度,v<sub>i1</sub>、v<sub>i2</sub>为第i组正交运动下的两次平移运动获取的消隐点对的像素坐标,w=K<sup>‑T</sup>K<sup>‑1</sup>是绝对二次曲线的图像,其具体形式如下:<maths num="0002"><math><![CDATA[<mrow><mi>w</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mn>1</mn><mo>/</mo><msubsup><mi>f</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mtd><mtd><mn>0</mn></mtd><mtd><mrow><mo>-</mo><msub><mi>u</mi><mn>0</mn></msub><mo>/</mo><msubsup><mi>f</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mn>1</mn><mo>/</mo><msubsup><mi>f</mi><mi>y</mi><mn>2</mn></msubsup></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>/</mo><msubsup><mi>f</mi><mi>y</mi><mn>2</mn></msubsup></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msub><mi>u</mi><mn>0</mn></msub><mo>/</mo><msubsup><mi>f</mi><mi>x</mi><mn>2</mn></msubsup></mrow></mtd><mtd><mrow><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>/</mo><msubsup><mi>f</mi><mi>y</mi><mn>2</mn></msubsup></mrow></mtd><mtd><mrow><msubsup><mi>u</mi><mn>0</mn><mn>2</mn></msubsup><mo>/</mo><msubsup><mi>f</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>v</mi><mn>0</mn><mn>2</mn></msubsup><mo>/</mo><msubsup><mi>f</mi><mi>y</mi><mn>2</mn></msubsup><mo>+</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001135871660000021.GIF" wi="1046" he="231" /></maths>通过上述两个式子以及八个消隐点坐标得到包含四个非线性方程的二元方程组,用非线性最小二乘法迭代求解出尺度因子f<sub>x</sub>和f<sub>y</sub>,其中,u<sub>0</sub>,v<sub>0</sub>为摄像机主点坐标;步骤4:标定外参数;一组对应像点的齐次坐标即可获得一个线性方程,如下:uu′f<sub>11</sub>+vu′f<sub>12</sub>+u′f<sub>13</sub>+uv′f<sub>21</sub>+vv′f<sub>22</sub>+v′f<sub>23</sub>+uf<sub>31</sub>+vf<sub>32</sub>+f<sub>33</sub>=0 (3)其中,f<sub>ij</sub>为基本矩阵的第i行、第j列上的元素;[uu′,vu′,u′,uv′,vv′,v′,u,v,1]f=0 (4)其中,f=[f<sub>11</sub>,f<sub>12</sub>,f<sub>13</sub>,f<sub>21</sub>,f<sub>22</sub>,f<sub>23</sub>,f<sub>31</sub>,f<sub>32</sub>,f<sub>33</sub>]<sup>T</sup>为9维列向量;如果给出n组对应点{m<sub>i</sub>,m′<sub>i</sub>},则有如下齐次线性方程组<maths num="0003"><math><![CDATA[<mrow><mi>A</mi><mi>f</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><msub><mi>u</mi><mn>1</mn></msub><msubsup><mi>u</mi><mn>1</mn><mo>&prime;</mo></msubsup></mrow></mtd><mtd><mrow><msub><mi>v</mi><mn>1</mn></msub><msubsup><mi>u</mi><mn>1</mn><mo>&prime;</mo></msubsup></mrow></mtd><mtd><msubsup><mi>u</mi><mn>1</mn><mo>&prime;</mo></msubsup></mtd><mtd><mrow><msub><mi>u</mi><mn>1</mn></msub><msubsup><mi>v</mi><mn>1</mn><mo>&prime;</mo></msubsup></mrow></mtd><mtd><mrow><msub><mi>v</mi><mn>1</mn></msub><msubsup><mi>v</mi><mn>1</mn><mo>&prime;</mo></msubsup></mrow></mtd><mtd><msubsup><mi>v</mi><mn>1</mn><mo>&prime;</mo></msubsup></mtd><mtd><msub><mi>u</mi><mn>1</mn></msub></mtd><mtd><msub><mi>v</mi><mn>1</mn></msub></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mrow><msub><mi>u</mi><mi>n</mi></msub><msubsup><mi>u</mi><mi>n</mi><mo>&prime;</mo></msubsup></mrow></mtd><mtd><mrow><msub><mi>v</mi><mi>n</mi></msub><msubsup><mi>u</mi><mi>n</mi><mo>&prime;</mo></msubsup></mrow></mtd><mtd><msubsup><mi>u</mi><mi>n</mi><mo>&prime;</mo></msubsup></mtd><mtd><mrow><msub><mi>u</mi><mi>n</mi></msub><msubsup><mi>v</mi><mi>n</mi><mo>&prime;</mo></msubsup></mrow></mtd><mtd><mrow><msub><mi>v</mi><mi>n</mi></msub><msubsup><mi>v</mi><mi>n</mi><mo>&prime;</mo></msubsup></mrow></mtd><mtd><msubsup><mi>v</mi><mi>n</mi><mo>&prime;</mo></msubsup></mtd><mtd><msub><mi>u</mi><mi>n</mi></msub></mtd><mtd><msub><mi>v</mi><mi>n</mi></msub></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mi>f</mi><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001135871660000022.GIF" wi="1302" he="222" /></maths>给定至少8组对应像点即可在相差一未知非零比例因子下线性求解上述齐次线性方程组,在此施加约束||f||=1,当对应像点精确匹配时,仅利用8组对应点求取基本矩阵;进而获取相差一个非零比例因子下的本质矩阵<img file="FDA0001135871660000023.GIF" wi="234" he="71" />其中<img file="FDA0001135871660000024.GIF" wi="141" he="55" />s为相差的比例因子;利用E<sup>*</sup>=Udiag((σ<sub>1</sub>+σ<sub>2</sub>)/2,(σ<sub>1</sub>+σ<sub>2</sub>)/2,0)V<sup>T</sup>作为<img file="FDA0001135871660000031.GIF" wi="515" he="71" />的近似,其中特征值σ<sub>1</sub>≥σ<sub>2</sub>≥σ<sub>3</sub>;从本质矩阵E<sup>*</sup>估计出摄像机间的外部参数<img file="FDA0001135871660000032.GIF" wi="121" he="62" />具有下述4种可能的组合形式:[UWV<sup>T</sup>|u<sub>3</sub>];[UWV<sup>T</sup>|‑u<sub>3</sub>];[UW<sup>T</sup>V<sup>T</sup>|u<sub>3</sub>];[UW<sup>T</sup>V<sup>T</sup>|‑u<sub>3</sub>]. (6)其中,<img file="FDA0001135871660000033.GIF" wi="374" he="222" />u<sub>3</sub>为U的第3列向量,<img file="FDA0001135871660000034.GIF" wi="151" he="76" />通过Chirality约束,依据偏航角为正锐角和平移向量<img file="FDA0001135871660000035.GIF" wi="27" he="55" />的X方向分量<img file="FDA0001135871660000036.GIF" wi="40" he="63" />为负值来判断正确参数组合;利用获得的摄像机内部参数和相差一比例因子下的摄像机间外部参数重建距离信息精确已知的两个特征点,重建距离与实际距离的比值即为比例因子s。
地址 116024 辽宁省大连市甘井子区凌工路2号