如何从 Python 中的 3D 点的 Delaunay 三角剖分生成四面体？

Question

我需要对一组 3D 点进行 Delaunay 三角剖分。我为它写了一个脚本（如下），但似乎输出中没有四面体。请给我一些意见/想法。我正在使用 Python3。非常感谢。

from scipy.spatial import Delaunay
import matplotlib.pyplot as plt
import numpy as np
points= np.array([[1,2,2],[1,3,6],[4,3,4],[5,3,2]])
tri= Delaunay(points)
fig= plt.figure()
ax= fig.gca(projection= '3d')
ax.plot_trisurf(points[:,0],points[:,1],points[:,2],triangles= tri.simplices)
plt.plot(points[:,0],points[:,1],points[:,2],'+')
plt.show()

Answer 1

四面体在tri.simplices 成员中给出，该成员包含一个n x 4 索引数组（n 是四面体的数量）。四面体以一组四个索引的形式给出，它们对应于points 数组中四面体的四个点的索引。

例如以下代码将绘制第一个四面体的线框：

tr = tri.simplices[0]  # indices of first tetrahedron
pts = points[tr, :]  # pts is a 4x3 array of the tetrahedron coordinates

# plotting the six edges of the tetrahedron
for ij in [[0, 1], [0, 2], [0, 3], [1, 2], [1, 3], [2, 3]]:
    ax.plot3D(pts[ij, 0], pts[ij, 1], pts[ij, 2])

更多示例代码请参见我之前的回答 here、here 和 here。