flood2 {hydflood} | R Documentation |
SpatRaster
along the German federal waterways Elbe and Rhine using the 1d water level
algorithm hyd1d::waterLevelFlood2()
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
waterLevelFlood2
provided by package hyd1d in
analogy to the INFORM 3 module 'Flut2'.
flood2(x, seq, filename = "", ...)
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 |
For every time step provided in seq
, flood2()
computes
a 1d water level using waterLevelFlood2
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.
SpatRaster
object with flood duration in the range of
[0, length(seq)]
.
Rosenzweig S, Giebel H, Schleuter M (2011). “Ökologische Modellierungen für die Wasser- und Schifffahrtsverwaltung – Das integrierte Flussauenmodell INFORM in seiner neuesten Fassung (Version 3). Bundesanstalt für Gewässerkunde, Koblenz, Germany.” doi:10.5675/bfg-1667.
df.gauging_data
,
getGaugingDataW
,
waterLevelFlood2
,
writeRaster
,
terraOptions
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 <- flood2(x = x, seq = seq)