import topotoolbox
import matplotlib.pyplot as plt
dem = topotoolbox.load_dem('bigtujunga')
flats, sills = dem.identifyflats()
_= flats.plot(cmap='terrain')
plt.show()