对于我的研究，我需要将 3D 多边形（面）和 3D 多面体（体积）都转换为体素（等距立方体），因此有限差分 PDE 求解器方案可以处理几何。我已经在 python 中有一个工作设置，但它目前非常缓慢且不稳定：
1）为了将多边形转换为体素，我在多边形上放置了一个网格，并将每个网格点分配给一个体素。为了防止“洞”，当我增加体素分辨率时，我必须显着增加网格分辨率。
2）为了将（凸）多面体转换为体素，我在网格的一个隔间上运行，其中包含它，并为所有网格点尝试 sage.all.Polyhedron.contains （如 sagemath）。此例程相当缓慢且有缺陷
AttributeError: 'Polyhedron' object has no attribute '_Vrepresentation'
错误，随机无法为某些边列表构建多面体。
我正在寻找一些pythonic方式
1）通过检查多边形与体素的多边形面的交集来正确地对多边形进行体素化
2) 有一个快速可靠的 polyhedron.contains() 例程。
到目前为止，我发现的唯一内容是 CGAL bindings: polyhedron facet intersection example。有没有人有他们的经验？