Get map unit keys (mukeys) from gSSURGO

These functions query the NRCS gSSURGO Web Feature Service to retrieve map unit keys based on different spatial filters.

Usage

ssurgo_mukeys_bbox(bbox)

ssurgo_mukeys_point(point, distance)

ssurgo_mukeys_bigbbox(bbox)

Arguments

bbox

Numeric vector of length 4: c(xmin, ymin, xmax, ymax) in WGS84 (EPSG:4326). Features that intersect the bounding box are returned.

point

Numeric vector of length 2: c(lon, lat) in WGS84 (EPSG:4326).

distance

Numeric. Distance in meters from the point. Use 0 for exact point intersection.

Value

Character vector of unique map unit keys (mukeys).

Details

These functions use the NRCS SDM Data Access Web Feature Service: https://sdmdataaccess.nrcs.usda.gov/SpatialFilterHelp.htm

The total extent cannot exceed 10,100,000,000 square meters (~3,900 square miles). Use `ssurgo_mukeys_bigbbox()` for large bounding boxes.

Examples

if (FALSE) { # \dontrun{
# Bounding box query
mukeys <- ssurgo_mukeys_bbox(bbox = c(-114.006, 32.1823, -113.806, 32.2823))

# Point with distance (600m radius)
mukeys <- ssurgo_mukeys_point(point = c(-91.22, 38.46), distance = 600)

# Large bounding box
mukeys <- ssurgo_mukeys_bigbbox(bbox = c(-120, 35, -110, 45))
} # }