gps数据坐标转换.docx

《gps数据坐标转换.docx》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、//笛卡尔坐标系typedefstructtagCRDCARTESIAN{doublex;doubley;doublez;}CRDCARTESIAN,*PCRDCARTESIAN;//typedefCRDCARTESIAN*PCRDCARTESIAN;//大地坐标系typedefstructtagCRDGEODETIC{doublelongitude;//经度doublelatitude;//纬度doubleheight;//大地高,可设为0}CRDGEODETIC;typedefCRDGEODETIC*PCRDGEODETIC;voidCoordCovert::Geodet

2、icToCartesian(PCRDCARTESIANpcc,PCRDGEODETICpcg,doubledSemiMajorAxis,doubledFlattening){doubleB;//纬度度数doubleL;//经度度数doubleL0;//中央经线度数doublel;//L-L0doublet;//tanBdoublem;//ltanBdoubleN;//卯酉圈曲率半径doubleq2;doublex;//高斯平面纵坐标doubley;//高斯平面横坐标doubles;//赤道至纬度B的经线弧长doublef;//参考椭球体扁率doublee2;//椭球第一偏心率

3、doublea;//参考椭球体长半轴//doubleb;//参考椭球体短半轴doublea1;doublea2;doublea3;doublea4;doubleb1;doubleb2;doubleb3;doubleb4;doublec0;doublec1;doublec2;doublec3;intDatum=84;//投影基准面类型：北京54基准面为54，西安80基准面为80，WGS84基准面为84intprjno=0;//投影带号intzonewide=3;doubleIPI=0.;//3.98/180.0B=pcg->latitude;//纬度L=pcg->longitu

4、de;//经度if(zonewide==6){prjno=(int)(L/zonewide)+1;L0=prjno*zonewide-3;}else{prjno=(int)((L-1.5)/3)+1;L0=prjno*3;}if(Datum==54){a=;f=1/298.3;}elseif(Datum==84){a=;f=1/298.;}L0=L0*IPI;L=L*IPI;B=B*IPI;e2=2*f-f*f;//(a*a-b*b)/(a*a);l=L-L0;t=tan(B);m=l*cos(B);N=a/sqrt(1-e2*sin(B)*sin(B));q2=e2/(1-

5、e2)*cos(B)*cos(B);a1=1+(double)3/4*e2+(double)45/64*e2*e2+(double)175/256*e2*e2*e2+(double)11025/16384*e2*e2*e2*e2+(double)43659/65536*e2*e2*e2*e2*e2;a2=(double)3/4*e2+(double)15/16*e2*e2+(double)525/512*e2*e2*e2+(double)2205/2048*e2*e2*e2*e2+(double)72765/65536*e2*e2*e2*e2*e2;a3=(double)15

6、/64*e2*e2+(double)105/256*e2*e2*e2+(double)2205/4096*e2*e2*e2*e2+(double)10359/16384*e2*e2*e2*e2*e2;a4=(double)35/512*e2*e2*e2+(double)315/2048*e2*e2*e2*e2+(double)31185/13072*e2*e2*e2*e2*e2;b1=a1*a*(1-e2);b2=(double)-1/2*a2*a*(1-e2);b3=(double)1/4*a3*a*(1-e2);b4=(double)-1/6*a4*a*(1-e2);c0

7、=b1;c1=2*b2+4*b3+6*b4;c2=-(8*b3+32*b4);c3=32*b4;s=c0*B+cos(B)*(c1*sin(B)+c2*sin(B)*sin(B)*sin(B)+c3*sin(B)*sin(B)*sin(B)*sin(B)*sin(B));x=s+(double)1/2*N*t*m*m+(double)1/24*(5-t*t+9*q2+4*q2*q2)*N*t*m*m*m*m+(double)1/720*(61-58*t*t+t*t*t*t)*N*t*m*m*m*m*m*