Add as.mesh3d.boundingbox

* permits fancier bounding box plotting with rgl

NAMESPACE

@@ -34,6 +34,7 @@ S3method(as.im3d,im3d)
 S3method(as.im3d,matrix)
 S3method(as.list,neuronlistfh)
 S3method(as.mesh3d,ashape3d)
+S3method(as.mesh3d,boundingbox)
 S3method(as.mesh3d,hxsurf)
 S3method(as.neuron,data.frame)
 S3method(as.neuron,default)

R/hxsurf.R

@@ -291,6 +291,27 @@ as.mesh3d.hxsurf<-function(x, Regions=NULL, material=NULL, drop=TRUE, ...){
   tmesh3d(vertices=verts, indices=inds, homogeneous = FALSE, material = material, ...)
 }
 
+#' @description \code{as.mesh3d.boundingbox} converts a nat
+#'   \code{\link{boundingbox}} object into an rgl compatible \code{mesh3d}
+#'   object.
+#' @rdname as.mesh3d
+#' @export
+#' @examples 
+#' bb=boundingbox(kcs20)
+#' mbb=as.mesh3d(bb)
+#' \donttest{
+#' plot3d(kcs20)
+#' # simple plot
+#' plot3d(bb)
+#' shade3d(mbb, col='red', alpha=0.3)
+#' 
+#' }
+as.mesh3d.boundingbox <- function(x, ...) {
+  centroid=colMeans(x)
+  size=diff(x)/2
+  mat=scaleMatrix(size[1], size[2], size[3])%*%translationMatrix(centroid[1], centroid[2], centroid[3])
+  cube3d(mat)
+}
 
 #' Convert an object to a nat hxsurf object
 #' 

tests/testthat/test-hxsurf.R

@@ -104,7 +104,13 @@ test_that("we can convert hxsurf to rgl::mesh3d",{
   expect_equal(as.hxsurf(tet.mesh3d, region = 'Inside'), tet.hxsurf)
 
   expect_equal(as.mesh3d(surf, Regions=c("LH_L","LH_R")),
-               as.mesh3d(subset(surf, c("LH_L","LH_R"), drop=TRUE)))  
+               as.mesh3d(subset(surf, c("LH_L","LH_R"), drop=TRUE)))
+})
+
+test_that("we can convert boundingbox to rgl::mesh3d",{
+  corners=rbind(rep(1,3), rep(-1,3))
+  bb=boundingbox(corners)
+  expect_equal(as.mesh3d(bb), rgl::cube3d())
 })
 
 test_that("we can convert ashape3d to rgl::mesh3d",{