空间直角坐标系与大地坐标系转换程序.doc
空间直角坐标系与大地坐标系转换程序includeiostream includecmathincludeiomanip using namespace std; define PI 2.0*asin1.0void main double a,b,c,d1,d2,f1,f2,m1,m2,B,L,H,X,Y,Z,W,N,e;cout请分别输入椭球的长半轴、短半轴(国际单位)endl; cinab; a6378137; 以WGS84为例b6356752.3142;esqrta*a-b*b/a;ca*a/b;int x;cout请输入0或1，0大地坐标系到空间直角坐标系；1空间直角坐标系到大地坐标系endl;cinx; switchx case 0 cout请分别输入该点大地纬度、经度、大地高(国际单位,纬度经度请按度分秒，分别输入)endl; cind1f1m1d2f2m2H; BPI*d1f1/60m1/3600/180; LPI*d2f2/60m2/3600/180; Wsqrt1-e*e*sinB*sinB; Na/W; XNH*cosB*cosL; YNH*cosB*sinL; ZN*1-e*eH*sinB; cout空间直角坐标系中X,Y,Z，坐标值(国际单位)分别为fixedsetprecision6X fixedsetprecision6Y fixedsetprecision6Zendl;break; case 1 cout请分别输入空间直角坐标系中X,Y,Z的值(国际单位)endl; cinXYZ; double t,m,n, P,k,B0; mZ/sqrtX*XY*Y; t0 B0atanm; 初值 nZ/sqrtX*XY*Y; Pc*e*e/sqrtX*XY*Y; k1a*a-b*b/b*b; tmP*n/sqrtkn*n; 现在为t1，之后代替t2,t3. Batant; Wsqrt1-e*e*sinB*sinB; Na/W; HZ/sinB - N*1-e*e; int i; fori1;fabsB-B010E-10;i每一次新的B与上一次计算的B比较，误差小于10E-10 rad B0B; nt; tmP*n/sqrtkn*n;迭代 Batant; Wsqrt1-e*e*sinB*sinB; Na/W; ifX0Y0 LatanY/XPI; ifX0Y0 LatanY/XPI; ifX0Y0 L2*PI-atanY/X; Latan2Y,X; HsqrtX*XY*Y/cosB-N; int Bd,Bf,Ld,Lf; double Bm,Lm; B180*B/PI;B转化为度做单位 BdB; BfB-Bd*60; BmB-Bd*60-Bf*60; L180*L/PI;L转化为度做单位 LdL; LfL-Ld*60; LmL-Ld*60-Lf*60; cout大地坐标系中纬度，经度，大地高(国际单位)分别为Bd Bf fixedsetprecision6BmendlLd Lf fixedsetprecision6Lmendlfixedsetprecision6Hendl; break;运行结果
