Numeric, to set the elevation scale relative to x and y
drape
RasterLayer, to 'drape' colors representing the values of
this layer on the 3D representation of layer x. In this
casex typically has elevation data
col
A color palette generating function such as rainbow,
heat.colors, and topo.colors, or one or your own making
at
A numeric variable of breakpoints defining intervals along
the range of x or a number definining the number of
intervals the range of x will be divided into.
rev
Logical. If TRUE, the color palette values are
reversed in order
useLegend
Logical. If TRUE (default) the content of the
slot x@legend@colortable is used instead of col and
at.
adjust
Logical. If TRUE, the x and y axes are scaled
relative to the cell (z) values
...
Any argument that can be passed to
surface3d
Note
Previous versions opened a new device with each call of
plot3D. This behaviour has been modified, and now a new device is
opened only if none is active. Thus, you should use rgl::open3d()
if you need to display a new scene without modifying the current one.
Author(s)
Robert J. Hijmans and Oscar Perpi<c3><83><c2><b1><c3><83><c2><a1>n
Examples
## Not run:
## rgl is needed to use plot3D
library(rgl)
data(volcano)
r <- raster(volcano)
extent(r) <- c(0, 610, 0, 870)
## level plot as reference
levelplot(r, col.regions=terrain.colors)
plot3D(r)
## Use different colors with a predefined function
plot3D(r, col=rainbow)
## or with a custom function using colorRampPalette
myPal <- colorRampPalette(brewer.pal(11, 'PuOr'))
plot3D(r, col=myPal)
## With at you can define an homogeneus color table for different Rasters
r2 <- r + 100
r3 <- r + 200
s <- stack(r, r2, r3)
maxVal <- max(maxValue(s))
minVal <- min(minValue(s))
N <- 40
breaks <- seq(minVal, maxVal, length=N)
plot3D(r, at=breaks)
plot3D(r2, at=breaks)
plot3D(r3, at=breaks)
## Default: x-axis and y-axis are adjusted with z-values. Therefore,
## labels with decorate3d() are useless
plot3D(r, adjust=TRUE)
decorate3d()
## Compare the graphic limits
par3d('bbox')
## with the extent of the Raster
extent(r)
## Set adjust=FALSE to fix it
plot3D(r, adjust=FALSE)
decorate3d()
## Once again, compare the graphic limits
par3d('bbox')
## with the extent of the Raster
extent(r)
## zfac controls the z values so z-axis will be distorted
plot3D(r, adjust=FALSE, zfac=2)
decorate3d()
par3d('bbox')
## With drape you can disconnect the z-axis from the colors
drape <- cut(r^4, 4)
plot3D(r, drape=drape)
## Compare with:
plot3D(r, at=4)
## End(Not run)