三角形三维等高线图

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

使用的数据与trisurf3d_demo2的第二个图中的数据相同。 tricontourf3d_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. n_angles = 48
  9. n_radii = 8
  10. min_radius = 0.25
  12. # Create the mesh in polar coordinates and compute x, y, z.
  13. radii = np.linspace(min_radius, 0.95, n_radii)
  14. angles = np.linspace(0, 2*np.pi, n_angles, endpoint=False)
  15. angles = np.repeat(angles[..., np.newaxis], n_radii, axis=1)
  16. angles[:, 1::2] += np.pi/n_angles
  18. x = (radii*np.cos(angles)).flatten()
  19. y = (radii*np.sin(angles)).flatten()
  20. z = (np.cos(radii)*np.cos(3*angles)).flatten()
  22. # Create a custom triangulation.
  23. triang = tri.Triangulation(x, y)
  25. # Mask off unwanted triangles.
  26. triang.set_mask(np.hypot(x[triang.triangles].mean(axis=1),
  27.                          y[triang.triangles].mean(axis=1))
  28.                 < min_radius)
  30. fig = plt.figure()
  31. ax = fig.gca(projection='3d')
  32. ax.tricontour(triang, z, cmap=plt.cm.CMRmap)
  34. # Customize the view angle so it's easier to understand the plot.
  35. ax.view_init(elev=45.)
  37. plt.show()

下载这个示例