calc_bioclim_raster: Bioclimatic index calculation raster

In sitscholl/rebecka_package: Utility functions used in the REBECKA project

Description

This function allows you to calculate bioclimatic indices based on the overlay of raster datasets.

Usage

calc_bioclim_raster(
  l,
  index.calc = "winkler",
  frml = NULL,
  agg.func = NULL,
  period = NULL,
  set.zero = NULL,
  project = T,
  crs = "+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0",
  res
)

Arguments

l	list, a list of RasterStacks or RasterBricks. The elements of the list have to be named according to the elements in 'frml' and must have the same names() attribute that corresponds to the day the raster values refer to.
index.calc	string, Name of the index that you want to calculate. See function get_index_params for further details.
frml	formula, A formula object that describes how to calculate the daily index values. The left side of the formula will be used to name the resulting RasterLayer and the right side to calculate the daily rasters.
agg.func	function, A function that will be used to aggregate the daily rasters calculated by 'frml' for every year
period	integer, A sequence of numbers that indicates the day of years that should be used to aggregate the index
set.zero	logical, Are values below 0 possible or should they be converted to 0?
project	logical, Should the resulting raster files be projected to another crs? Defaults to TRUE.
crs	character or object of class 'CRS', output coordinate system reference. Only used if project is TRUE. Defaults to WGS84
res	single or vector of two numerics, output resolution. Only used if project is TRUE. Defaults to the resolution of the input rasters.

Examples

calc_bioclim_raster(l = list('T_max' = tx.stack, 'T_max' = tn.stack),
  frml = winkler ~ ((T_max + T_min) / 2) - 10,
  agg.func = sum,
  period = 91:304,
  set.zero = T)

sitscholl/rebecka_package documentation built on Aug. 25, 2020, 4:20 a.m.