Skip to contents

This indicator allows the extraction of zonal statistics for resource layers previously downloaded from SoilGrids, thus in total supporting the calculation of zonal statistics for 10 different soil properties at 6 different depths for a total of 4 different model outputs (stat). Zonal statistics will be calculated for all SoilGrid layers that have been previously made available vie get_resources(). The required resource for this indicator is:

Format

A tibble with a column for the SoilGrid layer, the depth and the model output statistic as well as additional columns for all zonal statistics specified via stats_soil

Details

The following arguments can be set:

stats_soil

Function to be applied to compute statistics for polygons either single or multiple inputs as character. Supported statistics are: "mean", "median", "sd", "min", "max", "sum" "var".

engine

The preferred processing functions from either one of "zonal", "extract" or "exactextract" as character.

Examples

library(sf)
library(mapme.biodiversity)

temp_loc <- file.path(tempdir(), "mapme.biodiversity")
if (!file.exists(temp_loc)) {
  dir.create(temp_loc)
  resource_dir <- system.file("res", package = "mapme.biodiversity")
  file.copy(resource_dir, temp_loc, recursive = TRUE)
}

(try(aoi <- system.file("extdata", "sierra_de_neiba_478140_2.gpkg",
  package = "mapme.biodiversity"
) %>%
  read_sf() %>%
  init_portfolio(
    years = 2022,
    outdir = file.path(temp_loc, "res"),
    tmpdir = tempdir(),
    add_resources = FALSE,
    verbose = FALSE
  ) %>%
  get_resources("soilgrids",
    layers = c("clay", "silt"), depths = c("0-5cm", "5-15cm"), stats = "mean"
  ) %>%
  calc_indicators("soilproperties", stats_soil = c("mean", "median"), engine = "extract") %>%
  tidyr::unnest(soilproperties)))
#> Simple feature collection with 4 features and 10 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: 4 × 11
#>   WDPAID NAME            DESIG_ENG  ISO3  assetid layer depth stat   mean median
#>    <dbl> <chr>           <chr>      <chr>   <int> <chr> <chr> <chr> <dbl>  <dbl>
#> 1 478140 Sierra de Neiba National … DOM         1 clay  0-5cm mean   35.0   34.8
#> 2 478140 Sierra de Neiba National … DOM         1 clay  5-15… mean   35.4   35.3
#> 3 478140 Sierra de Neiba National … DOM         1 silt  0-5cm mean   33.8   33.7
#> 4 478140 Sierra de Neiba National … DOM         1 silt  5-15… mean   33.9   33.8
#> # ℹ 1 more variable: geom <POLYGON [°]>