readDWD.multia | R Documentation |
read multi_annual dwd data.
Intended to be called via readDWD()
.
All other observational data at dwdbase
can be read
with readDWD.data()
, except for the multi_annual and
subdaily/standard_format data.
readDWD.multia(
file,
fileEncoding = "latin1",
tryenc = TRUE,
comment.char = "\032",
quiet = rdwdquiet(),
...
)
file |
Name of file on harddrive, like e.g. DWDdata/multi_annual_mean_81-10_Temperatur_1981-2010_aktStandort.txt or DWDdata/multi_annual_mean_81-10_Temperatur_1981-2010_Stationsliste_aktStandort.txt |
fileEncoding |
|
tryenc |
Logical. If reading fails, try with encoding ="". This was added in version 1.8.5 (2023-09-27) because of non-reproducible issues in runLocalTests. DEFAULT: TRUE |
comment.char |
|
quiet |
Ignored.
DEFAULT: FALSE through |
... |
Further arguments passed to |
data.frame
Berry Boessenkool, berry-b@gmx.de, Feb 2019
readDWD()
## Not run: # Excluded from CRAN checks, but run in localtests
# Temperature aggregates (2019-04 the 9th file, 2022-05 the 8th):
durl <- selectDWD(res="multi_annual", per="mean_81-10")[8]
murl <- selectDWD(res="multi_annual", per="mean_81-10", meta=TRUE)[8]
# encoding issue not tested enough to be in the source code:
ma_temp <- dataDWD(durl)
ma_meta <- dataDWD(murl)
head(ma_temp)
head(ma_meta)
ma <- merge(ma_meta, ma_temp, all=TRUE)
berryFunctions::linReg(ma$Stationshoehe, ma$Jahr, main="annual average ~ elevation")
op <- par(mfrow=c(3,4), mar=c(0.1,2,2,0), mgp=c(3,0.6,0))
for(m in colnames(ma)[8:19])
{
berryFunctions::linReg(ma$Stationshoehe, ma[,m], xaxt="n", xlab="", ylab="", main=m)
abline(h=0)
}
par(op)
par(bg=8)
berryFunctions::colPoints(ma$geogr..Laenge, ma$geogr..Breite, ma$Jahr, add=F, asp=1.4)
DEU <- terra::vect(system.file("extdata/DEU.gpkg", package="rdwd"))
pdf("MultiAnn.pdf", width=8, height=10)
par(bg="grey90")
for(m in colnames(ma)[8:19])
{
terra::plot(DEU, border="grey40")
berryFunctions::colPoints(ma[-262,]$geogr..Laenge, ma[-262,]$geogr..Breite, ma[-262,m],
asp=1.4, # Range=range(ma[-262,8:19]),
col=berryFunctions::divPal(200, rev=TRUE), zlab=m, add=T)
}
dev.off()
berryFunctions::openFile("MultiAnn.pdf")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.