update repo

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)
 ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@