3D空间球面三角形顶点夹角计算问题求助

问题描述 #

我参考方案写了一个JavaScript函数ang3D2，输入三个向量，以第二个参数为顶点计算另外两个向量的夹角。三次不同顺序调用后得到三个角度，但其和始终是180°。实际输入的向量是球面上的非共面点，球面三角形内角和应该在180°到540°之间。我试过用atan2(crossesMag, dots)替代acos(dots)，但第一个角度计算有误。附上函数和调用代码，求修正方向。

原函数及调用代码 #

function ang3D2(v1, v2, v3) {
      // 复制并归一化输入向量
      let V1 = v1.copy().normalize();
      let V2 = v2.copy().normalize();
      let V3 = v3.copy().normalize();
  
      let Va = p5.Vector.sub(V1,V2).normalize();
      let Vb = p5.Vector.sub(V3,V2).normalize();
  
      let dots = p5.Vector.dot(Va,Vb);
  
      let crosses = p5.Vector.cross(Va,Vb);
      let crossesMag = crosses.mag();// / V2.z; ?
  
      let angle = (Math.acos(dots) * 180.0) / Math.PI; 
      //let angle2 = (Math.atan2(crossesMag,dots) * 180.0) / Math.PI;
  
      return Math.round(angle * 1000) / 1000; // 切换返回acos/atan2测试，atan2结果接近但第一个角度错误，acos仅适用于共面情况

    }//ang3d2()  

    //& 调用逻辑:

    function angs() {
      if (cyl1V && cyl2V && cyl3V) {
        ang1 = ang3D2(cyl2V, cyl1V, cyl3V);
        ang2 = ang3D2(cyl1V, cyl2V, cyl3V);
        ang3 = ang3D2(cyl1V, cyl3V, cyl2V);
        angS = ang1 + ang2 + ang3;
        console.log(angS);
      } 
    } //angs  

问题根源 #

你当前计算的是平面三角形的内角，而非球面三角形的内角。代码中Va = V1-V2、Vb = V3-V2是从顶点V2指向另外两点的向量，它们的夹角是平面上的角，三个角加起来自然是180°，和球面没有关系。

球面三角形的内角定义是：顶点处两个相邻球面边的切线夹角，等价于两个边所在平面（过球心）的二面角。

修正方案 #

正确计算逻辑 #

1. 输入的三个向量本身是球面上的点（从球心出发的向量），无需做减法操作。
2. 以顶点向量V2为基准，计算两个边所在平面的法向量：
   • 平面（球心、V2、V1）的法向量：n1 = cross(V2, V1)
   • 平面（球心、V2、V3）的法向量：n2 = cross(V2, V3)
3. 两个法向量的夹角即为球面三角形在该顶点的内角（注意调整角度范围到0-180°）。
4. 使用atan2(crossMag, dot)替代acos(dot)，避免精度丢失问题。

修正后的函数代码 #

function ang3D2(v1, v2, v3) {
  // 复制并归一化输入向量（确保是单位球面上的点）
  let V1 = v1.copy().normalize();
  let V2 = v2.copy().normalize();
  let V3 = v3.copy().normalize();

  // 计算两个边所在平面的法向量
  let n1 = p5.Vector.cross(V2, V1).normalize();
  let n2 = p5.Vector.cross(V2, V3).normalize();

  // 计算法向量的点积和叉积模长
  let dot = p5.Vector.dot(n1, n2);
  let crossMag = p5.Vector.cross(n1, n2).mag();

  // 计算二面角（即球面内角），转换为角度制
  let angleRad = Math.atan2(crossMag, dot);
  let angleDeg = (angleRad * 180.0) / Math.PI;
  
  // 确保角度在0-180°范围内（球面内角的有效范围）
  angleDeg = angleDeg > 180 ? 360 - angleDeg : angleDeg;

  return Math.round(angleDeg * 1000) / 1000;
}

// 调用逻辑保持不变
function angs() {
  if (cyl1V && cyl2V && cyl3V) {
    ang1 = ang3D2(cyl2V, cyl1V, cyl3V);
    ang2 = ang3D2(cyl1V, cyl2V, cyl3V);
    ang3 = ang3D2(cyl1V, cyl3V, cyl2V);
    angS = ang1 + ang2 + ang3;
    console.log(angS); // 现在结果应在180°-540°之间
  }
}

关键说明 #

• 球面内角的核心是二面角，不是平面向量夹角，这是你之前的核心误区。
• cross(V2, V1)得到的平面法向量，其夹角直接对应球面内角。
• atan2相比acos在角度接近0°或180°时精度更高，避免数值不稳定问题。

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