(完整word版)单像空间后方交会.doc

1、单像空间后方交会测绘学院成晓倩1 概述1.1定义利用一定数量的地面控制点和对应像点坐标求解单张像片外方位元素的方法称为空间后方交会。1.2所需控制点个数与分布共线条件方程的一般形式为：x x0f a1 ( X X S ) b1 (Y YS ) c1 ( Z Z S )a3 ( X X S ) b3 (Y YS ) c3 ( Z Z S )（ 1）y y0f a2 ( X X S ) b2 (Y YS ) c2 (Z Z S )a3 ( X X S ) b3 (Y YS ) c3 (Z Z S )式中包含有六个外方位元素，即 X S、 YS、 Z S、 、 、 ，只有确定了这六个外方位元素的值，

2、才能利用共线条件方程真正确定一张像片的任一像点与对应地面点的坐标关系。个数： 对任一控制点，我们已知其地面坐标( X i、 Yi 、 Zi ) 和对应像点坐标 ( xi、 yi ) ，代入共线条件方程可以列出两个方程式，因此，只少需要 3 个控制点才能解算出六个外方位元素。在实际应用中，为了避免粗差，应有多余检查点，因此，一般需要4 6 个控制点。分布： 为了最有效地控制整张像片，控制点应均匀分布于像片边缘，如下图所示。分布合理分布合理分布不合理由于共线条件方程是非线性的， 直接答解十分困难， 所以首先将共线方程改化为线性形式，然后再答解最为简单的线性方程组。2 空间后方交会的基本思路2.1共

3、线条件方程线性化的基本思路在共线条件方程中，令X a1 ( X X S ) b1 (Y YS ) c1 (Z Z S )Ya2 ( X X S ) b2 (YYS )c2 ( ZZ S )（ 2）Za3 ( X X S ) b3 (YYS )c3 ( ZZ S )则共线方程变为xx0XfZ（ 3）yy0YfZ对上式两侧同乘Z ，并移至方程同侧，则有f X( xx0 )Z 0（ 4）f Y( yy0 )Z 0令Fxf X( x x0 )Z（ 5）Fyf Y( y y0 )Z由于上式是共线方程的变形，因此，Fx、 Fy 是 X S、 YS、 Z S、 、 、的函数。对 Fx、 Fy 分别按泰劳级数

4、展开，并且只保留一次项，得Fx( Fx) 0FxX SFxYSFxZ SFxFxFxX SYSZ S（ 6）Fy( Fy) 0FyX SFyYSFyZ SFyFyFyX SYSZ S式中， (Fx )0、 ( Fy ) 0分别是 Fx 和 Fy 的初值；Fx 、 Fy 分别是 Fx 和 Fy 对各个外方位元素的偏导数；X S、 YS、 ZS、分别是 X S、YS、ZS、 、 、初值的增量。为了明确（ 6）式中常数项的意义，对（6）式两侧同乘以1，则Z1 Fx1 (Fx)0(1 )FxXZZZX S(1 )Fx(1 )FxZZ1 Fy1 (Fy) 0(1 )FyXZZZX S(1 )Fy(1 )

5、FyZZSS( 1 ) Fx YS( 1 ) Fx Z SZYSZZ S(1 )FxZ（ 7）( 1 ) Fy YS( 1 ) Fy Z SZYSZZ S(1 )FyZ考查（ 7）式中的常数项，有1( Fx) 01 f X (xx0 ) Z( xx0 ) ( f X )（ 8）ZZZ( x x0 ) ( x计x0 )( xx计 )式中 x 是像点坐标的观测值；x计 是由相应地面坐标和外方位元素初值计算出的像点坐标。这样（ 7）式中的常数项就有明确的意义，即为像点观测值和计算值之差。同样也可以得到，1 (Fy )01 f Y ( yy0 ) Z ( yy0 ) (f Y )（ 9）ZZZ( y

6、y0 ) ( y计y0 )( yy计 )现将（ 7）式改写为vxa11 X Sa12 YSa13 ZSa14a15a16l x（ 10）v ya21 X Sa22 YSa23 ZSa24a25a26l y式中， vx、v y 为残差； aij为系数；X S、 YS、 ZS、是待求值， l x、 l y是像点观测值和计算值之差。与（7）式相比较，显然有a111Fxa2 11FyZX SZX Sa121Fxa221FyZYSZYSa131Fxa231FyZZ SZZ Sa141Fxa241FyZZa151Fxa251Fy（10a）ZZFxFya161a261ZZl xxx计l yyy计式（ 10）

7、就是以外方位元素增量为待求值的共线条件方程线性化公式，也称误差方程式。要得到完整的线性化形式，关键是求各个系数aij，而求 aij 的关键是求出Fx、 Fy 对各个外方位元素的偏导数。如何求偏导数，将在共线方程线性化部分介绍。2.2答解外方位元素的基本过程每个控制点都可以按（10）式列出两个误差方程式，n 个控制点可列出2n个方程，用矩阵形式可表示为：VAX L（ 11）式中V vx1v1yvx2v y2vxnv yn T ；a111a121a161a211a221a261Aa112a122a162222；a21a22a26a11na12na16na21na22na26nXX SYSZ ST；

8、L l x1l y1l x2l y2l xnl yn T 。如果能答解这2n 个方程构成的方程组，则可得到外方位元素的增量。具体的求解过程应是一个迭代过程：（ 1）给出外方位元素的初值，X S0、YS0、ZS0、 0、 0、 0；（ 2）对每个控制点计算误差方程式系数aij 和 l x、 l y ，从而按（10）式组成误差方程式；（3）答解线性方程组，得到每个外方位元素的增量XS、YS、ZS、；（ 4）将增量和初值相加，得到新的外方位元素值；（ 5）各个增量是否小于规定的限差？若是，则停止迭代运算；若不是，则将新外方位元素值作为初值重复（ 2）（ 5）。2.3误差方程组的答解方法（最小二乘原理

9、）式（ 11）是一个由 2n 个方程组成的误差方程组，且方程个数多于待求值的个数，对这样的方程组应如何答解呢？在摄影测量中一般按最小二乘原理进行答解。按最小二乘原理 ，求出的待求参数的最佳估计值应使各误差方程式的残差平方和为最小，即满足V T V min（ 12）这样就转化为V T V 对待求值的求极值问题。下面以式（ 11）为例，说明求极值后误差方程式的变化。将对X S、YS、Z S、求极值，即令VTV分别V T V0X SV T V0YSV T V0Z SV T V0V T V（ 13）0V T V0这样将得到六个新的线性方程式，方程式的个数与待求值的个数相同。解这个方程组， 则可得到 X

10、 S、 YS、 ZS、的最佳估计值。在测量平差中把由式（11）变为式（13）的过程称为误差方程式的法化，法化后的方程式称为法方程式。显然，法方程式的系数和常数项将与误差方程式不同。究竟法方程式的系数、 常数项和原误差方程式有什么变化，又有什么关系呢？这可以通过较复杂的推导过程来找到。在这里， 我们略去推导过程， 只按矩阵方式给出结论。由于V TV (AX L)T(AX L)XTATAX 2XTATL LTL则V T V2AT AX 2ATL 0X整理后有ATAXATL令N A T A 即为法方程式的系数阵。两边同乘以 N 1 ，则可求出 X ，即X N 1ATL（ 14）该式即为X S、YS、

11、ZS、的解。3 共线条件方程的线性化在“共线条件方程线性化的基本思路”中，我们知道：共线条件方程线性化的关健是求各个偏导数（ Fx 和 Fy ），下面我们分别求取线元素和角元素的偏导数。3.1线元素的偏导数已知Fxf X(xx0 ) ZX a1 ( X X S ) b1 (Y YS ) c1 ( Z Z S )和Y a2 ( X X S ) b2 (Y YS ) c2 ( Z Z S )Fyf Y( yy0 )ZZ a3 ( X X S ) b3 (Y YS ) c3 ( Z Z S )则FxfX(x x0 )Za1fa3 ( x x0 )X SX SX SFxfX( x x0 )Zb1 fb

12、3 ( x x0 )YSYSYSFxfX(x x0 )Zc1 f c3 (x x0 )Z SZ SZS（ 15）FyYZf( y y0 )a2 fa3 ( y y0 )X SX SX SFyfY( y y0 ) Zb2 f b3 ( y y0 )YSYSYSFyfY( y y0 )Zc2fc3 ( y y0 )Z SZ SZ S如果把内方位元素也作为未知数进行答解，则FxZFx0FxXx0y0fFy0FyZFyYx0y0f3.2角元素的偏导数Fx 和 Fy 是角元素、的复合函数， 为了推导的方便，我们将对角元素求导数的过程分为三个步骤。第一步：求各个方向余弦ai、 bi 、 ci 对、的偏导数

13、，共有27 个；第二步：求X、Y、Z 对、的偏导数，共有9 个；第三步：求 Fx 、 Fy 对角元素、的偏导数，共6 个。【第一步】已知a1coscos- sinsinsina2-cos sin- sinsincosa3-sincosb1cossinb2coscosb3-sinc1sincoscossinsinc2-sinsincossincosc3coscos则有a1c1a1a3 sin,a1a2,a2a2a3 cos,a2a1,c2 ,a3a3sin sin,a30,c3 ,b10,b1b3 sin,b1b2,b20,b2b3 cos,b2b1,b30,b3cos,b30,c1a1,c1c

14、3 sin,c1c2,c2a2c2c3 cos,c2c1,c3a3 ,c3cossin,c30,【第二步】在求 X、Y、Z 对、的偏导数时，只需将各方向余弦的偏导数代入即可。但，为了公式的简单与工整，应利用各个方向余弦间的关系，将各个导数尽量化为由X、Y、Z 和xx0X1fai、 bi 、 ci 表示的形式。这样在乘以后，考虑Z ，各个系数很容易化成Zyy0YfZ以 x、y 表示的形式。下面仅以求X 为例，说明X、Y、Z 对、偏导数的推导过程。由()()() 可得：X X Sb1 Y YSc1 Z Z SX a1Xa1 ( X X S )b1 (Y YS )c1 (Z Z S )由于a1c1

15、、b10 、c1a1 ，代入上式有Xc1 ( X X S ) a1 ( Z Z S )这个式子显然不是用X、Y、Z 表示的形式， 如何将其化为以 X、Y、Z 表示的形式呢？前面我们讲旋转矩阵的性质时，已知道旋转矩阵的每个元素等于其代数余子式，因此有c1( 1)i j a2a3， a1 ( 1)i j b2b3b2b3c2，这里 i、 j 是元素的行、列数。代入上式，有c3X(a2 b3a3b2 )( X X S ) (b2 c3 b3 c2 )( Z Z S )a2b3 ( X X S ) a3 b2 ( X X S ) b2c3 ( Z Z S ) b3c2 ( Z Z S )b3 a2 (

16、 X X S ) c2 ( Z ZS ) b2 a3 ( X X S ) c3 (Z Z S )b3 a2 ( X X S ) c2 ( Z ZS ) b2 a3 ( X X S ) c3 (Z Z S ) b3 b2 (Y YS ) b3b2 (Y YS )b3 a2 ( X X S ) b2 (Y YS ) c2 (Z Z S ) b2 a3 ( X X S ) b3 (Y YS ) c3 (Z Z S )b3 Yb2 Z对其余的求导也作类似处理，可等到全部9 个偏导数如下：XXZ sinXb3 Yb2 ZYYXYZ cosYb3b1 ZXZXZX sinY cosZb2b1 YX【第三步

17、】直接对 Fx 、Fy 求、的偏导数，并将上式分别代入即可得到6 个角元素偏导数，即FxfX(x x0 ) Zf ( b3 Y b2 Z ) (x x0 )( b2 X b1Y )FxfX(xx0 ) ZZ f sin( xx0 )( X sinY cos)FxfX(x x0 ) Zf YFyY( y y0 ) Z（ 16）ff ( b3 Y b1 Z ) ( y y0 )( b2 X b1Y)FyfY( yy0 ) ZZ f cos( yy0 )( X sinY cos)FyfY( y y0 ) Zf X3.3误差方程的各个系数将（ 15）和（ 16）式代入（ 10a）中，即可得到误差方程式

18、的各个系数，即a111a2 11 a1 f a3 ( x x0 ) a 2f a (3y y )0ZZa121a2 21 b1 f b3 (x x0 )b 2f b (3y y )0ZZa131c3 (x x0 )a2 31 c1 fc 2f c (3y y )0ZZ( y y0 )2a14( x x0 )( y y0 )a2 4 ffb1fb1 f(xx0 ) 2 b2y y0 b3(xx0 )( y y0 ) b2 x x0 b3ffa15( x x0 )2a25(x x0 )( y y0 )fsinfsin(xx0 )( yy0 ) cosf sin( yy0 )2cosf cosffa

19、16y y0a26( x x0 )axx0 , a1,a 0ayy0 ,a0, a29117f181927f28(17)3.4误差方程的近似形式在近似垂直摄影情况下，各个角元素都很小，此时，0，且 ZH ，代入(17) 式，则误差方程的各个系数变为：a 11fa210Hfa 120a22H( xx 0 )a 13a23( y y0 )H2H( xx 0)(18)( x x0 )( y y0 )a 14 fa24ffa 15( xx 0 )( yy 0 )a25( y y0 ) 2f ffa 16yy 0a26( x x0 )这样，可得到误差方程的近似形式为vxfx x0ZS( x x0 ) 2

20、(x x0 )( y y0 )( y y0 )l xX SH ffHfv yfy y0Z S( x x0 )( y y0 ) f( y y0 )2( x x0 )l yYSHfHf(19)4 空间后方交会的计算过程(1) 获 取 已 知 数 据 。 包 括 ： n 个 控 制 点 的 地 面 坐 标 ( X i、 Yi 、 Z i ) ； 内 方 位 元 素x0、 y0、 f ；摄影航高H ；像片比例尺m 。(2)量测 n 个控制点对应的像点坐标( xi 、 yi ) ，并进行必要的系统误差改正。(3)确定外方位元素的初值X S 0 、 Y S 0、 Z S 0 、0、0 、0 。在近似垂直摄影情况下，各个初值可按如下方法确定：01X SnnX iYS01n