作业4--空间后方交会

《作业4--空间后方交会》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、实用文案作业报告空间后方交会专业：测绘工程班级：2008级(1)班姓名：陈闻亚指导教师：陈强2010年4月16日标准文档实用文案1作业任务------------------------------------------------------------------------------------32作业思想---------------------------------------------------------------------------------------33作业条件及数据-------------------

2、-------------------------------------------------34作业过程---------------------------------------------------------------------------35源程序-----------------------------------------------------------------------------46计算结果--------------------------------------------------------

3、-------------------177心得体会与建议-----------------------------------------------------------------------------17标准文档实用文案1作业任务计算近似垂直摄影情况下后方交会解。即利用摄影测量空间后方交会的方法，获取相片的6个外方位元素。限差为0.1。2作业思想利用摄影测量空间后方交会的方法求解。该方法的基本思想是利用至少三个一直地面控制点的坐标A（XA，YA，ZA）、B（XB，YB，ZB）C（XC，YC，ZC），与其影像上对应的三个像点的影像坐

4、标a（xa，ya）、b（xb，yb）、c（xc，yc），根据共线方程，反求该相片的外方位元素XS、YS、ZS、φ、ω、κ。3作业条件及数据已知摄影机主距f=153.24mm，四对点的像点坐标与相应的地面坐标列入下表：表1点号像点坐标地面坐标x（mm）y(mm)X(m)Y(m)Z(m)1-86.15-68.9936589.4125273.322195.172-53.40-82.2137631.0831324.51.728.69314.7876.6339100.9724934.982386.50410.4664.4340426.5430319.8

5、1757.314作业过程4．1获取已知数据相片比例尺1/m=1：10000，内方位元素f=153.24mm，x0，y0；获取控制点的地面测量坐标Xt、Yt、Zt。4．2量测控制点的像点坐标：本次作业中为已知。见表1。4．3确定未知数的初始值：在近似垂直摄影情况下，胶原素的初始值为0，即φ0=ω0=κ0=0；线元素中，ZS0=H=mf=1532.4m，XS0、YS0的取值可用四个控制点坐标的平均值，即：XS0==38437.00YS0==89106.624．4计算旋转矩阵R：利用胶原素的近似值计算方向余弦值，组成R阵。4．5逐点计算像点坐标的近

6、似值：利用未知数的近似值按共线方程式计算控制点像点坐标的近似值（x）（y）。4．6组成误差方程：标准文档实用文案逐点计算误差方程式的系数和常数项。4．7组成法方程式：计算法方程的系数矩阵ATA与常数项ATL。4．8求解外方位元素：根据法方程，由式X=（AtA）-1ATL解求外方位元素改正数，并与相应的近似值求和，得到外方位元素新的近似值。4．9求解外方位元素：将求得的外方位元素的改正数与规定的限差（0.1）比较，小于限差则计算终止，否则用新的近似值重复第4.4至4.8步骤的计算，知道满足要求为止。5源程序#include#i

7、nclude#includeconstdoublePRECISION=1e-5;typedefdoubleDOUBLE[5];intInputData(int&Num,DOUBLE*&Data,double&m,double&f);intResection(constint&Num,constDOUBLE*&Data,constdouble&m,constdouble&f);intInverseMatrix(double*matrix,constint&row);intmain(intargc,char*a

8、rgv[]){DOUBLE*Data=NULL;intNum;doublef(0),m(0);if(InputData(Num,Data,m,f)){if(Data