贝塞尔大地坐标解算

1、贝塞尔大地坐标解算(正算)#include#mclude#define PI 3.1415926535897323hitmt q.p;double h,i,j,L,B,Al,S;double a,g,c,d,e,f;void n(double b.double Ldouble al,double s,double a,double d,double e,double f);double z(double d.double m,double s);void i(double dd);pnntfCiW选择:nl.我国1954北京坐标系n2.我国1980西安坐标系n3.我国2000国家大 地坐标系1

2、1”)；scanf(d”.&p);lf(p=l)/Ia=6378245.0;g=6356863.0188;c=6399698.9018;d=0.00335232986926;e=0.00669342162297;M.00673852541468;else if(p=2)a=6378140.0;g=6356755.2882;c=6399596.6520;d=0.0033528131779;e=0.00669438499959;M.00673950181947;else if(p=3)a=6378137.0;g=6356752.3141;c=6399593.6259;d=l/298.2572221

3、01;e=0.0066943800229;A0.00673949677547;pnntf(”请输入注：输入格式为角度(例如：30,40501心=”);scanf(”lf%lf%lf”,&h,&i,&j);L=z(h,i,j);pnntf(”请输入注：输入格式为角度(例如：30,405(r)iiBl=”);scanf(”lf%lf%lf”,&h,&i,&j);B=z(h,i,j);pnntf(”请输入注：输入格式为角度(例如：304050加Al=”);scanf(”lf%lf%lf”,&h,&i,&j);Al=z(h,i,j);print.” 请输入 nS=”)；scanf(”lf,&S);n(

4、B,L,Al,S,a,d,e,f);scanf(”d”,&q);double z(double d.double m,double s)double e;double sign=(d0.0)?-l .0:1.0;iRd=0)sign=(m0.0)?-l .0:1.0;iRm=0)sign=(s0.0)?-1.0:1.0;)if(d0)d=d*(-L0);if(s0)s=s*(1.0);e=sign*(d*3600+m*60+s)*PV(3600*180);return e;void r(double dd)int de;int d.m;double s;double sign=(dd0.0)?

5、-l.0:1.0:if(dd0)dd=fabs(dd);dd=dd*3600*180PI;de=iiit(dd/3600);d=sign*de;dd=dd-de*3600;m=iiit(dd/60);s=dd-m*60;printf( ”d%ddnu s);void n(double b,double Ldouble aLdouble s,double a.double d.double eJouble f) double u Lm,M,A2,u2,ss,ss l,k,aa,bbxcjr,ril,rf2,B2X2;double aa2,bb2,cc2,k2l;u 1 =atan(sqit( l

6、-e)*tan(b);m=asiii(cos(u 1) *siii(a 1);if(m0.0)m=m+2*PI;M=ataii(tan(u l)/cos(al);if(M(2.8*PI/180*pow(10.0,-7.0);)ssl=ss;ss=a a*s+bb*siii(ssl)*cos(2*M+ssl)+cc*siii(2*ssl)*cos(4*M+2*ssl);A2=atan(tan(m)/cos(M+ss);if(A20.0)A2=A2+PI;if(al=PI)A2=A2+PI;u2=atan(-cos(A2)*tan(M+ss);itI =atan(siii(u 1 )*tan(a

7、1);if(nl=PI)nl=rrl+PI;n2=atan(siii(u2)*tan(A2);if(nl=PI)n2=rr2+PI;elseiR(M+ss)v=Pl)n2=rr2+PI;n-rr2-rrl;B2=atan(sqrt( 1 +f)*tan(u2);k2=e*cos(m)*cos(m);aa2=(e/2+e*e/8+e*e*e/16)-e*(l+e)*k2/16+3*e*k2*k2/128;bb2=e* (1+e)*k2/ 16-e *k2 *k2/3 2;cc2=e*k2*k2/256;H=n:-siii(m)*(aa2*ss+bb2*sm(ss)*cos(2*M+ss)+cc2

8、*sin(2*ss)*cos(4*M+2*ss);L2=l+ll;pnntf(”L2=”)；r(L2)；printf(”uiB2=)；r(B2);printf(MnA2=H);r(A2)；反算#iiiclude#iiiclude#define PI 3.1415926535897323hitmt q.p;double h,i,j,L,B,L2.B2;double a,g,c,d,e,f;void n(double b.double Ldouble 12,double b2,double ajouble d.double ejouble f);double z(double d.double m

9、,double s);void i(double dd);pnntf(”请选择:nl.我国1954北京坐标系n2.我国1980西安坐标系n3.我国2000国家大 地坐标系11”)；scanf(”d”.&p);if(p=l)a=6378245.0;g=6356863.0188;c=6399698.9018;d=0.00335232986926;e=0.00669342162297;仁0.00673852541468;else iRp=2)a=6378140.0;g=6356755.2882;c=6399596.6520;d=0.0033528131779;e=0.00669438499959;仁

10、0.00673950181947;else iRp=3)a=6378137.0;g=6356752.3141;c=6399593.6259;d=l/298.257222101;e=0.0066943800229;仁0.00673949677547;pnntf(”请输入注：输入格式为角度(例如：304050加Ll=”);L=z(h,iJ);pnntf(”请输入 注：输入格式为角度(例如：30,40卯)何1=”);B=z(h,i,j);pnntf(”请输入注：输入格式为角度(例如：3040*50)nL2=);L2=z(h,ij);pnntf(”请输入注：输入格式为角度(例如：3040,50)nB2

11、=);B2=z(h,i,j);n(B,L,L2.B2,a.d,e,f);scanf(”s”.&q);double z(double d.double m,double s)double e;double sign=(d0.0)?-l .0:1.0;if(d=O)sign=(mO.O)?-l .0:1.0;if(m=0)(sign=(s0.0)?-1.0:1.0;)if(d0)d=d*(-L0);iRm0)if(s0)s=s*(1.0);e=sign*(d*3600+m*60+s)*PI/(3600*180);return e;void r(double dd)int de;int dan;do

12、uble s;double sign=(dd0.0)?-l .0:1.0;if(dd0)dd=fabs(dd);dd=dd*3600*180/PI;de=iiit(dd/3600);d=sign*de;dd=dd-de*3600;m=iiit(dd/60);s=dd-m*60;printf( ”d%d，lf dm, s);void n(double b.double Ldouble 12,double b2,double a,double djouble e,double f) double i】Lm.NLA2,u2,ss,ssLk,aa,bb,cc,ii, 口QB2.L2,m0,A10,Al

13、,S;double aa2,bb2,cc2,k2l;u 1 =atan(sqit( l-e)*tan(b);u2=ataii(sqrt( l-e)*tan(b2);ss=acos(sin(u l)*sin(u2)+cos(u 1 )*cos(u2)*cos(12-l);m0=asin(cos(ul)*cos(u2)*sin(12-l)/sin(ss);n 0=12-1+0.0033 51 *ss*sin(m0);k2=e*cos(m0)*cos(m0);aa2=(e/2+e*e/8+e*e*e/16)-e*(l+e)*k2/16+3*e*k2*k2/128;ssl=ss+sm(m0)*aa2*

14、ss*sm(m0);m=asin(cos(u 1 )*cos(u2)*sin(rr0)/siii(ss 1);A10=ataii(sin(ji0)/(cos(u 1 )*tan(u2)-sin(u 1 )*cos(rrO);if(A10=0.0)A10=A10+PI;if(m=0.0)A10=A10+PI;M=ataii(siii(u 1) * tan( A10)/s iii(m);if(M=0.0)M=M+PI;zk2=e*cos(m)*cos(m);aa2=(e/2+e*e/8+e*e*e/16)-e*(l+e)*k2/16+3*e*k2*k2/128;bb2=e*(l+e)*k2/16-e*k2*k2/32;n-12 -l+sin(m) * (aa2 *ss 1 +bb2 *sin(ss 1 )*cos(2 *M+ss 1);ss=acos(sin(u l)*sin(u2)+cos(u 1 )*cos(u2)*cos(rr);A1 =atan(sin(iT)/(cos(u 1 )*tan(u2)-sin(u 1 )*cos(rr);if(Al=0.0)A1=A1+PI;if(m=0.0)A1=A1+PI;A2=atan(sin(iT)/(-cos(u2) *tan(u 1 )+sin(u2)*cos(iT);if(A2=0.0)