diff --git a/CodeR/03_Compute_angles_cone.r b/CodeR/03_Compute_angles_cone.r index af3b3f9515088dc6a70ea7673f98878df1dfbc64..e408c7562706c61d5936649588b54da3ea45005b 100644 --- a/CodeR/03_Compute_angles_cone.r +++ b/CodeR/03_Compute_angles_cone.r @@ -15,8 +15,11 @@ library("data.table") ## DataR and CodeR ## DataR <- paste0("/home/nibio/REPOS/holger_tower_pov/DataR") ## DataR <- paste0("/home/pepito/Documents/holger_tower_pov/DataR") -DataR <- paste0("/home/nica/Documents/Holger_tower_pov/DataR") -CodeR <- paste0("/home/nica/Documents/Holger_tower_pov/CodeR") +## DataR <- paste0("/home/nica/Documents/Holger_tower_pov/DataR") +## CodeR <- paste0("/home/nica/Documents/Holger_tower_pov/CodeR") + +DataR <- paste0("/home/nibio/REPOS/holger_tower_pov/DataR") +CodeR <- paste0("/home/nibio/REPOS/holger_tower_pov/CodeR") ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## Functions @@ -25,8 +28,8 @@ source(paste0(CodeR, "/Functions/Functions.r")) ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## plot data - full path DataPaths <- list.files(DataR, - pattern = "plot_", - ## pattern = "Plot_", + ## pattern = "plot_", + pattern = "Plot_", ## pattern = "merg", ## pattern = ".rds", ## pattern = "Tower", @@ -36,7 +39,7 @@ Data <- read.las(DataPaths) Data[, 1:3] ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ -Data <- Data[sample(1:nrow(Data), nrow(Data)*0.20), ] +Data <- Data[sample(1:nrow(Data), nrow(Data)*0.10), ] Data[, X := Scale(Data$X)] Data[, Y := Scale(Data$Y)] Data[, Z := Scale(Data$Z)] @@ -79,12 +82,17 @@ Data[, H := dist3d2(X, Y, Z, Data[, Angle := CO/H] Data[, Angle := (asin(Angle))*180/pi] -plot3d(Data[Angle > 10, ], aspect = "iso", size = 0.5, col = "grey21") +OpeningAngl <- 20 + +plot3d(Data[Angle < OpeningAngl, ], aspect = "iso", size = 1.5, col = "red") spheres3d(CornA, radius = 1, col = "yellow") spheres3d(CornB, radius = 1, col = "red") -plot3d(Data[Angle < 28, ], aspect = "iso", size = 0.5, col = "red", add = T) -pan3d(3) +plot3d(Data[Angle > OpeningAngl, ], aspect = "iso", size = 0.5, col = "grey21") +spheres3d(CornA, radius = 1, col = "yellow") +spheres3d(CornB, radius = 1, col = "red") +plot3d(Data[Angle < OpeningAngl, ], aspect = "iso", size = 0.5, col = "red", add = T) +pan3d(3) ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@