摄影测量空间后方交会外方位元素的解算程序

航摄相片的外方位元素表示的是摄影摄影瞬间相片上的点对于地面上的点之间的关系的一些参数，在测绘工作中，如果求出了一张航拍相片的外方位元素，那么就可以根据像素点的坐标计算出对应的地面点的坐标，而这些解算过程计算量是非常大的。所以，编写一个能够计算航空摄影测量外方位元素的程序就很有必要了。今天，小博根据所学专业知识完成了该程序的编写，下面是该程序的详细算法，希望对大家有所帮助。

为了方便大家在测绘技术上的交流，小博创建了一个QQ群，欢迎大家加入测绘技术交流QQ群234035436

程序图片

程序计算原理

单张影像的空间后方交会：利用已知地面控制点数据及相应像点坐标根据共线方程反求影像的外方位元素。  

数学模型为共线条件方程式：

求解过程： 

（1）获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影测量坐标。 

（2）量测控制点的像点坐标并做系统改正。 

（3）确定未知数的初始值。在竖直摄影且地面控制点大致分布均匀的情况下，按如下方法确定初始值，即：

             q=w=k=0 

             式中:m为摄影比例尺分母；n为控制点个数。  

（4）用三个角元素的初始值，计算个方向余弦，组成旋转矩阵R。 

（5）逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面

坐标代入共线方程式，逐点计算像点坐标的近似值（x）、（y）。 

（6）逐点计算误差方程式的系数和常数项，组成误差方程式。 

（7）计算法方程的系数矩阵A和常数项L，组成法方程式。 

（8）解法方程，求得外方位元素的改正数dXs，dYs，dZs，dq，dw，dk。 

（9）用前次迭代取得的近似值，