发明名称 一种摄像机分步自标定方法
摘要 本发明提供一种摄像机分步自标定方法,包括:第1阶段,对同一场景拍摄两张在不同焦距下的图像,对图像进行特征提取和匹配并剔除误匹配点,利用匹配点对求解摄像机的主点坐标;第2阶段,从不同角度对同一场景进行拍摄获得三幅可供匹配的图像,对图像进行特征提取和匹配并剔除误匹配点,然后基于Kruppa方程,将已获得的主点坐标代入方程实现倾斜因子以及摄像机在成像平面u、v轴方向上的尺度因子三个参数的求解。通过本发明的方法,摄像机的主点坐标的标定精度较高,并减小了Kruppa方程系数矩阵规模,降低了求解运算量,具有实时性。该方法适用于视觉系统的摄像机标定,可在三维测量、三维重建、机器导航以及增强现实等领域使用。
申请公布号 CN102750704B 申请公布日期 2015.04.22
申请号 CN201210219658.2 申请日期 2012.06.29
申请人 吉林大学 发明人 王欣;高焕玉
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 代理人
主权项 一种摄像机分步自标定方法,所述摄像机采取的成像矩阵为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>K</mi><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>f</mi><mi>u</mi></msub></mtd><mtd><mi>s</mi></mtd><mtd><msub><mi>u</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>f</mi><mi>v</mi></msub></mtd><mtd><msub><mi>v</mi><mn>0</mn></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FSB0000134362570000011.GIF" wi="363" he="182" /></maths>其中,(u<sub>0</sub>,v<sub>0</sub>)是摄像机主点坐标,s是倾斜因子,f<sub>u</sub>、f<sub>v</sub>分别为摄像机在成像平面u轴、v轴方向上的尺度因子;其特征在于,所述方法分为2个阶段进行:第1阶段,标定摄像机的主点坐标(u<sub>0</sub>,v<sub>0</sub>);第2阶段,结合第1阶段获得的主点坐标,标定摄像机其他三个参数,即倾斜因子s以及摄像机在成像平面u轴、v轴方向上的尺度因子f<sub>u</sub>、f<sub>v</sub>;其中在第1阶段中,通过调节摄像机的焦距,同时保持摄像机的位姿状态不变,对同一场景拍摄两张在不同焦距下的图像,对图像进行特征提取和匹配,并剔除误匹配点,最后利用匹配点对求解摄像机的主点坐标;在第2阶段中,控制摄像机的内参数不变,从不同角度对同一场景进行拍摄获得三幅的可供匹配的图像,对图像进行特征提取和匹配,并剔除误匹配点,然后基于Kruppa方程,确定Kruppa方程的比例因子初始值并利用遗传算法对其进行优化,通过优化后的比例因子以及第1阶段获得的主点坐标代入Kruppa方程实现摄像机其他三个参数的求解;其中所述第1阶段的具体流程包括:步骤1.1:保持摄像机的位姿状态不变,调整摄像机的焦距,获得摄像机对同一场景分别在焦距f<sub>1</sub>和f<sub>2</sub>条件下的两张图像;步骤1.2:利用SURF算法对获取的两张图像进行特征提取,然后对两图像所含的特征点进行初匹配操作,得到两图像特征点的匹配对集合;步骤1.3:采用对极几何约束的RANSAC算法对步骤1.2中得到的匹配对集合进行误匹配剔除操作,将误匹配点对从匹配点对集合中去除,然后将剩余的匹配点对集合保存起来;步骤1.4:判断步骤1.3中所获得的匹配点对数P是否大于第一阈值T1,若大于T1,则进行步骤1.5操作,否则转入步骤1.1,其中T1≥8;步骤1.5:将步骤1.4所获得的P个匹配特征点对写成如下的方程组形式:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mrow><mo>(</mo><msup><msub><mi>y</mi><mn>1</mn></msub><mo>&prime;</mo></msup><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>u</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msup><msub><mi>x</mi><mn>1</mn></msub><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>v</mi><mn>0</mn></msub><mo>=</mo><msub><mi>x</mi><mn>1</mn></msub><msup><msub><mi>y</mi><mn>1</mn></msub><mo>&prime;</mo></msup><mo>-</mo><msup><msub><mi>x</mi><mn>1</mn></msub><mo>&prime;</mo></msup><msub><mi>y</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msup><msub><mi>y</mi><mn>2</mn></msub><mo>&prime;</mo></msup><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>u</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msup><msub><mi>x</mi><mn>2</mn></msub><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>v</mi><mn>0</mn></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><msup><msub><mi>y</mi><mn>2</mn></msub><mo>&prime;</mo></msup><mo>-</mo><msup><msub><mi>x</mi><mn>2</mn></msub><mo>&prime;</mo></msup><msub><mi>y</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><mrow><mo>(</mo><msup><msub><mi>y</mi><mi>P</mi></msub><mo>&prime;</mo></msup><mo>-</mo><msub><mi>y</mi><mi>P</mi></msub><mo>)</mo></mrow><msub><mi>u</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>x</mi><mi>P</mi></msub><mo>-</mo><msup><msub><mi>x</mi><mi>P</mi></msub><mo>&prime;</mo></msup><mo>)</mo></mrow><msub><mi>v</mi><mn>0</mn></msub><mo>=</mo><msub><mi>x</mi><mi>P</mi></msub><msup><msub><mi>y</mi><mi>P</mi></msub><mo>&prime;</mo></msup><mo>-</mo><msup><msub><mi>x</mi><mi>P</mi></msub><mo>&prime;</mo></msup><msub><mi>y</mi><mi>P</mi></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FSB0000134362570000021.GIF" wi="970" he="372" /></maths>利用最小二乘法求解该方程组,得到的解(u<sub>0</sub>,v<sub>0</sub>)便是摄像机的主点坐标,其中,(x<sub>i</sub>,y<sub>i</sub>)和(x′<sub>i</sub>,y′<sub>i</sub>)是两图像的匹配特征点对集合中第i对特征点坐标,i=1,2,…,P。
地址 130012 吉林省长春市前进大街2699号