This function extracts areal statistics for the drivers of deforestation based on the data source produced by Fritz et al (2022). The required resource is:
Format
A tibble with 3 columns indicating the class of a deforestation driver, the absolute area in ha, and the percentage in relation to the total area of forest loss as indicated by the Fritz et al. (2022) resource.
Examples
# \dontshow{
mapme.biodiversity:::.copy_resource_dir(file.path(tempdir(), "mapme-data"))
# }
# \dontrun{
library(sf)
library(mapme.biodiversity)
outdir <- file.path(tempdir(), "mapme-data")
dir.create(outdir, showWarnings = FALSE)
aoi <- system.file("extdata", "sierra_de_neiba_478140_2.gpkg",
package = "mapme.biodiversity"
) %>%
read_sf() %>%
init_portfolio(
years = 2021,
outdir = outdir,
tmpdir = tempdir(),
add_resources = FALSE,
verbose = FALSE
) %>%
get_resources("fritz_et_al", res_drivers = 100) %>%
calc_indicators("deforestation_drivers") %>%
tidyr::unnest(deforestation_drivers)
aoi
#> Simple feature collection with 10 features and 8 fields
#> Geometry type: POLYGON
#> Dimension: XY
#> Bounding box: xmin: -71.80933 ymin: 18.57668 xmax: -71.33201 ymax: 18.69931
#> Geodetic CRS: WGS 84
#> # A tibble: 10 × 9
#> WDPAID NAME DESIG_ENG ISO3 assetid class area percent
#> <dbl> <chr> <chr> <chr> <int> <chr> <dbl> <dbl>
#> 1 478140 Sierra de Neiba National Park DOM 1 commercial… 0 0
#> 2 478140 Sierra de Neiba National Park DOM 1 commercial… 0 0
#> 3 478140 Sierra de Neiba National Park DOM 1 managed fo… 0 0
#> 4 478140 Sierra de Neiba National Park DOM 1 mining 0 0
#> 5 478140 Sierra de Neiba National Park DOM 1 natural di… 0 0
#> 6 478140 Sierra de Neiba National Park DOM 1 pasture 0 0
#> 7 478140 Sierra de Neiba National Park DOM 1 roads 0 0
#> 8 478140 Sierra de Neiba National Park DOM 1 wildfire 0 0
#> 9 478140 Sierra de Neiba National Park DOM 1 other subs… 16115. 1
#> 10 478140 Sierra de Neiba National Park DOM 1 shifting c… 0 0
#> # ℹ 1 more variable: geom <POLYGON [°]>
# }