单张影像空间后方交会程序.doc

摄影测量实习报告 -----单张影像空间后方交会程序 实习时间 2013.5.20-2013.5.24 学生班级 10测绘（1）班 学生姓名 房新明 学生学号 1072143138 所在院系 矿业工程学院 指导老师 张会战 邵亚琴 一．实习目的 1．深入理解单张影像空间后方交会的原理，体会在有多余观测情况下，用最小二乘平差方法编程实现解求影像外方位元素的过程。 2．利用Visual C++编写一个完整的单张影像空间后方交会程序，通过对提供的试验数据进行计算，输出像片的外方位元素并进行评定精度。 3．通过编写程序实现单张影像空间后方交会计算，掌握非线性方程线性化的过程、相应数据读入与存储的方法以及迭代计算的特点，巩固各类基础课程及计算机课程的学习内容，培养上机调试程序的动手能力，通过对实验结果的分析，增强综合运用所学知识解决专业实际问题的能力。 4、掌握空间后方交会的定义和实现算法 (1) 定义：空间后方交会是以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素Xs，Ys，Zs，φ，ω，κ。 (2) 算法：由于每一对像方和物方共轭点可列出2个方程，因此若有3个已知地面坐标的控制点，则可列出6个方程，解求6个外方位元素的改正数△Xs，△Ys，△Zs，△φ，△ω，△κ。实际应用中为了提高解算精度，常有多余观测方程，通常是在影像的四个角上选取4个或均匀地选择更多的地面控制点，因而要用最小二乘平差方法进行计算。 5、了解空间后方交会的基本过程 (1) 获取已知数据。从摄影资料中查取影像比例尺1/m，平均摄影距离（航空摄影的航高）、内方位元素x0，y0，f；获取控制点的空间坐标Xt，Yt，Zt。 (2) 量测控制点的像点坐标并进行必要的影像坐标系统误差改正，得到像点坐标。 (3) 确定未知数的初始值。单像空间后方交会必须给出待定参数的初始值，在竖直航空摄影且地面控制点大体对称分布的情况下，Xs0和Ys0为均值，Zs0为航高，φ、ω、κ的初值都设为0。或者κ的初值可在航迹图上找出或根据控制点坐标通过坐标正反变换求出。 (4) 计算旋转矩阵R。利用角元素近似值计算方向余弦值，组成R阵。 (5) 逐点计算像点坐标的近似值。利用未知数的近似值按共线条件式计算控制点像点坐标的近似值（x），(y)。 (6) 逐点计算误差方程式的系数和常数项，组成误差方程式。 (7) 计算法方程的系数矩阵ATA与常数项ATL，组成法方程式。 (8) 解求外方位元素。根据法方程，解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。 (9) 检查计算是否收敛。将所求得的外方位元素的改正数与规定的限差比较，通常对φ，ω，κ的改正数△φ，△ω，△κ给予限差，通常为0.1′，当3个改正数均小于0.1′时，迭代结束。否则用新的近似值重复(4)～（8）步骤的计算，直到满足要求为止。 二、实习要求 1、认真复习单张影像空间后方交会的原理及有关内容； 2、编写程序，计算左片号23和右片号24的外方位元素并进行精度评定； 3、验证数据部分外方位元素和旋转矩阵R已经给出，请验证自己编写的程序是否正确； 三、实习环境 1、硬件环境：windows操作系统 2、软件环境：VC++ 四、实习原理 1、以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，求解该影象在航空摄影时刻的像片外方位元素 ， ， ， ，ω，κ共线条件方程如下：x-x0=-f*[a1(X-Xs)+b1(Y-Ys)+c1(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)]y-y0=-f*[a2(X-Xs)+b2(Y-Ys)+c2(Z-Zs)]/[a3(X-Xs)+b3(Y-Ys)+c3(Z-Zs)]其中：x,y为像点的像平面坐标； Xs, ,f为影像的外方位元素；，为摄站点的物方空间坐标；X,Y,Z为物方点的物方空间坐标；旋转矩阵R ; 由于此共线条件方程是非线性方程，先对其进行线性化，像点观测值一般视为等权，即P=I；矩阵形式：V=AX-L，P=I；通过间接平差，为提高精度，增加多余观测方程，根据最小二乘平差原理，可计算出外方位元素的改正数。经过迭代计算，每次迭代用未知数的近似值与上次迭代计算的改正数之和作为新的近似值，重复计算，求出新的改正数，这样反复趋近，直到改正数小于某个限值为止。Mi=m0*√Qii 其中m0=±√【VV】/（2N-