发明名称 补偿摄影测量影像畸变的数字畸变模型的生成方法
摘要 补偿摄影测量影像畸变的数字畸变模型的生成方法,包括以下步骤:①建立控制场;②拍摄控制场并测定影像外方位元素;③量测各控制标志在影像上的实际位置;④选择像幅四个角隅标志,令其畸变为零;在控制场上,对此四个角隅标志引入投影差;对像幅内所有控制标志引入投影差,实现控制场至真二维控制场的变换;⑤计算像幅内各控制标志在影像上的理论位置;⑥通过上述实际位置和理论位置,确定影像上各控制标志的总畸变;⑦依据影像上与待测像素最临近的四个标志的总畸变,内插获取影像上每个像素的总畸变;从而构成该影像的数字畸变模型。本发明不仅可以方便地对不同成像机理的相机的影像畸变补偿,而且可获得相应的内方位元素,以实现继后的三维测量。
申请公布号 CN1645054A 申请公布日期 2005.07.27
申请号 CN200510018151.0 申请日期 2005.01.13
申请人 武汉大学 发明人 冯文灏;陈大为;商浩亮;侯文广;李欣
分类号 G01C11/00;G01B11/00;G03B37/06 主分类号 G01C11/00
代理机构 武汉天力专利事务所 代理人 程祥;冯卫平
主权项 1.补偿摄影测量影像畸变的数字畸变模型的生成方法,其特征是:包括以下步骤:①建立控制场;②用检定相机拍摄控制场;测定影像外方位元素;③量测各控制标志在影像上的实际位置;④选择像幅四个角隅标志,令其畸变为零;在控制场上,确定通过此四个角隅标志的拟合平面,并对此四个角隅标志引入投影差;对像幅内所有控制标志引入投影差,实现控制场至真二维控制场的变换;⑤依据上述四个角隅标志,按二维直接线性变换原理,解求8个变换系数,计算像幅内各控制标志在影像上的理论位置;⑥通过各控制标志在影像上的实际位置和理论位置,确定影像上各控制标志的总畸变(f1,f2,f3,f4);⑦获取影像上每个像素的总畸变:依影像上与待测像素最临近的四个标志的总畸变,按下述任意分布数据点关系式,内插获取该像素的总畸变f(x,y),直至获取影像上每个像素的总畸变;从而构成该影像的数字畸变模型;<math> <mrow> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> </mrow> </math> <math> <mrow> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mfrac> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> </mrow> </math> <math> <mrow> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mfrac> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mo>)</mo> </mrow> </mrow> </math> <math> <mrow> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>=</mo> <msub> <mi>f</mi> <mn>3</mn> </msub> <mo>+</mo> <mfrac> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mfrac> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mn>4</mn> </msub> <mo>-</mo> <msub> <mi>f</mi> <mn>3</mn> </msub> <mo>)</mo> </mrow> </mrow> </math> <math> <mrow> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </msub> <mo>=</mo> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>+</mo> <mfrac> <msqrt> <msup> <mrow> <mo>(</mo> <mi>x</mi> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mi>y</mi> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <msqrt> <msup> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>x</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>y</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mfrac> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mn>3,4</mn> <mo>)</mo> </mrow> </msub> <mo>-</mo> <msub> <mi>f</mi> <mrow> <mo>(</mo> <mn>1,2</mn> <mo>)</mo> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </math> 式中x、y为待测像素的横、纵坐标;(x1、y1)、(x2、y2)、(x3、y3)、(x4、y4)分别为待测像素最临近的四个标志的横、纵坐标;x(1,2)、y(1,2)为(x1、y1)和(x2、y2)的连线上的点的横、纵坐标;x(3,4)、y(3,4)为(x3、y3)和(x4、y4)的连线上的点的横、纵坐标;其中x(1,2)=x、x(3,4)=x。
地址 430072湖北省武汉市武昌珞珈山
您可能感兴趣的专利