发明名称 基于任意已知运动的单目摄像机标定方法
摘要 本发明基于任意已知运动的单目摄像机标定方法属于计算机视觉检测和图像处理领域,涉及大型锻件尺寸测量系统中单目摄像机的标定方法,具体是指一种基于任意已知运动的单目摄像机标定方法。标定方法利用四组已知运动获取12张图像,然后根据12张图像求取摄像机参数,包括主点坐标和尺度因子;首先,利用三坐标测量仪的实现靶标的已知运动获取靶标图像;再将畸变校正后的12张图像分组进行消隐点估计;采用最优化方法获取消隐点和“拟合直线簇”,标定单目摄像机。本发明基于极大似然估计的消隐点鲁棒性算法,提高消隐点估计的稳定性,是一种简单易实现的标定方法,特别适用于快速准确实现实验现场的标定。
申请公布号 CN104751452A 申请公布日期 2015.07.01
申请号 CN201510101366.2 申请日期 2015.03.08
申请人 大连理工大学 发明人 贾振元;刘巍;樊超楠;杨景豪;王灵丽;刘阳
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 关慧贞
主权项 一种基于任意已知运动的单目摄像机标定方法,其特征是,标定方法利用四组已知运动获取12张图像,然后根据12张图像求取摄像机参数,包括主点坐标(u<sub>0</sub>,v<sub>0</sub>)、尺度因子f<sub>x</sub>、f<sub>y</sub>;具体步骤如下:步骤1:利用三坐标测量仪的实现靶标的已知运动获取靶标图像控制三坐标测量仪带动靶标的做四组运动,每组运动含有X方向和Y方向两次纯平移,纯平移运动方向位移为D<sub>0</sub>,摄像机在每次运动中停留三个位置,其中摄像机光心位置分别为O<sub>0</sub>,O<sub>1</sub>,O<sub>2</sub>,进而得到三张图像I<sub>0</sub>,I<sub>1</sub>,I<sub>2</sub>,而且在确保靶标在视场内的前提下,每组正交运动之间摄像机均有不小于5度的俯仰角度或扫视角度,通过四组的两次纯平移运动得到共12张靶标图像;根据三坐标测量仪运动位置的三维坐标计算每组运动中两次纯平移运动的夹角θ<sub>i</sub>,i=1,2,3,4;步骤2:消隐点的鲁棒估计对纯平移运动前、后摄取的两幅图像中的特征点的像点进行畸变校正后,采取由交于同一点的“拟合直线簇”来确定消隐点和对应像点连线的最大似然估计;设定消隐点坐标为v=(u<sub>x</sub>,v<sub>x</sub>),则经过消隐点的“拟合直线簇”l<sub>i</sub>的方程为:y‑v<sub>x</sub>=k<sub>j</sub>(x‑u<sub>x</sub>)                    (1)其中,k<sub>j</sub>为对应像点A<sub>j</sub>=(u<sub>j</sub>,v<sub>j</sub>)、A′<sub>j</sub>=(u′<sub>j</sub>,v′<sub>j</sub>)连线的最大似然估计(即“拟合直线”)所对应的斜率,则对应像点到“拟合直线”的距离之和表示如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>d</mi><mrow><mo>(</mo><msub><mi>A</mi><mi>j</mi></msub><mo>,</mo><msub><mi>l</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>+</mo><mi>d</mi><mrow><mo>(</mo><msubsup><mi>A</mi><mi>j</mi><mo>&prime;</mo></msubsup><mo>,</mo><msub><mi>l</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mo>|</mo><msub><mi>k</mi><mi>j</mi></msub><msub><mi>u</mi><mi>j</mi></msub><mo>-</mo><msub><mi>v</mi><mi>j</mi></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>k</mi><mi>j</mi></msub><msub><mi>u</mi><mi>x</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><msqrt><msubsup><mi>k</mi><mi>j</mi><mn>2</mn></msubsup><mo>+</mo><mn>1</mn></msqrt></mfrac><mo>+</mo><mfrac><mrow><mo>|</mo><msub><mi>k</mi><mi>j</mi></msub><msubsup><mi>u</mi><mi>j</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>v</mi><mi>j</mi><mo>&prime;</mo></msubsup><mo>+</mo><mrow><mo>(</mo><msub><mi>v</mi><mi>x</mi></msub><mo>-</mo><msub><mi>k</mi><mi>j</mi></msub><msub><mi>u</mi><mi>x</mi></msub><mo>)</mo></mrow><mo>|</mo></mrow><msqrt><msubsup><mi>k</mi><mi>j</mi><mn>2</mn></msubsup><mo>+</mo><mn>1</mn></msqrt></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000678880910000011.GIF" wi="1651" he="188" /></maths>采用最优化方法获取消隐点和“拟合直线簇”,最优化目标为最小化所有像点到相应“拟合直线”距离之和,见式(3);优化变量为消隐点坐标分量以及各“拟合直线”的斜率;<img file="FDA0000678880910000021.GIF" wi="1358" he="142" />其中,n为“拟合直线簇”中直线的数目;采用该方法分别计算4组运动的消隐点;步骤3:标定单目摄像机采用步骤2的方法,将畸变校正后的12张图像分组进行消隐点估计,每组3张图像中对应点连线得到2个消隐点,4组共8个消隐点;将得到的8个消隐点坐标、主点坐标以及两平移运动度带入下式中:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>cos</mi><msub><mi>&theta;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>T</mi></msubsup><mi>w</mi><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></mrow><mrow><msqrt><msubsup><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow><mi>T</mi></msubsup><mi>w</mi><msub><mi>v</mi><mrow><mi>i</mi><mn>1</mn></mrow></msub></msqrt><msqrt><msubsup><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow><mi>T</mi></msubsup><mi>w</mi><msub><mi>v</mi><mrow><mi>i</mi><mn>2</mn></mrow></msub></msqrt></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3,4</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000678880910000022.GIF" wi="1491" he="199" /></maths>其中,θ为第i组运动中两平移运动之间的角度,v<sub>i1</sub>、v<sub>i2</sub>为第i组正交运动下的两次平移运动获取的消隐点对的像素坐标(齐次形式),w=K<sup>‑T</sup>K<sup>‑1</sup>是绝对二次曲线的图像,其具体形式如下:<img file="FDA0000678880910000023.GIF" wi="1474" he="234" />通过公式(4)、(5)以及八个消隐点坐标得到包含四个非线性方程的二元方程组,用非线性最小二乘法迭代求解出尺度因子f<sub>x</sub>和f<sub>y</sub>和主点坐标(u<sub>0</sub>,v<sub>0</sub>);因此,得到单目摄像机参数主点坐标(u<sub>0</sub>,v<sub>0</sub>)、尺度因子f<sub>x</sub>、f<sub>y</sub>,最终完成标定过程。
地址 116024 辽宁省大连市甘井子区凌工路2号