单像空间后方交会姓名：学号：时间：Echodidthisforyou.2013/4/252目录一、作业任务.................................................................................................................-3-二、计算原理.................................................................................................................-3-三、算法流程.................................................................................................................-7-四、源程序.....................................................................................................................-8-五、计算结果.................................................................................................................-8-六、结果分析.................................................................................................................-8-七、心得与体会.............................................................................................................-8-八、附页.........................................................................................................................-8-1.c++程序.............................................................................................................-8-2.C++程序截图...................................................................................................-15-3.matlb程序......................................................................................................-16-一、作业任务3已知条件：摄影机主距f=153.24mm，x0=0，y0=0,像片比例尺为1:40000，有四对点的像点坐标与相应的地面坐标如下表。以单像空间后方交会方法，求解该像片的外方位元素。二、计算原理1.获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄测坐标。2.测量控制点的像点坐标并作系统误差改正。3.确定未知数的初始值。在竖直摄影且地面控制点大体对称分布的情况下，按如下方法确定初始值，即式中：为摄影比例尺分母；为控制点个数4.用三个角元素的初始值按下式计算各方向余弦值，组成旋转矩阵矩阵中各元素的计算公式如下：4点号像点坐标地面坐标x(mm)y(mm)X(m)Y(m)Z(m)1-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.4664.4340426.5430319.81757.315.逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面坐标，带入以下共线方程式，逐点计算像点坐标的近似值、6.逐点计算误差方程式的系数和常数项，组成误差方程式。5由常数项计算公式：得到常数项矩阵计算式为：7.计算法方程的是系数矩阵和常数项，组成法方程式。8.解法方程，求得外方位元素的改正数。9.用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。6式中：代表迭代次数。10.将求得的外方位元素改正数与规定的限差比较，若小于限差，则迭代结束。否则用新的近似值重复4~9满足要求为止。11.由误差方程的计算式：以及单位权中误差的计算式（式中：表示多余观测数）以及平差中6个参数的协因数阵最终得到参数的中误差为7否是否是输入原始数据像点坐标计算，系统误差改正确定外方位元素初始值组成旋转矩阵R逐点组成误差方程式并法化所有像点完否计算改正后的外方位元素计算中误差，输出成果，结束解法方程，求外方位元素改正数外方位元素改正数是否小于限差是结束并提示错误信息迭代次数小于n三、算法流程8四、源程序源程序代码请见附页五、计算结果在经过三次迭代之后得到的最终成果如表1所示表1，最终计算结果Xs（米）Ys（米）Zs（米）ψ（弧度）ω（弧度）κ（弧度）39795.45225827476.4623857572.685988-0.0039870.0021140.067578六、结果分析由计算结果可知在拍摄照片瞬间，摄影中心在地面摄影测量坐标系中的坐标为（397...
VIP
