height_regression {lidaRtRee}

Regression of detected heights VS reference heights

Description

Computes a linear regression model between the reference heights and the detected heights of matched pairs.

Usage

height_regression(lr, ld, matched, plot = TRUE, species = NULL, ...)

Arguments

lr

data.frame or matrix. 3D coordinates (X Y Height) of reference positions

ld

data.frame or matrix. 3D coordinates (X Y Height) of detected positions

matched

data.frame. contains pair indices, typically returned by tree_matching

plot

boolean. indicates whether results should be plotted

species

vector of strings. species for standardized color use by call to species_color

...

arguments to be passed to methods, as in plot

Value

A list with two elements. First one is the linear regression model, second one is a list with stats (root mean square error, bias and standard deviation of detected heights compared to reference heights).

See Also

tree_matching

Examples

# create tree locations and heights
ref_trees <- cbind(c(1, 4, 3, 4, 2), c(1, 1, 2, 3, 4), c(15, 18, 20, 10, 11))
def_trees <- cbind(c(2, 2, 4, 4), c(1, 3, 4, 1), c(16, 19, 9, 15))

# tree matching
match1 <- tree_matching(ref_trees, def_trees)

# height regression
reg <- height_regression(ref_trees, def_trees, match1,
  species = c("ABAL", "ABAL", "FASY", "FASY", "ABAL"),
  asp = 1, xlim = c(0, 21), ylim = c(0, 21)
)
summary(reg$lm)
reg$stats