空間直角坐標系與大地坐標系轉(zhuǎn)換程序

1、空間直角坐標系與大地坐標系轉(zhuǎn)換程序#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<<"請分別輸入該點大地緯度、經(jīng)度、大地高（國際單位別輸入）"<<e ndl;，緯度經(jīng)度請按度分秒，分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); /現(xià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年數(shù)注(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<<"大地坐標系中緯度，經(jīng)度，大地高（國際單位）分別為"<<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;運行結(jié)果請輸入0或I,吐犬地坐標系到空間a角坐標系i i；空間直角坐標系到大地坐標系請分別輸入空間直角坐標系中的值（國際單位）-2686431-0253934292457.82254338 &4951 768393尢地坐林系中緯度，經(jīng)度J大地高（國際單位）分別為37 32 18.540000122 2 25.389009-0.006005Pvess any key to con±inue請輸/.0或1,0吠她坐標奈到空閭直爭坐標務(wù)1;空間fife坐標系到+她坐標系a請分別輸入該點經(jīng)度、大地高（國際單也緯度經(jīng)度請按度分秒，分別輸入）3732