|
@@ -0,0 +1,163 @@
|
|
|
+//转换常数
|
|
|
+const x_pi = (3.14159265358979324 * 3000.0) / 180.0;
|
|
|
+const pi = 3.14159265358979324;
|
|
|
+const a = 6378245.0;
|
|
|
+const ee = 0.00669342162296594323;
|
|
|
+
|
|
|
+function transformLon(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 * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0) / 3.0;
|
|
|
+ ret += ((20.0 * Math.sin(x * pi) + 40.0 * Math.sin((x / 3.0) * pi)) * 2.0) / 3.0;
|
|
|
+ ret += ((150.0 * Math.sin((x / 12.0) * pi) + 300.0 * Math.sin((x / 30.0) * pi)) * 2.0) / 3.0;
|
|
|
+ return ret;
|
|
|
+}
|
|
|
+
|
|
|
+function transformLat(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 * pi) + 20.0 * Math.sin(2.0 * x * pi)) * 2.0) / 3.0;
|
|
|
+ ret += ((20.0 * Math.sin(y * pi) + 40.0 * Math.sin((y / 3.0) * pi)) * 2.0) / 3.0;
|
|
|
+ ret += ((160.0 * Math.sin((y / 12.0) * pi) + 320 * Math.sin((y * pi) / 30.0)) * 2.0) / 3.0;
|
|
|
+ return ret;
|
|
|
+}
|
|
|
+// 判断是否在国内
|
|
|
+function outOfChina(lat, lon) {
|
|
|
+ if (lon < 72.004 || lon > 137.8347) return true;
|
|
|
+ if (lat < 0.8293 || lat > 55.8271) return true;
|
|
|
+ return false;
|
|
|
+}
|
|
|
+
|
|
|
+/*
|
|
|
+ WGS-84转换GCJ-02
|
|
|
+(即 天地图转高德、腾讯地图)
|
|
|
+[{lat: 21, lng: 110}]
|
|
|
+ */
|
|
|
+export const wgs_gcj_encrypts = (latlngs) => {
|
|
|
+ const point = [];
|
|
|
+ for (const latlng of latlngs) {
|
|
|
+ if (outOfChina(latlng.lat, latlng.lng)) {
|
|
|
+ point.push({
|
|
|
+ lat: latlng.lat,
|
|
|
+ lng: latlng.lng
|
|
|
+ });
|
|
|
+ return point;
|
|
|
+ }
|
|
|
+ let dLat = transformLat(latlng.lng - 105.0, latlng.lat - 35.0);
|
|
|
+ let dLon = transformLon(latlng.lng - 105.0, latlng.lat - 35.0);
|
|
|
+ const radLat = (latlng.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);
|
|
|
+ dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
|
|
|
+ const lat = latlng.lat + dLat;
|
|
|
+ const lng = latlng.lng + dLon;
|
|
|
+ point.push({
|
|
|
+ lat: lat,
|
|
|
+ lng: lng
|
|
|
+ });
|
|
|
+ }
|
|
|
+ return point;
|
|
|
+};
|
|
|
+
|
|
|
+/*
|
|
|
+ WGS-84转换GCJ-02
|
|
|
+(即 天地图转高德、腾讯地图)
|
|
|
+[[21,110]]
|
|
|
+ */
|
|
|
+export const wgs_gcj_encrypts2 = (latlngs) => {
|
|
|
+ const point = [];
|
|
|
+ for (const latlng of latlngs) {
|
|
|
+ // if (outOfChina(latlng[0], latlng[1])) {
|
|
|
+ // point.push({
|
|
|
+ // lat: latlng[0],
|
|
|
+ // lng: latlng[1]
|
|
|
+ // });
|
|
|
+ // return point;
|
|
|
+ // }
|
|
|
+ let dLat = transformLat(latlng[1] - 105.0, latlng[0] - 35.0);
|
|
|
+ let dLon = transformLon(latlng[1] - 105.0, latlng[0] - 35.0);
|
|
|
+ const radLat = (latlng[0] / 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);
|
|
|
+ dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
|
|
|
+ const lat = latlng[0] + dLat;
|
|
|
+ const lng = latlng[1] + dLon;
|
|
|
+ point.push([lat, lng]);
|
|
|
+ }
|
|
|
+ return point;
|
|
|
+};
|
|
|
+/*
|
|
|
+ BD-09转换GCJ-02
|
|
|
+ (即 百度转高德、腾讯地图)
|
|
|
+ */
|
|
|
+export const bd_google_encrypt = (latlngs) => {
|
|
|
+ const point = [];
|
|
|
+ for (const latlng of latlngs) {
|
|
|
+ const x = latlng.lng - 0.0065;
|
|
|
+ const y = latlng.lat - 0.006;
|
|
|
+ const z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * x_pi);
|
|
|
+ const theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * x_pi);
|
|
|
+ const gg_lon = z * Math.cos(theta);
|
|
|
+ const gg_lat = z * Math.sin(theta);
|
|
|
+ point.push({
|
|
|
+ lat: gg_lon,
|
|
|
+ lng: gg_lat
|
|
|
+ });
|
|
|
+ }
|
|
|
+ return point;
|
|
|
+};
|
|
|
+
|
|
|
+/*
|
|
|
+ GCJ-02转换BD-09
|
|
|
+ (即 高德、腾讯转百度地图)
|
|
|
+ */
|
|
|
+export const google_bd_encrypt = (latlngs) => {
|
|
|
+ const point = [];
|
|
|
+ for (const latlng of latlngs) {
|
|
|
+ const x = latlng.lng;
|
|
|
+ const y = latlng.lat;
|
|
|
+ const z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * x_pi);
|
|
|
+ const theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * x_pi);
|
|
|
+ const bd_lon = z * Math.cos(theta) + 0.0065;
|
|
|
+ const bd_lat = z * Math.sin(theta) + 0.006;
|
|
|
+ point.push({
|
|
|
+ lat: bd_lat,
|
|
|
+ lng: bd_lon
|
|
|
+ });
|
|
|
+ }
|
|
|
+ return point;
|
|
|
+};
|
|
|
+
|
|
|
+/*
|
|
|
+ GCJ-02 到 WGS-84 的转换
|
|
|
+ (即 高德、腾讯转天地图)
|
|
|
+ */
|
|
|
+export const gcj_wgs_encrypts = (latlngs) => {
|
|
|
+ const point = [];
|
|
|
+ for (const latlng of latlngs) {
|
|
|
+ if (outOfChina(latlng.lat, latlng.lng)) {
|
|
|
+ point.push({
|
|
|
+ lat: latlng.lat,
|
|
|
+ lng: latlng.lng
|
|
|
+ });
|
|
|
+ return point;
|
|
|
+ }
|
|
|
+ let dLat = transformLat(latlng.lng - 105.0, latlng.lat - 35.0);
|
|
|
+ let dLon = transformLon(latlng.lng - 105.0, latlng.lat - 35.0);
|
|
|
+ const radLat = (latlng.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);
|
|
|
+ dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * pi);
|
|
|
+ const lat = dLat;
|
|
|
+ const lng = dLon;
|
|
|
+ point.push({
|
|
|
+ lat: lat,
|
|
|
+ lng: lng
|
|
|
+ });
|
|
|
+ }
|
|
|
+ return point;
|
|
|
+};
|