Closed markwang0 closed 5 months ago
Thanks Mark, this does look strange. Will look into this as soon as I get a chance. Note that there have been some updates to the depression filling algorithm on the main branch of the github repo (#243) that may not yet be reflected on the distributed packages. It may be worth cloning from github directly and seeing if this changes the results.
Thank you @mdbartos. Here is the result of running the above code with a pysheds installation from GitHub:
@markwang0 I think the issue here might be filling depressions. Is TauDEM filling depressions? The depression filling step fills in a large portion of the DEM. On my machine, if I skip the filling depressions step, pysheds' routing produces better looking output.
Pysheds D8 with just resolve flats Pysheds dinf with just resolving flats
That makes sense @groutr -- when I skip fill_depressions and just do resolve_flats, I get results that look like yours. I did not fill depressions when running TauDEM
Describe the bug
D8 and Dinf flow directions calculated with pysheds do not appear as expected when run on a particular DEM. The flow direction rasters have diagonal streaks of values.
To Reproduce
Download the DEM that causes the problem:
Use pysheds to calculate flow directions on the DEM:
These are the resulting rasters:
Expected behavior
I would expect the results to look like these flow direction rasters calculated with TauDEM:
Environment
M1 MacBook Pro Python 3.12.3 arm64