| fetchRaCA {soilDB} | R Documentation |
Get Rapid Carbon Assessment (RaCA) data via state, geographic bounding-box, RaCA site ID, or series query from the SoilWeb system.
fetchRaCA(series = NULL, bbox = NULL, state = NULL, rcasiteid = NULL, get.vnir = FALSE)
series |
a soil series name, case insensitive |
bbox |
a bounding box in WGS84 geographic coordinates e.g. |
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) |
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.
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
D.E. Beaudette, USDA-NRCS staff
http://www.nrcs.usda.gov/wps/portal/nrcs/detail/soils/survey/?cid=nrcs142p2_054164 fetchRaCA() Tutorial
## 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)