csdms-contrib / fwdet

Calculate floodwater depth based on an inundation polygon (e.g. from remote sensing) and a DEM. See: http://csdms.colorado.edu/wiki/Model:FwDET
GNU General Public License v2.0
56 stars 28 forks source link

Water Depth issue #19

Closed jeremyEudaric closed 3 months ago

jeremyEudaric commented 4 months ago

Dear Team thank you for your help and time , its finally working (l hope) l resized my DEM thanks to QGIS (Warp re projected) and l used a DEM of Tiff (32-float type) from Copernicus. However l have an issue with the results l have a water depth between 0 and 0.12 that doesn´t make sens. Could you give me you point of view if you have time.

results

Also looks like fwdet does not take all the inundation map as you can see on this screen. Do you have an idea as well how l could fix this issue ? :)

C

Thank you very much for your kindness and time

cefect commented 4 months ago

This can happen if there is a poor match between your DEM and the polygon. I suggest starting with a smaller area to try and understand the issue. If this doesn't work, please upload the data and I can take a look.

Regards

jeremyEudaric commented 4 months ago

Thank you for your reply and time unfortunate l try with a smaller area and l still have the same issue fwdet does not take all the inundation map and l have a water depth between 0 and 0.12 that doesn´t make sens. Do you have an idea why ? (l am using a DEM from Copernicus 90 m)

You will fin d attached the link to download the data.

https://drive.google.com/file/d/1hE4NjL_sDL8gQ-zibVEm3iGU9c1MsPXg/view?usp=drive_link

Thank you very much for your time

jeremyEudaric commented 4 months ago

Screenshot from 2024-05-17 12-27-08 l tried on a small portion but its still the same issue l guess the issue comes from my polygon. Are they a specific rules for the polygon ( which format to save it etc ) ?

l have this error l think its the main issue

"Unable to execute algorithm Could not load source layer for INPUT_B: /tmp/processing_rHdWLB/973a830eb6a14868965b64d13f5e6623/OUTPUT.tif not found Unable to execute algorithm Could not load source layer for INPUT_B: /tmp/processing_rHdWLB/973a830eb6a14868965b64d13f5e6623/OUTPUT.tif not found Execution failed after 0.32 seconds"

cefect commented 3 months ago

Please provide the log file also to help me debug

jeremyEudaric commented 3 months ago

Thank you for our reply and time you will find attached the log files

https://drive.google.com/drive/folders/18km7boOaIIRpoDeJ7S1IlfVWItbRtolq?usp=drive_link

cefect commented 3 months ago

your DEM is not really a DEM. fwdet expects a single band raster with elevation values (typically in meters). As yours only ranges from 0-1, so too do your water depths only range from 0-1. I'll update the documentation to make this more clear.

I'm not sure where you obtained the DEM from, but perhaps there is some scaling value you could apply to convert it to meters. Alternatively, I've had luck using MERIT DEM or some other global DEM (e.g., FABDEM). However, I suspect there is something more accurate/detailed available for Greece.

image

jeremyEudaric commented 3 months ago

Thank you for your reply l downloaded the DEM from Copernicus. l Tried with a new DEM from https://www.usgs.gov/ with a single band (l added the file to the link its called new DEM) but l still have the issue with the poor match as explained above or sometime l have no results at all like this: pppp

l will try with MERIT DEM l will let you know soon as possible , thanks a lot

cefect commented 3 months ago

Send me the new inputs and output. I can take a look. Maybe try forcing the minimum and maximum on the layer styles of your results layer, sometimes if there are very few real values QGiS will fail to automatically apply a style.

jeremyEudaric commented 3 months ago

https://drive.google.com/drive/folders/18km7boOaIIRpoDeJ7S1IlfVWItbRtolq?usp=drive_link The new polygon and DEM, thanks a lot (DEM single band)

cefect commented 3 months ago

DEM_new.tif has had zeros burned from the inundation polygon. I'm not sure why you've done this, but you certainly don't want to be burning out these regions. Run with the un-burned DEM and you should get something reasonable.

image