secondtrailsopening: secondtrailsopening

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

Description

secondtrailsopening

Usage

secondtrailsopening(
  DTM,
  plots,
  verticalcreekheight,
  treeselectionoutputs,
  CostMatrix = list(list(list(Slope = 3, Cost = 3), list(Slope = 5, Cost = 5),
    list(Slope = 12, Cost = 20), list(Slope = 22, Cost = 60), list(Slope = 27, Cost =
    600), list(Slope = Inf, Cost = 1000)), list(list(CostType = "Initial", CostValue =
    1000), list(CostType = "Access", CostValue = 1000), list(CostType = "BigTrees",
    CostValue = 500), list(CostType = "Reserves", CostValue = 500), list(CostType =
    "Futures", CostValue = 50), list(CostType = "maintrails", CostValue = 1e-04),
    list(CostType = "SecondTrails", CostValue = 0.1))),
  scenarios = NULL,
  fact = 3,
  advancedloggingparameters = loggingparameters()
)

Arguments

DTM	Digital terrain model (raster)
plots	Studied plots (SpatialPolygonsDataFrame)
verticalcreekheight	Relative elevation from nearest channel network (Large RasterLayer)
treeselectionoutputs	A list with: your inventory with: "DistCrit", "Slope", "SlopeCrit", "LoggingStatus", "Selected", "Up", "VolumeCumSum", "ProbedHollowProba", "ProbedHollow" new columns (see the outputs metadata in the vignette). the objective volume with or without a bonus (if hollow trees exploitation) 6 sets of spatial points: harvestable, selected, future and reserve, hollow and energy wood trees
CostMatrix	List of list defining conditional weight over binned slopes values
scenarios	The chosen scenario from scenariosparameters function
fact	Aggregation factor of cost raster resolution to initial DTM
advancedloggingparameters	Advanced parameters of the logging simulator

Value

Secondary trails polygons according to scenario conditions.

Examples

data(DTMParacou)
data(Paracou6_2016)
data(HarvestablePolygons)
data(MainTrails)
data(Plots)
data(PlotSlope)
data("SpeciesCriteria")
data(VerticalCreekHeight)

scenarios <- scenariosparameters(scenario = "RIL3")

inventory <- ONFGuyafortaxojoin(addtreedim(inventorycheckformat(Paracou6_2016),
 volumeparameters = ForestZoneVolumeParametersTable),SpeciesCriteria)

AccessPolygons <- FilterAccesExplArea(harvestablepolygons = HarvestablePolygons,
maintrails = MainTrails,
advancedloggingparameters = loggingparameters())

treeselectionoutputs <- treeselection(inventory,
topography = DTMParacou, plotslope = PlotSlope, maintrails = MainTrails,
harvestablepolygons = HarvestablePolygons,
speciescriteria = SpeciesCriteria, objective = scenarios$objective,
scenario = "RIL3",
objectivelax = TRUE,
advancedloggingparameters = loggingparameters())

secondtrails <- secondtrailsopening(DTM = DTMParacou,
   plots = Plots,

   treeselectionoutputs = treeselectionoutputs,
   verticalcreekheight = VerticalCreekHeight,
   CostMatrix = list(list(list(Slope = 3, Cost = 3),
                           list(Slope = 5, Cost = 5),
                           list(Slope = 12, Cost = 20),
                           list(Slope = 22, Cost = 60),
                           list(Slope = 27, Cost = 600),
                           list(Slope = Inf, Cost = 1000)),
                      list(list(CostType = "Initial", CostValue = 1000),
                           list(CostType = "Access", CostValue = 10000),
                           list(CostType = "BigTrees", CostValue = 500),
                           list(CostType = "Reserves", CostValue = 500),
                           list(CostType = "Futures", CostValue = 50),
                           list(CostType = "maintrails", CostValue = 1E-4),
                           list(CostType = "SecondTrails", CostValue = 0.1))),
 scenarios = scenarios,
 fact = 3,
 advancedloggingparameters = loggingparameters())

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