发明名称 一种透镜畸变参数的标定方法
摘要 本发明属于测量技术领域,涉及对透镜畸变参数标定方法的改进。本发明的标定步骤是:放置靶标-(2)拍摄靶标图象(3)获取特征点的坐标-计算特征点的交比值CR-获取x<SUB>ia</SUB>,x<SUB>ib</SUB>,x<SUB>ic</SUB>,x<SUB>id</SUB>,y<SUB>ia</SUB>,y<SUB>ib</SUB>,y<SUB>ic</SUB>,y<SUB>id</SUB>的表达式—求取畸变参数k<SUB>1</SUB>,k<SUB>2</SUB>。本发明方法算法简单,需要的已知条件少,容易实现,标定精度高,校正效果好。
申请公布号 CN1508527A 申请公布日期 2004.06.30
申请号 CN02156701.8 申请日期 2002.12.17
申请人 北京航空航天大学 发明人 张广军;贺俊吉;杨宪铭
分类号 G01M11/02 主分类号 G01M11/02
代理机构 中国航空专利中心 代理人 梁瑞林
主权项 1、一种透镜畸变参数的标定方法,其特征在于,利用透视投影成像中交比不变原理来进行镜头畸变参数的标定,具体的步骤是:(1)放置靶标,在被标定的镜头前放置一个靶标,所使用的靶标为平面形靶标,要求靶标平面上的图案中至少有一条直线,该直线上至少有4个已知其间距的特征点;(2)拍摄靶标图象,由被标定镜头拍摄靶标的一幅图像,尽量使靶标图像充满画面;(3)获取特征点的坐标,在靶标图象中选择位于一条直线的至少4个特征点,使用图像处理算法,索贝尔(Sobel)算子,提取图像中有关的横向和纵向直线,然后计算这些横线和纵线的交叉点的坐标,最后得到位于同一条直线的至少4个特征点的坐标;(4)计算特征点的交比值CR,根据获取的位于同一直线上的4个特征点A、B、C、D的位置,和靶标图案的已知参数,利用下式:<maths num="001"><![CDATA[ <math><mrow><mrow><mi>CR</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>D</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>SR</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo>)</mo></mrow></mrow><mrow><mi>SR</mi><mrow><mo>(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>D</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mi>AC</mi><mi>BC</mi></mfrac><mfrac><mi>BD</mi><mi>AD</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo></mrow><mo>&lsqb;</mo><mn>1</mn><mo>&rsqb;</mo></mrow></math>]]></maths>计算得到上述4个特征点的交比值CR;(5)获取x<sub>ia</sub>,x<sub>ib</sub>,x<sub>ic</sub>,x<sub>id</sub>,y<sub>ia</sub>,y<sub>ib</sub>,y<sub>ic</sub>,y<sub>id</sub>的表达式,利用前面获得的4个特征点A、B、C、D的图像坐标,分别记为A(x<sub>a</sub>,y<sub>a</sub>),B(x<sub>b</sub>,y<sub>b</sub>),C(x<sub>c</sub>,y<sub>c</sub>),D(x<sub>d</sub>,y<sub>d</sub>),代入下式:<img file="A0215670100022.GIF" wi="985" he="594" />得到x<sub>ia</sub>,x<sub>ib</sub>,x<sub>ic</sub>,x<sub>id</sub>,y<sub>ia</sub>,y<sub>ib</sub>,y<sub>ic</sub>,y<sub>id</sub>的表达式,(x<sub>ia</sub>,x<sub>ib</sub>,x<sub>ic</sub>,x<sub>id</sub>,y<sub>ia</sub>,y<sub>ib</sub>,y<sub>ic</sub>,y<sub>id</sub>分别表示A、B、C、D四点的未经畸变的图像坐标(横坐标和纵坐标));(6)将上面得到的x<sub>ia</sub>,x<sub>ib</sub>,x<sub>ic</sub>,x<sub>id</sub>,y<sub>ia</sub>,y<sub>ib</sub>,y<sub>ic</sub>,y<sub>id</sub>和CR代入下式:<img file="A0215670100031.GIF" wi="754" he="261" />得到两个关于畸变参数的k<sub>1</sub>,k<sub>2</sub>的一元二次方程,求解这两个方程,便得到畸变参数k<sub>1</sub>,k<sub>2</sub>的数值。
地址 100083北京市海淀区学院路37号