Cannot replicate the example using 'elevation.tiff'

[jddratli]

I am just getting started with pysheds and I am trying to replicate the example that is posted, but I get different results for the flow direction, accumulation, catchment delineation and river network extraction. I have essentially copied and pasted the script into my own Jupyter notebook, so I am confused as to why. I have attached images of the results that I get. The code is the same as the example code on the ReadMe page. If someone can help me out, I would be greatly appreciative.

-John

[philippkraft]

Looks like you forgot to fill pits and depressions and resolve flats. If I leave them out I get the same picture as you.

# Condition DEM
# ----------------------
# Fill pits in DEM
pit_filled_dem = grid.fill_pits(dem)

# Fill depressions in DEM
flooded_dem = grid.fill_depressions(pit_filled_dem)
    
# Resolve flats in DEM
inflated_dem = grid.resolve_flats(flooded_dem)

[mdbartos]

Thanks @philippkraft!

[jddratli]

Okay, thanks! Sounds like that's probably it so I'll give it another try. John

…

________________________________ From: Philipp Kraft ***@***.***> Sent: Thursday, February 9, 2023 12:42 PM To: mdbartos/pysheds ***@***.***> Cc: Ratcliff, John David ***@***.***>; Author ***@***.***> Subject: Re: [mdbartos/pysheds] Cannot replicate the example using 'elevation.tiff' (Issue #211) You don't often get email from ***@***.*** Learn why this is important<https://aka.ms/LearnAboutSenderIdentification> Looks like you forgot to fill pits and depressions and resolve flats. If I leave them out I get the same picture as you. # Condition DEM # ---------------------- # Fill pits in DEM pit_filled_dem = grid.fill_pits(dem) # Fill depressions in DEM flooded_dem = grid.fill_depressions(pit_filled_dem) # Resolve flats in DEM inflated_dem = grid.resolve_flats(flooded_dem) — Reply to this email directly, view it on GitHub<#211 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ALLVP4BOIQ5ZSPFDR4FRHZLWWUT63ANCNFSM6AAAAAAUATPNTE>. You are receiving this because you authored the thread.Message ID: ***@***.***>