# NOT RUN {
if(require(aqp)) {
# path to local PedonPC back-end DB
dsn <- "S:/Service_Center/NRCS/pedon/pedon.accdb"
if(file.exists(dsn)) {
# get routinely used soil data SoilProfileCollection object
f <- fetchPedonPC(dsn)
# determine which profiles have densic contacts
idx <- which(f$densic.contact)
# plot only those profiles with densic contact
if(length(idx))
plot(f[idx, ], name='hzname')
}
}
# }
Run the code above in your browser using DataLab