如何通过matplotlib绘制四面体网格？

我想matplotlib绘制一个四面体网格，下面是一个简单的四面体网格：如何通过matplotlib绘制四面体网格？

xyz = np.array([ 
    [-1,-1,-1], 
    [ 1,-1,-1], 
    [ 1, 1,-1], 
    [-1, 1,-1], 
    [-1,-1, 1], 
    [ 1,-1, 1], 
    [ 1, 1, 1], 
    [-1, 1, 1]], dtype=np.float) 

tets = np.array([ 
    [0,1,2,6], 
    [0,5,1,6], 
    [0,4,5,6], 
    [0,7,4,6], 
    [0,3,7,6], 
    [0,2,3,6]], dtype=np.int)

当然，在实际应用中，四面体网格中的数量可能很大。我无法在Google中找到任何有用的帮助信息。那么通过matplotlib绘制四面体网格的更好方法是什么？

此外，我可以得到网格的所有三角形面。

tri = np.array([ 
[0 2 1] 
[0 1 5] 
[0 6 1] 
[0 3 2] 
[0 2 6] 
[0 6 3] 
[0 7 3] 
[0 5 4] 
[0 6 4] 
[0 4 7] 
[0 6 5] 
[0 6 7] 
[1 2 6] 
[5 1 6] 
[2 3 6] 
[3 7 6] 
[4 5 6] 
[7 4 6]],dtype=np.int)

来源

2016-12-30 Huayi Wei

阅读：https://mathema.tician.de/software/meshpy/ – eyllanesc

和这个：https：//github.com/inducer/meshpy/blob/master/examples/test_tri_quadratic.py – eyllanesc

matplotlib或许是任务的错误的工具。 3D绘图很难，例如有ParaView可以试用。您可以使用meshio（我的项目）到您的目写入适当的数据类型：

import meshio 
meshio.write('out.vtu', xyz, {'tetra': tets})

来源

2018-01-20 08:58:09

这是一个很好的工作。 'meshio'可以读取和写入obj网格文件吗？ –

从来没有听说过这种格式。 –

以下是wiki https://en.wikipedia.org/wiki/Wavefront_.obj_file。 –

一个可以使用的mpl_toolkits.mplot3d.art3d.Poly3DCollection：

import mpl_toolkits.mplot3d as a3 
axes = a3.Axes3D(pl.figure()) 
vts = xyz[tri, :] 
tri = a3.art3d.Poly3DCollection(vts) 
tri.set_alpha(0.2) 
tri.set_color('grey') 
axes.add_collection3d(tri) 
axes.plot(point[:,0], point[:,1], point[:,2], 'ko') 
axes.set_axis_off() 
axes.set_aspect('equal') 
pl.show()

来源

2017-01-12 23:19:16

如何通过matplotlib绘制四面体网格？

回答

相关问题