futurereserve: Future & reserve trees designation

In VincyaneBadouard/LoggingLab: An R Package to Simulate Forest Logging

Future & reserve trees designation

Description

designates which trees will be harvested in future seasons: "future trees", and which should be protected as seed trees: "reserve trees".

Usage

futurereserve(
  inventory,
  speciescriteria,
  advancedloggingparameters = loggingparameters()
)

Arguments

inventory	Input inventory (see the inputs formats and metadata in the vignette) (data.frame)
speciescriteria	Table of species exploitability criteria : species names, economic interest level, minimum and maximum felling diameter, in the same format as SpeciesCriteria (data.frame)
advancedloggingparameters	Other parameters of the logging simulator loggingparameters (list)

Details

Future trees are all trees satisfying the following conditions:

• species of 1st economic rank

• DBH between 35cm ('FutureTreesMinDiameter') and the species MinFD or UpMinFD if it has been raised for its species.

• in the harvestable zones

Reserve trees are randomly chosen among future trees so that the number of reserve trees is equal to the number of harvested trees.

Value

Input inventory with selected future and reserve trees (LoggingStatus = "future"/"reserve"), for your logging criteria (SpeciesCriteria)

See Also

Paracou6_2016, SpeciesCriteria, loggingparameters

Examples

data(Paracou6_2016)
data(ParamCrownDiameterAllometry)
data(ForestZoneVolumeParametersTable) # The volume parameters data in the global environment
data(SpeciesCriteria)
data(DTMParacou)
data(MainTrails)
data(HarvestableAreaOutputsCable)

scenario <- "RIL2"

inventory <- addtreedim(cleaninventory(Paracou6_2016, PlotMask),
volumeparameters = ForestZoneVolumeParametersTable)

inventory <- commercialcriteriajoin(inventory, SpeciesCriteria)

harvestableOutputs <- harvestable(inventory,
 topography = DTMParacou,
 diversification = TRUE,
 specieslax = FALSE,
 plotslope = HarvestableAreaOutputsCable$PlotSlope,
 maintrails = MainTrails,
 harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons,
 scenario = scenario,
 advancedloggingparameters = loggingparameters())

inventory <- harvestableOutputs$inventory
HVinit <- harvestableOutputs$HVinit

inventory <- selected(inventory, scenario = "manual", fuel = "0",
diversification = TRUE, specieslax = FALSE, objectivelax = TRUE,
topography = DTMParacou,
advancedloggingparameters = loggingparameters(), VO = 125,
 HVinit = HVinit)$inventory

futurereserve(inventory, SpeciesCriteria)

VincyaneBadouard/LoggingLab documentation built on Oct. 16, 2024, 9:42 p.m.