js大地坐标与经纬度坐标互转

js大地坐标与经纬度坐标互转

//高斯投影坐标反算成经纬度
function GaussToBL(X,Y){
   
    let ProjNo;
    let ZoneWide; 带宽
    let output = new Array(2);
    let longitude1,latitude1, longitude0, X0,Y0, xval,yval;//latitude0,
    let e1,e2,f,a, ee, NN, T,C, M, D,R,u,fai, iPI;
    iPI = 3.14159265358979324/180.0; 3.1415926535898/180.0;
    // a = 6378245.0; f = 1.0/298.3; //54年北京坐标系参数
    a=6378140.0; f=1.0/298.257; //80年西安坐标系参数
    ZoneWide = 6; 6度带宽
    ProjNo = parseInt(X/1000000) ; //查找带号
    longitude0 = (ProjNo-1) * ZoneWide + ZoneWide / 2;
    longitude0 = longitude0 * iPI ; //中央经线


    X0 = ProjNo*1000000+500000;
    Y0 = 0;
    xval = X-X0; yval = Y-Y0; //带内大地坐标
    e2 = 2*f-f*f;
    e1 = (1.0-Math.sqrt(1-e2))/(1.0+Math.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 = u+(<