diff --git a/CodeR/01_test_data_structure.r b/CodeR/01_test_data_structure.r index 7b82fc426ef2a7d02297396396bbf148945a3447..b7a5cd37c1ea0534f8bc0f69b0b2d0c58fc4f4c9 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) ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@