Frontpage 3D示例

此示例再现Frontpage 3D示例。

Frontpage 3D示例

  1. # This import registers the 3D projection, but is otherwise unused.
  2. from mpl_toolkits.mplot3d import Axes3D  # noqa: F401 unused import
  4. from matplotlib import cbook
  5. from matplotlib import cm
  6. from matplotlib.colors import LightSource
  7. import matplotlib.pyplot as plt
  8. import numpy as np
  10. filename = cbook.get_sample_data('jacksboro_fault_dem.npz', asfileobj=False)
  11. with np.load(filename) as dem:
  12.     z = dem['elevation']
  13.     nrows, ncols = z.shape
  14.     x = np.linspace(dem['xmin'], dem['xmax'], ncols)
  15.     y = np.linspace(dem['ymin'], dem['ymax'], nrows)
  16.     x, y = np.meshgrid(x, y)
  18. region = np.s_[5:50, 5:50]
  19. x, y, z = x[region], y[region], z[region]
  21. fig, ax = plt.subplots(subplot_kw=dict(projection='3d'))
  23. ls = LightSource(270, 45)
  24. # To use a custom hillshading mode, override the built-in shading and pass
  25. # in the rgb colors of the shaded surface calculated from "shade".
  26. rgb = ls.shade(z, cmap=cm.gist_earth, vert_exag=0.1, blend_mode='soft')
  27. surf = ax.plot_surface(x, y, z, rstride=1, cstride=1, facecolors=rgb,
  28.                        linewidth=0, antialiased=False, shade=False)
  29. ax.set_xticks([])
  30. ax.set_yticks([])
  31. ax.set_zticks([])
  32. fig.savefig("surface3d_frontpage.png", dpi=25)  # results in 160x120 px image

下载这个示例