单像空间后方交会实验报告(c 版).docxVIP

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