经纬度圆心+米级半径的三圆交集计算问题求助

问题分析 #

需要计算三个圆心为经纬度、半径为米的圆的交集，现有平面坐标测试代码正常，但结合经纬度转笛卡尔坐标后返回false，谷歌地图验证参数确实存在交集。核心问题有两个：测试代码的参数错误，以及经纬度转笛卡尔坐标的方式不适合小范围平面计算。

错误点修正 #

1. 测试代码的参数错误 #

你的调用代码中，第二个圆的y坐标错误传入了x坐标值：

// 错误：第二个圆的y坐标用了[0]（x值）
get_cartesian(50, 4.501)[0]
// 正确应该是[1]（y值）
get_cartesian(50, 4.501)[1]

这个错误会导致第二个圆的位置完全偏离，直接触发两圆距离大于半径和的判断，返回false。

2. 经纬度转坐标的优化 #

地球是球面，直接转3D笛卡尔坐标后用平面几何计算交集会引入误差。对于小范围（半径1000米）场景，更适合用局部平面近似，将经纬度转换为以某个点为原点的局部直角坐标系，单位为米，完全适配平面交集计算函数。

修正后的代码 #

局部平面坐标转换函数 #

// 将目标经纬度转换为以原点经纬度为基准的局部平面坐标（单位：米）
function latLonToLocal(lat, lon, originLat, originLon) {
    const earthRadius = 6367000;
    // 转换为弧度
    const latRad = lat * Math.PI / 180;
    const lonRad = lon * Math.PI / 180;
    const originLatRad = originLat * Math.PI / 180;
    const originLonRad = originLon * Math.PI / 180;

    // 计算东西方向（x）和南北方向（y）的偏移量
    const x = earthRadius * (lonRad - originLonRad) * Math.cos(originLatRad);
    const y = earthRadius * (latRad - originLatRad);

    return [x, y];
}

修正后的三圆交集函数（保留原逻辑，修正变量重复声明） #

function calculateThreeCircleIntersection(x0, y0, r0,
                                        x1, y1, r1,
                                        x2, y2, r2){

    var a, horizontal_d_1_2_circle, vertical_d_1_2_circle, distance_1_2_circle, h, rx, ry;
    const EPSILON = 0.0000001;

    horizontal_d_1_2_circle = x1 - x0;
    vertical_d_1_2_circle = y1 - y0;

    distance_1_2_circle = Math.sqrt((vertical_d_1_2_circle*vertical_d_1_2_circle) + (horizontal_d_1_2_circle*horizontal_d_1_2_circle));

    if (distance_1_2_circle > (r0 + r1))
    {
        console.log("Circle1 and Circle2 are separate");
        return false;
    }
    if (distance_1_2_circle < Math.abs(r0 - r1))
    {
        console.log("One circle is inside the other");
        return false;
    }


    a = ((r0*r0) - (r1*r1) + (distance_1_2_circle*distance_1_2_circle)) / (2.0 * distance_1_2_circle) ;

    const point2_x = x0 + (horizontal_d_1_2_circle * a/distance_1_2_circle);
    const point2_y = y0 + (vertical_d_1_2_circle * a/distance_1_2_circle);

    h = Math.sqrt((r0*r0) - (a*a));

    rx = -vertical_d_1_2_circle * (h/distance_1_2_circle);
    ry = horizontal_d_1_2_circle * (h/distance_1_2_circle);

    const intersectionPoint1_x = point2_x + rx;
    const intersectionPoint2_x = point2_x - rx;
    const intersectionPoint1_y = point2_y + ry;
    const intersectionPoint2_y = point2_y - ry;

    console.log("INTERSECTION Circle1 AND Circle2:", "(" + intersectionPoint1_x + "," + intersectionPoint1_y + ")" + " AND (" + intersectionPoint2_x + "," + intersectionPoint2_y + ")");

    const d1 = Math.sqrt(Math.pow(intersectionPoint1_x - x2, 2) + Math.pow(intersectionPoint1_y - y2, 2));
    const d2 = Math.sqrt(Math.pow(intersectionPoint2_x - x2, 2) + Math.pow(intersectionPoint2_y - y2, 2));

    if(Math.abs(d1 - r2) < EPSILON) {
        console.log("INTERSECTION Circle1 AND Circle2 AND Circle3:", "(" + intersectionPoint1_x + "," + intersectionPoint1_y + ")");
        return [intersectionPoint1_x, intersectionPoint1_y];
    }
    else if(Math.abs(d2 - r2) < EPSILON) {
        console.log("INTERSECTION Circle1 AND Circle2 AND Circle3:", "(" + intersectionPoint2_x + "," + intersectionPoint2_y + ")");
        return [intersectionPoint2_x, intersectionPoint2_y];
    }
    else {
        console.log("INTERSECTION Circle1 AND Circle2 AND Circle3:", "NONE");
        return null;
    }
}

正确的测试调用 #

// 以第一个点为原点，转换所有坐标
const originLat = 50;
const originLon = 4.5;

const [x0, y0] = [0, 0];
const [x1, y1] = latLonToLocal(50, 4.501, originLat, originLon);
const [x2, y2] = latLonToLocal(50.001, 4.501, originLat, originLon);

// 调用交集计算函数
calculateThreeCircleIntersection(x0, y0, 1000, x1, y1, 1000, x2, y2, 1000);

说明 #

• 局部平面近似在小范围（几公里内）的误差可以忽略，完全满足需求
• 如果需要处理大范围场景，才需要考虑球面几何的交集计算，但逻辑会复杂很多
• 修正后的函数会返回具体的交点坐标（如果存在），或者返回null

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