kkyong77/soilDB

soilDB: Simplified Access to NCSS Soil Databases

Installation

Get the stable version from CRAN:

install.packages('soilDB', dep=TRUE)

Get the development version from Github:

remotes::install_github("ncss-tech/soilDB", dependencies=FALSE, upgrade=FALSE, build=FALSE)

Website

http://ncss-tech.github.io/AQP/

Important changes in soilDB-2.0

NASIS

• low-level functions return empty data.frame objects when local database (or selected set) is empty
• fetchNASIS() is now a wrapper around pedon and component "fetch" functions
• uncode() is now used in all queries to local NASIS database
• many (coded) column names in pedon @site and @horizons have changed
• lower cases is used for all un-coded data

Database APIs

• SDA

  • fetchSDA
  • fetchSDA_spatial
  • SDA_query
  • SDA_query_features

• SSURGO/KSSL via SoilWeb

  • fetchKSSL
  • fetchOSD
  • siblings
  • OSDquery
  • seriesExtent

• SCAN/SNOTEL

  • fetchSCAN
  • SCAN_SNOTEL_metadata

• Henry Mount Soil and Water Database

  • fetchHenry

• NASIS local database

  • fetchNASIS

• NASIS WWW interface

  • parseWebReport
  • fetchNASISWebReport

Utility Functions

• estimateSTR
• STRplot
• KSSL_VG_model
• simplfyFragmentData
• simplifyColorData
• uncode
• code

Related Documentation

• fetchKSSL
• SDA_query
• fetchOSD
• SCAN/SNOTEL Data

Related Packages

• aqp
• sharpshootR

Examples

library(soilDB)
library(sharpshootR)
library(latticeExtra)

# get morphology + extended summaries
soils <- c('cecil', 'altavista', 'lloyd', 'wickham', 'wilkes',  'chewacla', 'congaree')
s <- fetchOSD(soils, extended = TRUE)

# experimental viz of hillslope position, from SSURGO component records
res <- vizHillslopePosition(s$hillpos)
print(res$fig)

# profile sketches
par(mar=c(0,1,0,4), xpd=NA)
plot(s$SPC, plot.order=res$order, cex.names=1, axis.line.offset = -0.1, width=0.2)


# siblings
s <- 'Amador'
amador <- siblings(s, only.major = FALSE, component.data = TRUE)

# limit to named soil series
sib.data <- subset(amador$sib.data, subset= ! compkind %in% c('Miscellaneous area', 'Family', 'Taxon above family'))

# get parsed OSD records
sibs <- fetchOSD(c(s, unique(amador$sib$sibling)), extended = TRUE)

# order by subgroup taxonomy
# invert colors
par(mar=c(0,0,0,0), fg='white', bg='black')
SoilTaxonomyDendrogram(sibs$SPC, dend.width = 1.5, y.offset = 0.4, scaling.factor = 0.02, width=0.2, cex.taxon.labels = 1, cex.names = 1)

Dependency Graph