Skip to contents

Read SLGA soil data from Google Earth Engine

Source: R/read_slga.R
read_slga.Rd

Reads CSIRO Soil and Landscape Grid of Australia (SLGA) data from Google Earth Engine. Supports 7 soil attributes, 6 depth layers, and 3 statistics. This is a static dataset (no date required).

Usage

read_slga(
  region,
  collection = c("CLY", "SND", "SLT", "AWC", "BDW", "PHC", "NTO"),
  depth = "0-5",
  stat = c("mean", "ci_lower", "ci_upper"),
  backend = c("rest", "rgee"),
  cache = TRUE,
  max_tries = 3L,
  initial_delay = 1
)

Arguments

region

An sf::sf, sf::st_sfc(), or terra::ext() object defining the spatial extent. Required.

collection

Character. Soil attribute to extract. One of:

  • "CLY" — Clay content (percent)

  • "SND" — Sand content (percent)

  • "SLT" — Silt content (percent)

  • "AWC" — Available water capacity (mm)

  • "BDW" — Bulk density (g/cm3)

  • "PHC" — pH (CaCl2)

  • "NTO" — Total nitrogen (percent)

Default "CLY".

depth

Character. Depth layer. One of "0-5", "5-15", "15-30", "30-60", "60-100", "100-200". Default "0-5".

stat

Character. Statistic. "mean" (default), "ci_lower" (5th percentile), or "ci_upper" (95th percentile).

backend

Character. "rest" (default) or "rgee".

cache

Logical. Use disk cache? Default TRUE.

max_tries

Integer. Retry attempts. Default 3L.

initial_delay

Numeric. Initial retry delay in seconds. Default 1.

Value

A terra::rast() SpatRaster. CRS: EPSG:4326. Values depend on the collection (e.g., percent for clay, pH units for PHC).

Data availability

Static dataset. Australia only. 90 m spatial resolution. 6 depth layers: 0-5, 5-15, 15-30, 30-60, 60-100, 100-200 cm.

References

Viscarra Rossel, R.A. et al. (2015). A digital soil map of Australia. Soil Research, 53(7), 745-757. doi:10.1071/SR15158

See also

read_gee(), collect_gee_data()

Other GEE readers: read_chirps(), read_era5(), read_gee(), read_landsat(), read_modis_lst(), read_modis_ndvi(), read_sentinel2(), read_srtm(), read_worldclim()

Examples

if (FALSE) { # interactive()
# Clay content, 0-5cm, mean
clay <- read_slga(region = terra::ext(138, 140, -36, -34))

# pH at 30-60cm depth
ph <- read_slga(region = terra::ext(138, 140, -36, -34),
                collection = "PHC", depth = "30-60")

# Sand content with confidence interval
sand_lo <- read_slga(region = terra::ext(138, 140, -36, -34),
                     collection = "SND", stat = "ci_lower")
}