主权项 |
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。 |