Dear all,
I have several sets of x-y-z points and I need to estimate the volume that
encompass all my points.
Recently I got some adivice to show the "convex hull" of my points using
geometry package (see code below).
But now I need to calculate the volume of my set of points.
Any advice are wellcome.
Thanks in advance
Miltinho Astronauta
Brazil
-----
require(geometry)
library(rgl)
x<-rnorm(100)
k<-runif(100)*20+5
y<-rnorm(100)
z<-sin(x)+cos(y)^2
my.data.3d<-data.frame(cbind(x,y,z))
my.data.3d = sqrt(3) * my.data.3d / drop(sqrt((my.data.3d^2) %*% rep(1,3)))
data.ts.surf = t( convhulln(my.data.3d,"QJ") ) # see the qhull
documentations for the options
## Not run:
rgl.triangles(my.data.3d[data.ts.surf,1],my.data.3d[data.ts.surf,2],my.data.3d[data.ts.surf,3],col="blue",alpha=.2)
for(i in 1:(8*360)) rgl.viewpoint(i/8)
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
