This function will convert USDA-NCSS horizon boundary topography codes into a vertical offset, suitable for use in plotSPC. Default values are reasonable starting points for encoding smooth, wavy, irregular, or broken style horizon boundary topography as defined in Field Book for Describing and Sampling Soils, version 3.0.

hzTopographyCodeToOffset(
  x,
  codes = c("smooth", "wavy", "irregular", "broken"),
  offset = c(0, 4, 8, 12)
)

Arguments

x

vector of boundary topography codes to be converted

codes

character vector of topography terms ('smooth') or codes ('S'), case insensitive, see details

offset

vertical offset (depth units) used to create "chevron" effect

Value

vector of vertical offsets with same length as x

Details

Additional examples are available in the Visualization of Horizon Boundaries tutorial. Missing data in x (NA) or codes that are not defined in codes are returned with an offset of 0.

Either format (or mixture) are accepted, case insensitive:

  • terms: c('smooth', 'wavy', 'irregular', 'broken')

  • coded values: c('s', 'w', 'i', 'b')

See also

Author

D.E. Beaudette