资源描述：
《大地坐标与空间直角坐标的转换程序代码.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、#include"stdio.h"#include"math.h"#include"stdlib.h"#include"iostream"#definePI3.1415926535897323doublea,b,c,e2,ep2;intmain(){intm,n,t;doubleRAD(doubled,doublef,doublem);voidRBD(doublehd);voidBLH_XYZ();voidXYZ_BLH();voidB_ZS();voidB_FS();voidGUS_ZS();voidGUS_FS();printf("大地测量学");sp1:printf("请
2、选择功能:");printf("1.大地坐标系到大地空间直角坐标的转换");printf("2.大地空间直角坐标到大地坐标系的转换");printf("3.贝塞尔大地问题正算");printf("4.贝塞尔大地问题反算");printf("5.高斯投影正算");printf("6.高斯投影反算");printf("0.退出程序");scanf("%d",&m);if(m==0)exit(0);sp2:printf("请选择椭球参数(输入椭球序号):");printf("1.克拉索夫斯基椭球参数");printf("2.IUGG_1975椭球参
3、数");printf("3.CGCS_2000椭球参数");printf("0.其他椭球参数(自行输入)");scanf("%d",&n);switch(n){case1:a=6378245.0;b=6356863.0188;c=6399698.9018;e2=0.00669342162297;ep2=0.00673852541468;break;case2:a=6378140.0;b=6356755.2882;c=6399596.6520;e2=0.00669438499959;ep2=0.00673950181947;break;case3:a=6378137.0;
4、b=6356752.3141;c=6399593.6259;e2=0.00669438002290;ep2=0.00673949677547;break;case0:{printf("请输入椭球参数：");printf("长半径a=");scanf("%lf",&a);printf("短半径b=");scanf("%lf",&b);c=a*a/b;ep2=(a*a-b*b)/(b*b);e2=(a*a-b*b)/(a*a);break;}default:printf("输入错误!请重新输入!");gotosp2;}while(1){switch(m){cas
5、e1:BLH_XYZ();break;case2:XYZ_BLH();break;case3:B_ZS();break;case4:B_FS();break;case5:GUS_ZS();break;case6:GUS_FS();break;default:printf("输入错误!请重新输入!");gotosp1;}printf("是否继续进行此功能计算?");printf("（若继续进行此功能计算，则输入1；若选择其他功能进行计算，则输入2；若退出，则输入0.）");scanf("%d",&t);switch(t){case1:break
6、;case2:gotosp1;case0:exit(0);}}}doubleRAD(doubled,doublef,doublem){doublee;doublesign=(d<0.0)?-1.0:1.0;if(d==0){sign=(f<0.0)?-1.0:1.0;if(f==0){sign=(m<0.0)?-1.0:1.0;}}if(d<0)d=d*(-1.0);if(f<0)f=f*(-1.0);if(m<0)m=m*(-1.0);e=sign*(d*3600+f*60+m)*PI/(3600*180);returne;}voidRBD(doublehd){intt;intd
7、,f;doublem;doublesign=(hd<0.0)?-1.0:1.0;if(hd<0)hd=fabs(hd);hd=hd*3600*180/PI;t=int(hd/3600);d=sign*t;hd=hd-t*3600;f=int(hd/60);m=hd-f*60;printf("%d'%d'%lf'",d,f,m);}voidBLH_XYZ(){doubleB,L,H,N,W;doubled,f,m;doubleX,Y,Z;printf("请输入大地坐
