Closed rogerzzl closed 7 months ago
Thanks @rogerzzl, my guess is that if you are using D8, the flats did not resolve properly. Try the fix listed here by using a smaller epsilon value: https://github.com/mdbartos/pysheds/issues/224#issuecomment-1684378212
I tried to use small epsilon value like resolve_flats(flooded_dem, eps=1e-7), but the acc result is still the same as the previous output
oops, I tried to look back the codes in sgrid.py, the param max_iter might cause this issue, so once I tentatively added resolve_flats(flooded_dem,eps=1e-7, max_iter=5000), the right result is there!!
Thank you Dr. Matthew Bartos for your useful suggestions and efforts with the awesome tool Pysheds, I'll close this issue then.
All the best!
Firstly many thanks to the contributors for the pysheds tool. I met an issue about generating accumulation and didnot figure it out yet.
my codes:
there was a fracture in the main channel as shown in the red box in my accumulation figure (
While I used ArcPy to test that the channels at this red box position should be linked together, and this is also the reality
.
).