// 坐标系的转换

var pi_180 = Math.PI / 180;
var _180_pi = 180 / Math.PI;

var projectionTypes = {};
projectionTypes.bj54 = {
  a: 6378245,//长半轴
  e2: 0.006693421622966//第一偏心率平方
};
projectionTypes.wgs84 = {
  a: 6378137,//长半轴
  e2: 0.00669437999013//第一偏心率平方
};

module.exports = {
  // 7坐标系统转换方法
  // sourceType:源投影类型 targetType：目标投影类型 params： 7参数 blh：源大地坐标
  coordinateTransfUse7params(sourceType, targetType, params, blh) {
    sourceType = sourceType.toLowerCase();
    targetType = targetType.toLowerCase();
    var XYZ1 = _BLH2XYZ(sourceType, blh);
    var XYZ2 = _transXYZBy7params(XYZ1, params);
    var blh2 = _XYZ2BLH(targetType, XYZ2);
    console.log(blh2);
  },

  _BLH2XYZ(projectionType, blhObj) { //将大地坐标喜欢换为空间直角坐标系
    var e2 = projectionTypes[projectionType].e2; //第一偏心率平方值
    var a = projectionTypes[projectionType].a; //长半轴
    var N = a / Math.sqrt(1 - e2 * Math.sin(blhObj.b * pi_180) * Math.sin(blhObj.b * pi_180));
    var X = (N + blhObj.h) * Math.cos(blhObj.b * pi_180) * Math.cos(blhObj.l * pi_180);
    var Y = (N + blhObj.h) * Math.cos(blhObj.b * pi_180) * Math.sin(blhObj.l * pi_180);
    var Z = [N * (1 - e2) + blhObj.h] * Math.sin(blhObj.b * pi_180);
    var resultObj = {};
    resultObj.X = X;
    resultObj.Y = Y;
    resultObj.Z = Z;
    return resultObj;
  },

  _transXYZBy7params(xyz0, params) { //根据7参数对空间直角坐标系进行运算
    var x0 = xyz0.X;
    var y0 = xyz0.Y;
    var z0 = xyz0.Z;
    var dx = params.dx;
    var dy = params.dy;
    var dz = params.dz;
    var rx = params.rx;
    var ry = params.ry;
    var rz = params.rz;
    var m = params.m;
    var x1 = dx + (1 + m) * x0 + Math.sin(rz / 60 / 60 * pi_180) * y0 - Math.sin(ry / 60 / 60 * pi_180) * z0;
    var y1 = dy + (1 + m) * y0 - Math.sin(rz / 60 / 60 * pi_180) * x0 + Math.sin(rx / 60 / 60 * pi_180) * z0;
    var z1 = dz + (1 + m) * z0 + Math.sin(ry / 60 / 60 * pi_180) * x0 - Math.sin(rx / 60 / 60 * pi_180) * y0;
    var resultObj = {};
    resultObj.X = x1;
    resultObj.Y = y1;
    resultObj.Z = z1;
    return resultObj;
  },

  _XYZ2BLH(projectionType, xyz1) { 
    var a = projectionTypes[projectionType].a;
    var e2 = projectionTypes[projectionType].e2; //第一偏心率平方值 bj54

    // console.log('哈哈哈哈1-----',xyz1.X,'----', xyz1.Y, '----', xyz1.Z);

    var X = xyz1.X;
    var Y = xyz1.Y;
    var Z = xyz1.Z;

    var L
    if(X >= 0){
      L = Math.atan(Y/X) * _180_pi;
    } else {
      L = Math.atan(Y/X) * _180_pi;
      if(L >= 0){
        L = L - 180;
      } else {
        L = L + 180; 
      }
    }    

    // console.log('哈哈哈哈2-----', L);

    var B = Math.atan(Z / Math.sqrt(X * X + Y * Y));
    while (true) { //迭代计算
      var N = a / Math.sqrt(1 - e2 * Math.sin(B) * Math.sin(B));
      var tempB = Math.atan((N * e2 * Math.sin(B) + Z) / Math.sqrt(X * X + Y * Y));
      B = Math.atan((Z + N * e2 * Math.sin(tempB)) / Math.sqrt(X * X + Y * Y));
      if (tempB == B) {        
        break;
      }
    }
    var H = Z / Math.sin(B) - N * (1 - e2);
    B = B * _180_pi;
    var resultObj = {};
    resultObj.L = L;
    resultObj.B = B;
    resultObj.H = H;
    // if (resultObj.B < 0) {
    //   resultObj.B = resultObj.B + 180;
    // }
    // if (resultObj.L < 0) {
    //   resultObj.L = resultObj.L + 180;
    // }
    // console.log('哈哈哈哈3-----',resultObj.L,'----', resultObj.B)
    return resultObj;
  }
}

记录一下arcgis将国家2000转化为经纬度的过程

1. 将Excel数据导入arcgis，右键选择displayXY，定义好坐标系。
   
2. 将数据导出为shape 格式
   
3. 选择投影转化工具

data management tools -> projections and transformations ->

输入之前导出的shape数据，选择好要转换的投影坐标系，我选择输出的投影坐标为地理坐标系，所以输出的坐标中就变成的地理坐标。

这个意思是会在原有的数据后面增加两列

注意：input coordinate format 要选择shape格式，不要选dd_2，否则会处理失败！

完!

下图3号处的中央子午线度数依据坐标大致定位，本例中坐标为浙江某地，经度大约120


下图中num一列可有可无


注意XY与经纬度对应的顺序