Map floods with Sentinel-1 radar images. We replicate in this package the work of Bauer-Marschallinger et al. (2022)1 on the TU Wien Bayesian-based flood mapping algorithm. This implementation is entirely based on dask and data access via STAC with odc-stac. The algorithm requires three pre-processed input datasets stored and accessible via STAC at the Earth Observation Data Centre For Water Resources Monitoring (EODC). It is foreseen that future implementations can also use data from other STAC catalogues. This notebook explains how microwave backscattering can be used to map the extent of a flood. The workflow detailed in this notebook forms the backbone of this package. For a short overview of the Bayesian decision method for flood mapping see this ProjectPythia book.
To install the package, do the following:
pip install git+https://github.com/TUW-GEO/dask-flood-mapper
Storm Babet hit the Denmark and Northern coast of Germany at the 20th of October 2023 Wikipedia. Here an area around Zingst at the Baltic coast of Northern Germany is selected as the study area.
Define the time range and geographic region in which the event occurred.
time_range = "2022-10-11/2022-10-25"
bbox = [12.3, 54.3, 13.1, 54.6]Use the flood module and calculate the flood extent with the Bayesian decision method applied tp Sentinel-1 radar images. The object returned is a xarray with lazy loaded Dask arrays. To get the data in memory use the compute method on the returned object.
from dask_flood_mapper import flood
flood.decision(bbox=bbox, datetime=time_range).compute()It is also possible to remotely process the data at the EODC Dask Gateway with the added benefit that we can then process close to the data source without requiring rate-limiting file transfers over the internet.
For ease of usage of the Dask Gateway install the eodc package besides the dask-gateway package. Also, see the EODC documentation.
pip install dask-gateway eodc
# or use pipenv sync -dConnect to the gateway (this requires an EODC account).
from eodc.dask import EODCDaskGateway
from rich.prompt import Prompt
your_username = Prompt.ask(prompt="Enter your Username")
gateway = EODCDaskGateway(username=your_username)Create a cluster.
cluster = gateway.new_cluster(cluster_options)
client = cluster.get_client()
cluster.adapt(minimum=2, maximum=5)
Map the flood the same way as we have done when processing locally.
flood.decision(bbox=bbox, datetime=time_range).compute()Please find the contributing guielines in the specific file here.
Credits go to EODC (https://eodc.eu) for developing the infrastructure and the management of the data required for this workflow. This work has been supported as part of the interTwin project (https://www.intertwin.eu). The interTwin project is funded by the European Union Horizon Europe Programme - Grant Agreement number 101058386.
Views and opinions expressed are however those of the authors only and do not necessarily reflect those of the European Union Horizon Europe/Horizon 2020 Programmes. Neither the European Union nor the granting authorities can be held responsible for them.
- Bauer-Marschallinger, Bernhard, Senmao Cao, Mark Edwin Tupas, Florian Roth, Claudio Navacchi, Thomas Melzer, Vahid Freeman, and Wolfgang Wagner. “Satellite-Based Flood Mapping through Bayesian Inference from a Sentinel-1 SAR Datacube. Remote Sensing 14, no. 15 (January 2022): 3673. https://doi.org/10.3390/rs14153673.
This repository is covered under the MIT License.