FORTLS (version 1.0.1)

estimation.plot.size: Assess Consistency of Metrics for Simulated TLS Plots

Description

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).

Usage

estimation.plot.size(tree.list.tls,
                     plot.parameters = list(radius.max = 25,
                                            k.tree.max = 50,
                                            BAF.max = 4),
                     average = FALSE, all.plot.designs = FALSE)

Arguments

tree.list.tls

Data frame with information of trees detected from TLS point cloud data in the same format as tree.detection ‘Value’.

plot.parameters

Optional list containing parameters for circular fixed area, k-tree and angle-count plot designs. The parameters are as follows:

average

Logical; if TRUE, plot means values and standard deviation of estimations will be represented. By default, it will be set as FALSE.

all.plot.designs

Logical; if TRUE, charts for each plot design are drawn together. By default, it will be set as FALSE.

Value

Invisible NULL

Details

If there are strata in the tree.list.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.list.tls argument, they will be considered equal.

The outputs of this function are inspired by Fig. 3 of Brunner and Gizachew (2014).

References

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.

See Also

tree.detection, tree.detection.multiple

Examples

# NOT RUN {
# Loading dataset with trees detected from TLS single-scans

data("Rioja.data")

tree.list.tls <- Rioja.data$tree.list.tls


# Without strata and plot parameters by default

estimation.plot.size(tree.list.tls)

estimation.plot.size(tree.list.tls, average = TRUE)

estimation.plot.size(tree.list.tls, all.plot.designs = TRUE)


# Considering two strata and different parameters

# Loading dataset with strata

plot.attributes <- Rioja.data$plot.attributes


# Merging the plot.attributes data set with strata information

tree.list.tls <- merge(tree.list.tls, plot.attributes, by = "id")


estimation.plot.size(tree.list.tls,
                     plot.parameters = list(radius.max = 10, k.tree.max = 10, BAF.max = 2))

estimation.plot.size(tree.list.tls, ,
                     plot.parameters = list(radius.max = 15, k.tree.max = 20, BAF.max = 4),
                     average = TRUE)

estimation.plot.size(tree.list.tls,
                     plot.parameters = list(radius.max = 20, k.tree.max = 30, BAF.max = 2),
                     all.plot.designs = TRUE)


# }