| estimateSoilDepth {aqp} | R Documentation |
Estimate the soil depth of a single profile within a SoilProfileCollection object.
estimateSoilDepth(f, name = "hzname", top = "hzdept", bottom = "hzdepb", p = 'Cr|R|Cd', no.contact.depth = NULL, no.contact.assigned = NULL)
f |
A SoilProfileCollection object of length 1, e.g. a single profile |
name |
the name of the column that contains horizon designations |
top |
the name of the column that contains horizon top depths |
bottom |
the name of the column that contains horizon bottom depths |
p |
a REGEX pattern for determining contact with bedrock |
no.contact.depth |
in the absense of contact with bedrock, a depth at which we can assume a standard depth |
no.contact.assigned |
assumed standard depth |
a single integer, the soil depth
D.E. Beaudette and J.M. Skovlin
getSoilDepthClass, profileApply
data(sp1)
depths(sp1) <- id ~ top + bottom
# apply to each profile in a collection, and save as site-level attribute
sp1$depth <- profileApply(sp1, estimateSoilDepth, name='name', top='top', bottom='bottom')
## Not run:
##D # sample data
##D data(gopheridge, package='soilDB')
##D
##D # run on a single profile
##D estimateSoilDepth(gopheridge[1, ])
##D
##D # apply to an entire collection
##D profileApply(gopheridge, estimateSoilDepth)
## End(Not run)