python - 使用 scipy.spatial.Delaunay 代替 matplotlib.tri.Triangulation 的内置版本

Question

似乎 matplotlib.tri.Triangulation 使用了一个有缺陷且可能不正确的 Delaunay 三角剖分实现，该三角剖分将被 qHull 取代.

我正在尝试使用 mpl_toolkits.mplot3d.plot_trisurf() 绘制 trisurf 并遇到一堆无用的异常(IndexError 和 主要是 KeyError，没有指出到底出了什么问题)。

因为 scipy.spatial.Delaunay 已经使用了 qHull，我想知道是否有办法构建一个 matplotlib.tri.Triangulation 对象以与 一起使用mpl_toolkits.mplot3d.plot_trisurf() 使用 scipy 的 Delaunay 三角剖分实现。

我尝试通过 triangles 参数将 delaunay.points 直接传递给 matplotlib.tri.Triangulate，但这会导致ValueError:三角形最小元素超出范围。

Answer 1

http://docs.scipy.org/doc/scipy-0.13.0/reference/generated/scipy.spatial.Delaunay.html http://matplotlib.org/dev/api/tri_api.html

因此您需要将点和三角形从 qhull 传递到 Triangulation 构造函数:

import numpy as np
import scipy.spatial
import matplotlib
import math

import matplotlib.tri as mtri
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

# First create the x and y coordinates of the points.
n_angles = 20
n_radii = 10
min_radius = 0.15
radii = np.linspace(min_radius, 0.95, n_radii)
angles = np.linspace(0, 2*math.pi, n_angles, endpoint=False)
angles = np.repeat(angles[..., np.newaxis], n_radii, axis=1)
angles[:, 1::2] += math.pi/n_angles
x = (radii*np.cos(angles)).flatten()
y = (radii*np.sin(angles)).flatten()

# Create the Delaunay tessalation using scipy.spatial
pts = np.vstack([x, y]).T
tess = scipy.spatial.Delaunay(pts)

# Create the matplotlib Triangulation object
x = tess.points[:, 0]
y = tess.points[:, 1]
tri = tess.vertices # or tess.simplices depending on scipy version
triang = mtri.Triangulation(x=pts[:, 0], y=pts[:, 1], triangles=tri)

# Plotting
z = x*x + y*y
fig = plt.figure()
ax = fig.gca(projection='3d')
ax.plot_trisurf(triang, z)
plt.show()

输出(使用 matplotlib current master):