LossCalc: Area of Forest Loss

In mapme-initiative/mapme.forest: Package to Support Project Evaluation and Planning in the Forest Sector

Area of Forest Loss

Description

Calculation of the area where forest loss occured in a tree cover time series. It uses information on the occurence of forest loss and a forest mask to focus the analysis to pixles where lost occured. The area is then summed up for every polygon specified in the studysite object.

Usage

LossCalc(
  inputForestMap = NULL,
  inputLossMap = NULL,
  studysite = NULL,
  latlon = NULL,
  polyName = NULL,
  ncores = 1,
  saveCSV = FALSE,
  years = 2001:2018
)

Arguments

inputForestMap	A RasterLayer object with forest cover represented in binary (e.g 0 represents no forest; 1 represents forest). Each layer is expected to represent one year in the years object in a consecutive order.
inputLossMap	A RasterLayer object with values indicating the year tree cover was lost. The annual area of tree cover lost is calculated based on the values found in this data set.
studysite	An object of type sf with a given number of polygons defining the areas of interest. Forest area statistics will be calculated for features in the studysite object.
latlon	logical: Indicates whether or not the inputForestMap object is based on a geographic coordinate system or is projected to a planar coordinate system. In the former case, area is approximated by raster::area.
polyName	charachter of length 1. Indicates the column in the data frame of the studysite object to uniquely identify the features of interest. The function will fail if there is no unique identification of the polygons
ncores	The number of cores to use, i.e. at most how many child processes will be run simultaneously.
saveCSV	Default is FALSE. You can specify a directory on your local machine where the results area saved in .csv format. The features will be identified by the column specified in polyName.
years	A vector of type numeric indicating the years which are represented by pixels. For GFW data and the default these are 2001 to 2018.

Value

The studysite object with it's attribute table being ammended by the results of the calculation. Additionally, a csv file containing the polyName attribute as well as the results of the calculation can be saved to disk.

Note

This function relies heavily on parallization, indicating the importance of both, a high number of CPUs and large enough RAM.

Author(s)

Darius Görgen (MapTailor Geospatial Consulting GbR) [e-mail]
Maintainer: MAPME-Initiative [e-mail]
Contact Person: Dr. Johannes Schielein
Copyright: MAPME-Initiative
License: GPL-3

Examples

library(sf)
library(raster)
library(mapme.forest)

aoi = st_read(system.file("extdata", "aoi_polys.gpkg", package = "mapme.forest"))
yearlyRaster = stack(system.file("extdata", "pkgTest_yearlyCover.tif",
                                  package = "mapme.forest"))
lossYear = raster(system.file("extdata", "pkgTest_lossyear.tif",
                              package = "mapme.forest"))

result = LossCalc(inputForestMap = yearlyRaster,
                 inputLossMap = lossYear,
                 studysite = aoi[1,],
                 latlon = TRUE,
                 polyName = "id",
                 ncores = 1,
                 saveCSV = FALSE,
                 years = 2000:2018)
str(result)

mapme-initiative/mapme.forest documentation built on April 15, 2022, 4:50 a.m.