baseline_lev: baseline_lev

View source: R/03_baseline_lev.R

baseline_levR Documentation

baseline_lev

Description

baseline_lev() derives land productivity at the beginning of the time series on study, resulting in a 3-class SpatRaster object with (1) low, (2) medium and (3) high productivity

Usage

baseline_lev(
  obj2process = NULL,
  yearsBaseline = 3,
  drylandProp = 0.4,
  highprodProp = 0.1,
  cores2use = 1,
  filename = ""
)

Arguments

obj2process

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

yearsBaseline

Numeric. Number of years to be averaged and used as baseline. Optional. Default is 3

drylandProp

Numeric. Proportion of drylands over total land, either expressed as a fraction of unity or percentage. Optional. Default is 0.4

highprodProp

Numeric. Proportion of land classified as 'highly productive' over total land, either expressed as a fraction of unity or percentage. Optional. Default is 0.1

cores2use

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

filename

Character. Output filename. Optional

Details

baseline_lev() uses the proportion of drylands over the total land ('drylandProp') to classify the level of productivity into low level. UNPD declares that 40 percent of the World’s land resources are drylands (Middleton et al., 2011) and, therefore, 40 percent of pixels at the global level can be classified as low productivity land. This assumption is the default, but it should be adjusted for local and regional studies. In addition, baseline_lev() classifies by default 10 percent of pixels as high level of land productivity and the rest (100 - ('drylandProp' + 10)) as medium level. Proportion of pixels classified as 'high' can be also modified by passing the argument 'highprodProp'

Value

SpatRaster object

Author(s)

Xavier Rotllan-Puig

References

Middleton, N., L. Stringer, A. Goudie, and D. Thomas. 2011. “The Forgotten Billion. MDG Achievement in the Drylands.” New York, NY, 10017, USA: United Nations Development Programme.

Examples


sb <- terra::rast(paste0(system.file(package='LPDynR'), "/extdata/sb_cat.tif"))
baseline_lev(obj2process = sb,
             yearsBaseline = 3,
             drylandProp = 0.4,
             cores2use = 2)



LPDynR documentation built on Oct. 16, 2023, 5:06 p.m.