mortalityunderpolygon: Trees removed under polygons

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

Trees removed under polygons

Description

Trees removed under polygons

Usage

mortalityunderpolygon(inventory, polygons, deathcause)

Arguments

inventory	Input inventory (see the inputs formats and metadata in the vignette) (data.frame) with at least the following columns: "idTree": unique ID of the tree "Xutm": X coordinate of the tree "Yutm": Y coordinate of the tree "Selected": == "1" for the trees selected for logging.
polygons	Polygon(s) (trails, landings) under which trees are removed. (POLYGON or MULTIPOLYGON with a crs in UTM)
deathcause	The cause of the tree death under the polygon(s) among: "maintrail", "2ndtrail", "treefall2nd", "landing" (character)

Value

Input inventory (data.frame) with a DeathCause column filled with the given deathcause.

Examples

# inventory
data(Paracou6_2016)
data(DTMParacou)
data(SpeciesCriteria)
data(MainTrails)
data(HarvestableAreaOutputsCable)

Paracou6_2016 <- tibble::add_column(Paracou6_2016, DBH= NA) # add a DBH column
Paracou6_2016$DBH = Paracou6_2016$CircCorr/pi # and calculate it
Paracou6_2016 <- dplyr::filter(Paracou6_2016, DBH >= 10)

inventory <- addtreedim(Paracou6_2016,
                        volumeparameters = ForestZoneVolumeParametersTable)

inventory <- suppressMessages(treeselection(inventory,
  topography = DTMParacou,
  speciescriteria = SpeciesCriteria,
  scenario = "RIL2",
  specieslax = FALSE, objectivelax = TRUE,
  harvestablearea = HarvestableAreaOutputsCable$HarvestableArea,
  plotslope = HarvestableAreaOutputsCable$PlotSlope,
  maintrails = MainTrails,
  harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons
  )$inventory)

# polygon
data(SecondaryTrails)
polygons <- SecondaryTrails$SmoothedTrails
deathcause <- "2ndtrail"

Rslt <- mortalityunderpolygon(inventory,
                              polygons,
                              deathcause)
library(ggplot2)
library(sf)
inventory <- st_set_crs(st_as_sf(Rslt, coords = c("Xutm", "Yutm")),
                                 st_crs(polygons))

DeadTrees <- sf::st_as_sf(
dplyr::filter(inventory, DeathCause == deathcause),
              coords = c("Xutm", "Yutm"))
ggplot() +
  geom_sf(data = inventory) +
  geom_sf(data = DeadTrees,
  aes(colour = "DeadTrees"), show.legend = "point") +
  geom_sf(data = polygons,
  alpha = 0.5, fill = "red") +
  scale_colour_manual(values = c("DeadTrees" = "red")) +
  labs(color = "Status")

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