单像空间后方交会.doc

单像空间后方交会测绘学院 成晓倩概述定义利用一定数量的地面控制点和对应像点坐标求解单张像片外方位元素的方法称为空间后方交会。所需控制点个数与分布共线条件方程的一般形式为:x x0fa1(Xa3(Xa2(XXS)XS)XS)b1(Yb3(Yb2(YYS)YS)YS)c1(Zc3(Zc2(ZZS)ZS)ZS)(1)y y0fa3(XXS)b3(YYS)c3(ZZS)式中包含有六个外方位元素, 即XS、YS、ZS、、、,只有确定了这六个外方位元素的值,才能利用共线条件方程真正确定一张像片的任一像点与对应地面点的坐标关系。个数:对任一控制点,我们已知其地面坐标(Xi、Yi、Zi)和对应像点坐标(xi、yi),代入共线条件方程可以列出两个方程式, 因此,只少需要3个控制点才能解算出六个外方位元素。在实际应用中,为了避免粗差,应有多余检查点,因此,一般需要 4~6个控制点。分布:为了最有效地控制整张像片,控制点应均匀分布于像片边缘,如下图所示。分布合理 分布合理 分布不合理由于共线条件方程是非线性的, 直接答解十分困难,所以首先将共线方程改化为线性形式,然后再答解最为简单的线性方程组。空间后方交会的基本思路共线条件方程线性化的基本思路Xa1(XXS)b1(YYS)c1(ZZS)Ya2(XXS)b2(YYS)c2(ZZS)Za3(XXS)b3(YYS)c3(ZZS)则共线方程变为在共线条件方程中,令(2)Xx0 fZYy0 fZ(3)对上式两侧同乘 Z,并移至方程同侧,则有令fX (xfY (yx0)Z 0y0)Z 0(4)Fx fXFy fY(x x0)Z(y y0)Z(5)由于上式是共线方程的变形,因此,Fx、Fy是XS、YS、ZS、、、的函数。对Fx、Fy分别按泰劳级数展开,并且只保留一次项,得0Fx (Fx)Fy (Fy)0Fx FxXS YSXS YSSSFy X Fy YXS YSFx Fx FxZSZSSFy Z Fy FyZSFxFy (6)式中,(Fx)0、(Fy)0分别是Fx和Fy的初值; Fx、Fy分别是Fx和Fy对各个外方位元素的偏导数;XS、YS、ZS、 、 、分别是XS、YS、ZS、、、 初值的增量。1为了明确(6)式中常数项的意义,对( 6)式两侧同乘以 ,则Z1FxZ1(Fx)0Z( 1)FxZ XS1)Fx YZ YS( 1)Fx ZZ ZSSS( 1)FxZ( 1)Z( 1)FxZSS(7)1FyZ1(Fy)0Z( 1)FyZ XS1)Fy YZ YS( 1)Fy ZZ ZS1 Fy 1( ) ( )Z Z1 Fy( )ZXFxS(XFyS(1(Fx)01[fX(xx0)Z][(xx0)(fX)]ZZZ考查(7)式中的常数项,有[(xx0)(x计x0)]=-(xx计)(8)式中x是像点坐标的观测值;x计是由相应地面坐标和外方位元素初值计算出的像点坐标。1(Fy)01[fY(yy0)Z][(yy0) (fY)]ZZZ这样(7)式中的常数项就有明确的意义,即为像点观测值和计算值之差。同样也可以得到,现将(7)式改写为[(yy0)(y计y0)]=-(yy计)(9)vxa11XSa12YSa13ZSa14a15a16lxvya21XSa22YSa23ZSa24a25a26ly(10)式中,vx、vy为残差;aij为系数;XS、YS、ZS、 、 、是待求值,lx、ly是像点观测值和计算值之差。与( 7)式相比较,显然有a11a121 FxZ XS1 FxZ YSa21a221 FyZXS1 FyZYSa13=1 FxZ ZSa23=1 FyZZSa= 1 Fx a= 1Fy14 24Z Z1 Fxa=1Fya=(10a)15 25Z Za= 1 Fx a= 1Fy16Zlx x x计26Zly y y计式(10)就是以外方位元素增量为待求值的共线条件方程线性化公式, 也称误差方程式。要得到完整的线性化形式,关键是求各个系数aij,而求aij的关键是求出Fx、Fy对各个外方位元素的偏导数。如何求偏导数,将在共线方程线性化部分介绍。答解外方位元素的基本过程每个控制点都可以按( 10)式列出两个误差方程式, n个控制点可列出 2n个方程,用矩阵形式可表示为:式中vvvV1 1 2x y xV AX Lvvvn nTy x y ;(11)aaa1 1 111 12 16aaa1 1 121 22 26aaa2 2 211 12 16aaaA 2 2 2 ;21 22 26aaan n n11 12 16aaan n n21 22 26TX XS YS ZS ;llllL1