// 坐标系的转换
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;
}
}