Passer au contenu

This function generates a standardised raster of managed forests, optionally including disturbed forests and a separate layer for productive managed forests (based on NPP data). The input raster should contain LULC classification values which are reclassified according to defined managed forest class codes.

Utilisation

make_managed_forests(
  raster_in,
  pus,
  iso3,
  lulc_product = c("esri_10m", "dynamic_world", "esa_worldcover", "local"),
  forest_classes = NULL,
  include_disturbed_forest = FALSE,
  make_productive = FALSE,
  raster_npp = NULL,
  name_out,
  output_path = NULL
)

Arguments

raster_in

A SpatRaster representing forest management classifications.

pus

A SpatRaster used as a template for extent, resolution, and alignment.

iso3

A 3-letter ISO country code string.

lulc_product

Character. LULC product used for class value lookups: "esri_10m" (default), "dynamic_world", "esa_worldcover", or "local". When "local", forest_classes must be explicitly provided.

forest_classes

Integer vector or NULL. Class codes representing managed forests. If NULL, automatically determined from lulc_product. Default is NULL.

include_disturbed_forest

Logical. If TRUE, includes disturbed forest class (code 20) in forest_classes.

make_productive

Logical. If TRUE, also generates a productive managed forest raster.

raster_npp

A SpatRaster of NPP data used for generating the productive forest layer. Required if make_productive = TRUE.

name_out

Character. Name stem for output file(s) (optional, not currently used).

output_path

Optional character. Directory to save output raster(s) if provided.

Valeur de retour

A SpatRaster with either one layer (managed forests) or two layers (managed forests and productive managed forests).

Détails

If make_productive = TRUE, the function multiplies the managed forest raster by a normalised raster of Net Primary Productivity (NPP) to generate a layer representing productive managed forests.

Exemples

if (FALSE) { # \dontrun{
managed_forests <- make_managed_forests(
  raster_in = forest_lulc,
  pus = planning_units,
  iso3 = "NPL",
  make_productive = TRUE,
  raster_npp = npp_raster
)
} # }