Extracts the major soil order at the GPS points provided.

extract_daas_soil_order(x, cache = FALSE)



List. An object with names containing vector pairings of longitude and latitude values expressed as decimal degree values in that order, each individual vectors' items should be named "x" (longitude) and "y" (latitude), respectively. The list item names should be descriptive of the individual vectors and will be included in a "location" column of the output.


Boolean. Store soil data locally for later use? If FALSE, the downloaded files are removed when R session is closed. To take advantage of cached files in future sessions, use cache = TRUE. Defaults to FALSE. Value is optional. All future requests will use the cached data unless remove_cache() is used to remove the cached file.


A data.table::data.table with the provided GPS

coordinates and the respective Digital Atlas of Australian Soils (DAAS soil order), "Spatial Data Conversion of the Atlas of Australian Soils to the Australian Soil Classification v01".


The first run will take additional time to download and extract the soils data. If cache = TRUE, any use after this will be much faster due to the locally cached geospatial soils data.


if (FALSE) { # interactive()
locs <- list(
  "Merredin" = c(x = 118.28, y = -31.48),
  "Corrigin" = c(x = 117.87, y = -32.33),
  "Tamworth" = c(x = 150.84, y = -31.07)

extract_daas_soil_order(x = locs, cache = FALSE)