资源描述：
《一个老师给高斯投影正、反算c++源码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、//高斯投影正、反算//////6度带宽 54年北京坐标系//高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres)voidGaussProjCal(doublelongitude,doublelatitude,double*X,double*Y){ intProjNo=0; intZoneWide; ////带宽 double longitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval; double a,f,e2,ee,
2、NN,T,C,A,M,iPI; iPI=0.0174532925199433; ////3.1415926535898/180.0; ZoneWide=6; ////6度带宽 a=6378245.0; f=1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ProjNo=(int)(longitude/ZoneWide); longitude0=ProjNo*ZoneWide+ZoneWide/2; l
3、ongitude0=longitude0*iPI; latitude0=0; longitude1=longitude*iPI; //经度转换为弧度 latitude1=latitude*iPI; //纬度转换为弧度 e2=2*f-f*f; ee=e2*(1.0-e2); NN=a/sqrt(1.0-e2*sin(latitude1)*sin(latitude1)); T=tan(latitude1)*tan(latitude1); C=ee*cos(latitude1)*cos(latitude1)
4、; A=(longitude1-longitude0)*cos(latitude1); M=a*((1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)*latitude1-(3*e2/8+3*e2*e2/32+45*e2*e2*e2/1024)*sin(2*latitude1) +(15*e2*e2/256+45*e2*e2*e2/1024)*sin(4*latitude1)-(35*e2*e2*e2/3072)*sin(6*latitude1)); xval=NN*(A+(1-T+C)*A*A
5、*A/6+(5-18*T+T*T+72*C-58*ee)*A*A*A*A*A/120); yval=M+NN*tan(latitude1)*(A*A/2+(5-T+9*C+4*C*C)*A*A*A*A/24 +(61-58*T+T*T+600*C-330*ee)*A*A*A*A*A*A/720); X0=1000000L*(ProjNo+1)+500000L; Y0=0; xval=xval+X0; yval=yval+Y0; *X=xval; *Y=yval;}//高斯投影由大地坐标(Unit:M
6、etres)反算经纬度(Unit:DD)voidGaussProjInvCal(doubleX,doubleY,double*longitude,double*latitude){ intProjNo; intZoneWide; ////带宽 doublelongitude1,latitude1,longitude0,latitude0,X0,Y0,xval,yval; doublee1,e2,f,a,ee,NN,T,C,M,D,R,u,fai,iPI; iPI=0.0174532925199433; /
7、///3.1415926535898/180.0; a=6378245.0; f=1.0/298.3; //54年北京坐标系参数 ////a=6378140.0; f=1/298.257; //80年西安坐标系参数 ZoneWide=6; ////6度带宽 ProjNo=(int)(X/1000000L); //查找带号 longitude0=(ProjNo-1)*ZoneWide+ZoneWide/2; longitude0=longitude0*iPI; //中央经线 X0=P
8、rojNo*1000000L+500000L; Y0=0; xval=X-X0; yval=Y-Y0; //带内大地坐标 e2=2*f-f*f; e1=(1.0-sqrt(1-e2))/(1.0+sqrt(1-e2)); ee=e2/(1-e2); M=yval; u=M/(a*(1-e2/4-3*e2*e2/64-5*e2*e2*e2/256)); fai
