摄影测量后方交会程序

《摄影测量后方交会程序》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、摄影测量后方交会程序（c/c++）输入数据截图：结果截图：程序源代码（其中的矩阵求逆在前面已经有了，链接）：#include#include#includeconstdoublePRECISION=1e-5;typedefdoubleDOUBLE[5];intInputData(int&Num,DOUBLE*&Data,double&m,double&f);intResection(constint&Num,constDOUBLE*&Data,constdouble&m,constdoub

2、le&f);intInverseMatrix(double*matrix,constint&row);intmain(intargc,char*argv[]){DOUBLE*Data=NULL;intNum;doublef(0),m(0);if(InputData(Num,Data,m,f)){if(Data!=NULL){delete[]Data;}return1;}if(Resection(Num,Data,m,f)){if(Data!=NULL){delete[]Data;}return1;}if(Data!=NULL){delete[

3、]Data;}printf("解算完毕...");do{printf("计算结果保存于"结果.txt"文件中""请选择操作（输入P打开结果数据，R打开原始数据，其它退出程序）：");fflush(stdin);//刷新输入流charorder=getchar();if('P'==order

4、

5、'p'==order){system("结果.txt");}elseif('R'==order

6、

7、'r'==order){system("data.txt");}elsebreak;system("cls");}while(1);system

8、("PAUSE");return0;}/***********************************************函数名：InputData*函数介绍：从文件（data.txt）中读取数据,*文件格式如下：*点数m(未知写作0)*内方位元素(fx0y0)*编号xyXYZ*下面是一个实例：40153.24001-86.15-68.9936589.4125273.322195.172-53.4082.2137631.0831324.51728.693-14.78-76.6339100.9724934.982386.50410.

9、4664.4340426.5430319.81757.31*参数：(in/out)Num(点数)，*(in/out)Data(存放数据),m,f,x0,y0*返回值：int，0成功，1文件打开失败，2控制点个*数不足，3文件格式错误*作者：vcrs*完成时间：09-10-4**********************************************/intInputData(int&Num,DOUBLE*&Data,double&m,double&f){doublex0,y0;FILE*fp_input;if(!(fp_inp

10、ut=fopen("data.txt","r"))){return1;}fscanf(fp_input,"%d%lf",&Num,&m);if(Num<4){return2;}fscanf(fp_input,"%lf%lf%lf",&f,&x0,&y0);f/=1000;if(m<0

11、

12、f<0){return3;}Data=newDOUBLE[Num];double*temp=newdouble[Num-1];doublescale=0;inti;for(i=0;i

13、,"%*d%lf%lf%lf%lf%lf",&Data[i][0],&Data[i][1],&Data[i][2],&Data[i][3],&Data[i][4])!=5){return3;}//单位换算成mData[i][0]/=1000.0;Data[i][1]/=1000.0;}//如果m未知则归算其值if(0==m){for(i=0;i

14、ta[i][1]-Data[i+1][1]);scale+=temp[i]/2.0;}m=scale/(Num-1);}fclose(fp_input);delet