单像空间后方交会习题.doc

单像空间后方交会<遥感07-1、2学生适用>测绘学院　王双亭1 概述1.1　定义　　　利用一定数量的地面控制点和对应像点坐标求解单张像片外方位元素的方法称为空间后方交会.1.2　所需控制点个数与分布共线条件方程的一般形式为：　　〔1〕式中包含有六个外方位元素,即,只有确定了这六个外方位元素的值,才能利用共线条件方程真正确定一张像片的任一像点与对应地面点的坐标关系.个数：对任一控制点,我们已知其地面坐标和对应像点坐标,代入共线条件方程可以列出两个方程式,因此,只少需要3个控制点才能解算出六个外方位元素.在实际应用中,为了避免粗差,应有多余检查点,因此,一般需要4～6个控制点.分布：为了最有效地控制整张像片,控制点应均匀分布于像片边缘,如下图所示.分布合理分布合理分布不合理　　由于共线条件方程是非线性的,直接答解十分困难,所以首先将共线方程改化为线性形式,然后再答解最为简单的线性方程组.2 空间后方交会的基本思路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个.[第一步]已知则有[第二步]在求对的偏导数时,只需将各方向余弦的偏导数代入即可.但,为了公式的简单与工整,应利用各个方向余弦间的关系,将各个导数尽量化为由和表示的形式.这样在乘以后,考虑,各个系数很容易化成以表示的形式.下面仅以求为例,说明对偏导数的推导过程.由可得：由于、、,代入上式有这个式子显然不是用表示的形式,如何将其化为以表示的形式呢？前面我们讲旋转矩阵的性质时,已知道旋转矩阵的每个元素等于其代数余子式,因此有,,这里是元素的行、列数.代入上式,有　　对其余的求导也作类似处理,可等到全部9个偏导数如下：[第三步]直接对、求的偏导数,并将上式分别代入即可得到6个角元素偏导数,即　　〔16〕3.3误差方程的各个系数将〔15〕和〔16〕式代入〔10a〕中,即可得到误差方程式的各个系数,即 <17>3.4误差方程的近似形式在近似垂直摄影情况下,各个角元素都很小,此时,,且,代入<17>式,则误差方程的各个系数变为： <18>这样,可得到误差方程的近似形式为<19>4　空间后方交会的计算过程<1>获取已知数据.包括：个控制点的地面坐标；内方位元素；摄影航高；像片比例尺.<2>量测个控制点对应的像点坐标,并进行必要的系统误差改正.<3>确定外方位元素的初值.在近似垂直摄影情况下,各个初值可按如下方法确定： <4>计算各个方向余弦,组成旋转矩阵R.<5>逐点计算像点坐标值,常数项和误差方程式系数,即逐点组建误差方程式.<6>计算法方程式的系数阵和常数阵,组成法方程式.<7>按<14>式解各个外方位元素的增量<或改正数>,并与相应初值求和,得到外方位元素的新初值.<8>检查计算是否收敛.将所求的外方位元素改正数与规定的限差比较,通常只对角元素的改正数设定限差<一般为>.当三个角改正数都小于限差时,迭代结束；否则,用新的初值重复<4>～<8>步.<9>外方位元素的精度估算.5　空间后方交会的理论精度根据平差原理,平差后的单位权中误差为 <20>式中是单位权中误差；是第个方程的残差；是控制点个数；是未知数个数,在此为6. 根据误差传播定律, <21>式中是法方程系数阵的逆阵.若则第个未知数的中误差为 <22>作业：1、 何为单像空间后方交会？2、 试分析空间后方交会所需控制点个数和分布？3、 简述共线条件方程线性化的基本思路.4、 证明,进而证明5、 请给出组建误差方程的程序框图.6、 试说明后方交会的计算过程.10 / 10。