Browse Source

天地图便宜纠正

CzRger 1 month ago
parent
commit
33405423e8
2 changed files with 151 additions and 3 deletions
  1. 148 0
      src/components/easyMap/gcj02Mecator.js
  2. 3 3
      src/components/easyMap/initMapInfo.ts

+ 148 - 0
src/components/easyMap/gcj02Mecator.js

@@ -0,0 +1,148 @@
+// 导入proj控件,使用其方法注入gcj02坐标系
+import * as proj from 'ol/proj';
+
+var forEachPoint = function (func) {
+  return function (input, opt_output, opt_dimension) {
+    var len = input.length;
+    var dimension = opt_dimension ? opt_dimension : 2;
+    var output;
+    if (opt_output) {
+      output = opt_output;
+    } else {
+      if (dimension !== 2) {
+        output = input.slice();
+      } else {
+        output = new Array(len);
+      }
+    }
+    for (var offset = 0; offset < len; offset += dimension) {
+      func(input, output, offset);
+    }
+    return output;
+  };
+};
+var gcj02 = {};
+var i = 0;
+var PI = Math.PI;
+var AXIS = 6378245.0;
+var OFFSET = 0.00669342162296594323; // (a^2 - b^2) / a^2
+
+function delta(wgLon, wgLat) {
+  var dLat = transformLat(wgLon - 105.0, wgLat - 35.0);
+  var dLon = transformLon(wgLon - 105.0, wgLat - 35.0);
+  var radLat = (wgLat / 180.0) * PI;
+  var magic = Math.sin(radLat);
+  magic = 1 - OFFSET * magic * magic;
+  var sqrtMagic = Math.sqrt(magic);
+  dLat = (dLat * 180.0) / (((AXIS * (1 - OFFSET)) / (magic * sqrtMagic)) * PI);
+  dLon = (dLon * 180.0) / ((AXIS / sqrtMagic) * Math.cos(radLat) * PI);
+  return [dLon, dLat];
+}
+
+function outOfChina(lon, lat) {
+  if (lon < 72.004 || lon > 137.8347) {
+    return true;
+  }
+  if (lat < 0.8293 || lat > 55.8271) {
+    return true;
+  }
+  return false;
+}
+
+function transformLat(x, y) {
+  var 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 transformLon(x, y) {
+  var 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;
+}
+
+gcj02.toWGS84 = forEachPoint(function (input, output, offset) {
+  var lng = input[offset];
+  var lat = input[offset + 1];
+  if (!outOfChina(lng, lat)) {
+    var deltaD = delta(lng, lat);
+    lng = lng - deltaD[0];
+    lat = lat - deltaD[1];
+  }
+  output[offset] = lng;
+  output[offset + 1] = lat;
+});
+
+gcj02.fromWGS84 = forEachPoint(function (input, output, offset) {
+  var lng = input[offset];
+  var lat = input[offset + 1];
+  if (!outOfChina(lng, lat)) {
+    var deltaD = delta(lng, lat);
+    lng = lng + deltaD[0];
+    lat = lat + deltaD[1];
+  }
+  output[offset] = lng;
+  output[offset + 1] = lat;
+});
+
+var sphericalMercator = {};
+
+var RADIUS = 6378137;
+var MAX_LATITUDE = 85.0511287798;
+var RAD_PER_DEG = Math.PI / 180;
+
+sphericalMercator.forward = forEachPoint(function (input, output, offset) {
+  var lat = Math.max(Math.min(MAX_LATITUDE, input[offset + 1]), -MAX_LATITUDE);
+  var sin = Math.sin(lat * RAD_PER_DEG);
+
+  output[offset] = RADIUS * input[offset] * RAD_PER_DEG;
+  output[offset + 1] = (RADIUS * Math.log((1 + sin) / (1 - sin))) / 2;
+});
+
+sphericalMercator.inverse = forEachPoint(function (input, output, offset) {
+  output[offset] = input[offset] / RADIUS / RAD_PER_DEG;
+  output[offset + 1] = (2 * Math.atan(Math.exp(input[offset + 1] / RADIUS)) - Math.PI / 2) / RAD_PER_DEG;
+});
+
+var projzh = {};
+projzh.ll2gmerc = function (input, opt_output, opt_dimension) {
+  let output = gcj02.fromWGS84(input, opt_output, opt_dimension);
+  return projzh.ll2smerc(output, output, opt_dimension);
+};
+projzh.gmerc2ll = function (input, opt_output, opt_dimension) {
+  let output = projzh.smerc2ll(input, input, opt_dimension);
+  return gcj02.toWGS84(output, opt_output, opt_dimension);
+};
+projzh.smerc2gmerc = function (input, opt_output, opt_dimension) {
+  let output = projzh.smerc2ll(input, input, opt_dimension);
+  output = gcj02.fromWGS84(output, output, opt_dimension);
+  return projzh.ll2smerc(output, output, opt_dimension);
+};
+projzh.gmerc2smerc = function (input, opt_output, opt_dimension) {
+  let output = projzh.smerc2ll(input, input, opt_dimension);
+  output = gcj02.toWGS84(output, output, opt_dimension);
+  return projzh.ll2smerc(output, output, opt_dimension);
+};
+
+projzh.ll2smerc = sphericalMercator.forward;
+projzh.smerc2ll = sphericalMercator.inverse;
+
+// 定义GCJ02
+const gcj02Extent = [-20037508.342789244, -20037508.342789244, 20037508.342789244, 20037508.342789244];
+const gcj02Mecator = new proj.Projection({
+  code: 'GCJ-02',
+  extent: gcj02Extent,
+  units: 'm',
+});
+proj.addProjection(gcj02Mecator);
+// 将4326/3857转为gcj02坐标的方法定义
+proj.addCoordinateTransforms('EPSG:4326', gcj02Mecator, projzh.ll2gmerc, projzh.gmerc2ll);
+proj.addCoordinateTransforms('EPSG:3857', gcj02Mecator, projzh.smerc2gmerc, projzh.gmerc2smerc);
+
+// 我使用的react,所以这里需要导出定义的gcj02Mecator,提供给外部使用
+export default gcj02Mecator;
+

+ 3 - 3
src/components/easyMap/initMapInfo.ts

@@ -1,20 +1,20 @@
 import * as layer from 'ol/layer'
 import * as source from 'ol/source'
+import gcj02Mecator from './gcj02Mecator'
 // @ts-ignore
 import LutuImg from './images/bg-land.png'
-
 const isInternet = window.location.origin.includes('74.10') ? false : true
 
 const baseMapView = {
   center: [109.6915958479584, 19.111636735969318],
-  projection: "EPSG:4326",
+  projection: 'EPSG:4326', // 设置地图的投影
   zoom: 9
   // extent: [120.8953306326286,31.3667480047968,121.37735577911297,31.692561298253832]
 }
 const initInternet = (obj) => {
-
   const _layer = isInternet ? new layer.Tile({
     source: new source.XYZ({
+      projection: gcj02Mecator,
       url: 'http://wprd0{1-4}.is.autonavi.com/appmaptile?lang=zh_cn&size=1&style=7&x={x}&y={y}&z={z}'
     }),
     visible: obj.visible,