PyVista 中 voxelize 方法的边界处理优化

在 PyVista 项目中，voxelize() 和 voxelize_volume() 方法用于将三维网格模型转换为体素表示。近期社区发现这两个方法在处理网格边界时存在一个值得优化的行为：默认情况下会裁剪掉网格的末端边界。

问题描述

当对一个简单的立方体网格（边界范围 [-0.5, 0.5]）进行体素化时，使用默认参数会发现输出结果丢失了部分边界。例如，使用密度参数 0.2 进行体素化后，输出网格的边界变为 [-0.5, 0.3]，而不是预期的 [-0.5, 0.5]。

这种行为源于方法内部使用了 numpy.arange() 来生成体素网格，而 arange() 函数默认不包含终点值。这种实现方式会导致网格的末端边界被裁剪掉，在某些应用场景下可能不符合用户预期。

技术分析

在计算机图形学和科学计算中，体素化过程通常需要考虑以下几种边界处理方式：

1. 包含边界：确保原始网格的所有顶点都包含在体素化结果中
2. 排除边界：当前 PyVista 的默认行为
3. 完全包围：通过适当扩展边界确保原始网格被完全包含

PyVista 目前通过 enclosed=True 参数支持第三种方式，但缺乏对第一种方式的直接支持。

解决方案

社区提出了两种改进方案：

1. 短期方案：添加 include_end_bound 参数，允许用户选择是否包含末端边界

   voxelize(mesh, density=None, check_surface=True, enclosed=False, include_end_bound=False)
   

2. 长期方案：设计更灵活的边界控制接口，例如：

   voxelize(mesh, density=None, check_surface=True, range=["include", "exclude"])
   

   这种设计可以支持：

   • 显式指定边界值
   • 选择包含/排除边界
   • 自动扩展边界确保完全包围

实现进展

目前已经提交了两个相关实现：

1. 为 voxelize() 方法添加了 include_end_bound 参数
2. 为 voxelize_volume() 方法添加了相同的功能

这些改进将首先作为可选功能引入，待社区验证效果后，可能会考虑将其设为默认行为。

应用建议

对于需要精确控制体素化边界的应用场景，建议：

1. 如果要求保留原始网格的所有几何信息，使用 include_end_bound=True
2. 如果要求确保原始网格完全被体素网格包围，使用 enclosed=True
3. 对于需要特殊边界处理的情况，可以等待未来更灵活的边界控制接口

这一改进将使得 PyVista 的体素化功能更加灵活和精确，能够满足更广泛的科学计算和可视化需求。