import topotoolbox
import matplotlib.pyplot as plt
dem = topotoolbox.load_dem('perfectworld')
excess = dem.excesstopography(threshold=0.3, method='fsm2d')
_= excess.plot(cmap='terrain')
plt.show()