get_ssurgo_soil_profile: Retrieve soil profile data and convert it to an object of...

In apsimx: Inspect, Read, Edit and Run 'APSIM' "Next Generation" and 'APSIM' Classic

Retrieve soil profile data and convert it to an object of class ‘soil_profile’

Description

Generate a synthetic soil profile based on the information in SSURGO database

Usage

get_ssurgo_soil_profile(
  lonlat,
  shift = -1,
  nmapunit = 1,
  nsoil = 1,
  xout = NULL,
  soil.bottom = 200,
  method = c("constant", "linear"),
  nlayers = 10,
  verbose = FALSE
)

Arguments

lonlat	Longitude and latitude vector (e.g. c(-93, 42))
shift	simple mechanism for creating an area of interest by displacing the point indicated in lonlat by some amount of distance (e.g. 300 - in meters)
nmapunit	number of mapunits to select (see ssurgo2sp)
nsoil	number of soils to select (see ssurgo2sp)
xout	see ssurgo2sp
soil.bottom	see ssurgo2sp
method	interpolation method see ssurgo2sp
nlayers	number for layer for the new soil profile
verbose	default FALSE. Whether to print messages.

Details

Data source is USDA-NRCS Soil Data Access. See package soilDB for more details

Value

this function will always return a list. Each element of the list will be an object of class ‘soil_profile’

Examples

## Not run: 
require(soilDB)
require(sp)
require(sf)
require(spData)
## Soil inforation for a single point
sp <- get_ssurgo_soil_profile(lonlat = c(-93, 42))
plot(sp[[1]])
plot(sp[[1]], property = "water")

## End(Not run)

apsimx documentation built on March 18, 2022, 7:52 p.m.