steadiness: steadiness

View source: R/02_steadiness.R

steadinessR Documentation

steadiness

Description

steadiness derives the Steadiness Index from a land productivity variable

Usage

steadiness(obj2process = NULL, cores2use = 1, filename = "")

Arguments

obj2process

SpatRaster object (or its file name). If time series, each layer is one year

cores2use

Numeric. Number of cores to use for parallelization. Optional. Default is 1 (no parallelization)

filename

Character. Output filename. Optional

Details

The Steadiness Index is based on the combination of two metrics calculated per pixel: (1) the slope derived from a linear regression of the different years of the time series and (2) the net change on the same period. It results in a 4-class SpatRaster object ranging from (1) strong negative to (4) strong positive ecosystem dynamics. See Ivits et al. (2013) for further explanations.

Values = 0 in the final map indicates that there is a scarcity of data in the productivity variable (i.e. only 1 year with data), so that the indicator cannot be calculated

Value

SpatRaster object

Author(s)

Xavier Rotllan-Puig

References

Ivits, E., M. Cherlet, W. Mehl, and S. Sommer. 2013. “Ecosystem Functional Units Characterized by Satellite Observed Phenology and Productivity Gradients: A Case Study for Europe.” Ecological Indicators 27: 17–28. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1016/j.ecolind.2012.11.010")}

Examples


sb <- terra::rast(paste0(system.file(package='LPDynR'), "/extdata/sb_cat.tif"))
steadiness(obj2process = sb)



LPDynR documentation built on Sept. 13, 2024, 1:10 a.m.