从Skimage高程图（2D numpy.ndarray）生成3D表面图

Question

In the skimage Segmentation tutorial , a 3D surface plot of the elevation map generated from the sobel function was plotted.

>>> from skimage.filters import sobel
>>> elevation_map = sobel(coins)

Question: elevation_map appears to be a 2D numpy.ndarray . How do we generate the 3D map shown using this?

Answer 1

This is likely produced using Paraview/VTK;

Try to play around the following:

from skimage import data
from skimage.filters import sobel
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
from matplotlib import cm
from scipy.ndimage import zoom
coins = data.coins()
coins = zoom(coins, 10)

elevation_map = sobel(coins)

fig = plt.figure()
ax = fig.add_subplot(111, projection='3d')

m, n=elevation_map.shape
X, Y = np.meshgrid(np.arange(n), np.arange(m))

ax.plot_surface(X, Y, elevation_map, cmap=cm.viridis, antialiased=False)
ax.axis("off")
ax.set_facecolor('black')
plt.show()