aqp (version 1.18)

get.increase.depths: Return the horizon top depths from a call to get.increase.matrix()

Description

get.increase.depths performs the conversion of the square matrix output of get.increase.matrix back to horizon top depth for where criteria were met.

Note that the `threshold.fun` result is allowed to contain NA, but that will result in no output for affected cells.

Usage

get.increase.depths(p, attr, threshold.fun, vertical.distance)

Arguments

p

a SoilProfileCollection, containing a single profile

attr

horizon attribute name to get the "increase" of

threshold.fun

a function that returns the threshold (as a function of attr); may return a constant single value

vertical.distance

the vertical distance (determined from difference SPC top depth variable) within which increase must be met

Value

Returns a numeric vector of depths where the increase retuirement is met. Typically the first is the one of interst to taxonomy.

get.increase.depths converts to horizon dop depth by using above matrix output to determine depths where increase is met.

See Also

getArgillicBounds, crit.clay.argillic

Examples

# NOT RUN {
data(sp1, package = 'aqp')
depths(sp1) <- id ~ top + bottom
site(sp1) <- ~ group

p <- sp1[1]
attr <- 'prop' # clay contents 
foo <- get.increase.depths(p, threshold.fun = crit.clay.argillic, 
                           attr = attr, vertical.distance = 30)
foo
# }