mukey) grid from SoilWeb Web Coverage Service (WCS)
Download chunks of the gNATSGO or gSSURGO map unit key grid via bounding-box from the SoilWeb WCS.
mukey.wcs(aoi, db = c("gnatsgo", "gssurgo"), res = 30, quiet = FALSE)
area of interest (AOI) defined using either a
name of the gridded map unit key grid to access, should be either 'gnatsgo' or 'gssurgo'
grid resolution, units of meters. The native resolution of gNATSGO and gSSURGO (this WCS) is 30m.
logical, passed to
raster object containing indexed map unit keys and associated raster attribute table or a try-error if request fails
aoi should be specified as one of:
bbox object, or a
bounding-box specified as (xmin, ymin, xmax, ymax) e.g. c(-114.16, 47.65, -114.08, 47.68)
coordinate reference system of BBOX, e.g. '+init=epsg:4326'
The WCS query is parameterized using
raster::extent derived from the above AOI specification, after conversion to the native CRS (EPSG:5070) of the gNATSGO / gSSURGO grid.
Databases available from this WCS can be queried using
WCS_details(wcs = 'mukey').
The gNATSGO grid includes raster soil survey map unit keys which are not in SDA.
D.E. Beaudette and A.G. Brown