Skip to content
Snippets Groups Projects
Commit 51cec653 authored by Nicolas Cattaneo's avatar Nicolas Cattaneo
Browse files

old

parent f8a859b1
No related branches found
No related tags found
No related merge requests found
...@@ -70,8 +70,8 @@ CornA2 <- copy(CornA) ...@@ -70,8 +70,8 @@ CornA2 <- copy(CornA)
CornA2[, X := X - CornA$X] CornA2[, X := X - CornA$X]
CornA2[, Y := Y - CornA$Y] CornA2[, Y := Y - CornA$Y]
## plot3d(Data2, aspect = "iso", size = 0.2, col= "red") plot3d(Data2, aspect = "iso", size = 0.2, col= "red")
## spheres3d(c(0, 0, 40), radius = 1, col = "yellow") spheres3d(c(0, 0, 40), radius = 1, col = "yellow")
## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
## compute angles on the plane. ## compute angles on the plane.
...@@ -96,8 +96,8 @@ Data2 <- rbind(Data3, Data2[angle_x %inrange% list(270.1, 359.9), ]) ...@@ -96,8 +96,8 @@ Data2 <- rbind(Data3, Data2[angle_x %inrange% list(270.1, 359.9), ])
rm(Data3) rm(Data3)
gc() gc()
## plot3d(Data2, aspect = "iso", size = 0.2, col= "grey21") plot3d(Data2, aspect = "iso", size = 0.2, col= "grey21")
## spheres3d(c(0, 0, 40), radius = 1, col = "yellow") spheres3d(c(0, 0, 40), radius = 1, col = "yellow")
## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
## Compute angles around the Z axis ## Compute angles around the Z axis
...@@ -112,16 +112,16 @@ Data2[, angle_x := VoxR::axis_angle(Data2, axis = "X", ...@@ -112,16 +112,16 @@ Data2[, angle_x := VoxR::axis_angle(Data2, axis = "X",
## FOV ## FOV
VerticalOpen <- 23 VerticalOpen <- 30
## VerticalOpen <- 10 ## VerticalOpen <- 10
OrizontalOpen <- 23 OrizontalOpen <- 30
## ##
VertcalAngle <- 50 VertcalAngle <- 50
HorizAngle <- 90 HorizAngle <- 90
plot3d(Data2, aspect = "iso", size = 0.2, plot3d(Data2, aspect = "iso", size = 0.2,
col= "grey21") col= "grey71")
spheres3d(c(0, 0, 40), radius = 1, col = "yellow") spheres3d(c(0, 0, 40), radius = 1, col = "yellow")
plot3d(Data2[angle_z < (90-VertcalAngle) + VerticalOpen & plot3d(Data2[angle_z < (90-VertcalAngle) + VerticalOpen &
...@@ -129,7 +129,7 @@ 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 + OrizontalOpen &
angle_x > HorizAngle, ], angle_x > HorizAngle, ],
aspect = "iso", size = 0.7, aspect = "iso", size = 0.7,
col= "red", add = T) col= "gry21", add = T)
pan3d(2) pan3d(2)
## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ## @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment