data-raw/biodiversity_intactness_index.R

In mapme-initiative/mapme.biodiversity: Efficient Monitoring of Global Biodiversity Portfolios

library(mapme.biodiversity)
library(terra)

mapme_options(outdir = NULL)

x <- read_sf(
  system.file("extdata", "shell_beach_protected_area_41057_B.gpkg",
              package = "mapme.biodiversity")) %>%
get_resources(get_biodiversity_intactness_index("lbii.asc"))

bii <- prep_resources(x)
bii <- bii [[1]]

outdir <- "inst/resources/biodiversity_intactness_index"
dir.create(outdir, recursive = TRUE, showWarnings = FALSE)

bii <- mask(bii, x)
writeRaster(bii, file.path(outdir, "lbii.asc"))
file.remove(c(file.path(outdir, "lbii.prj"),
              file.path(outdir, "lbii.asc.aux.xml")))