addtreedim: Compute tree dimensions (tree and crown height and diameter,...

In thomasgaquiere/Maria: Simulate Different Forms of Logging on a Forest

Description

Compute tree dimensions (tree and crown height and diameter, harvestable volume, wood density, and AGB)

Usage

addtreedim(
  inventory,
  volumeparameters,
  crowndiameterparameters = ParamCrownDiameterAllometry,
  advancedloggingparameters = loggingparameters()
)

Arguments

inventory	Input inventory (see the inputs formats and metadata in the vignette) (data.frame)
volumeparameters	Parameters of the tree volume tables (in the same format of ForestZoneVolumeParametersTable) used to compute the harvestable volume of each tree, depending on its geographic zone if several locations (data.frame)
crowndiameterparameters	Parameters of crown diameter allometry (in the same format of ParamCrownDiameterAllometry) used to compute the crown diameter of each tree, depending on its DBH (Diameter at Breast Height) and its species, genus or family. (data.frame)
advancedloggingparameters	Other parameters of the logging simulator loggingparameters (list)

Details

'addtreedim' compute some tree dimensions with the 'advancedloggingparameters' argument:

• Tree height (in m) ('TreeHeightAllometry')

• Tree Harvestable Volume (m^3) ('TreeHarvestableVolumeAllometry')

• Trunk height (in m) ('TrunkHeightAllometry')

• Crown height (in m)

• Crown diameter (in m) ('CrownDiameterAllometry')

• Wood density (g/cm^3) (by BIOMASS package)

• Tree above-ground biomass (AGB) (in Mg) (by BIOMASS package)

Value

The inventory (data.frame) with additional variables: TreeHeight (m), TreeHarvestableVolume (m^3), TrunkHeight (m), CrownHeight (m), CrownDiameter (m), estimated 'WoodDensity' (g/cm^3) and tree above-ground biomass ('AGB') (in Mg)

Examples

data(Paracou6_2016)
data(ForestZoneVolumeParametersTable) # The volume parameters data in the global environment
data(ParamCrownDiameterAllometry)

if (!("DBH" %in% names(Paracou6_2016))) {
tibble::add_column(Paracou6_2016, DBH = NA) #if DBH doesn't exist create it
Paracou6_2016$DBH = Paracou6_2016$CircCorr/pi} # and calculate it
Paracou6_2016 <- dplyr::filter(Paracou6_2016, DBH >= 10)

addtreedim(Paracou6_2016, volumeparameters = ForestZoneVolumeParametersTable)

thomasgaquiere/Maria documentation built on Dec. 24, 2021, 1:20 a.m.