# NOT RUN {
#Create points with paths
n <- 20
orig <- data.table(X = runif(n, min = -5, max = 5),
Y = runif(n, min = -5, max = 5),
Z = runif(n, min = -5, max = 5))
end <- data.table(X = runif(n, min = -5, max = 5),
Y = runif(n, min = -5, max = 5),
Z = runif(n, min = -5, max = 5))
#Create a potential AABB
AABBs <- data.table(X = 0, Y = 0, Z = 0)
edge_length <- c(2, 2, 2)
#Plot
# }
# NOT RUN {
cube <- rgl::cube3d()
cube <- rgl::scale3d(cube,
edge_length[1]/2,
edge_length[2]/2,
edge_length[3]/2)
box <- rgl::translate3d(cube, AABBs[[1]], AABBs[[2]], AABBs[[3]])
rgl::shade3d(box, col= "green", alpha = 0.6)
rgl::points3d(orig, size = 5, col = "black")
rgl::points3d(end, size = 5, col = "red")
for(i in 1:nrow(orig)) {
rgl::lines3d(c(orig[[1]][i], end[[1]][i]),
c(orig[[2]][i], end[[2]][i]),
c(orig[[3]][i], end[[3]][i]), col = "grey")
}
# }
# NOT RUN {
#Estimation
lines_interception(orig, end, AABBs, edge_length, progress = FALSE)
# }
Run the code above in your browser using DataLab