使用Trimesh计算两个凸网格交集为何返回空结果？

我尝试用Trimesh计算两个凸网格的交集体积，但调用trimesh.Trimesh.intersection()方法时，就算是明显重叠的基础形状，也返回空交集。

举个简单例子，我计算两个偏移放置的立方体的交集：

# 创建偏移变换矩阵
T = np.eye(4)
T[:3, 3] = np.array([0.5, 0.5, 0.5])

# 创建基础立方体网格
mesh1 = trimesh.creation.box(extents=[1.0, 1.0, 1.0])
mesh2 = trimesh.creation.box(extents=[1.0, 1.0, 1.0], transform=T)

# 计算交集
test_intersection = mesh1.intersection(mesh2)

print(test_intersection)
print(test_intersection.volume)

# 可视化
scene = trimesh.Scene([mesh1, mesh2])
scene.show()

运行结果显示两个形状的交集为空：

<trimesh.Trimesh(vertices.shape=(0, 3), faces.shape=(0, 3))>
0.0

但可视化时能明显看到立方体重叠。

为什么网格交集为空？我是不是误解了.intersection()方法的用法？
（环境：Python 3.12.7，pip安装的Trimesh 4.5.2）

问题原因及解决方法 #

1. 核心原因：依赖后端库缺失 #

Trimesh的intersection功能并非完全内置，它依赖外部几何计算库（如openvdb、blender或pygalmesh）。默认情况下如果未安装这些依赖，Trimesh会使用一个简易 fallback 实现，对边界接触的凸形状极易计算出错，返回空集。

2. 具体解决步骤 #

• 安装完整依赖：推荐直接安装Trimesh的全依赖包，或者单独安装适配体积运算的openvdb：

  # 安装Trimesh全依赖
  pip install trimesh[all]
  # 或单独安装openvdb
  pip install pyopenvdb
  

• 显式指定计算后端：安装依赖后，调用intersection时指定后端，确保使用可靠算法：

  test_intersection = mesh1.intersection(mesh2, engine='openvdb')
  

• 验证结果：完成上述操作后，重新运行代码，会得到正确的交集体积（此处应为0.125）。

3. 额外注意事项 #

• 检查网格有效性：确保网格是闭合流形的，可以用mesh.is_watertight属性验证，损坏的网格也会导致交集计算失败。
• 替代计算方式：针对凸形状，也可以使用trimesh.proximity.boolean_intersection，稳定性更强：

  from trimesh.proximity import boolean_intersection
  test_intersection = boolean_intersection([mesh1, mesh2], engine='openvdb')
  

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