From 51cec653b6ed71748d3a8eef0e24c6907d483c9f Mon Sep 17 00:00:00 2001 From: nicsoattc <nicolas.cattaneo@nibio.no> Date: Fri, 1 Sep 2023 10:00:41 +0200 Subject: [PATCH] old --- CodeR/01_test_data_structure.r | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) diff --git a/CodeR/01_test_data_structure.r b/CodeR/01_test_data_structure.r index 7b82fc4..b7a5cd3 100644 --- a/CodeR/01_test_data_structure.r +++ b/CodeR/01_test_data_structure.r @@ -70,8 +70,8 @@ CornA2 <- copy(CornA) CornA2[, X := X - CornA$X] CornA2[, Y := Y - CornA$Y] -## plot3d(Data2, aspect = "iso", size = 0.2, col= "red") -## spheres3d(c(0, 0, 40), radius = 1, col = "yellow") +plot3d(Data2, aspect = "iso", size = 0.2, col= "red") +spheres3d(c(0, 0, 40), radius = 1, col = "yellow") ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## compute angles on the plane. @@ -96,8 +96,8 @@ Data2 <- rbind(Data3, Data2[angle_x %inrange% list(270.1, 359.9), ]) rm(Data3) gc() -## plot3d(Data2, aspect = "iso", size = 0.2, col= "grey21") -## spheres3d(c(0, 0, 40), radius = 1, col = "yellow") +plot3d(Data2, aspect = "iso", size = 0.2, col= "grey21") +spheres3d(c(0, 0, 40), radius = 1, col = "yellow") ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## Compute angles around the Z axis @@ -112,16 +112,16 @@ Data2[, angle_x := VoxR::axis_angle(Data2, axis = "X", ## FOV -VerticalOpen <- 23 +VerticalOpen <- 30 ## VerticalOpen <- 10 -OrizontalOpen <- 23 +OrizontalOpen <- 30 ## VertcalAngle <- 50 HorizAngle <- 90 plot3d(Data2, aspect = "iso", size = 0.2, - col= "grey21") + col= "grey71") spheres3d(c(0, 0, 40), radius = 1, col = "yellow") plot3d(Data2[angle_z < (90-VertcalAngle) + VerticalOpen & @@ -129,7 +129,7 @@ plot3d(Data2[angle_z < (90-VertcalAngle) + VerticalOpen & angle_x < HorizAngle + OrizontalOpen & angle_x > HorizAngle, ], aspect = "iso", size = 0.7, - col= "red", add = T) + col= "gry21", add = T) pan3d(2) ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ -- GitLab