如何计算3D空间中垂直对齐共面凸多边形的重叠区域？

核心思路 #

因为两个多边形共面，完全可以将3D坐标转换为平面内的2D局部坐标系，把问题简化为2D凸多边形的交集面积计算，避免复杂的3D凸包操作。

具体步骤 #

1. 3D转2D局部坐标系

   • 确定多边形所在平面的局部正交坐标系：
     • 取平面内任意一点（比如其中一个多边形的顶点）作为原点。
     • 在平面内选取两个正交的单位向量作为u轴和v轴：
       1. 取多边形的一条边向量，归一化后作为u轴方向。
       2. 计算平面的法向量（可通过多边形的两条边向量叉乘得到），再用法向量叉乘u轴向量，得到v轴方向并归一化。
     • 将所有3D顶点投影到这个u-v坐标系，得到对应的2D坐标。

2. 计算2D凸多边形的交集
   由于原多边形都是凸的，交集也必然是凸多边形（若存在），两种可行方法：

   • Sutherland-Hodgman裁剪算法：用其中一个多边形的每条边依次裁剪另一个多边形，最终得到交集多边形。该算法逻辑直观，适合凸多边形场景。
   • 顶点+交点收集法：
     1. 找出所有位于对方多边形内部的顶点。
     2. 计算两个多边形所有边对的交点（仅保留线段间的有效交点）。
     3. 对收集到的点集求2D凸包（用Graham扫描法或Andrew算法即可，实现难度远低于3D凸包），得到有序的顶点序列。

3. 用鞋带公式计算面积
   对排序后的2D顶点序列，直接使用鞋带公式计算面积，该面积即为原3D重叠区域的面积（因投影是正交单位投影，无缩放误差）。

针对示例的验证 #

将示例中的3D点转换到2D坐标系（以(0,0,0)为原点，u轴为(√2/2,√2/2,0)，v轴为(0,0,1)）：

• 多边形1的2D坐标：(0,0), (√2,0), (√2,1), (0,1)
• 多边形2的2D坐标：(0.5√2,0), (1.5√2,0), (1.5√2,1), (0.5√2,1)
• 交集的2D坐标：(0.5√2,0), (√2,0), (√2,1), (0.5√2,1)
• 鞋带公式计算面积：(√2 - 0.5√2) * 1 = 0.5√2，与预期形状的面积一致。

关键细节提示 #

• 点在凸多边形内的判断：用同侧法——检查点是否在凸多边形所有边的同一侧（需保证顶点按顺时针/逆时针有序）。
• 线段交点计算：在2D坐标系下用参数法求解，仅保留两条线段之间的有效交点（而非直线交点）。

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