单像空间后方交会实习报告.doc

遥感信息工程学院 201130259XXXX XX 单像空间后方交会实习报告 （遥感信息工程学院2011级） 班 级： XXXXX 姓 名： XX 学 号： 201130259XXXX 指导教师： 袁修孝 1、 实习目的 1、掌握空间后方交会的定义和实现算法 2、了解摄影测量平差的基本过程 3、熟练VC程序编写2、 实习原理 空间后方交会以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素Xs，Ys，Zs，。由于每一对像方和物方共轭点可列出2个方程，因此若有3个已知地面坐标的控制点，则可列出6个方程，解求6个外方位元素的改正数Xs，Ys，Zs，。在本次编程实习中，为了提高解算精度，引入多余观测方程，特在影像上选取了4个地面控制点，用最小二乘平差方法进行计算。3、 实习方法 在老师的指导下，运用VC6.0上机编写程序。四、实习过程 1、弄清楚平差计算基本过程（即编程基本思路）： 获取已知数据。包括影像比例尺1/m，平均摄影距离（航空摄影的航高）H，内方位元素x0、y0、f，控制点的空间坐标X、Y、Z。 量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。 确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，Xs0和Ys0为均值，Zs0为航高，、的初值都设为0。或者的初值可在航迹图上找出或根据控制点坐标通过坐标正反变换求出。 计算旋转矩阵R。利用角元素近似值计算方向余弦值，组成R阵。 逐点计算像点坐标的近似值。利用未知数的近似值按共线条件式计算控制点像点坐标的近似值（x），(y)。输入原始数据归算像点坐标x，y计算和确定初值Xs0，Ys0，Zs0,0，0，0组成旋转矩阵R计算(x),(y)和lx，ly逐点组成误差方程式并法化所有点完否？计算改正后的外方位元素解法方程，求未知数改正数未知数改正数限差否？整理并输出计算结果正常结束否迭代次数小于限差否？否是否输出中间结果和出错信息非正常结束 图1 逐点计算误差方程式的系数和常数项，组成误差方程式。 计算法方程的系数矩阵ATA与常数项ATL，组成法方程式。 解求外方位元素。根据法方程，解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。 检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较，通常对，的改正数，给予限差，通常为0.000001弧度，当3个改正数均小于0.000001弧度时，迭代结束。否则用新的近似值重复步骤的计算，直到满足要求为止。 计算单位权中误差以及各个未知数的中误差。 2、在纸上绘出计算机程序框图如图1所示； 3、上机编写、调试程序，输出结果如图2所示。 图2 4、 实习小结通过这次编程实习，我深入理解了单像空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。在理论知识的基础上，上机编程积累了实践经验，增强了综合运用所学知识解决实际问题的能力。在跟老师讲解编程思路的时候，老师提出了一些概念性