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,
  dividebyten = TRUE,
  quiet = rdwdquiet(),
  ...
)

Value

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

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(), see readDWD.raster(). DEFAULT: NULL
toraster: Logical: convert output (list of matrixes + meta informations) to a list with data (terra::rast) + meta (list from the first subfile, but with vector of dates)? DEFAULT: TRUE
dividebyten: Logical: Divide the numerical values by 10? See readDWD. toraster??? DEFAULT: TRUE
quiet: Ignored. DEFAULT: FALSE through rdwdquiet()
...: Further arguments passed to dwdradar::readRadarFile(), i.e. na and clutter

Author

Berry Boessenkool, berry-b@gmx.de, Aug 2019. Significant input for the underlying dwdradar::readRadarFile() came from Henning Rust & Christoph Ritschel at FU Berlin.

Examples

Run this code

if (FALSE)  # 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)
plotRadar(r$dat, main=paste("mm in 24 hours preceding", r$meta$date))