The following examples demonstrates how to use the tree_metrics
function of rTLS
package to estimate: i) tree height, ii) crown area, and iii) diameter at breast height (DBH).
First, we need to read the point clouds in R. This can be done using the function fread
of data.table
for a fast reading or any other traditional approach for reading tables.
For this example specifically, we will use a point cloud already embedded in rTLS
, pc_tree
. This point cloud was created from a group of scans of a tree. It has a coarse resolution of 0.05 m for TLS standards in order to integrate it as an example into the package. The file can be load following:
library(rTLS)
data(pc_tree)
#Visualize the point cloud
::plot3d(pc_tree) rgl
tree_metrics
The next step is run tree_metrics
. This can be done following:
tree_metrics(pc_tree)
In general, the tree height is estimated based on the range of the Z axis, the crown area is calculated applying a convex hull on the point cloud, and the DBH is calculated extracting the area of the convex hull on the subset of points between a given region, and then estimating the diameter of a circle. Since the estimation of the DBH is applied to a subset of points between 1.25 and 1.35 m, it tends to be affected by presence of other points that do not below to the main trunk. The selection of the region for DBH could be further modified using region.diameter
.