R绘图系统16.doc

16.1library(rgl)tetra - function() t1 - tetrahedron3d()t2vb - t1$vbt2vb1, - -3t2 - tmesh3d(t2vb, t1$it)plane - qmesh3d(rbind(rep(-3.01, 4), c(-2, -2, 2, 2), c(-3, 3, 3, -3), rep(1, 4), matrix(1:4, ncol=1)open3d(windowRect=c(0, 0, 600, 600)# clear3d()shade3d(plane, color=white, specular=black)wire3d(plane)wire3d(t1, lwd=3)wire3d(t2, lwd=3)segments3d(rbind(t2$vb1, t2$it, t1$vb1, t1$it), rbind(t2$vb2, t2$it, t1$vb2, t1$it), rbind(t2$vb3, t2$it, t1$vb3, t1$it), col=gray, lwd=3)view3d(40, -30)t1 - cube3d()t1tube - t1t1tube$ib - t1tube$ib, -(3:4)t2vb - t1$vbt2vb1, - -5t2 - qmesh3d(t2vb, t1$ib)plane - qmesh3d(rbind(rep(-5.01, 4), c(-2, -2, 2, 2), c(-3, 3, 3, -3), rep(1, 4), matrix(1:4, ncol=1)open3d(windowRect=c(0, 0, 600, 600)# clear3d()shade3d(plane, color=white, ambient=white, specular=white, emission=white)wire3d(plane)shade3d(t1tube, color=white, specular=black)wire3d(t1, lwd=3)wire3d(t2, lwd=3)segments3d(rbind(t2$vb1, t2$ib, 4, t1$vb1, t1$ib, 4), rbind(t2$vb2, t2$ib, 4, t1$vb2, t1$ib, 4), rbind(t2$vb3, t2$ib, 4, t1$vb3, t1$ib, 4), col=gray, lwd=3)view3d(55, -20, fov=0)rgl.postscript(Figures/threed-3dproj.eps)system(epstopdf Figures/threed-3dproj.eps)system(convert Figures/threed-3dproj.pdf Web/threed-3dproj.png)16.2library(rgl)tetra - function() t2vb - t1$vbt2vb1, - -3t2vb2, c(1, 4) - t2vb2, c(1, 4)*.8t2vb3, c(1, 4) - t2vb3, c(1, 4)*.8t2vb2, 2:3 - t2vb2, 2:3*.6t2vb3, 2:3 - t2vb3, 2:3*.6t2 - tmesh3d(t2vb, t1$it)t3vb - t1$vbt3vb1, - -10t3vb2, - 0t3vb3, - 0t3 - tmesh3d(t3vb, t1$it)open3d(windowRect=c(0, 0, 600, 600)# clear3d()shade3d(plane, color=white, specular=black)wire3d(plane)wire3d(t1, lwd=3)wire3d(t2, lwd=3)segments3d(rbind(t3$vb1, t3$it, t1$vb1, t1$it), rbind(t3$vb2, t3$it, t1$vb2, t1$it), rbind(t3$vb3, t3$it, t1$vb3, t1$it), col=gray, lwd=3)shade3d(translate3d(scale3d(cube3d(), .1, .1, .1), -10, 0, 0)view3d(50, -20, zoom=.8)t1 - cube3d()t1tube - t1t1tube$ib - t1tube$ib, -(3:4)t2vb - t1$vbt2vb2, t1$ib, 4 - t2vb2, t1$ib, 4*.4t2vb3, t1$ib, 4 - t2vb3, t1$ib, 4*.4t2vb2, t1$ib, 3 - t2vb2, t1$ib, 3*.6t2vb3, t1$ib, 3 - t2vb3, t1$ib, 3*.6t2vb1, - -4.5t2 - qmesh3d(t2vb, t1$ib)t2tube - t2t2tube$ib - t2tube$ib, -(3:4)t3 - translate3d(scale3d(cube3d(), .1, .1, .1), -10, 0, 0)plane - qmesh3d(rbind(rep(-4.51, 4), c(-2, -2, 2, 2), c(-3, 3, 3, -3), rep(1, 4), matrix(1:4, ncol=1)open3d(windowRect=c(0, 0, 600, 300)# clear3d()shade3d(plane, color=white, ambient=white, specular=white, emission=white)wire3d(plane)shade3d(t1tube, color=white, specular=black)wire3d(t1, lwd=3)shade3d(t2tube, color=white, specular=black)wire3d(t2, lwd=3)segments3d(rbind(-10, t1$vb1, t1$ib, 4), rbind(0, t1$vb2, t1$ib, 4), rbind(0, t1$vb3, t1$ib, 4), col=gray, lwd=3)shade3d(t3)view3d(50, -15, fov=0, zoom=.7)rgl.postscript(Figures/threed-3dvp.eps)system(epstopdf Figures/threed-3dvp.eps)system(convert Figures/threed-3dvp.pdf Web/threed-3dvp.png)16.3quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)par(mar=rep(0, 4)persp(quakeDens)16.4quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)par(mar=rep(0, 4)persp(quakeDens, scale=FALSE, expand=0.02, theta=60, d=.1, r=.1, xlab=longitude, ylab=latitude, zlab=)16.5quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)par(mar=rep(0, 4)zinterp - with(quakeDens, z-1, -1 + z-1, -ncol(z) + z-nrow(z), -1 + z-nrow(z), -ncol(z)persp(quakeDens, scale=FALSE, expand=0.02, theta=60, d=.1, r=.1, axes=FALSE, box=FALSE, col=gray(.4 + 1:6/10)cut(zinterp, 6)16.6quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)shallowCantyQuakes - subset(cantyQuakes, DEPTH 20)trellis.device(pdf, color=FALSE, file=Figures/threed-cloud%d.pdf, onefile=FALSE)for (i in seq(40, 80, 20) print(cloud(-DEPTH LONG + LAT, shallowCantyQuakes, xlim=c(171, 173), ylim=c(-44.5, -42.5), pch=16, col=rgb(0, 0, 0, .5), screen=list(z=i, x=-70)dev.off()16.7quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)shallowCantyQuakes - subset(cantyQuakes, DEPTH 20)library(scatterplot3d)par(lab=c(3, 3, 0)s3d - with(shallowCantyQuakes, scatterplot3d(-DEPTH LONG + LAT, angle=30, scale.y=0.45, type=n, pch=16, color=rgb(0, 0, 0, .5), x.ticklabs=pretty(LONG, 3), grid=FALSE, zlim=c(-20, 0)quakeDensXY - kde2d(shallowCantyQuakes$LONG, shallowCantyQuakes$LAT, n=30)lapply(contourLines(quakeDensXY, nlevels=8), function(cl) polygon(s3d$xyz.convert(cl$x, cl$y, rep(-20, length(cl$x), lwd=.5, col=gray(.8 - cl$level/20), border=NA) )quakeDensXZ - kde2d(shallowCantyQuakes$LONG, -shallowCantyQuakes$DEPTH, n=30)lapply(contourLines(quakeDensXZ, nlevels=8), function(cl) polygon(s3d$xyz.convert(cl$x, rep(-43.2, length(cl$x), cl$y), lwd=.5, col=gray(.8 - 2*cl$level), border=NA) )quakeDensYZ - kde2d(shallowCantyQuakes$LAT, -shallowCantyQuakes$DEPTH, n=30)lapply(contourLines(quakeDensYZ, nlevels=8), function(cl) polygon(s3d$xyz.convert(rep(171.5, length(cl$x), cl$x, cl$y), lwd=.5, col=gray(.8 - cl$level), border=NA) )lapply(contourLines(quakeDensYZ, nlevels=8), function(cl) polygon(s3d$xyz.convert(rep(173.5, length(cl$x), cl$x, cl$y), lwd=.5, col=gray(.8 - cl$level), border=NA) )with(shallowCantyQuakes, s3d$points3d(-DEPTH LONG + LAT, pch=16, col=rgb(0, 0, 0, .3)s3d$box()16.8library(rgl)quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)open3d(windowRect=c(0, 0, 900, 450)# clear3d(all)persp3d(quakeDens$x, quakeDens$y, quakeDens$z, aspect=c(1, 0.55, .2), col=white, box=FALSE, axes=FALSE, xlab=, ylab=, zlab=)par3d(userMatrix=rotationMatrix(-80/180*pi, 1, 0, 0)%*% rotationMatrix(-65/180*pi, 0, 0, 1), zoom=.5)snapshot3d(Figures/threed-rglpersp.png)system(cp Figures/threed-rglpersp.png Web/)16.9library(rgl)library(misc3d)quakes - read.csv(Quake/earthquakes.csv)NZquakes - quakesc(LAT, LONG, MAG, DEPTH)cantyQuakes - quakesquakes$LAT -44 & quakes$LONG 171 & quakes$LONG 173.5, library(MASS)quakeDens - kde2d(cantyQuakes$LONG, cantyQuakes$LAT, n=30)shallowCantyQuakes - subset(cantyQuakes, DEPTH 20)d - with(shallowCantyQuakes, kde3d(LONG, LAT, -DEPTH, h=c(.1, .1, 2), n = 30) )open3d(windowRect=c(0, 0, 900, 900)with(shallowCantyQuakes, plot3d(LONG, LAT, -DEPTH, aspect=c(1, 0.55, 1), axes=TRUE, box=FALSE, xlab=, ylab=, zlab=) contour3d(d$d, c(.4, .1), d$x, d$y, d$z, color=rep(gray80, 2), color2=gray, specular=black, engine=rgl, add=TRUE, alpha=.5) )par3d(userMatrix=rotationMatrix(-60/180*pi, 1, 0, 0)%*% rotationMatrix(-40/180*pi, 0, 0, 1), zoom=.9)snapshot3d(Figures/contour3d.png)system(cp Figures/contour3d.png Web/threed-contour3d.png)16.10library(rgl)open3d(windowRect=c(0, 0, 600, 600)clear3d(all)light3d()material3d(shininess=100, specular=black)# Headradius - function(d) pchisq(d2, 3)shade3d(ellipse3d(diag(3), level=radius(1), centre=c(0, 0, 1), color=yellow)# Neck# logo is 100x76png(rlogoExtended.png, width=500, height=250)library(grid)grid.rect(gp=gpar(col=NA, fill=yellow)library(png)rlogo - as.raster(readPNG(system.file(extra, Rlogo.png, package=RGraphics)rlogorlogo = #FFFFFF - yellowgrid.raster(rlogo, x=.6, y=.01, width=.08, just=c(bottom)dev.off()shade3d(cylinder3d(cbind(0, 0, c(-1.4, 1), e1=cbind(0, 0, 1), e2=cbind(1, 0, 0), sides=100), color=yellow, texture=rlogoExtended.png, texcoords=list(x=rep(seq(1, 0, length.out=101), each=4)-c(1:2, 403:404), y=rep(c(0, 1, 1, 0), 50)old - function() shade3d(cylinder3d(cbind(0, 0, c(-1.3, 1), e1=cbind(0, 0, 1), e2=cbind(1, 0, 0), sides=100), color=yellow)# Eyeseyeball - ellipse3d(diag(3), level=radius(.4)shade3d(translate3d(eyeball, .8, .35, .7), color=white)shade3d(translate3d(eyeball, .8, -.35, .7), color=white)# Translate radius of eye, rotate, translate position of eyepupil - rotate3d(translate3d(ellipse3d(diag(3), level=radius(.05), .4, 0, 0), 30/180*pi, 0, 0, 1)shade3d(translate3d(pupil, .8, .35, .7), color=black)shade3d(translate3d(pupil, .8, -.35, .7), color=black)# points3d(1.21, c(-.35, .35), .7, cex=3)# Noseshade3d(cylinder3d(cbind(c(1, 1.3), 0, .3), radius=.2, e1=cbind(1, 0, 0), e2=cbind(0, 1, 0), sides=100), color=yellow)shade3d(ellipse3d(diag(3), level=radius(.2), centre=c(1.3, 0, .3), color=yellow)# Mouthshade3d(elli