Arrange OSD Data Graphically

# need these
library(soilDB)
library(sharpshootR)

# soils of interest
s.list <- c('hornitos', 'perkins', 'argonaut', 'inks', 'mokelumne', 'dunstone', 'auburn', 'pentz', 'pardee', 'peters', 'amador', 'laniger')

# fetch and convert data into an SPC
s <- fetchOSD(s.list)
# estimate soil depth based on horizon designation
s$soil.depth <- profileApply(s, estimateSoilDepth, name='hzname', top='top', bottom='bottom')

# plot profiles, order by soil depth
par(mar=c(0,0,0,0))
plot(s, name='hzname', plot.order=order(s$soil.depth), cex.names=0.85, axis.line.offset=-5)

plot of chunk check-via-plotting

# plot dendrogram + profiles
SoilTaxonomyDendrogram(s, cex.taxon.labels=0.8)

plot of chunk plot-dendrogram


This document is based on aqp version 1.7-4 and soilDB version 1.3-4.