import topotoolbox
import matplotlib.pyplot as plt
dem = topotoolbox.load_dem('perfectworld')
fd = topotoolbox.FlowObject(dem)
acc = fd.flow_accumulation()
_= acc.plot(cmap='Blues',norm="log")
plt.show()