2025.9.4

riku

2025-09-04 7f6661cca40e3530111d628222fa25462022ec78

 src/utils/map/calculate.js

@@ -94,20 +94,72 @@
}

/**
 * 高德地图坐标转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) {
    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);
    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);
  }
@@ -130,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;
  },

  /**
@@ -180,22 +283,43 @@
  /**
   * 高德地图坐标转GPS坐标算法
   */
  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];
    }
  gcj02towgs84,

  //从GPS转高德
  wgs84_To_Gcj02,

  /**
   * 计算一组经纬度坐标的中心点
   * @param {Array} coordinates - 经纬度数组，格式: [{lng: number, lat: number}, ...]
   * @returns {Object} 中心点坐标 {lng: number, lat: number}
   */
  calculateCenterCoordinates(coordinates) {
    if (coordinates.length === 0) return { lng: 0, lat: 0 };
    if (coordinates.length === 1) return coordinates[0];

    let x = 0,
      y = 0,
      z = 0;
    coordinates.forEach((p) => {
      const lng = (p.lng * Math.PI) / 180; // 经度转弧度
      const lat = (p.lat * Math.PI) / 180; // 纬度转弧度

      // 转换为三维坐标
      x += Math.cos(lat) * Math.cos(lng);
      y += Math.cos(lat) * Math.sin(lng);
      z += Math.sin(lat);
    });

    // 取平均值
    const avgX = x / coordinates.length;
    const avgY = y / coordinates.length;
    const avgZ = z / coordinates.length;

    // 转换回经纬度
    const lng = (Math.atan2(avgY, avgX) * 180) / Math.PI;
    const hypotenuse = Math.sqrt(avgX ** 2 + avgY ** 2);
    const lat = (Math.atan2(avgZ, hypotenuse) * 180) / Math.PI;

    return { lng, lat };
  }
};