着色示例

显示如何制作阴影浮雕图的示例，如Mathematica （http://reference.wolfram.com/mathematica/ref/ReliefPlot.html）或通用映射工具（https://gmt.soest.hawaii.edu/）

  1. import numpy as np
  2. import matplotlib.pyplot as plt
  3. from matplotlib.colors import LightSource
  4. from matplotlib.cbook import get_sample_data
  7. def main():
  8.     # Test data
  9.     x, y = np.mgrid[-5:5:0.05, -5:5:0.05]
  10.     z = 5 * (np.sqrt(x**2 + y**2) + np.sin(x**2 + y**2))
  12.     filename = get_sample_data('jacksboro_fault_dem.npz', asfileobj=False)
  13.     with np.load(filename) as dem:
  14.         elev = dem['elevation']
  16.     fig = compare(z, plt.cm.copper)
  17.     fig.suptitle('HSV Blending Looks Best with Smooth Surfaces', y=0.95)
  19.     fig = compare(elev, plt.cm.gist_earth, ve=0.05)
  20.     fig.suptitle('Overlay Blending Looks Best with Rough Surfaces', y=0.95)
  22.     plt.show()
  25. def compare(z, cmap, ve=1):
  26.     # Create subplots and hide ticks
  27.     fig, axs = plt.subplots(ncols=2, nrows=2)
  28.     for ax in axs.flat:
  29.         ax.set(xticks=[], yticks=[])
  31.     # Illuminate the scene from the northwest
  32.     ls = LightSource(azdeg=315, altdeg=45)
  34.     axs[0, 0].imshow(z, cmap=cmap)
  35.     axs[0, 0].set(xlabel='Colormapped Data')
  37.     axs[0, 1].imshow(ls.hillshade(z, vert_exag=ve), cmap='gray')
  38.     axs[0, 1].set(xlabel='Illumination Intensity')
  40.     rgb = ls.shade(z, cmap=cmap, vert_exag=ve, blend_mode='hsv')
  41.     axs[1, 0].imshow(rgb)
  42.     axs[1, 0].set(xlabel='Blend Mode: "hsv" (default)')
  44.     rgb = ls.shade(z, cmap=cmap, vert_exag=ve, blend_mode='overlay')
  45.     axs[1, 1].imshow(rgb)
  46.     axs[1, 1].set(xlabel='Blend Mode: "overlay"')
  48.     return fig
  51. if __name__ == '__main__':
  52.     main()

着色示例

着色示例2

参考

本例中显示了以下函数、方法和类的使用：

  1. import matplotlib
  2. matplotlib.colors.LightSource
  3. matplotlib.axes.Axes.imshow
  4. matplotlib.pyplot.imshow

下载这个示例