Thank you all for the solutions.. I now have 3 ways of plotting a 'perfect' sphere :)
1> rgl 2> scatterplot3d 3> persp (with polygon) as below.
Personally I like the transformation solution as in 3. Thanks a mill.
Derick Schoonbee
From: David Brahm <[EMAIL PROTECTED]> Reply-To: [EMAIL PROTECTED] To: [EMAIL PROTECTED] CC: [EMAIL PROTECTED] Subject: Re: [R] Plot a sphere Date: Fri, 26 Dec 2003 10:29:08 -0500
Derick Schoonbee <[EMAIL PROTECTED]> wrote: > Would somebody please be so kind as to direct me in plotting a 3D sphere?
Here's one way. I generate an empty 3D plot with "persp", then fill it with
polygons transformed with "trans3d" (as found in the help for "persp"). I
didn't do hidden surface removal (you didn't mention whether you wanted it),
but if you do, just reorder the polygons from "back" to "front" and paint them
a solid color (e.g. col="red"), so hidden ones get painted over.
pmat <- persp(0:1, 0:1, matrix(,2,2), xlim=c(-1,1), ylim=c(-1,1), zlim=c(-1,1),
theta=25, phi=30, expand=.9, xlab="X", ylab="Y", zlab="Z")
trans3d <- function(x,y,z, pmat) { # From the help for "persp"
tr <- cbind(x,y,z,1) %*% pmat
list(x = tr[,1]/tr[,4], y= tr[,2]/tr[,4])
}
theta <- seq(0, 2*pi, length=51) phi <- seq(0, pi, length=26) x <- cos(theta) %o% sin(phi) y <- sin(theta) %o% sin(phi) z <- rep(1, length(theta)) %o% cos(phi)
for (j in seq(phi)[-1]) for (i in seq(theta)[-1]) {
idx <- rbind(c(i-1,j-1), c(i,j-1), c(i,j), c(i-1,j))
polygon(trans3d(x[idx], y[idx], z[idx], pmat))
}
--
-- David Brahm ([EMAIL PROTECTED])______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
