如何求解3D盒体与三角形的相交截面或相交面积？

核心通用算法选择 #

不需要枚举所有相交场景，基于凸多边形裁剪的通用逻辑即可完成全场景计算，推荐方案如下：

• 第一步先做快速相交预校验，用分离轴定理（SAT）过滤完全不相交的场景，仅需要13个轴的投影校验即可完成相交性判断，时间复杂度极低
• 第二步确认相交后，采用Sutherland-Hodgman多边形裁剪算法，以三角形为原始待裁剪多边形，用盒体的6个面作为裁剪面逐次执行裁剪：每次裁剪仅需要判断多边形顶点在裁剪面的内侧/外侧，按规则补充交点生成新的多边形顶点，6个面全部裁剪完成后剩余的多边形就是二者的相交截面
• 第三步计算相交面积，盒体和三角形都是凸几何体，凸凸相交得到的截面一定是共面的凸多边形，直接用凸多边形鞋带公式即可计算出面积，仅需要遍历所有截面顶点坐标即可完成计算

实现注意事项 #

所有坐标计算建议统一转换到盒体的局部坐标系下执行，盒体局部坐标系下6个裁剪面的方程可以简化为x=±x半长、y=±y半长、z=±z半长，能大幅降低面方程判断、交点计算的复杂度，减少计算误差。

核心逻辑伪代码示例 #

// 输入：三角形顶点数组tri[3]，轴对齐/定向盒体参数box（中心点、半长数组、旋转矩阵）
// 输出：相交面积，无相交返回0
float calculate_tri_box_intersection_area(Triangle tri, Box box) {
  // 坐标统一转换到盒体局部坐标系
  Polygon current_poly = tri.convert_to_box_local(box);
  // SAT快速排除不相交情况
  if (!sat_check_tri_box_intersect(current_poly, box.half_extent)) return 0.0f;
  // 逐次用6个盒体面裁剪多边形
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_X_POS, box.half_extent[0]);
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_X_NEG, -box.half_extent[0]);
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_Y_POS, box.half_extent[1]);
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_Y_NEG, -box.half_extent[1]);
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_Z_POS, box.half_extent[2]);
  current_poly = sutherland_hodgman_clip(current_poly, PLANE_Z_NEG, -box.half_extent[2]);
  // 少于3个顶点无有效面积
  if (current_poly.vertex_count < 3) return 0.0f;
  // 鞋带公式计算凸多边形面积
  return convex_polygon_area_by_shoelace(current_poly);
}

内容的提问来源于stack exchange，提问作者Кирилл Ровнер