fetchRaCA {soilDB}R Documentation

Fetch KSSL Data (EXPERIMENTAL)

Description

Get Rapid Carbon Assessment (RaCA) data via state, geographic bounding-box, RaCA site ID, or series query from the SoilWeb system.

Usage

fetchRaCA(series = NULL, bbox = NULL, state = NULL, rcasiteid = NULL, get.vnir = FALSE)

Arguments

series

a soil series name, case insensitive

bbox

a bounding box in WGS84 geographic coordinates e.g. c(-120, 37, -122, 38), constrained to a 5-degree block

state

a two-letter US state abbreviation, case insensitive

rcasiteid

an RaCA site id (e.g. 'C1609C01')

get.vnir

boolean, should associated VNIR spectra be downloaded? (see details)

Details

The VNIR spectra associated with RaCA data are quite large [each gzip-compressed VNIR spectra record is about 6.6kb], so requests for these data are disabled by default. Note that VNIR spectra can only be queried by soil series or geographic BBOX.

Value

pedons:

a SoilProfileCollection object containing site/pedon/horizon data

trees:

a data.frame object containing tree DBH and height

veg:

a data.frame object containing plant species

stock:

a data.frame object containing carbon quantities (stocks) at standardized depths

sample:

a data.frame object containing sample-level bulk density and soil organic carbon values

spectra:

a numeric matrix containing VNIR reflectance spectra from 350–2500 nm

Author(s)

D.E. Beaudette, USDA-NRCS staff

References

http://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/?cid=nrcs142p2_054164 fetchRaCA() Tutorial

See Also

fetchOSD

Examples

## Not run: 
##D   # search by series name
##D   s <- fetchRaCA(series='auburn')
##D   
##D   # search by bounding-box
##D   # s <- fetchRaCA(bbox=c(-120, 37, -122, 38))
##D   
##D   # check structure
##D   str(s, 1)
##D   
##D   # extract pedons
##D   p <- s$pedons
##D   
##D   # how many pedons
##D   length(p)
##D   
##D   # plot 
##D   par(mar=c(0,0,0,0))
##D   plot(p, name='hzn_desgn', max.depth=150)
## End(Not run)

[Package soilDB version 1.6.3 Index]