taxaExtent {soilDB}R Documentation

Get SoilWeb 800m Major Component Soil Taxonomy Grids

Description

This function downloads a generalized representation of the geographic extent of any single taxon from the top 4 levels of Soil Taxonomy, or taxa matching a given formative element used in Great Group or subgroup taxa. Data are provided by SoilWeb, ultimately sourced from the current SSURGO snapshot. Data are returned as raster objects representing area proportion falling within 800m cells. Currently area proportions are based on major components only. Data are only available in CONUS and returned using an Albers Equal Area / NAD83(2011) coordinate reference system (EPSG: 5070).

Usage

taxaExtent(
  x,
  level = c("order", "suborder", "greatgroup", "subgroup"),
  formativeElement = FALSE,
  timeout = 60,
  as_Spatial = getOption("soilDB.return_Spatial", default = FALSE)
)

Arguments

x

single taxon label (e.g. haploxeralfs) or formative element (e.g. pale), case-insensitive

level

the taxonomic level within the top 4 tiers of Soil Taxonomy, one of 'order', 'suborder', 'greatgroup', 'subgroup'

formativeElement

logical, search using formative elements instead of taxon label

timeout

time that we are willing to wait for a response, in seconds

as_Spatial

Return raster (RasterLayer) classes? Default: FALSE.

Details

See the Geographic Extent of Soil Taxa tutorial for more detailed examples.

Taxon Queries

Taxon labels can be conveniently extracted from the "ST_unique_list" sample data, provided by the SoilTaxonomy package.

Formative Element Queries

Greatgroup:

The following labels are used to access taxa containing the following formative elements (in parentheses)

Subgroup:

The following labels are used to access taxa containing the following formative elements (in parenthesis).

Value

a SpatRaster object (or RasterLayer when as_Spatial=TRUE)

Author(s)

D.E. Beaudette and A.G. Brown

Examples



## Not run: 
  library(terra)
  
  # soil order
  taxa <- 'vertisols'
  x <- taxaExtent(taxa, level = 'order')
  
  # suborder
  taxa <- 'ustalfs'
  x <- taxaExtent(taxa, level = 'suborder')
  
  # greatgroup
  taxa <- 'haplohumults'
  x <- taxaExtent(taxa, level = 'greatgroup')
  
  # subgroup
  taxa <- 'Typic Haploxerepts'
  x <- taxaExtent(taxa, level = 'subgroup')
  
  # greatgroup formative element
  taxa <- 'psamm'
  x <- taxaExtent(taxa, level = 'greatgroup', formativeElement = TRUE)
  
  # subgroup formative element
  taxa <- 'abruptic'
  x <- taxaExtent(taxa, level = 'subgroup', formativeElement = TRUE)
  
  # coarsen for faster plotting
  a <- terra::aggregate(x, fact = 5, na.rm = TRUE)
  
  # quick evaluation of the result
  terra::plot(a, axes = FALSE)

## End(Not run)

[Package soilDB version 2.8.3 Index]