求助：地图两点间距指定点精确距离的中间点计算方案

你的核心需求是在两点的大圆航线上，找到距离指定点精确距离的点，现有代码的中间点计算逻辑存在错误，导致无法保证精度。以下是纯JavaScript实现的方案，不依赖任何地图API，前后端通用。

第一步：先实现准确的两点距离计算（Haversine公式） #

首先需要确保能精确计算两点间的大圆距离，这是后续计算的基础：

// 计算两点间的大圆距离（单位：米）
function calculateDistance(lat1, lon1, lat2, lon2) {
  const earthRadius = 6371000; // 地球平均半径，单位米
  const dLat = (lat2 - lat1) * Math.PI / 180;
  const dLon = (lon2 - lon1) * Math.PI / 180;
  
  const a = 
    Math.sin(dLat/2) * Math.sin(dLat/2) +
    Math.cos(lat1 * Math.PI / 180) * Math.cos(lat2 * Math.PI / 180) * 
    Math.sin(dLon/2) * Math.sin(dLon/2); 
  const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
  return earthRadius * c; // 返回距离（米）
}

第二步：实现精确的坐标位移函数 #

方法1：从起点向终点方向移动指定距离 #

这个方法直接计算从起点出发，沿着到终点的方向移动指定距离后的坐标，不管目标点是否在两点之间（超出的话就是延长线）：

// 计算从起点向终点方向移动指定距离后的坐标
function movePointAlongPath(startLat, startLng, endLat, endLng, distance) {
  const earthRadius = 6371000;
  // 转换为弧度
  const lat1 = startLat * Math.PI / 180;
  const lon1 = startLng * Math.PI / 180;
  const lat2 = endLat * Math.PI / 180;
  const lon2 = endLng * Math.PI / 180;

  // 计算从起点到终点的方位角
  const dLon = lon2 - lon1;
  const y = Math.sin(dLon) * Math.cos(lat2);
  const x = Math.cos(lat1)*Math.sin(lat2) - Math.sin(lat1)*Math.cos(lat2)*Math.cos(dLon);
  const bearing = Math.atan2(y, x);

  // 计算新坐标
  const angularDistance = distance / earthRadius;
  const lat3 = Math.asin( Math.sin(lat1)*Math.cos(angularDistance) + 
                          Math.cos(lat1)*Math.sin(angularDistance)*Math.cos(bearing) );
  const lon3 = lon1 + Math.atan2( Math.sin(bearing)*Math.sin(angularDistance)*Math.cos(lat1), 
                                  Math.cos(angularDistance) - Math.sin(lat1)*Math.sin(lat3) );

  // 转换回角度
  return {
    lat: lat3 * 180 / Math.PI,
    lng: lon3 * 180 / Math.PI
  };
}

方法2：两点间按比例插值（仅适用于点在两点之间的情况） #

如果确定要找的点在两点的大圆线段上，可以直接按距离比例插值，计算更高效：

// 两点间按距离比例插值，返回距离起点指定距离的点（必须满足distance ≤ 两点间总距离）
function getPointAtDistance(startLat, startLng, endLat, endLng, distance) {
  const earthRadius = 6371000;
  // 转换为弧度
  const lat1 = startLat * Math.PI / 180;
  const lon1 = startLng * Math.PI / 180;
  const lat2 = endLat * Math.PI / 180;
  const lon2 = endLng * Math.PI / 180;

  // 计算两点间的大圆距离
  const totalDistance = calculateDistance(startLat, startLng, endLat, endLng);
  if (distance > totalDistance) {
    throw new Error("指定距离超过两点间总距离");
  }
  
  // 计算角度比例
  const delta = distance / earthRadius;
  const deltaTotal = totalDistance / earthRadius;

  // 大圆插值公式计算中间点
  const a = Math.sin((deltaTotal - delta) / 2) / Math.sin(deltaTotal / 2);
  const b = Math.sin(delta / 2) / Math.sin(deltaTotal / 2);
  
  const x = a * Math.cos(lat1) * Math.cos(lon1) + b * Math.cos(lat2) * Math.cos(lon2);
  const y = a * Math.cos(lat1) * Math.sin(lon1) + b * Math.cos(lat2) * Math.sin(lon2);
  const z = a * Math.sin(lat1) + b * Math.sin(lat2);
  
  const lat3 = Math.atan2(z, Math.sqrt(x*x + y*y));
  const lon3 = Math.atan2(y, x);

  return {
    lat: lat3 * 180 / Math.PI,
    lng: lon3 * 180 / Math.PI
  };
}

使用示例 #

用你的测试坐标验证：

// 测试坐标
const startLat = 31.5387761;
const startLng = 72.9696664;
const endLat = 31.569235;
const endLng = 72.983557;

// 计算从起点向终点移动1000米后的点
const targetPoint = movePointAlongPath(startLat, startLng, endLat, endLng, 1000);
console.log(targetPoint);

// 验证距离：检查新点到起点的实际距离
const checkDistance = calculateDistance(startLat, startLng, targetPoint.lat, targetPoint.lng);
console.log("实际距离：", checkDistance.toFixed(2), "米"); // 结果应接近1000米

注意事项 #

• 以上方法基于地球是正球体的假设，使用平均半径6371km，精度足够大部分民用地图项目（误差在1%以内）；如果需要测绘级高精度，可改用WGS84椭球参数，但计算逻辑会更复杂。
• 角度与弧度的转换是常见错误点，需确保转换逻辑正确。

内容的提问来源于stack exchange，提问作者Ahmad Tahir