发明名称 基于摄影测量技术的数码相机投影中心位置精确标定方法
摘要 本发明涉及基于摄影测量技术的数码相机投影中心位置精确标定方法,有效解决解算精度高、可操作性强的问题,方法是,布设回光反射测量标志、控制点标志、编码标志、定向靶和基准尺,构成标定场;距离标定场3米处设置电子经纬仪和待检定相机,相机镜头上粘贴回光反射测量标志;相机在不同水平位置对标定场拍照,得到多张照片,用待检定相机在标定场前的不同位置以不同姿态对标定场拍照,与多张照片提取控制点标志的数据,形成标校数据;对控制点标志中心像点的坐标进行改正,解算得到3个投影中心位置,本发明可一体化实现控制点标志图像中心量测、高精度相机标校与投影中心精确测定三项工作,投影中心的位置标定精度优于5mm。
申请公布号 CN105424058A 申请公布日期 2016.03.23
申请号 CN201510749536.8 申请日期 2015.11.06
申请人 中国人民解放军信息工程大学 发明人 丛佃伟;许其凤;吕志伟;李军正;张同成
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 郑州天阳专利事务所(普通合伙) 41113 代理人 聂孟民
主权项 一种基于摄影测量技术的数码相机投影中心位置精确标定方法,其特征在于,包括以下步骤:(1)、在测量平面上布设至少200个回光反射测量标志,同时布置至少12个控制点标志和12个编码标志,在平面前设置定向靶和基准尺,构成相机投影中心标定场;(2)、距离标定场3米处,设置两台高精度电子经纬仪和待检定相机,两台电子经纬仪平行于标定场,相机镜头处粘贴回光反射标志;(3)、在标定场前不同水平位置利用待检定相机对标志场进行拍照,得到不同水平状态下的多张照片,同时在每个水平固定位置使用两台电子经纬仪测定相机镜头上粘贴的回光反射标志中心的坐标,用灰度值加权质心算法对获取的不同位置的回光反射标志中心的坐标进行计算,得到标志中心像点的坐标(x,y);(4)、用待检定相机在标定场前的不同位置,待检定相机至少9台,均匀分布在标志场前方,以不同姿态对标定场进行拍照,在每个位置至少拍摄三张照片,与步骤(3)中获取的多张照片一起,提取控制点标志的数据,将这些数据采用光束法自标定的方法对相机参数进行标校,得到标校参数,形成标校数据;(5)、根据步骤(4)中获取的标校数据,对步骤(3)中获取的照片中的控制点标志中心像点的坐标(x,y)进行改正,改正数(Δx′,Δy′)的计算公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open = '' close = '}'><mtable><mtr><mtd><mrow><msup><mi>&Delta;x</mi><mo>&prime;</mo></msup><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><mfrac><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>f</mi></mrow><mi>f</mi></mfrac><mo>+</mo><msub><mi>K</mi><mi>x</mi></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mn>4</mn></msup><mo>+</mo><msub><mi>K</mi><mi>3</mi></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mi>6</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msub><mi>P</mi><mn>1</mn></msub><mrow><mo>(</mo><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>P</mi><mn>2</mn></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mover><mi>y</mi><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>&Delta;y</mi><mo>&prime;</mo></msup><mo>=</mo><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>-</mo><mfrac><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo><mi>&Delta;</mi><mi>f</mi></mrow><mi>f</mi></mfrac><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mn>4</mn></msup><mo>+</mo><msub><mi>K</mi><mi>3</mi></msub><mover><mi>y</mi><mo>&OverBar;</mo></mover><msup><mi>r</mi><mi>6</mi></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>+</mo><msub><mi>P</mi><mn>2</mn></msub><mrow><mo>(</mo><msup><mi>r</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><msup><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>P</mi><mn>2</mn></msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>&CenterDot;</mo><mover><mi>y</mi><mo>&OverBar;</mo></mover></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub></mrow></mtd><mtd><mrow><mover><mi>y</mi><mo>&OverBar;</mo></mover><mo>=</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub></mrow></mtd><mtd><mrow><msup><mi>r</mi><mn>2</mn></msup><mo>=</mo><msup><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup><mo>+</mo><msup><mover><mi>y</mi><mo>&OverBar;</mo></mover><mn>2</mn></msup></mrow></mtd></mtr></mtable></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000840464070000011.GIF" wi="886" he="494" /></maths>改正后的控制点标志中心像点的坐标(x′,y′)为<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = '' close = '}'><mtable><mtr><mtd><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mo>=</mo><mi>x</mi><mo>+</mo><msup><mi>&Delta;x</mi><mo>&prime;</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>y</mi><mo>&prime;</mo></msup><mo>=</mo><mi>y</mi><mo>+</mo><msup><mi>&Delta;y</mi><mo>&prime;</mo></msup></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000840464070000012.GIF" wi="283" he="174" /></maths>根据共线条件方程方法解算得到3个投影中心位置坐标(X<sub>S</sub>,Y<sub>S</sub>,Z<sub>S</sub>),共线条件方程为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open = '' close = '}'><mtable><mtr><mtd><mrow><msup><mi>x</mi><mo>&prime;</mo></msup><mo>=</mo><mo>-</mo><mi>f</mi><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>y</mi><mo>&prime;</mo></msup><mo>=</mo><mo>-</mo><mi>f</mi><mfrac><mrow><msub><mi>a</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>a</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>X</mi><mo>-</mo><msub><mi>X</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>b</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Y</mi><mo>-</mo><msub><mi>Y</mi><mi>S</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>Z</mi><mo>-</mo><msub><mi>Z</mi><mi>S</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000840464070000013.GIF" wi="819" he="295" /></maths>其中,f、x<sub>0</sub>、y<sub>0</sub>为相机内方位元素;K<sub>1</sub>、K<sub>2</sub>、K<sub>3</sub>为相机的径向畸变;P<sub>1</sub>、P<sub>2</sub>为相机的偏心畸变;b<sub>1</sub>、b<sub>2</sub>为相机的像平面畸变;a<sub>1</sub>、a<sub>2</sub>、a<sub>3</sub>、b<sub>1</sub>、b<sub>2</sub>、b<sub>3</sub>、c<sub>1</sub>、c<sub>2</sub>、c<sub>3</sub>为共线条件方程中旋转矩阵的元素;X、Y、Z为地面点在以摄站为原点的地辅系中的坐标;解算得到3个投影中心位置后,将步骤(3)中两台电子经纬仪获取的相机镜头上的回光反射标志中心的坐标归心到镜头中心后,与解算得到的投影中心位置进行比较,计算得到实际相机投影中心与镜头中心位置的差值,得到投影中心与镜头中心的位置关系,实现对数码相机中心位置的精确标定。
地址 450052 河南省郑州市陇海中路66号