AreaCalc: Area of forest change

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

Area of forest change

Description

Calculates the changes in forest area based on a rasterObject forest mask and a SpatialPolygonsDataFrame containing areas of interest.

Usage

AreaCalc(
  inputForestMap = 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.
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 its data frame appended columnwise for every single entry in the years object representing the total area of forest in a given layer of the inputForestMap object. If latlon=TRUE the returned area is in km², otherwise in the squared unit of the input projection (most commonly in meters).

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"))
result = AreaCalc(inputForestMap = yearlyRaster,
                  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.