Set projection and extent for DWD raster data. Optionally (and per default) also reprojects to latlon data. 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 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.
projectRasterDWD(r, proj = "radolan", extent = "radolan",
latlon = TRUE, quiet = FALSE)Raster object
Desired projection.
Use NULL to not set proj+extent but still consider latlon.
Can be a raster::crs output,
a projection character string (will be passed to crs),
or a special charstring for internal defaults, namely:
"radolan" (readDWD.binary + .asc + .radar), "seasonal" (.raster) or "nc" (.nc).
DEFAULT: "radolan"
Desired extent. Can be an extent object,
a vector with 4 numbers, or "radolan" / "rw" / "seasonal" / "nc"
with internal defaults.
DEFAULT: "radolan"
Logical: reproject r to lat-lon crs? DEFAULT: TRUE
Logical: suppress progress messages? DEFAULT: FALSE
Raster object with projection and extent, invisible
raster::crs,
raster::projection,
raster::extent,
raster::projectRaster,
readDWD.binary, readDWD.raster, readDWD.asc
# NOT RUN {
# To be used after readDWD.binary, readDWD.raster, readDWD.asc
# }
Run the code above in your browser using DataLab