if (require(rgl)) {
r <- raster(system.file("external/test.grd", package="raster"))
x = r
x[] = 1:ncell(x)
x <- mask(x, r)
plot3D(r, drape=x)
decorate3d(xlab = "x", ylab = "y", zlab = "z", axes=TRUE)
}
Run the code above in your browser using DataLab