## Not run:
# 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")
# ## End(Not run)
Run the code above in your browser using DataCamp Workspace