aqp (version 1.18)

glom: Return a "ragged" group of horizons selected from SoilProfileCollection

Description

glom() returns a "clod" of horizons from a (often single profile) SoilProfileCollection that have a common attribute. You "glom" SPC horizons into a "clod".

A clod is a "ragged" group of horizons. In this case, "ragged" means that number of horizons, horizon depths, distinctness and topography vary from pedon to pedon.

All horizons included within the specified interval are returned in their entirety (not just the portion within the interval). A 'clod' references the original pedon data -- it is not resampled like slice or slab.

Horizon intersection is based on unique ID hzidname(spc) and attribute of interest.

Currently, just depth/depth intervals are supported as common attribute for creating a clod i.e. all horizons overlapping interval 25 to 100 cm -- but many other attributes are possible and will be implemented perhaps via formula interface similar to slice.

Usage

glom(p, z1, z2 = NA, as.data.frame = FALSE)

Arguments

p

A single-profile SoilProfileCollection; e.g. glom is called via profileApply()

z1

Top depth (required) - depth to intersect horizon at; if `z2` specified, top depth of intersect interval.

z2

OPTIONAL: Bottom depth - bottom depth of intersection interval

as.data.frame

OPTIONAL: return a data.frame, by intersection with horizons(p)? default: FALSE

Value

data.frame or SoilProfileCllection representation of the "clod" of horizons that have a common depth range

See Also

slice, slab

Examples

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

p <- sp1[1]

foo <- glom(p, 25, 100)

# there are 4 horizons in the clod glommed from depths 25 to 100 on profile 1 in sp1
nrow(foo) 
# }