主权项 |
一种基于摄影测量技术的数码相机投影中心位置精确标定方法,其特征在于,包括以下步骤:(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>Δx</mi><mo>′</mo></msup><mo>=</mo><mo>-</mo><msub><mi>x</mi><mn>0</mn></msub><mo>-</mo><mfrac><mrow><mover><mi>x</mi><mo>‾</mo></mover><mo>·</mo><mi>Δ</mi><mi>f</mi></mrow><mi>f</mi></mfrac><mo>+</mo><msub><mi>K</mi><mi>x</mi></msub><mover><mi>x</mi><mo>‾</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>‾</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>‾</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>‾</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>‾</mo></mover><mo>·</mo><mover><mi>y</mi><mo>‾</mo></mover><mo>+</mo><msub><mi>b</mi><mn>1</mn></msub><mover><mi>x</mi><mo>‾</mo></mover><mo>+</mo><msub><mi>b</mi><mn>2</mn></msub><mover><mi>y</mi><mo>‾</mo></mover></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>Δy</mi><mo>′</mo></msup><mo>=</mo><mo>-</mo><msub><mi>y</mi><mn>0</mn></msub><mo>-</mo><mfrac><mrow><mover><mi>y</mi><mo>‾</mo></mover><mo>·</mo><mi>Δ</mi><mi>f</mi></mrow><mi>f</mi></mfrac><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mover><mi>y</mi><mo>‾</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>‾</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>‾</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>‾</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>‾</mo></mover><mo>·</mo><mover><mi>y</mi><mo>‾</mo></mover></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mrow><mover><mi>x</mi><mo>‾</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>‾</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>‾</mo></mover><mn>2</mn></msup><mo>+</mo><msup><mover><mi>y</mi><mo>‾</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>′</mo></msup><mo>=</mo><mi>x</mi><mo>+</mo><msup><mi>Δx</mi><mo>′</mo></msup></mrow></mtd></mtr><mtr><mtd><mrow><msup><mi>y</mi><mo>′</mo></msup><mo>=</mo><mi>y</mi><mo>+</mo><msup><mi>Δy</mi><mo>′</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>′</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>′</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)中两台电子经纬仪获取的相机镜头上的回光反射标志中心的坐标归心到镜头中心后,与解算得到的投影中心位置进行比较,计算得到实际相机投影中心与镜头中心位置的差值,得到投影中心与镜头中心的位置关系,实现对数码相机中心位置的精确标定。 |