高斯投影正、反算 代码

 

 
//高斯投影正、反算 
//6度带宽 54年北京坐标系 
//高斯投影由经纬度(Unit:DD)反算大地坐标(含带号，Unit:Metres) 

void GaussProjCal(double longitude, double latitude, double *X, double *Y) 
{
 int ProjNo=0; int ZoneWide; 带宽 
 double longitude1,latitude1, longitude0,latitude0, X0,Y0, xval,yval;
 double a,f, e2,ee, 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; 
 longitude0 = 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