Calculate waterdepths for 3Di results. For results of type 'raw' the variable 's1' is used as waterlevel. For results of type 'aggregate', the variable 's1_max' is used as waterlevel.
Calculate concentrations for 3Di water quality results.
- Interpolated or gridcell-constant waterlevels or concentrations
- Interfaces with threediresults via threedigrid
- Progress indicator support
- Low memory consumption
For the interpolated mode, the 'lizard'-method is used. For a detailed description, read the docstring for the LizardLevelCalculator.
For the maximum waterlevel calculation, the maximum waterlevel for each point is taken before the interpolation is applied. This can lead to situations where the highest waterlevel for a pixel for a certain timestep is higher than the maximum waterlevel for the pixel.
Make sure GDAL is available as (from osgeo import gdal)
$ pip install threedidepth
From the cli:
$ threedidepth gridadmin.h5 results_3di.nc dem.tif waterdepth.tif $ threediwq gridadmin.h5 water_quality_results_3di.nc substance1 <xmin ymin xmax ymax> concentration.tif
Or python:
>>> threedidepth.calculate_waterdepth(...) >>> threedidepth.calculate_water_quality(...)
For development, clone the repository and use a docker compose setup:
$ docker compose build --build-arg uid=`id -u` --build-arg gid=`id -g` lib $ docker compose up --no-start $ docker compose start $ docker compose exec lib bash
Create a virtualenv:
# note that Dockerfile prepends .venv/bin to $PATH (docker)$ python3 -m venv .venv --system-site-packages
Install dependencies & package and run tests:
(docker)$ pip install -r requirements.txt (docker)$ pip install -e .[test] (docker)$ pytest
Update packages:
(docker)$ rm -rf .venv (docker)$ python3 -m venv .venv --system-site-packages (docker)$ pip install -e . (docker)$ pip freeze | grep -v threedidepth > requirements.txt