大地坐标与空间坐标转换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空间直角坐标换算到大地坐标输入1 ；大地坐标换算到空间直角坐标输入2 Mendl;cmm;if(m=l)kd();elsedk()；elsecout输入错误，重新输入,endl;retuin 0;void kd()double x,y,z,bLb2,b,h;cout输入空间直角坐标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 大地经度 L=,rfixedsetprecision(8)lendl;coutM大地高 H=Hfixedsetpiecision(3)hendl;void dk()double B、L,H,x,y,z;coutM输入大地坐标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空间坐标系 X=nfixedsetprecision(3)xendl;/fixed 为小数点后三位 coutM空间坐标系 Y=nyendl;coutM空间坐标系 Z=Hzendl;.G:C+ 空间坐标转换DebugDADLexb隹间直角坐标换算到大地坐标输入1 ；大地坐标换算到空间直扇坐标输入2输入空间直角坐标X,Y,Z 4898979.4862828427.1255656854.24?大地纬度 B=45.15313630犬地登度 L=30 . 00000000 大地高 H =1632465.589CD 回Press any key to continue.G:C+空间坐标转换 DebugDADLexK空间直角坐标换算到大地坐标输入1 ；大地坐标换算到空间直角坐标输入2 2输入大地坐标B, L, H 45.15313630301632465589空间坐底系X=4