get_taper: Calculates stem diameter or height based on taper equation...

In sgezan/Nothopack: G&Y Forest Simulator for Roble-Rauli-Coigue Stands in Chile

Description

It reads input of tree-level information (SPECIE, DBH, HT), stand-level (dominant specie, zone), and one taper component (di or hi) to calculate the other element (hi or di) based on a selection of available taper equation models.

Usage

get_taper(
  SPECIES = NA,
  ZONE = NA,
  DBH = NA,
  HT = NA,
  di = NA,
  hi = NA,
  model = 2
)

Arguments

SPECIES	Dominant specie (1:Rauli, 2:Roble, 3:Coigue, 4:Others or Mixed)
ZONE	Growth zone of the corresponding stand
DBH	Diameter at breast height (cm) of tree
HT	Total tree height (m)
di	Stem diameter (cm) at given stem height hi (m)
hi	Stem height (m) at given stem diameter di (cm)
model	Type of selected taper model (1:zone specific, 2:all zones) (default = 2)

Value

The missing component (di or hi) from the requested taper model. Or it returns the complete taper on vectors d and h (evaluated each centimeter up to total height).

Author(s)

S.A. Gezan and P. Moreno

References

Gezan, S.A. and Moreno, P. (2000). Informe Procedimientos y Resultados Modelos de Ahusamiento. Reporte Interno. Projecto FONDEF D97I1065. Chile

Examples

# Example 1: Unknown diameter inside bark stem diameter
(di <- get_taper(SPECIES=1, ZONE=2, DBH=8.1, HT=10.2, hi=4.3, model=1)$di)

# Example 2: Unknown stem height
(hi <- get_taper(SPECIES=1, ZONE=2, DBH=8.1, HT=10.2, di=9.8, model=1)$hi)

# Example 3: Plotting complete taper model for a given tree
tree <- get_taper(SPECIES=1, ZONE=2, DBH=8.1, HT=10.2, di=9.8, model=1)
plot(tree$h,tree$d, ylab='di (cm)', xlab='hi (m)', type='l', col='blue')

sgezan/Nothopack documentation built on April 25, 2021, 8:03 a.m.