Any ideas why the catchment delineated has holes? How can I filled the holes?

[staceyzhang-cbcl]

I was testing the functions with two DEM files. Both catchments I delineated with a randomly picked pour point had holes in them where I dont think there should be. How do I eliminate them? Did I do something wrong?

For the first DEM, I filled pits, depressions and resolved flats. The first two was good, with assert not depressions.any()/pits.any() returning nothing. But the flats seemed not resolved. Could it be the reason?

Thank you!

[mdbartos]

Greetings. When running resolve_flats, try setting a lower epsilon value, e.g.:

resolve_flats(..., eps=1e-7)

[staceyzhang-cbcl]

Thanks so much for your quick response!

This solved one of my two DEM problems. The second DEM still has one hole (but not two, I take that as a partial success). I wonder how low the epsilon value can be. I tried eps=1e-11 but it still didn't resolve all flats.

[mdbartos]

@Eric-Jenkins : Could be a good test case to look at.

@staceyzhang-cbcl : Are you able to send the DEM?

[staceyzhang-cbcl]

I tried. But the file is too large to upload. It's 125MB after compression.

[staceyzhang-cbcl]

OK, I uploaded it to my google drive.

The pour point is 576410.1845, 5048888.3975.