大地測量學(xué)坐標換算系統(tǒng)

1、void C坐標括換算系統(tǒng)Dlg:OnBnClickedOk()UpdateData(TRUE);CDialogEx:OnOK();b=(1-rfa)*a; e1=sqrt(a*a-b*b)/a;e2=sqrt(a*a-b*b)/b;B=(double)BSecond/3600+(double)BMinute/60+(double)BDegree;L=(double)LSecond/3600+(double)LMinute/60+(double)LDegree;BRadian=B/180*PI;LRadian=L/180*PI;if(number=3)L0=(int)(L/3+0.5)*3;e

2、lseL0=(int)(L/6)+1)*6-3;l=(L-L0)*3600.00;double m0=a*(1-e1*e1);double m2=3*e1*e1*m0/2;double m4=5*e1*e1*m2/4;double m6=7*e1*e1*m4/6;double m8=9*e1*e1*m6/8;double a0=m0+m2/2+3*m4/8+5*m6/16+35*m8/128;double a2=m2/2+m4/2+15*m6/32+7*m8/16;double a4=m4/8+3*m6/16+7*m8/32;double a6=m6/32+m8/16;double a8=m8

3、/128;double W=sqrt(1-e1*e1*sin(BRadian)*sin(BRadian);double N=a/W;double n=e2*cos(BRadian);double t=tan(BRadian);double p=180*3600/PI;Double X=a0*BRadian-a2*sin(2*BRadian)/2+a4*sin(4*BRadian)/4-a6*sin(6*BRadian)/6+a8*sin(8*BRadian)/8;double Item1=(N*sin(BRadian)*cos(BRadian)*l*l)/(2*p*p);Double Item

4、2=(N*sin(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*(5-t*t+9*n*n+4*n*n*n*n)*l*l*l*l)/(24*p*p*p*p);double Item3=(N*sin(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*(61-58*t*t+t*t*t*t)*l*l*l*l*l*l)/(720*p*p*p*p*p*p);double Item4=N*cos(BRadian)*l/p;double Item5=N*cos(B

5、Radian)*cos(BRadian)*cos(BRadian)*(1-t*t+n*n)*l*l*l/(6*p*p*p);double Item6=N*cos(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*cos(BRadian)*(5-18*t*t+t*t*t*t+14*n*n-58*n*n*t*t)*l*l*l*l*l/(120*p*p*p*p*p);double x=X+Item1+Item2+Item3;double y=Item4+Item5+Item6;Result result;result.X=x;result.Y=y;dou

6、ble Bf,Bf1,BfRadian,tf,Mf,Nf,nf;Bf=x/a0;BfRadian=Bf/180*PI;Bf1=(x+sin(2*BfRadian)*a2/2-sin(4*BfRadian)*a4/4+sin(6*BfRadian)*a6/6+sin(8*BfRadian)*a8/8)/a0;while(Bf1-Bf <0.00001)Bf=Bf1;BfRadian=Bf/180*PI;Bf1=(x+sin(2*BfRadian)*a2/2-sin(4*BfRadian)*a4/4+sin(6*BfRadian)*a6/6+sin(8*BfRadian)*a8/8)/a0;

7、Bf=Bf1;BfRadian=Bf/180*PI;tf=tan(BfRadian);double Wf=sqrt(1-e1*e1*sin(BfRadian)*sin(BfRadian);Nf=a/Wf;Mf=a*(1-e1*e1)/(Wf*Wf*Wf);nf=e2*cos(BfRadian);double item1=tf*y*y/(2*Mf*Nf);double item2=tf*(5+3*tf*tf+nf*nf-9*nf*nf*tf*tf)*y*y*y*y/(24*Mf*Nf*Nf*Nf);double item3=tf*(61+90*tf*tf+45*tf*tf*tf*tf)*y*y*

8、y*y*y*y/(720*Mf*Nf*Nf*Nf*Nf*Nf);FSB=Bf-item1+item2-item3;double item4=y/(Nf*cos(BfRadian);double item5=(1+2*tf*tf+nf*nf)*y*y*y/(6*Nf*Nf*Nf*cos(BfRadian);double item6=(5+28*tf*tf+24*tf*tf*tf*tf)*y*y*y*y*y/(120*Nf*Nf*Nf*Nf*Nf*cos(BfRadian);FSL=(item4-item5+item6)+L0;result.BDegree1=(int)FSB;result.BMinute1=(int)(FSB-result.BDegree1)*60);result.BSecond1=(FSB-result.BDegree1)*60-result.BMinute1)*60;result.LDegree=