if (FALSE) {
data(LIDAR)
mba.int <- mba.surf(LIDAR, 300, 300, extend=TRUE)$xyz.est
##Image plot
image(mba.int, xaxs="r", yaxs="r")
##Perspective plot
persp(mba.int, theta = 135, phi = 30, col = "green3", scale = FALSE,
ltheta = -120, shade = 0.75, expand = 10, border = NA, box = FALSE)
##For a good time I recommend using rgl
library(rgl)
ex <- 10
x <- mba.int[[1]]
y <- mba.int[[2]]
z <- ex*mba.int[[3]]
zlim <- range(z)
zlen <- zlim[2] - zlim[1] + 1
colorlut <- heat.colors(as.integer(zlen))
col <- colorlut[ z-zlim[1]+1 ]
open3d()
surface3d(x, y, z, color=col, back="lines")
}
Run the code above in your browser using DataLab