# load test data
data(sp2)
# promote to SPC
depths(sp2) <- id ~ top + bottom
# one profile
p <- sp2[1,]
# 23 variables in horizon data
length(horizonNames(sp2))
# remove all but essential ones
replaceHorizons(p) <- horizons(p)[,c(idname(p), hzidname(p), horizonDepths(p))]
# inspect result (a clean slate)
horizons(p)
Run the code above in your browser using DataLab