空间直角坐标系与大地坐标系转换程序

1、空间直角坐标系与大地坐标系转换程序#in clude<iostream>#in clude<cmath>#in clude<ioma nip>using n ames pace std;#defi ne PI (2.0*asi n(1.0)void mai n() double a,b,c,d1,d2,f1,f2,m1,m2,B,L,H,X,Y,Z,W,N,e;cout<<"请分别输入椭球的长半轴、短半轴（国际单位）"<<e ndl;/cin> >a>>b;a=6378137;/ 以 WGS

2、84为例b=6356752.3142;e=sqrt(a*a-b*b)/a;c=a*a/b;int x;cout<<"请输入0或1，0:大地坐标系到空间直角坐标系;标系"<<endl;1 ：空间直角坐标系到大地坐cin>>x;switch(x)case 0:cout<<"请分别输入该点大地纬度、经度、大地高（国际单位别输入）"<<e ndl;，纬度经度请按度分秒，分cin >>d1>>f1>>m1>>d2>>f2>>m2>

3、;>H;B=P l*(d1+f1/60+m1/3600)/180;L=P l*(d2+f2/60+m2/3600)/180;W=sqrt(1-e*e*si n( B)*si n(B);N=a/W;X=(N+H)*cos(B)*cos(L);Y=(N+H)*cos(B)*si n(L);Z=(N*(1-e*e)+H)*si n(B);cout<<"空间直角坐标系中 X,YZ ,坐标值（国际单位）分别为"<<fixed<<set precisi on( 6)<<X<<""<<fix

4、ed<<set precisi on( 6)<<Y<<""<<fixed<<set precisi on( 6)<<Z<<e ndl;break;case 1:cout<<"请分别输入空间直角坐标系中X,YZ的值（国际单位）"<<endl;double t,m, n, P,k,B0;m=Z/sqrt(X*X+Y*Y);/t0/初值t=m+P*n/sqrt(k+ n*n); /现在为t1，之后代替t2,t3B0=ata n(m);n=Z/sqrt(X

5、*X+Y*Y);P=c*e*e/sqrt(X*X+Y*Y);k=1+(a*a-b*b)/(b*b);B=ata n(t);W=sqrt(1-e*e*si n( B)*si n(B);N=a/W;H=Z/si n(B) - N*(1-e*e);int i;69*q8-09*(p8-8)he869*(p8-8)上 8Qus8/二 d/8*08LH8-ETE8 -qnop乏-(8)SOMA*A+x*x)tbsHHMx/A)ue苍丄d* 卡 1/、(OVA)空 oax)m/、=d+(x/A)uelfl =(OVA)空 ovx)=w=d+(x/A)IQ(OAA)空 ovx)=w-(8)u一s*(8)u一

6、s*卩半 LFbsHM年数注(u*u +呂七bs、u*d+E丄Qi巴P2OL-山OL 寸、糊兆 痙24 8 bH購衣总T片8 甜/M+XOL-山OLA(§-8)sq段H)04Lf=(L-Ld)*60;Lm=(L-Ld)*60-Lf)*60;"<<Bd<<" "<<Bf<<""<<Lf<<"cout<<"大地坐标系中纬度，经度，大地高（国际单位）分别为"<<fixed<<set precisi on(

7、 6)<<Bm<<e ndl<<Ld<<""<<fixed<<set precisi on( 6)<<Lm<<e ndl<<fixed<<set precisi on( 6)<<H<<e ndl;break;运行结果请输入0或I,吐犬地坐标系到空间a角坐标系i i；空间直角坐标系到大地坐标系请分别输入空间直角坐标系中的值（国际单位）-2686431-0253934292457.82254338 &4951 768393尢地坐林系中纬度，经度J大地高（国际单位）分别为37 32 18.540000122 2 25.389009-0.006005Pvess any key to con±inue请输/.0或1,0吠她坐标奈到空闾直争坐标务1;空间fife坐标系到+她坐标系a请分别输入该点经度、大地高（国际单也纬度经度请按度分秒，分别输入）3732