# NOT RUN {
library(soilDB)
# soils of interest
s.list <- c('musick', 'cecil', 'drummer', 'amador', 'pentz', 'reiff',
'san joaquin','montpellier','grangeville','pollasky','ramona')
# fetch and convert data into an SPC
s <- fetchOSD(s.list, extended=TRUE)
res <- vizHillslopePosition(s$hillpos)
print(res$fig)
# }
Run the code above in your browser using DataLab