Plots empirical linear charts of density (N, trees/ha) and basal area (G, \({m}^{2}/ha\)) estimates (derived from simulated TLS plots) as a function of plot size (estimation-size charts) for different plot designs (circular fixed area, k-tree and angle-count), through continuous size increments (radius, k and BAF respectively). Size increments are set at 0.1 m, 1 tree and 0.1 \({m}^{2}/ha\) for fixed area, k-tree and angle-count plot designs, respectively. These size-estimation line charts represent the consistency in predicting the stand variables across different values of radius, k and BAF. Size-estimation charts can be drawn for individual sample plots (including all plots together in the same charts) or for mean values (global mean computed for all the sample plots, or for group means if different strata are considered). Finally, different plot designs can be compared if specified in the arguments, producing one size-estimation chart per variable (N and G).
estimation.plot.size(tree.tls,
plot.parameters = data.frame(radius.max = 25,
k.max = 50,
BAF.max = 4),
dbh.min = 4,
average = FALSE, all.plot.designs = FALSE)
Invisible NULL
Data frame with information of trees detected from TLS point cloud data in the same format as tree.detection.single.scan
and tree.detection.multi.scan
‘Value’.
Optional data frame containing parameters for circular fixed area, k-tree and angle-count plot designs. The parameters are as follows:
radius.max
: maximum plot radius (m) considered for circular fixed area plots. If the radius.max
specified is larger than the farthest tree from the plot centre, the horizontal distance from the farthest tree will be considered the maximum radius
. By default, the radius.max
will be 25 m.
k.tree.max
: maximum number of trees considered for k-tree plots. If k.tree.max
specified is larger than the maximum number of trees of the densest plot, this number of trees will be considered the maximum k.tree.max
. By default, k.tree.max
is 50.
BAF.max
: maximum basal area factor (\({m}^{2}/ha\)) considered for angle-count plots. By default, BAF.max
is 4.
Optional minimum dbh (cm) considered for detecting trees. By default it will be set at 4 cm.
Logical; if TRUE
, plot means values and standard deviation of estimations will be represented. By default, it will be set as FALSE
.
Logical; if TRUE
, charts for each plot design are drawn together. By default, it will be set as FALSE
.
Juan Alberto Molina-Valero and Adela Martínez-Calvo.
If there are strata in the tree.tls
argument, they will be differentiated in charts with different colours. Strata must be specified in a numeric column named stratum
.
The all.plot.designs
argument only works for single strata, and therefore if there are additional strata in the tree.tls
argument, they will be considered equal.
The outputs of this function are inspired by Fig. 3 of Brunner and Gizachew (2014).
Brunner, A., & Gizachew, B. (2014). Rapid detection of stand density, tree positions, and tree diameter with a 2D terrestrial laser scanner. European Journal of Forest Research, 133(5), 819-831.
tree.detection.single.scan
, tree.detection.multi.scan
, tree.detection.several.plots
# \donttest{
# Loading dataset with trees detected from TLS single-scans
data("Rioja.data")
tree.tls <- Rioja.data$tree.tls
# Without strata and plot parameters by default
estimation.plot.size(tree.tls)
estimation.plot.size(tree.tls, average = TRUE)
estimation.plot.size(tree.tls, all.plot.designs = TRUE)
# }
Run the code above in your browser using DataLab