2D空间三点夹角计算结果随向量方向异常变化的问题排查

问题描述 #

需要计算2D空间中线段AB和BC形成的夹角，但现有代码的计算结果会随向量方向变化出现异常，无法正确反映实际的方向夹角。

现有实现代码 #

function getAngle2D(A, B, C) {
        const A2D = { x: A.x, y: A.y };
        const B2D = { x: B.x, y: B.y };
        const C2D = { x: C.x, y: C.y };

        // 向量AB和BC
        const AB = { x: B2D.x - A2D.x, y: B2D.y - A2D.y };
        const BC = { x: C2D.x - B2D.x, y: C2D.y - B2D.y };

        // 向量AB和BC的点积
        const dotProduct = AB.x * BC.x + AB.y * BC.y;

        // 向量AB和BC的模长
        const magnitudeAB = Math.sqrt(AB.x ** 2 + AB.y ** 2);
        const magnitudeBC = Math.sqrt(BC.x ** 2 + BC.y ** 2);

        // 计算夹角的余弦值
        const cosTheta = dotProduct / (magnitudeAB * magnitudeBC);

        // 计算弧度制夹角
        const angleRadians = Math.acos(cosTheta);

        // 转换为角度制
        const angleDegrees = angleRadians * (180 / Math.PI);

        return angleDegrees;
    }

测试点集 #

• 第一组点（绿色）

const points1 = [
    {x: 0.5067973136901855, y: 0.5480734705924988},
    {x: 0.5307877659797668, y: 0.5041729211807251},
    {x: 0.5793206095695496, y: 0.5759181380271912}
];

• 第二组点（蓝色）

const points2 = [
    {x: 0.3797118365764618, y: 0.7955537438392639},
    {x: 0.3800569176673889, y: 0.7407315969467163},
    {x: 0.4513204097747803, y: 0.7465168237686157}
];

夹角示意图 #

问题原因 #

现有代码仅用点积和Math.acos计算夹角，存在两个核心问题：

1. 向量选取错误：使用AB（从A到B）和BC（从B到C）向量，实际计算的是向量AB与BC的夹角，而非顶点B处∠ABC的夹角
2. 无法区分方向：Math.acos仅返回0-180度的最小夹角，无法判断夹角是顺时针还是逆时针方向，导致向量方向变化时结果不符合预期

修正后的代码 #

function getAngle2D(A, B, C) {
    // 提取顶点B出发的两个向量：BA和BC
    const BA = { x: A.x - B.x, y: A.y - B.y };
    const BC = { x: C.x - B.x, y: C.y - B.y };

    // 计算点积和叉积
    const dotProduct = BA.x * BC.x + BA.y * BC.y;
    const crossProduct = BA.x * BC.y - BA.y * BC.x;

    // 计算向量模长（用Math.hypot避免开平方精度问题）
    const magnitudeBA = Math.hypot(BA.x, BA.y);
    const magnitudeBC = Math.hypot(BC.x, BC.y);

    // 使用atan2结合叉积和点积，直接得到带方向的弧度值（范围-π到π）
    let angleRadians = Math.atan2(crossProduct, dotProduct);
    
    // 转换为0-360度的正角度
    let angleDegrees = angleRadians * (180 / Math.PI);
    if (angleDegrees < 0) {
        angleDegrees += 360;
    }

    return angleDegrees;
}

修正说明 #

• 向量调整：改用从顶点B出发的BA和BC向量，准确计算∠ABC的夹角
• 方向判断：通过叉积符号确定夹角方向，Math.atan2返回的结果能区分顺时针和逆时针角度
• 范围转换：将负角度转换为0-360度的正角度，符合常规的方向角需求

内容的提问来源于stack exchange，提问作者Stéphane M.