View source: R/02_steadiness.R
steadiness | R Documentation |
steadiness derives the Steadiness Index from a land productivity variable
steadiness(obj2process = NULL, cores2use = 1, filename = "")
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 |
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
SpatRaster object
Xavier Rotllan-Puig
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")}
sb <- terra::rast(paste0(system.file(package='LPDynR'), "/extdata/sb_cat.tif"))
steadiness(obj2process = sb)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.