如何在python中保存三角形网格

Question

我有一组顶点和三角形面，它们一起形成一个三角形网格：

import numpy as np

verts = [[0.1, 1.,  1. ]  [1.,  1.,  0.1]  [1.,  0.1, 1. ]  [1.,  1.,  1.9]  [1.,  1.9, 1. ]
 [1.9, 1.,  1. ] ]
faces = [[ 2,  1,  0]  [ 0,  3,  2]  [ 1,  4,  0]  [ 0,  4,  3]  [ 5,  1,  2]  [ 3,  5,  2]
 [ 5,  4,  1]  [ 4,  5,  3]]

我知道从这些三角形网格中创建的唯一方法是使用Poly3Dcollection

from mpl_toolkits.mplot3d.art3d import Poly3DCollection

myMesh = Poly3DCollection(verts[faces])

接下来，我想使用pygalmesh模块创建体积网格。它应该吸收一个表面网格并输出一个体积网格，如图here所示。

根据教程，我应该能够使用以下方法创建体积网格：

import pygalmesh

mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(
    "elephant.vtu",
    facet_angle=25.0,
    facet_size=0.15,
    facet_distance=0.008,
    cell_radius_edge_ratio=3.0,
    verbose=False
)

但是，当我跑步时：

mesh = pygalmesh.generate_volume_mesh_from_surface_mesh(myMesh)

我已返回：

TypeError: expected str, bytes or os.PathLike object, not Poly3DCollection

我认为这是因为我创建的网格不正确，或者是在使用pygalmesh.generate_volume_mesh_from_surface_mesh之前将其保存为某种形式。我不确定。我在模块github中提出了这个问题，但尚未收到良好的反馈。

Answer 1

[按照@Paddy Harrison的建议，我使用了meshio的write函数。具体来说，对于我的数据集，它是这样完成的：

import meshio 

points = np.array(verts)
cells = [("triangle", np.array(faces)]
meshio.write_points_cells('out.vtu',points,cells)