aRchi (version 2.1.3)

add_radius: Add the radius to a skeleton based on point distance to the skeleton

Description

Add the radius to a skeleton based on point distance to the skeleton

Usage

add_radius(aRchi, sec_length = 0.5, by_axis = TRUE, method = "median")
# S4 method for aRchi
add_radius(aRchi, sec_length = 0.5, by_axis = TRUE, method = "median")

Value

The aRchi file with an updated radius in the QSM slot.

Arguments

aRchi

an object of class aRchi containing at least a point cloud and a QSM.

sec_length

numeric. The length of the section to compute the radius. See details.

by_axis

logical. If TRUE the radius is calculated for each section within each axis, if FALSE the radius is calculated within section for all axes simultaneously.

method

character. The axis radius can be either the mean (method = "mean") or the median (method = "median") distance to the skeleton.

Details

The point distance to skeleton is likely to vary considerably over short distances (e.g. at a branching point). Therefore, the radius is computed by averaging the point distance to the skeleton over sections of user defined length.

Examples

# \donttest{
# import a point cloud
tls=system.file("extdata","Tree_2_point_cloud.las",package = "aRchi")
tls = lidR::readLAS(tls)

aRchi = aRchi::build_aRchi()
aRchi = aRchi::add_pointcloud(aRchi,point_cloud = tls)

# build a skeleton from the point cloud
aRchi = aRchi::skeletonize_pc(aRchi)

# smooth the skeleton
aRchi = aRchi::smooth_skeleton(aRchi)

# add the diameter to the skeleton
aRchi = aRchi::add_radius(aRchi)

# plot the QSM
plot(aRchi,skeleton = FALSE,show_point_cloud = FALSE)
# }