This function generates daily climatological maps.
daily_clima(b, dates = NULL, verbose = FALSE)RasterBrick/Stack containing the historical observations or a proxy
(typically a reanalysis) that is used to derive the climatological
information. This needs to contain daily layers for multiple years.
b@z should contain dates for comparison
(e.g. 2017.01.01).
Dates for which we need to calculate daily climatology. By default, this is a leap year.
Logical defining whether message should be printed on screen (FALSE by default).
The function returns a RasterBrick (if dates contains one date)
or a list of RasterBricks (if dates contains more than one date).
Extent, resolution and land-sea mask match those of b.
Values are the subset of b related to the given dates.
the climatology is calculated over the period spanning 4 days before and
if (FALSE) {
# Generate dummy RasterBrick
set.seed(0)
r <- raster(nrows = 2, ncols = 2,
xmn = 0, xmx = 360, ymn = -90, ymx = 90, vals = 30)
b <- raster::brick(lapply(1:(365 * 3),
function(i) raster::setValues(r, runif(n = raster::ncell(r),
min = 0, max = 100))))
raster::setZ(b) <- seq.Date(from = as.Date("1993-01-01"),
to = as.Date("1995-12-31"),
by = "day")
daily_clima(b, as.Date("1996-01-01"))
}