readRadarFile: read binary radolan radar file

Description

Read a single binary DWD Radolan file. To be used in rdwd, especially for proper [1/10 mm] unit correction in rdwd::readDWD.
If any files ar not read correctly, please let me know, referencing the Kompositformatbeschreibung at https://www.dwd.de/DE/leistungen/radolan/radolan.html.
The meta-info is extracted with readHeader (not exported, but documented).
Binary bits are converted to decimal numbers with Fortran routines, see https://github.com/brry/dwdradar/tree/master/src. They are called via bin2num (not exported, but documented).

Usage

readRadarFile(binfile, na = NA, clutter = NA)

Value

Invisible list with dat (matrix) and meta

(list with elements from header, see Kompositformatbeschreibung).

Arguments

binfile: Name of a single binary file
na: Value to be set for missing data (bit 14). DEFAULT: NA
clutter: Value to be set for clutter data (bit 16). DEFAULT: NA

Author

Maintained by Berry Boessenkool, berry-b@gmx.de, May + Oct 2019.
Original codebase by Henning Rust & Christoph Ritschel at FU Berlin

Examples

Run this code


f <- system.file("extdata/raa01_sf_2019-10-14_1950", package="dwdradar")
out <- readRadarFile(f)
out$meta

if(requireNamespace("terra", quietly=TRUE))
  terra::plot(terra::rast(out$dat))

# for more files, see the tests.
# for real-world usage, readDWD.binary / readDWD.radar in the rdwd package