green_crown_volume: Compute the green crown volume, area and return a mesh

In Blecigne/lidUrb: Terrestrial and Mobile Laser Scaning processing

Compute the green crown volume, area and return a mesh

Description

Compute the green crown volume, area and return a mesh

Usage

green_crown_volume(las, npts_in_clust = 200)

Arguments

las	a LAS file with a classification field for the wood class.
npts_in_clust	numeric. The average number of points in a cluster. Defines the size of each cluster and therefore the size of the convex hulls.

Value

a mesh object of the green crown and the green crown volume and area

References

Zhu, Z., Kleinn, C., & Nölke, N. (2020). Towards tree green crown volume: a methodological approach using terrestrial laser scanning. Remote Sensing, 12(11), 1841.

Examples

# import data
file = system.file("extdata", "tree_leaves.las", package="lidUrb")
las = lidR::readLAS(file,select = "xyz")

# filter noise
original = las
original@data[,original_index := 1:nrow(las@data)]

# filter noise and reduce point density
las = lidUrb::reduce_point_density(original,0.02)
las = lidUrb::filter_noise(las)

# segment leaves
las = lidUrb::LW_segmentation_dbscan(las)
las@data[,wood := as.numeric(p_wood >= 0.95)] # binary class (0 for leaves, 1 for wood)

# compute green crown volume
GCV = lidUrb::green_crown_volume(las)

# plot the green crown mesh over the original LAS
lidR::plot(las,color="wood",size=2,colorPalette = c("chartreuse4","cornsilk2"),clear_artifacts = FALSE)
rgl::shade3d(GCV$mesh,col = "chartreuse4",add=T)

# green crown volume and area
GCV$Green_crow_volume
GCV$Green_crown_area

Blecigne/lidUrb documentation built on Feb. 19, 2024, 9:12 a.m.