treecoverloss_emissions | R Documentation |
This functions allows to efficiently calculate emission statistics for areas of interest. For each year in the analysis timeframe, the forest losses from Hansen et al. (2013) are overlayed with the respective emission layer from Harris et al. (2021) and area-wise emission statistics are calculated for each year.
calc_treecoverloss_emissions(years = 2000:2023, min_size = 10, min_cover = 35)
years |
A numeric vector with the years for which to calculate emissions caused by treecover loss. |
min_size |
The minimum size of a forest patch in ha. |
min_cover |
The minimum threshold of stand density for a pixel to be considered forest in the year 2000. |
The required resources for this indicator are:
gfw_treecover
gfw_lossyear
gfw_emissions
A function that returns an indicator tibble with emissions as variable and emitted CO2 equivalent (in Mg) as value.
## Not run:
library(sf)
library(mapme.biodiversity)
outdir <- file.path(tempdir(), "mapme-data")
dir.create(outdir, showWarnings = FALSE)
mapme_options(
outdir = outdir,
verbose = FALSE
)
aoi <- system.file("extdata", "sierra_de_neiba_478140_2.gpkg",
package = "mapme.biodiversity"
) %>%
read_sf() %>%
get_resources(
get_gfw_treecover(version = "GFC-2023-v1.11"),
get_gfw_lossyear(version = "GFC-2023-v1.11"),
get_gfw_emissions()
) %>%
calc_indicators(
calc_treecoverloss_emissions(years = 2016:2017, min_size = 1, min_cover = 30)
) %>%
portfolio_long()
aoi
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.