examples/example_apa_seg.R

In JonasGlatthorn/APAtree: Computation of the 'Area Potentially Available' (APA) to Trees

library(APAtree)
data(tree_enrico, package = "APAtree")
data(plot_enrico, package = "APAtree")

tree_enrico$height_class <- tree_enrico$height > 20

# only calculating an apa_list for two plots and with a coarse resolution of 1 m
# to save time.
apa_list_enrico <- 
  apa_list(plot_dat = subset(plot_enrico, id_plot %in% c("5.2", "8.2")), 
           tree_dat = tree_enrico,
           buffer_column = "buffer_geometry",
           core_column = "border_geometry",
           plot_id_column = "id_plot",
           tree_id_column = "id_tree",
           weight_column = "crown_radius_95",
           agg_class_column = c("species", "height_class"),
           res = 10,
           apa_polygon = FALSE)

apa_list_enrico <- 
  apa_ndiv(apa_list_enrico, dis_trait_column = "species")

apa_list_enrico <-
  apa_seg(apa_list_enrico, nsim = 3, parallel = FALSE)

# stand-level segregation index of SpeciesNDiv
apa_list_enrico$plot_dat[, c("id_plot", "seg_species_ndiv")]