前方交会实验报告

《前方交会实验报告》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验报告实验名称：立体像对前方交会实验日期：2010.9专业班级：资环学院测绘2008013班学号：200801322姓名：周少良一．实验内容及要求掌握前方交会的远离，编写计算机程序，在计算机上调试，输出结果并对结果进行检查。通过上机调试程序加强动手能力的培养，通过对实验结果的分析，增强综合运用所学知识解决实际问题的能力。实验的具体要求如下：1.源程序名要遵循统一规定：以CH13**命名；2.输入像点坐标及外方位元素，要求一定以文件形式输入、输出数据；3.3编写前方交会程序，编程语言不限；4.编译及运行前方交会程序，输出计算结果；5.提交成果：程序框图，程序源代码，

2、计算结果文件。二．实验过程读入数据，采用文件方法，分别赋给相应的外方位元素要求用户输入像点坐标计算摄影基线分量，Bu,Bv,Bw组成旋转矩阵，R1，R2分别计算u1，v1，w1和u2，v2，w2计算N1，N2输出结果1.程序框图2.程序源代码：（采用vc6.0编译环境，c语言书写）#include#include#includedoubleXs[2],Ys[2],Zs[2],phi[2],omega[2],kappa[2],f;doublem1[3],m2[3];//同名像点坐标doubled1[3],d2[3]

3、;//模型点物方坐标//旋转矩阵voidRotationmatrix(doublephi,doubleomega,doublekappa,doublematrix[]){matrix[0]=cos(phi)*cos(kappa)-sin(phi)*sin(omega)*sin(kappa);matrix[1]=-cos(phi)*sin(kappa)-sin(phi)*sin(omega)*cos(kappa);matrix[2]=-sin(phi)*cos(omega);matrix[3]=cos(omega)*sin(kappa);matrix[4]=cos(o

4、mega)*cos(kappa);matrix[5]=-sin(omega);matrix[6]=sin(phi)*cos(kappa)+cos(phi)*sin(omega)*sin(kappa);matrix[7]=-sin(phi)*sin(kappa)+cos(phi)*sin(omega)*cos(kappa);matrix[8]=cos(phi)*cos(omega);}//矩阵相乘voidMultmatrix(double*A,double*B,double*Result,intm,intp,intn){inti;intj;intk;for(i=0;i

5、

6、*******************************""*************");fprintf(fp_out,"立体像对前方交会程序(C)资源与环境学院班级：""测绘工程2008013学号：200801322姓名：周少良");fprintf(fp_out,"**************************************""**************************************""*************");fprintf(fp_out,"读入的是外方位元素是：");in

7、ti=0;for(i=0;i<2;i++){fprintf(fp_out,"%10.7lft%10.7lft%10.7lft%10.7lft%10.7lft%10.7lft%10.7lf",Xs[i],Ys[i],Zs[i],phi[i],omega[i],kappa[i],f);}fprintf(fp_out,"**************************************""**************************************""**************");fprintf(fp_out,"输入的