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



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

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模块来创建体积网格。它应该接收一个曲面网格并输出一个体积网格,如图所示。

根据教程,我应该能够使用创建体积网格

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模块中提出了这个问题,但还没有收到好的反馈。

根据@Paddy Harrison的建议,我使用了meshiowrite函数。具体来说,对于我的数据集,它是这样做的:
import meshio 
points = np.array(verts)
cells = [("triangle", np.array(faces)]
meshio.write_points_cells('out.vtu',points,cells)

相关内容

  • 没有找到相关文章