rdwd (version 1.2.0)

readDWD.radar: read dwd gridded radolan radar data

Description

read gridded radolan radar data. Intended to be called via readDWD.

Usage

readDWD.radar(file, gargs = NULL, toraster = TRUE, ...)

Arguments

file

Name of file on harddrive, like e.g. DWDdata/hourly/radolan/recent/bin/ raa01-rw_10000-1802020250-dwd---bin.gz

gargs

Named list of arguments passed to R.utils::gunzip. The internal defaults are: remove=FALSE (recommended to keep this so file does not get deleted) and skip=TRUE (which reads previously unzipped files as is). If file has changed, you might want to use gargs=list(skip=FALSE, overwrite=TRUE) or alternatively gargs=list(temporary=TRUE). The gunzip default destname means that the unzipped file is stored at the same path as file. DEFAULT gargs: NULL

toraster

Logical: convert output (list of matrixes + meta informations) to a list with data (raster stack) + meta (list from the first subfile, but with vector of dates)? DEFAULT: TRUE

…

Further arguments passed to dwdradar::readRadarFile, i.e. na and clutter

Value

Invisible list with dat (matrix or raster, depending on toraster) and meta (list with elements from header)

See Also

readDWD, especially readDWD.binary https://wradlib.org for much more extensive radar analysis in Python Kompositformatbeschreibung at https://www.dwd.de/DE/leistungen/radolan/radolan.html for format description

Examples

# NOT RUN {
 # Excluded from CRAN checks, but run in localtests
# recent radar files 
rrf <- indexFTP("hourly/radolan/recent/bin", base=gridbase, dir=tempdir())
lrf <- dataDWD(rrf[773], base=gridbase, joinbf=TRUE, dir=tempdir(), read=FALSE)
r <- readDWD(lrf)

rp <- projectRasterDWD(r$dat)
raster::plot(rp, main=r$meta$date)
addBorders()
# }

Run the code above in your browser using DataCamp Workspace