调试自动生成网格融合图片功能（待完成）

riku

2025-09-03 b2392458ebf42594b9fc5390fda40d7a0a12f923

 src/utils/map/calculate.js

@@ -4,6 +4,167 @@

import { map } from './index_old';

function _prepare4convert(lnglats) {
  var coor = [];
  var maxLength = 1000;
  var start = 0;
  var end = start + maxLength;
  while (end <= lnglats.length) {
    coor.push(lnglats.slice(start, end));
    start += maxLength;
    end += maxLength;
  }
  if (start < lnglats.length) {
    coor.push(lnglats.slice(start));
  }
  return coor;
}

function parse2LngLat(lnglats) {
  // 创建包含4个节点的折线及文字标注
  var path = [];
  lnglats.forEach(function (value) {
    // eslint-disable-next-line no-undef
    path.push(new AMap.LngLat(value[0], value[1]));
  });
  return path;
}

//定义一些常量
const PI = 3.1415926535897932384626;
const a = 6378245.0; //长半轴
const ee = 0.00669342162296594323; //扁率/*** GCJ02 转换为 WGS84* @param lng* @param lat* @returns {*[]}*/

/**
 * 判断是否在国内，不在国内则不做偏移
 * @param lng
 * @param lat
 * @returns {boolean}
 */
function out_of_china(lng, lat) {
  // 纬度3.86~53.55,经度73.66~135.05
  return !(lng > 73.66 && lng < 135.05 && lat > 3.86 && lat < 53.55);
}

function transformlat(lng, lat) {
  // lat = +lat lng = +lng
  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 * PI) + 20.0 * Math.sin(2.0 * lng * PI)) *
      2.0) /
    3.0;
  ret +=
    ((20.0 * Math.sin(lat * PI) + 40.0 * Math.sin((lat / 3.0) * PI)) * 2.0) /
    3.0;
  ret +=
    ((160.0 * Math.sin((lat / 12.0) * PI) + 320 * Math.sin((lat * PI) / 30.0)) *
      2.0) /
    3.0;
  return ret;
}

function transformlng(lng, lat) {
  // lat = +latlng = +lng
  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 * PI) + 20.0 * Math.sin(2.0 * lng * PI)) *
      2.0) /
    3.0;
  ret +=
    ((20.0 * Math.sin(lng * PI) + 40.0 * Math.sin((lng / 3.0) * PI)) * 2.0) /
    3.0;
  ret +=
    ((150.0 * Math.sin((lng / 12.0) * PI) +
      300.0 * Math.sin((lng / 30.0) * PI)) *
      2.0) /
    3.0;
  return ret;
}

/**
 * 高德地图坐标转GPS坐标算法
 */
function gcj02towgs84(lng, lat) {
  // lat = +latlng = +lng
  if (out_of_china(lng, lat)) {
    return [lng, lat];
  } else {
    let dlat = transformlat(lng - 105.0, lat - 35.0);
    let dlng = transformlng(lng - 105.0, lat - 35.0);
    let radlat = (lat / 180.0) * PI;
    let magic = Math.sin(radlat);
    magic = 1 - ee * magic * magic;
    let sqrtmagic = Math.sqrt(magic);
    dlat = (dlat * 180.0) / (((a * (1 - ee)) / (magic * sqrtmagic)) * PI);
    dlng = (dlng * 180.0) / ((a / sqrtmagic) * Math.cos(radlat) * PI);
    let mglat = Math.round((lat * 2 - lat - dlat) * 1000000) / 1000000;
    let mglng = Math.round((lng * 2 - lng - dlng) * 1000000) / 1000000;
    return [mglng, mglat];
  }
}

//从GPS转高德
function wgs84_To_Gcj02(lon, lat) {
  if (out_of_china(lon, lat)) {
    return [lon, lat];
  } else {
    let dLat = transformlat(lon - 105.0, lat - 35.0);
    let dLon = transformlng(lon - 105.0, lat - 35.0);
    let radLat = (lat / 180.0) * PI;
    let magic = Math.sin(radLat);
    magic = 1 - ee * magic * magic;
    let sqrtMagic = Math.sqrt(magic);
    dLat = (dLat * 180.0) / (((a * (1 - ee)) / (magic * sqrtMagic)) * PI);
    dLon = (dLon * 180.0) / ((a / sqrtMagic) * Math.cos(radLat) * PI);
    let mgLat = lat + dLat;
    let mgLon = lon + dLon;
    return [mgLon, mgLat];
  }
}

/**
 * 将gps经纬度转换为高德地图经纬度
 * @param {*} lnglats
 * @param {*} callback
 */
function _convertLatlng(index, coor, lnglats, type = 'gps', callback) {
  if (index < coor.length) {
    coor[index].forEach((c) => {
      let r;
      if (type == 'gps') {
        r = wgs84_To_Gcj02(c[0], c[1]);
      } else {
        r = gcj02towgs84(c[0], c[1]);
      }
      lnglats.push({ lng: r[0], lat: r[1] });
    });
    _convertLatlng(index + 1, coor, lnglats, type, callback);

    // var path = parse2LngLat(coor[index]);
    // // eslint-disable-next-line no-undef
    // AMap.convertFrom(path, type, function (status, result) {
    //   if (result.info === 'ok') {
    //     lnglats.push.apply(lnglats, result.locations);
    //     _convertLatlng(index + 1, coor, lnglats, type, callback);
    //   }
    // });
  } else {
    callback(lnglats);
  }
}

export default {
  /**
   * 根据坐标点、距离和角度，得到另一个坐标点
@@ -21,6 +182,57 @@
    var lng = ((dx / ed + (pos[0] * Math.PI) / 180.0) * 180.0) / Math.PI;
    var lat = ((dy / ec + (pos[1] * Math.PI) / 180.0) * 180.0) / Math.PI;
    return [lng, lat];
  },

  /**
   * 获取两个经纬度之间的角度（0度-360度）
   */
  getAngle(lng_a, lat_a, lng_b, lat_b) {
    var a = ((90 - lat_b) * Math.PI) / 180;
    var b = ((90 - lat_a) * Math.PI) / 180;
    var AOC_BOC = ((lng_b - lng_a) * Math.PI) / 180;
    var cosc =
      Math.cos(a) * Math.cos(b) + Math.sin(a) * Math.sin(b) * Math.cos(AOC_BOC);
    var sinc = Math.sqrt(1 - cosc * cosc);
    var sinA = (Math.sin(a) * Math.sin(AOC_BOC)) / sinc;
    var A = (Math.asin(sinA) * 180) / Math.PI;
    var res = 0;
    if (lng_b > lng_a && lat_b > lat_a) res = A;
    else if (lng_b > lng_a && lat_b < lat_a) res = 180 - A;
    else if (lng_b < lng_a && lat_b < lat_a) res = 180 - A;
    else if (lng_b < lng_a && lat_b > lat_a) res = 360 + A;
    else if (lng_b > lng_a && lat_b == lat_a) res = 90;
    else if (lng_b < lng_a && lat_b == lat_a) res = 270;
    else if (lng_b == lng_a && lat_b > lat_a) res = 0;
    else if (lng_b == lng_a && lat_b < lat_a) res = 180;
    return res;
  },

  /**
   * 获取两经纬度间的距离
   */
  getDistance(lng1, lat1, lng2, lat2) {
    lat1 = lat1 || 0;
    lng1 = lng1 || 0;
    lat2 = lat2 || 0;
    lng2 = lng2 || 0;

    var rad1 = (lat1 * Math.PI) / 180.0;
    var rad2 = (lat2 * Math.PI) / 180.0;
    var a = rad1 - rad2;
    var b = (lng1 * Math.PI) / 180.0 - (lng2 * Math.PI) / 180.0;
    var r = 6378137;
    var distance =
      r *
      2 *
      Math.asin(
        Math.sqrt(
          Math.pow(Math.sin(a / 2), 2) +
            Math.cos(rad1) * Math.cos(rad2) * Math.pow(Math.sin(b / 2), 2)
        )
      );

    return distance;
  },

  /**
@@ -55,5 +267,24 @@
      coors_GD.push(p);
    }
    return coors_GD;
  }
  },

  convertFromGPS: function (gps, callback, type = 'gps') {
    var coor = _prepare4convert(gps);
    _convertLatlng(0, coor, [], type, function (result) {
      var gd = [];
      result.forEach((r) => {
        gd.push([r.lng, r.lat]);
      });
      callback(gd);
    });
  },

  /**
   * 高德地图坐标转GPS坐标算法
   */
  gcj02towgs84,

  //从GPS转高德
  wgs84_To_Gcj02
};