ForestElementsR (version 2.0.1)

h_age_gnfi3: Estimate Tree Height Growth With the 3rd German National Forest Inventory Growth Model (2012)

Description

Tree height growth model of the third German National Forest Inventory of 2012 bwi3_methods_2017ForestElementsR. Allows to estimate a tree's height at any age if its height is known at a given age.

Usage

h_age_gnfi3(species_id, age_yr, h_m_known, age_yr_known)

Value

A single height value or vector of height values corresponding to age_yr

Arguments

species_id

Vector of species id's preferably following the ger_nfi_2012 species coding. Ideally, these species_id's are provided as a fe_species_ger_nfi_2012 object. See Details for how other species codings are handled.

age_yr

Single numeric value or vector of ages (in years) for which the height is to be calculated

h_m_known

Vector of known height (m) values at age age_yr_known

age_yr_known

Vector of ages (years) for which the height h_m_known is known

Details

Originally, the function was parameterized for species and species groups corresponding to the national forest inventory's species coding (fe_species_ger_nfi_2012). We have attributed in addition these the original parameters also to the species codings fe_species_tum_wwk_short, and fe_species_bavrn_state_short. When called with a given species coding, the function will try to use the "nearest" of these three alternatives. Fallback option is the attempt to use fe_species_tum_wwk_short.

References

See Also

Other growth functions: age_d_gnfi3(), age_h_gnfi3(), d_age_gnfi3()

Examples

  # A European beech tree has a height of 25.2 m at age 75. Estimate
  # its height at age 83
  h_age_gnfi3(100, 83, 25.2, 75) # 100 is ger_nfi_2012 code for E. beech

  # Do the same backward in time, height at age 67
  h_age_gnfi3(100, 67, 25.2, 75)

  # Apply for more than one tree, different species, same age
  h_known   <- c(23.1, 16.2, 35.2, 19.3, 21.8)
  species   <- as_fe_species_tum_wwk_short(c(3, 3, 3, 6, 6))
  h_age_gnfi3(
    species, age_yr = 70, h_m_known = h_known, age_yr_known = 60
  )