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

单像空间后方交会姓名：学号：时间：Echo did this for you .2013/4/25作业任务已知条件：摄影机主距f=153.24mm，x0=0，y0=0, 像片比例尺为1:40000，有四对点的像点坐标与相应的地面坐标如下表。点号像点坐标地面坐标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.31以单像空间后方交会方法，求解该像片的外方位元素。计算原理获取已知数据。从航摄资料中查取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄测坐标。测量控制点的像点坐标并作系统误差改正。确定未知数的初始值。在竖直摄影且地面控制点大体对称分布的情况下，按如下方法确定初始值，即式中：为摄影比例尺分母；为控制点个数用三个角元素的初始值按下式计算各方向余弦值，组成旋转矩阵矩阵中各元素的计算公式如下：逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地面坐标，带入以下共线方程式，逐点计算像点坐标的近似值、逐点计算误差方程式的系数和常数项，组成误差方程式。由常数项计算公式：得到常数项矩阵计算式为：计算法方程的是系数矩阵和常数项，组成法方程式。解法方程，求得外方位元素的改正数。用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。式中：代表迭代次数。将求得的外方位元素改正数与规定的限差比较，若小于限差，则迭代结束。否则用新的近似值重复4~9满足要求为止。由误差方程的计算式：以及单位权中误差的计算式（式中：表示多余观测数）以及平差中6个参数的协因数阵最终得到参数的中误差为算法流程源程序源程序代码请见附页计算结果在经过三次迭代之后得到的最终成果如表1所示表1，最终计算结果Xs（米）Ys（米）Zs（米）ψ（弧度）ω（弧度）κ（弧度）39795.45225827476.4623857572.685988-0.0039870.0021140.067578结果分析由计算结果可知在拍摄照片瞬间，摄影中心在地面摄影测量坐标系中的坐标为（39795.452258，27476.462385，7572.685988）（单位：米），航向倾角ψ为-0.003987弧度，旁向倾角ω为0.002114弧度，像片旋角κ为0.067578弧度。表2，精度评定结果　参数XsYsZsψωκ中误差1.1073 1.2494 0.4881 0.0002 0.0002 0.0001 心得与体会通过本次实验，我对单张像片的空间后方交会的计算原理及实现过程有了很深的认识，并在牢记各种计算公式的推导过程基础上做到了熟练应用。在本次实验中，我遇到了很多困难，但是都一一得到了解决，在解决困难的过程中的编程能力得到了提高，对其所涉及到的知识的印象也得到了加深。这个程序本身还有一些不足，例如:在迭代过程中的矩阵的相乘、输出等，没有采用编写函数的形式，而是逐一进行计算，这是程序整体看起来较臃肿。在编程过程中为了确保计算的准确性，我还利用matlab同步编写了一个相同功能的计算单张相片空间后方交会的迭代程序，并以此来验算原c++程序的正确性。现在编写完这个程序，心中很自豪，也有一些遗憾，因为我想用c#来编写，因为c#在暑期实习用过之后，时隔近一年在没有复习，不免把学会的知识又还给了老师，本想通过这次实验对c#进行复习，但是由于时间原因以及一些个人因素，没能达到自己原本的目标。总之，感谢老师给了我这次锻炼自己、提升自己的机会。附页c++程序#includeiostreamusing namespace std;#includefstream#include iomanipconst int n=6;void inverse(double c[n][n]);templatetypename T1,typename T2void transpose (T1*mat1,T2*mat2,int a,int b);templatetypename T1,typename T2void multi(T1*mat1,T2*mat2,T2 * result,int a,int b,int c);int main(){double x[4][2]={-0.08615,-0.06899,-0.05340,0.08221,-0.01478,-0.07663,0.01046,0.06443};double X[4][3]={36589.41,25273.32,2195.17,37631.08,31324.51,728.69,39100