import topotoolbox
import matplotlib.pyplot as plt
from matplotlib.colors import ListedColormap
dem = topotoolbox.load_dem('bigtujunga')
prom, (x, y) = dem.prominence(tolerance=90)
_ = plt.subplot()
_ = dem.plot_hs(cmap=ListedColormap([0.9, 0.9, 0.9]), exaggerate=3)
_ = plt.scatter(x, y, prom, alpha=0.6, edgecolor='black')
plt.show()