ec-jrc / lisflood-utilities

LISFLOOD OS - Utilities
https://ec-jrc.github.io/lisflood
European Union Public License 1.2
22 stars 9 forks source link

Use stations.txt and ldd to produce catchment mask map #12

Closed domeniconappo closed 3 years ago

domeniconappo commented 4 years ago

In cutmaps tool, we can add another way to cut files, instead of directly providing maskmap or coordinates for cutting.

We may pass a stations.txt file and ldd of EFAS domain and follow the logic below.

https://github.com/ec-jrc/lisflood-calibration/blob/master/CAL_3_PREP_MISC.py#L75-L154

This is the script included in calibration tool: running the few pcraster commands from lines 75 to 154 gives you the maskmap and the outlets map (to produce tss files)

For doing that you could create a text file station.txt with coordinates in EFAS system (corrected , ask Juliana to put them in the center of cell), that's another thing, the station should be in our coordinates system (LEAE ).

Problem: currently the tool does not need PCRaster deps but now we are forced to. With latest PCRaster and conda environment is very easy to setup, though.