私の研究では、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だけです。誰もそれらを経験したことがありますか？