空间后方交会C程序代码.doc

《空间后方交会C程序代码.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

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,constdouble&f

2、);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!

3、=NULL){  delete[]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.tx

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

9、831324.51728.693-14.78-76.6339100.9724934.982386.50410.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,doubl

10、e&m,double&f){doublex0,y0;FILE*fp_input;if(!(fp_input=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;

13、inti;for(i=0;i

14、ata[i][2]-Data[i+1][2])/(Data[i][0]-Data[i+1][0])+