Delaunay细分生成单形：二维的三角形三维的四面体所以这里你得到了四面体。您可以提取它们的面（三角形），然后继续使用ax.plot_trisurf。

import matplotlib.pyplot as plt
from scipy.spatial import Delaunay
import numpy as np

# extract the four triangular faces of a tetrahedron
def simplex2triangles(sx):
    return np.array([
      [sx[0], sx[1], sx[2]],  
      [sx[0], sx[1], sx[3]],  
      [sx[0], sx[2], sx[3]],  
      [sx[1], sx[2], sx[3]]  
    ])

# cube vertices
points= np.array([
    [0,0,0],
    [4,0,0],
    [4,4,0],
    [0,4,0],
    [0,0,4],
    [4,0,4],
    [4,4,4],
    [0,4,4]
])

# apply Delaunay tetrahedralization
tri = Delaunay(points)
# extract tetrahedra
tetrahedra = tri.simplices
# extract all triangles
triangles = np.vstack(np.apply_along_axis(simplex2triangles, 1, tri.simplices))
# there are duplicated triangles, remove them (remark: it could remain some dups up to the order...)
triangles = np.unique(triangles, axis = 0)

# plot
x = points[:, 0]
y = points[:, 1]
z = points[:, 2]
ax = plt.axes(projection='3d')
ax.view_init(15, -30)
ax.plot_trisurf(x, y, z, triangles = triangles, 
    linewidth=0.5, linewidths=1, edgecolors='r', alpha=.70, antialiased=True)

我不知道为什么facecolors不工作。

还要注意我在Python中相当新手。