不规则空间数据的等高线图

不规则空间数据在规则网格上插值的等高线图与非结构三角形网格的三棱图的比较。

由于 contour 和 contourf 期望数据存在于规则网格上，因此绘制不规则间隔数据的等高线图需要不同的方法。这两个选项是：

  • 首先将数据插值到常规网格。这可以通过机载装置完成，例如，通过LinearTriInterpolator或使用外部功能，例如 通过scipy.interpolate.griddata。然后用常规的等高线绘制插值数据。
  • 直接使用tricontour或tricontourf，它将在内部进行三角测量。

此示例显示了两种方法。

  1. import matplotlib.pyplot as plt
  2. import matplotlib.tri as tri
  3. import numpy as np
  5. np.random.seed(19680801)
  6. npts = 200
  7. ngridx = 100
  8. ngridy = 200
  9. x = np.random.uniform(-2, 2, npts)
  10. y = np.random.uniform(-2, 2, npts)
  11. z = x * np.exp(-x**2 - y**2)
  13. fig, (ax1, ax2) = plt.subplots(nrows=2)
  15. # -----------------------
  16. # Interpolation on a grid
  17. # -----------------------
  18. # A contour plot of irregularly spaced data coordinates
  19. # via interpolation on a grid.
  21. # Create grid values first.
  22. xi = np.linspace(-2.1, 2.1, ngridx)
  23. yi = np.linspace(-2.1, 2.1, ngridy)
  25. # Perform linear interpolation of the data (x,y)
  26. # on a grid defined by (xi,yi)
  27. triang = tri.Triangulation(x, y)
  28. interpolator = tri.LinearTriInterpolator(triang, z)
  29. Xi, Yi = np.meshgrid(xi, yi)
  30. zi = interpolator(Xi, Yi)
  32. # Note that scipy.interpolate provides means to interpolate data on a grid
  33. # as well. The following would be an alternative to the four lines above:
  34. #from scipy.interpolate import griddata
  35. #zi = griddata((x, y), z, (xi[None,:], yi[:,None]), method='linear')
  38. ax1.contour(xi, yi, zi, levels=14, linewidths=0.5, colors='k')
  39. cntr1 = ax1.contourf(xi, yi, zi, levels=14, cmap="RdBu_r")
  41. fig.colorbar(cntr1, ax=ax1)
  42. ax1.plot(x, y, 'ko', ms=3)
  43. ax1.axis((-2, 2, -2, 2))
  44. ax1.set_title('grid and contour (%d points, %d grid points)' %
  45.               (npts, ngridx * ngridy))
  48. # ----------
  49. # Tricontour
  50. # ----------
  51. # Directly supply the unordered, irregularly spaced coordinates
  52. # to tricontour.
  54. ax2.tricontour(x, y, z, levels=14, linewidths=0.5, colors='k')
  55. cntr2 = ax2.tricontourf(x, y, z, levels=14, cmap="RdBu_r")
  57. fig.colorbar(cntr2, ax=ax2)
  58. ax2.plot(x, y, 'ko', ms=3)
  59. ax2.axis((-2, 2, -2, 2))
  60. ax2.set_title('tricontour (%d points)' % npts)
  62. plt.subplots_adjust(hspace=0.5)
  63. plt.show()

不规则空间数据的等高线图例

参考

此示例中显示了以下函数和方法的用法：

  1. import matplotlib
  2. matplotlib.axes.Axes.contour
  3. matplotlib.pyplot.contour
  4. matplotlib.axes.Axes.contourf
  5. matplotlib.pyplot.contourf
  6. matplotlib.axes.Axes.tricontour
  7. matplotlib.pyplot.tricontour
  8. matplotlib.axes.Axes.tricontourf
  9. matplotlib.pyplot.tricontourf

下载这个示例