map_grass_pathway: Map C4 and C3 proportion of grasses

In rebeccalpowell/grassmapr: Predict C3/C4 Grass Distribution

Description

Wrapper function to map proportion of grasses with C4 and C3 photosynthetic pathway - using climate layers and crossover temperature.

Usage

map_grass_pathway(temp.stack, precip.stack, C4.threshold, GS.threshold,
  precip.threshold, filename = "", ...)

Arguments

temp.stack	Raster* object. Each layer corresponds to a different temporal window (e.g., month) for a temperature climate variable. Object may be single or multi-layer.
precip.stack	Raster* object. Each layer corresponds to a different temporal window (e.g., month) for a precipitation climate variable. Object may be single or multi-layer.
C4.threshold	Numeric. Threshold value (lower-bound inclusive) for C4 climate temperature.
GS.threshold	Numeric. Threshold value (lower-bound inclusive) for Growing Season temperature.
precip.threshold	Numeric. Threshold value (lower-bound inclusive) for Growing Season precipitation.
filename	Character. Optional output root filename passed to writeRaster, default output file type is GeoTiff. If not specified, output is written to a temporary file.
...	Other arguments passed to writeRaster.

Value

RasterBrick. Layer[[1]] corresponds to percent of grasses with C4 photosynthetic pathway, Layer[[2]] corresponods to percent of grasses with C3 photosynthetic pathway.

See Also

mask_climate, calc_C4_ratio, calc_pft_cover.

Examples

# Load data for example
data(temp_NA)  # mean monthly temperature (deg. C)
data(prec_NA)  # mean monthly precipitation (mm)

# Specify climate thresholds
C4_temp <- 22    # mean monthly temperature >= 22 deg. C
GS_temp <- 5     # mean monthly temperature >= 5 deg. C
GS_prec <- 25    # mean monthly precipitation >= 25 mm

# Predict percent of grasses by photosynthetic pathway.
grass_map <- map_grass_pathway(temp.stack = temp_NA,
                               precip.stack = prec_NA,
                               C4.threshold = C4_temp,
                               GS.threshold = GS_temp,
                               precip.threshold = GS_prec)

# Plot C4 and C3 grass layers.
plot(grass_map)

rebeccalpowell/grassmapr documentation built on June 20, 2019, 8:07 p.m.