A explicação do script a seguir está no artigo R na Prática – Parte 31
#fonte rdocumentation.org
z <- 2 * volcano #exagero vertical
x <- 10 * (1:nrow(z))
y <- 10 * (1:ncol(z))
def.par=par(no.readonly=TRUE)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_22.jpeg",width=5,height=5,
#units="in",res=300)
par(bg = "white")
persp(x, y, z, theta = 135, phi = 30,
col = "green3", scale = FALSE,
ltheta = -120, shade = 0.75,
border = NA, box = FALSE)
par(def.par)
#dev.off()