基于Leaflet的不同半径地理圆交点计算问题求助

问题根源 #

1. 直接套用平面坐标系的两圆交点逻辑，但经纬度属于球面坐标，不同纬度的经度距离差异极大，修改位置或半径后，平面计算的误差会彻底打乱pf的补偿逻辑。
2. 代码59-89段的问题大概率是没做球面坐标到平面投影的转换，就直接对经纬度做几何运算，或者pf参数直接作用在经纬度上，完全不符合球面几何的距离规律。

修复步骤 #

1. 先转平面投影坐标再计算 #

Leaflet默认使用Web墨卡托（EPSG:3857）投影，先把圆心经纬度转成平面坐标（单位为米），所有几何操作都在这个坐标系下完成：

// 经纬度转平面坐标
const proj = L.CRS.EPSG3857;
const center1 = proj.project(L.latLng(lat1, lng1));
const center2 = proj.project(L.latLng(lat2, lng2));
// 半径保持米单位（和投影单位一致）
const r1 = radius1; 
const r2 = radius2;

2. 正确计算平面两圆交点 #

先计算两圆心的平面距离，再按平面几何公式求交点：

const d = Math.hypot(center2.x - center1.x, center2.y - center1.y);

// 先判断是否存在交点
if (d > r1 + r2 || d < Math.abs(r1 - r2) || d === 0) {
  // 处理无交点/重合情况，比如清空标记
  return;
}

// 计算交点的平面坐标
const a = (r1**2 - r2**2 + d**2) / (2 * d);
const h = Math.sqrt(r1**2 - a**2);
const midX = center1.x + (a * (center2.x - center1.x)) / d;
const midY = center1.y + (a * (center2.y - center1.y)) / d;

const int1 = {
  x: midX + (h * (center2.y - center1.y)) / d,
  y: midY - (h * (center2.x - center1.x)) / d
};
const int2 = {
  x: midX - (h * (center2.y - center1.y)) / d,
  y: midY + (h * (center2.x - center1.x)) / d
};

3. 让pf参数正常工作 #

pf作为权重因子，必须作用在平面坐标上，不能直接修改经纬度。比如要平衡两圆心的中间点：

// 用pf加权调整平面中点
const adjustedMid = {
  x: center1.x * pf + center2.x * (1 - pf),
  y: center1.y * pf + center2.y * (1 - pf)
};
// 如果要调整交点，同理在平面坐标上加权

调整完成后再转回经纬度：

const adjustedLatLng = proj.unproject(L.point(adjustedMid.x, adjustedMid.y));

4. 代码59-89段的排查重点 #

• 检查是否直接对经纬度做加减乘除：经纬度是角度单位，不是平面坐标，这类操作完全错误。
• 检查pf的应用位置：如果是在计算交点前调整圆心，必须先转平面坐标再加权；如果是调整交点，要在平面坐标计算完成后再操作。
• 检查半径单位：Leaflet的L.circle默认半径是米，如果代码里用了千米或其他单位，距离计算全错，pf必然失效。

关键提醒 #

所有几何计算必须在平面投影坐标系下完成，球面坐标的非线性特性会让任何直接的经纬度运算出现不可控误差，这也是pf参数失效的核心原因。

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