Width_canopy_function: Canopy width measured with LIDAR. In PROTOLIDAR: PRocess TOol LIdar DAta in R

Description

From the LIDAR dataset can be calculate the width of the grapevine plant.The function returns the average, minimum and maximum value of the width measured in meters.

Usage

 `1` ```Width_canopy_function(data_3D, distance_left, distance_right, min_canopy, max_canopy) ```

Arguments

 `data_3D` Here need use the output of the Extract_plant_3D_function `distance_left` the left distance of the plant, generally here we can write the half of the distance between plants, but is better the real distance from the center of the plant. Measured in meters. `distance_right` the right distance of the plant, generally here we can write the half of the distance between plants, but is better the real distance from the center of the plant. Measured in meters. `min_canopy` the minimum height of the canopy, measured in meters. `max_canopy` the maximum height of the canopy, measured in meters.

Details

Maximum and minimum values of height of canopy could be approximative values.

Author(s)

Monica Fernanda Rinaldi

Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42``` ```## Should be DIRECTLY executable !! ---- ## Data_3D is the output from the Extrac_plant_3D_function. ## For example: data (LIDAR_data) x <- LIDAR_data[,1] y <- LIDAR_data[,2] z <- LIDAR_data[,3] zdistance <- 190 # total LIDAR scan distance measured in cm. miny <- 0 # minimum height of the plant measured in cm. maxy <- 2000 # maximum height of the plant measured in cm. minx <- 450 # minimum width from where LIDAR starts to measure (cm). maxx <- 1470# maximum width from where LIDAR starts to measure (cm). minz <- 0 # the beginning of the LIDAR scan measured in cm. maxz <- 186 # the end of the LIDAR scan measured in cm (length of interest). out <- Extract_plant_grapevine_function(x,y,z,zdistance,miny,maxy,minx,maxx,minz,maxz) z_min <- 1.1 z_max <- 1.13 y_min <- 0.4 y_max <- 0.5 distance_left <- -0.6 distance_right <- 0.51 data_3D <- Extract_plant_3D_function(out,z_min,z_max,y_min,y_max,distance_left,distance_right) min_canopy <- 0.4 # is the minimum height of the canopy, approximately . Measured in meters. max_canopy <- 2 # is the maximum height of the canopy, approximately . Measured in meters. ## The function is currently defined as 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)) } width_canopy <- Width_canopy_function(data_3D,distance_left,distance_right,min_canopy,max_canopy)[,1] ```

PROTOLIDAR documentation built on May 29, 2017, 7:12 p.m.