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

单像空间后方交会

测绘学院成晓倩
1概述
1.1定义
利用一定数量的地面控制点和对应像点坐标求解单张像片外方位元素的方法称为空间后方交会。

1.2所需控制点个数与分布
共线条件方程的一般形式为：
（1）
式中包含有六个外方位元素，即，只有确定了这六个外方位元素的值，才能利用共线条件方程真正确定一张像片的任一像点与对应地面点的坐标关系。
个数：对任一控制点，我们已知其地面坐标和对应像点坐标，代入共线条件方程可以列出两个方程式，因此，只少需要3个控制点才能解算出六个外方位元素。
在实际应用中，为了避免粗差，应有多余检查点，因此，一般需要4～6个控制点。
分布：为了最有效地控制整张像片，控制点应均匀分布于像片边缘，如下图所示。

分布合理
分布合理
分布不合理

由于共线条件方程是非线性的，直接答解十分困难，所以首先将共线方程改化为线性形式，然后再答解最为简单的线性方程组。

空间后方交会的基本思路
2.1共线条件方程线性化的基本思路
在共线条件方程中，令
（2）
则共线方程变为
（3）
对上式两侧同乘，并移至方程同侧，则有
（4）
令
（5）
由于上式是共线方程的变形，因此，是的函数。
对分别按泰劳级数展开，并且只保留一次项，得
（6）
式中，、分别是和的初值；、分别是和对各个外方位元素的偏导数；分别是初值的增量。
为了明确（6）式中常数项的意义，对（6）式两侧同乘以，则
（7）
考查（7）式中的常数项，有
（8）
式中是像点坐标的观测值；是由相应地面坐标和外方位元素初值计算出的像点坐标。这样（7）式中的常数项就有明确的意义，即为像点观测值和计算值之差。
同样也可以得到，
（9）
现将（7）式改写为
（10）
式中，为残差；为系数；是待求值，是像点观测值和计算值之差。与（7）式相比较，显然有
（10a）
式（10）就是以外方位元素增量为待求值的共线条件方程线性化公式，也称误差方程式。要得到完整的线性化形式，关键是求各个系数，而求的关键是求出对各个外方位元素的偏导数。
如何求偏导数，将在共线方程线性化部分介绍。

2.2答解外方位元素的基本过程
每个控制点都可以按（10）式列出两个误差方程式，n个控制点可列出2n个方程，用矩阵形式可表示为：
（11）
式中
；
；
；
。
如果能答解这2n个方程构成的方程组，则可得到外方位元素的增量。

具体的求解过程应是一个迭代过程：
（1）给出外方位元素的初值，；
（2）对每个控制点计算误差方程式系数和，从而按（10）式组成误差方程式；
（3）答解线性方程组，得到每个外方位元素的增量；
（4）将增量和初值相加，得到新的外方位元素值；
（5）各个增量是否小于规定的限差？若是，则停止迭代运算；若不是，则将新外方位元素值作为初值重复（2）～（5）。

2.3误差方程组的答解方法（最小二乘原理）
式（11）是一个由2n个方程组成的误差方程组，且方程个数多于待求值的个数，对这样的方程组应如何答解呢？在摄影测量中一般按最小二乘原理进行答解。
按最小二乘原理，求出的待求参数的最佳估计值应使各误差方程式的残差平方和为最小，即满足
（12）
这样就转化为对待求值的求极值问题。
下面以式（11）为例，说明求极值后误差方程式的变化。将分别对求极值，即令
（13）
这样将得到六个新的线性方程式，方程式的个数与待求值的个数相同。解这个方程组，则可得到的最佳估计值。在测量平差中把由式（11）变为式（13）的过程称为误差方程式的法化，法化后的方程式称为法方程式。显然，法方程式的系数和常数项将与误差方程式不同。究竟法方程式的系数、常数项和原误差方程式有什么变化，又有什么关系呢？这可以通过较复杂的推导过程来找到。在这里，我们略去推导过程，只按矩阵方式给出结论。
由于

则

整理后有

令
即为法方程式的系数阵。
两边同乘以，则可求出，即
（14）
该式即为的解。

3共线条件方程的线性化
在“共线条件方程线性化的基本思路”中，我们知道：共线条件方程线性化的关健是求各个偏导数（和），下面我们分别求取线元素和角元素的偏导数。

3.1线元素的偏导数
已知
和
则
（15）
如果把内方位元素也作为未知数进行答解，则

3.2角元素的偏导数
和是角元素的复合函数，为了推导的方便，我们将对角元素求导数的过程分为三个步骤。
第一步：求各个方向余弦对的偏导数，共有27个；
第二步：求对的偏导数，共有9个；
第三步：求、对角元素的偏导数，共6个。
【第一步】
已知

则有

【第二步】
在求对的偏导数时，只需将各方向余弦的偏导数代入即可。但，为了公式的简单与工整，应利用各个方向余弦间的关系，将各个导数尽量化为由和表示的形式。这样在乘以后，考虑，各个系数很容易化成以表示的形式。