如何在R中将椭球转换为mesh3d？

Question

我使用cluster包中的ellipsoidhull方法从一组点中获取包含椭球（mvee）的最小体积。此方法返回类椭球的对象。我需要绘制生成的椭圆体。我尝试使用rgl包中的wire3d方法绘制椭球，但此方法将类mesh3d的对象作为输入参数。如何将椭圆体对象转换为mesh3d对象？

Answer 1

如果您实际上并不关心网格，并且只是想要绘制透明椭圆体，您可以使用它：

library(rgl)
library(cluster)
open3d()

ellipsoid3d <- function(cen, a = 1,b = 1,c = 1,n = 65, ...){
  f <- function(s,t){ 
    cbind(   a * cos(t)*cos(s) + cen[1],
             b *        sin(s) + cen[2],
             c * sin(t)*cos(s) + cen[3])
  }
  persp3d(f, slim = c(-pi/2,pi/2), tlim = c(0, 2*pi), n = n, add = T, ...)
}

set.seed(123)
n <- 6
for (i in 1:n){
   cen <- 3*runif(3)
   a <- runif(1)
   b <- runif(1)
   c <- runif(1)

   clr <- c("red","blue","green")[i %% 3 + 1 ]
   elpf <- ellipsoid3d(cen,a=a,b=b,c=c,col=clr,alpha=0.5)
}

产量：

我修改了cuttlefish44的有趣答案来获取此信息 - 请参阅此链接：enter link description here

还有一个来自dww的qmesh3d答案你可以用类似的方式修改以获得一个mesh3d，如果这是你真正想要的，但我认为这更优雅。

Answer 2

library(cluster)
xyz <- cbind(rnorm(10), rnorm(10), rnorm(10))
e <- ellipsoidhull(xyz)
A <- e$cov
center <- e$loc
r <- sqrt(e$d2)

library(Rvcg)
sphr <- vcgSphere()
library(rgl)
ell <- translate3d(
  scale3d(
    transform3d(sphr, chol(A)), 
    r, r, r),
  center[1], center[2], center[3])

shade3d(ell, color="red", alpha=0.3)
points3d(xyz)