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);
  }
@@ -231,41 +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(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];
    }
  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 };
  }
};