三角形三维填充等高线图

非结构化三角形网格的填充等高线图。

使用的数据与trisurf3d_demo2的第二个图中的数据相同。 tricontour3d_demo显示此示例的未填充版本。

三角形三维填充等高线图示例

  1. # This import registers the 3D projection, but is otherwise unused.
  2. from mpl_toolkits.mplot3d import Axes3D  # noqa: F401 unused import
  4. import matplotlib.pyplot as plt
  5. import matplotlib.tri as tri
  6. import numpy as np
  8. # First create the x, y, z coordinates of the points.
  9. n_angles = 48
  10. n_radii = 8
  11. min_radius = 0.25
  13. # Create the mesh in polar coordinates and compute x, y, z.
  14. radii = np.linspace(min_radius, 0.95, n_radii)
  15. angles = np.linspace(0, 2*np.pi, n_angles, endpoint=False)
  16. angles = np.repeat(angles[..., np.newaxis], n_radii, axis=1)
  17. angles[:, 1::2] += np.pi/n_angles
  19. x = (radii*np.cos(angles)).flatten()
  20. y = (radii*np.sin(angles)).flatten()
  21. z = (np.cos(radii)*np.cos(3*angles)).flatten()
  23. # Create a custom triangulation.
  24. triang = tri.Triangulation(x, y)
  26. # Mask off unwanted triangles.
  27. triang.set_mask(np.hypot(x[triang.triangles].mean(axis=1),
  28.                          y[triang.triangles].mean(axis=1))
  29.                 < min_radius)
  31. fig = plt.figure()
  32. ax = fig.gca(projection='3d')
  33. ax.tricontourf(triang, z, cmap=plt.cm.CMRmap)
  35. # Customize the view angle so it's easier to understand the plot.
  36. ax.view_init(elev=45.)
  38. plt.show()

下载这个示例