Download SoilGrids layers — get

A function to download SoilGrids dataset, created by Hengl et al. (2017)

Usage

get_soilgrids(
  aoi,
  layer = c("all", "sand", "silt", "clay", "soc", "phh2o", "bdod", "cfvo"),
  pred.quantile = "mean"
)

Arguments

aoi: SpatVector. A polygon layer with area of interest.
layer: character. A string indicating what layers should be downloaded. Either one of the following: 'all', 'sand', 'clay', 'silt', 'soc', bdod, cfvo or 'phh2o'
pred.quantile: character. A string indicating what predition quantiles should be downloaded. Either one of the following: 'mean' (default), 'Q0.05', 'Q0.95' or 'Q0.05' (see Details).

Value

SpatRaster

References

Hengl, Tomislav, Jorge Mendes de Jesus, Gerard B. M. Heuvelink, Maria Ruiperez Gonzalez, Milan Kilibarda, Aleksandar Blagotić, Wei Shangguan, et al. “SoilGrids250m: Global Gridded Soil Information Based on Machine Learning.” PLOS ONE 12, no. 2 (February 16, 2017): e0169748. https://doi.org/10.1371/journal.pone.0169748.

Examples

library(terra)
library(purrr)

f <- system.file("extdata/extent.shp", package="rusleR")
v <- vect(f)

sand <- get_soilgrids(v, layer = "sand")
#> Downloading Sand rasters...

plot(sand)