apply.months {climetrics}

Apply a function on rasters of each month over years

Description

Apply a user specified function over raster time series for each month separately.

Usage

  apply.months(x,FUN,...)

Arguments

x

A Raster or Raster Time Series object

FUN

A function to apply on rasters of each month over different years; The defaults if mean

...

additional argument: if x is a Raster or SpatRaster object, corresponding dates should be specified in the dates argument

Details

The output contains 12 raster layers corresponding to the 12 months (Jan-Dec.). For each of the available 12 months (from January to December), the corresponding rasters are selected and the FUN is applied. For example if the FUN is mean, the first layer in the output would be the mean of the rasters in January over different years, and the last later (12th layer) would be the mean of the rasters in December over different years.

Value

A Raster object (RasterStackBrick or SpatRaster depending on the input) contains 12 layers correspond to 12 months (From Jan to Dec)

Author(s)

Shirin Taheri; Babak Naimi

[e-mail]; [e-mail]

Examples

filePath <- system.file("external/", package="climetrics") # path to the dataset folder


pr <- rast(paste0(filePath,'/precip.tif'))

n <- readRDS(paste0(filePath,'/dates.rds')) # read corresoinding dates

head(n) # Dates corresponds to the layers in climate variables (pr, tmin, tmax, tmean)

####################

# use rts function in the rts package to make a raster time series:

pr.t <- rts(pr,n) 


###########################

p12 <- apply.months(pr.t,'mean')


p12

plot(p12)