单向后方交会实验

《单向后方交会实验》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、课程设计报告班级：测绘一班学号:20133279姓名：罗超日期：2015.4.26目录一、计算原理3二、算法流程4三、源程序5四、计算结果13五、结果分析13六、心得体会13一、计算原理已知条件摄影机主距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.63

2、39100.9724934.982386.50410.4664.4340426.5430319.81757.31以单像空间后方交会方法，求解该像片的外方位元素。二、算法流程（1）获取已知数据。从航摄资料中差取平均航高与摄影机主距；获取控制点的地面测量坐标并转换为地面摄影坐标。（2）量测控制点的像点坐标并作系统误差改正。（3）确定未知数的初始值。在竖直摄影且地面控制点大体对称分布的情况下，按如下方法确定初始值，即，，（4）用三个角元素的初始值按下式，计算各个方向余弦值，组成旋转矩阵R（5）逐点计算像点坐标的近似值。利用未知数的近似值和控制点的地

3、面坐标；带入共线方程式，逐点近似像点坐标的近似值（x）、（y）。（6）逐点计算误差方程式的系数和常数项，组成误差方程式。（7）计算法方程的系数矩阵和常数项，组成法方程式。（8）解法方程，求得外方位元素的改正数。（9）用前次迭代取得的近似值，加本次迭代的改正数，计算外方位元素的新值。（10）将求得的外方位元素改正数与规定的限差比较，若小于限差，则迭代结束。负责用新的近似值重复（4）-（9），直到满足要求为止。用共线方程进行空间后方交会的程序框如图所示。输入原始数据像点坐标计算，系统误差正确定外方位因素初始值组成旋转矩阵R逐点组成误差方程式并法化

4、所有像点完否是解法方程，求外方位元素改正数计算改正后的外方位元素外方位元素改正数是否小于限差是输出计算成果，计算并结束结束并显示错误信息迭代次数小于n是否否三、源程序#include”iostream"usingnamespacestd;#include”fstream"#includeconstintn=6;voidinverse(doublec[n][n]);templatevoidtranspose(TI*mat1,T2*mat2,inta,intb);template

5、ypenameTI,typenameT2>voidmulti(T1*mat1,T2*mat2,T2*result,inta,intb,intc);intmain(){doublex[4][2]={-0.08615,-0.06899,-0.05340,0.08221，-0.01478，-0.07663,0.01046,0.06443};doubleX[4][3]={36589.41,25273.32,2195.17,37631.08,31324.51,728.69,39100.97,24934.98,2386.50,40426.54,30319

6、.81,757.31};inti,j,num=0;//num为迭代次数doubleX0[6]={0};//设定未知数(XS,YS,ZS,三个角度)初始值doublef=0.15324;//摄影机主距f=153.24mmdoublea=1/40000.0;//像片比例尺为1:40000doubleR[3][3]={0};//初始化旋转知阵Rdoubleapprox_x[8]={0};//用于存放像点估计值doubleA[8][6]二{0};//定义了一个系数阵doubleAT[6][8]二{0};//定义了A的转置知阵doubleL[8]={0

7、}://定义常数项constdoublepi=3.1415926535897932;doubleAsum[6][6]={0};doubleresult2[6]={0};doubleresultl[6][8]={0};doublesumXYZ[3]={0};cout.precision(5);cout<<”已知像点坐标为:”;for(i=0;i<4;i++)for(j=0;j<2;j++){cout<

8、y"<