大地坐標(biāo)與空間坐標(biāo)轉(zhuǎn)換c++程序

1、#mclude#iiiclude#mclude#iiiclude #define E 0.006693422/E=e2define A 6378245/A=a#define PAI 3.141592654using namespace std;hitvoid kd();void dk();int m;cout空間直角坐標(biāo)換算到大地坐標(biāo)輸入1 ；大地坐標(biāo)換算到空間直角坐標(biāo)輸入2 Mendl;cmm;if(m=l)kd();elsedk()；elsecout輸入錯(cuò)誤，重新輸入,endl;retuin 0;void kd()double x,y,z,bLb2,b,h;cout輸入空間直角坐標(biāo)X,Y.

2、Z,endl;cmx;cmy;cmz;b 1 =atan(z/sqrt(x*x+y*y);b2=atan(A*E*sin(bl)/sqrt(lE*pow(sin(bl),2)+z)/(sqrt(x*x+y*y);wliile(fabs(b2-b l)pow( 10.-6) fbl=b2;b2=atan(A*E *sin(b l)/sqrt( 1-E *pow(sui(bl ),2)+z)/(sqrt(x*x+y*y);b=b2;l=atan(y/x);h=sqrt(x*x+y*y)/cos(b)A/sqrt(lE*pow(sin(b),2);b=b*180/PAI;l=l*180/PAI;co

3、utM大地緯度 B=Hfixedsetprecision(8)bendl;coutM 大地經(jīng)度 L=,rfixedsetprecision(8)lendl;coutM大地高 H=Hfixedsetpiecision(3)hendl;void dk()double B、L,H,x,y,z;coutM輸入大地坐標(biāo)B, L, HHendl;cmB;cmL;cmH;B=B*PAI/180;L=L*PAU180;x=(A/sqrt(lE*pow(sin(B),2)+H)*cos(B)*cos(L);y=(A/sqrt(lE*pow(sin(B),2)+H)*cos(B)*sin(L);z=(A*( l-

4、E)/sqrt( 1-E *pow(sin(B),2)+H)*siii(B);coutM空間坐標(biāo)系 X=nfixedsetprecision(3)xendl;/fixed 為小數(shù)點(diǎn)后三位 coutM空間坐標(biāo)系 Y=nyendl;coutM空間坐標(biāo)系 Z=Hzendl;.G:C+ 空間坐標(biāo)轉(zhuǎn)換DebugDADLexb隹間直角坐標(biāo)換算到大地坐標(biāo)輸入1 ；大地坐標(biāo)換算到空間直扇坐標(biāo)輸入2輸入空間直角坐標(biāo)X,Y,Z 4898979.4862828427.1255656854.24?大地緯度 B=45.15313630犬地登度 L=30 . 00000000 大地高 H =1632465.589CD 回Press any key to continue.G:C+空間坐標(biāo)轉(zhuǎn)換 DebugDADLexK空間直角坐標(biāo)換算到大地坐標(biāo)輸入1 ；大地坐標(biāo)換算到空間直角坐標(biāo)輸入2 2輸入大地坐標(biāo)B, L, H 45.15313630301632465589空間坐底系X=4