154 lines
3.8 KiB
JavaScript
154 lines
3.8 KiB
JavaScript
|
|
// 02转84坐标
|
||
|
|
export function gcj02ToWgs84([lng, lat]) {
|
||
|
|
const a = 6378245.0;
|
||
|
|
const ee = 0.00669342162296594323;
|
||
|
|
const pi = Math.PI;
|
||
|
|
|
||
|
|
if (outOfChina02(lng, lat)) {
|
||
|
|
return [lng, lat];
|
||
|
|
} else {
|
||
|
|
let dlat = transformLat02(lng - 105.0, lat - 35.0);
|
||
|
|
let dlng = transformLng02(lng - 105.0, lat - 35.0);
|
||
|
|
const radLat = (lat / 180.0) * pi;
|
||
|
|
let magic = Math.sin(radLat);
|
||
|
|
magic = 1 - ee * magic * magic;
|
||
|
|
const sqrtMagic = Math.sqrt(magic);
|
||
|
|
dlat = (dlat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * pi);
|
||
|
|
dlng = (dlng * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
|
||
|
|
const mgLat = lat + dlat;
|
||
|
|
const mgLng = lng + dlng;
|
||
|
|
return [lng * 2 - mgLng, lat * 2 - mgLat];
|
||
|
|
}
|
||
|
|
}
|
||
|
|
|
||
|
|
function outOfChina02(lng, lat) {
|
||
|
|
return lng < 72.004 || lng > 137.8347 || lat < 0.8293 || lat > 55.8271;
|
||
|
|
}
|
||
|
|
|
||
|
|
// 经纬度转换
|
||
|
|
function transformLat02(lng, lat) {
|
||
|
|
let ret =
|
||
|
|
-100.0 +
|
||
|
|
2.0 * lng +
|
||
|
|
3.0 * lat +
|
||
|
|
0.2 * lat * lat +
|
||
|
|
0.1 * lng * lat +
|
||
|
|
0.2 * Math.sqrt(Math.abs(lng));
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(6.0 * lng * Math.PI) +
|
||
|
|
20.0 * Math.sin(2.0 * lng * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(lat * Math.PI) + 40.0 * Math.sin((lat / 3.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((160.0 * Math.sin((lat / 12.0) * Math.PI) +
|
||
|
|
320 * Math.sin((lat * Math.PI) / 30.0)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
return ret;
|
||
|
|
}
|
||
|
|
|
||
|
|
function transformLng02(lng, lat) {
|
||
|
|
let ret =
|
||
|
|
300.0 +
|
||
|
|
lng +
|
||
|
|
2.0 * lat +
|
||
|
|
0.1 * lng * lng +
|
||
|
|
0.1 * lng * lat +
|
||
|
|
0.1 * Math.sqrt(Math.abs(lng));
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(6.0 * lng * Math.PI) +
|
||
|
|
20.0 * Math.sin(2.0 * lng * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(lng * Math.PI) + 40.0 * Math.sin((lng / 3.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((150.0 * Math.sin((lng / 12.0) * Math.PI) +
|
||
|
|
300.0 * Math.sin((lng / 30.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
return ret;
|
||
|
|
}
|
||
|
|
|
||
|
|
// 84转02坐标
|
||
|
|
export function wgs84ToGcj02([lng, lat]) {
|
||
|
|
if (outOfChina84(lng, lat)) {
|
||
|
|
return { lng, lat };
|
||
|
|
}
|
||
|
|
let dlat = transformLat84(lng - 105.0, lat - 35.0);
|
||
|
|
let dlng = transformLng84(lng - 105.0, lat - 35.0);
|
||
|
|
const radLat = (lat / 180.0) * Math.PI;
|
||
|
|
let magic = Math.sin(radLat);
|
||
|
|
magic = 1 - 0.00669342162296594323 * magic * magic;
|
||
|
|
const sqrtMagic = Math.sqrt(magic);
|
||
|
|
dlat =
|
||
|
|
(dlat * 180.0) /
|
||
|
|
(((6378245.0 * (1 - 0.00669342162296594323)) / (magic * sqrtMagic)) *
|
||
|
|
Math.PI);
|
||
|
|
dlng =
|
||
|
|
(dlng * 180.0) / ((6378245.0 / sqrtMagic) * Math.cos(radLat) * Math.PI);
|
||
|
|
const gcjLng = lng + dlng;
|
||
|
|
const gcjLat = lat + dlat;
|
||
|
|
return [gcjLng, gcjLat];
|
||
|
|
}
|
||
|
|
|
||
|
|
// 判断是否在中国境内
|
||
|
|
function outOfChina84(lng, lat) {
|
||
|
|
return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
|
||
|
|
}
|
||
|
|
|
||
|
|
// 转换公式所需的常量
|
||
|
|
function transformLat84(x, y) {
|
||
|
|
let ret =
|
||
|
|
-100.0 +
|
||
|
|
2.0 * x +
|
||
|
|
3.0 * y +
|
||
|
|
0.2 * y * y +
|
||
|
|
0.1 * x * y +
|
||
|
|
0.2 * Math.sqrt(Math.abs(x));
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(6.0 * x * Math.PI) + 20.0 * Math.sin(2.0 * x * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(y * Math.PI) + 40.0 * Math.sin((y / 3.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((160.0 * Math.sin((y / 12.0) * Math.PI) +
|
||
|
|
320 * Math.sin((y * Math.PI) / 30.0)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
return ret;
|
||
|
|
}
|
||
|
|
|
||
|
|
function transformLng84(x, y) {
|
||
|
|
let ret =
|
||
|
|
300.0 +
|
||
|
|
x +
|
||
|
|
2.0 * y +
|
||
|
|
0.1 * x * x +
|
||
|
|
0.1 * x * y +
|
||
|
|
0.1 * Math.sqrt(Math.abs(x));
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(6.0 * x * Math.PI) + 20.0 * Math.sin(2.0 * x * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((20.0 * Math.sin(x * Math.PI) + 40.0 * Math.sin((x / 3.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
ret +=
|
||
|
|
((150.0 * Math.sin((x / 12.0) * Math.PI) +
|
||
|
|
300.0 * Math.sin((x / 30.0) * Math.PI)) *
|
||
|
|
2.0) /
|
||
|
|
3.0;
|
||
|
|
return ret;
|
||
|
|
}
|