python - 使用 MeshPy 进行 Delaunay 三角剖分

Question

我一直在使用scipy.spatial.Delaunay()到目前为止，我们已经完成了点集的 Delaunay 三角剖分。然而，这通常会导致创建细小的三角形，这是我试图避免的。有人告诉我here那MeshPy为 Shewchuk 的 Triangle 提供包装器，它允许构建高质量、可调整的网格。

但是，我一直无法找到任何使用 Meshpy 创建 Delaunay 三角剖分的示例代码。 ，仅给出节点。 MeshPy网站以 3D 形式给出了这个基本示例:

from meshpy.tet import MeshInfo, build

mesh_info = MeshInfo()
mesh_info.set_points([
    (0,0,0), (2,0,0), (2,2,0), (0,2,0),
    (0,0,12), (2,0,12), (2,2,12), (0,2,12),
    ])
mesh_info.set_facets([
    [0,1,2,3],
    [4,5,6,7],
    [0,4,5,1],
    [1,5,6,2],
    [2,6,7,3],
    [3,7,4,0],
    ])
mesh = build(mesh_info)

但是手动提供面(或三角形)正是我不想做的。

注:Triangle的网站( http://www.cs.cmu.edu/~quake/triangle.delaunay.html ) 显示了 Delaunay 三角剖分的一些快速代码。但我该如何在 MeshPy 中做到这一点？ ？

Answer 1

如果您对三角形网格的质量不满意，也许 optimesh (我的一个项目)值得一看。安装

pip install optimesh

并在命令行上使用它，例如

optimesh in.vtk out.e

您还可以在 Python 中使用它。