高斯平面直角坐标系西安80与WGS84经纬度互转C#程序

注意区分3度带和6度带

六度带中央经线经度的计算：当地中央经线经度＝6°×当地带号－3°

三度带中央经线经度的计算:中央经线经度＝3°×当地带号

1) WGS84经纬度转高斯平面直角坐标系xy

    int ProjNo = 0;
    double longitude1, latitude1, longitude0, latitude0, X0, Y0;
    double e2, ee, NN, T, C, A, M, iPI;
    iPI = 0.0174532925199433; //3.1415926535898/180.0;
    int ZoneWide = 3; //按照西安80高斯投影直角坐标系，经度小于121.5°，3°带
    double _a = 6378137.0;  //80年西安坐标系参数
    double _f = 1 / 298.257223563; //80年西安坐标系参数
    ProjNo = (int)(longitude / ZoneWide);
    longitude0 = ZoneWide == 3 ? ProjNo * 3 : ProjNo * 6 - 3;
    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)