R/Width_canopy_function.R

In PROTOLIDAR: PRocess TOol LIdar DAta in R

Width_canopy_function <-
function(data_3D,distance_left,distance_right,min_canopy,max_canopy){

  x_plant <- y_plant <- z_plant <- NULL

  canopy <- subset(data_3D, data_3D$z_plant >= distance_left & data_3D$z_plant <= distance_right & data_3D$y_plant >= min_canopy & data_3D$y_plant <= max_canopy,select=c(x_plant,y_plant,z_plant))

  mean_width_canopy <- mean(abs(canopy[,1]))

  min_width_canopy  <- min(abs(canopy[,1]))

  max_width_canopy  <- max(abs(canopy[,1]))

  return(data.frame(mean_width_canopy,min_width_canopy,max_width_canopy))

}