flood3 {hydflood} | R Documentation |
Function to compute flood extent or flood duration SpatRaster
along the German federal waterways Elbe and Rhine using the 1d water level
algorithms hyd1d::waterLevel()
and hyd1d::waterLevelPegelonline()
Description
Computes flood extent, if length(seq)
equals 1, or flood
duration for the active floodplains along the German federal waterways Elbe
and Rhine based on 1d water levels computed by
waterLevel
or
waterLevelPegelonline
provided by package hyd1d.
Usage
flood3(x, seq, filename = "", ...)
Arguments
x |
has to by type |
seq |
has to be type |
filename |
supplies an optional output filename and has to be type
|
... |
additional arguments as for |
Details
For every time step provided in seq
, flood3()
computes
a 1d water level along the requested river section. This 1d water level is
transfered to a wl
(water level) raster layer, which is in fact a
copy of the csa
(cross section areas) layer, and then
compared to the dem
(digital elevation model) layer. Where the
wl
layer is higher than the dem
, layer flood duration is
increased by 1.
Value
SpatRaster
object with flood duration in the range of
[0, length(seq)]
.
References
Weber A (2022). “Flood durations and potential natural vegetation distribution on the floodplains of River Rhine and River Elbe, Germany.” doi:10.1594/PANGAEA.948042.
Weber A (2023). “Flood durations of the year 2022 on the floodplains of River Rhine and River Elbe, Germany.” doi:10.1594/PANGAEA.961117.
See Also
waterLevel
,
waterLevelPegelonline
,
writeRaster
,
terraOptions
Examples
options("hydflood.datadir" = tempdir())
library(hydflood)
# import the raster data and create a raster stack
c <- st_crs("EPSG:25833")
e <- ext(309000, 310000, 5749000, 5750000)
x <- hydSpatRaster(ext = e, crs = c)
# create a temporal sequence
seq <- seq(as.Date("2016-12-01"), as.Date("2016-12-31"), by = "day")
# compute a flood duration
fd <- flood3(x = x, seq = seq)