import topotoolbox
import matplotlib.pyplot as plt
dem = topotoolbox.load_dem('perfectworld')
eroded = dem.erode(size = (3,3))
_= eroded.plot(cmap='terrain')
plt.show()