rdwd (version 1.4.0)

projectRasterDWD: project DWD raster data

Description

Set projection and extent for DWD raster data. Optionally (and per default) also reprojects to latlon data. WARNING: reprojection to latlon changes values slightly. For the tested RX product, this change is significant, see: https://github.com/brry/rdwd/blob/master/misc/ExampleTests/Radartests.pdf In raster::plot, use zlim with the original range if needed.

Usage

projectRasterDWD(
  r,
  proj = "radolan",
  extent = "radolan",
  targetproj = "ll",
  quiet = rdwdquiet()
)

Arguments

r

Raster object

proj

Current projection to be given to r. Can be - a raster::crs() input (e.g. a projection character string), - NULL to not set proj+extent (but still consider targetproj), - or a special charstring for internal defaults, namely: "radolan" (readDWD.binary + .asc + .radar), "seasonal" (.raster) or "nc" (.nc). DEFAULT: "radolan"

extent

Current raster::extent() to be given to r. Ignored if proj=NULL. Can be an extent object, a vector with 4 numbers, or "radolan" / "rw" / "seasonal" / "nc" with internal defaults. DEFAULT: "radolan"

targetproj

r is reprojected to this raster::crs(). Use NULL to not reproject (i.e. only set proj and extent). DEFAULT: "ll" with internal default for lat-lon.

quiet

Logical: suppress progress messages? DEFAULT: FALSE through rdwdquiet()

Value

Raster object with projection and extent, invisible

Details

The internal defaults are extracted from the Kompositformatbeschreibung at https://www.dwd.de/DE/leistungen/radolan/radolan.html, as provided 2019-04 by Antonia Hengst. The nc extent was obtained by projecting Germanys bbox to EPSG 3034 (specified in the DWD documentation). Using that as a starting point, I then refined the extent to a visual match, see developmentNotes.R

See Also

raster::crs / extent / projectRaster readDWD.binary / raster / asc / radar / nc

Examples

# NOT RUN {
# To be used after readDWD.binary etc
# }