#--------------------------------------- # BY Jesus Eduardo CanulK # # If you like please upvote and comment # v=4.2k 4.3k 4.4k #-----—--------------------------------- z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) #10 meter spacing (S-N) y <- 10 * (1:ncol(z)) #10 meter spacing (E-W) z0 <- min(z) - 25 z <- rbind(z0, cbind(z0, z, z0), z0) x <- c(min(x) - 1e-10, x, max(x) + 1e-10) y <- c(min(y) - 1e-10, y, max(y) + 1e-10) t <- "transparent" Persp <- function(titulo,subtitulo,color){ persp(x, y, 2.5*z, theta = 135, phi = 30, scale = FALSE, ltheta = -120, shade = 0.45, border = NA, box = FALSE, main=paste("Renderizado de",titulo), col= color, col.main="white", col.sub="white", sub=subtitulo) } #------------------------ #Example 1 png(file="Avolcan.png",bg=t) Persp("color Blanco","","#ffffff") #----------------------------- #Example 2 fill <- matrix(nrow = nrow(z)-1, ncol = ncol(z)-1) fill[ , i2 <- c(1,ncol(fill))] <- "#000080" fill[i1 <- c(1,nrow(fill)) , ] <- "#000080" fcol <- fill png(file="Zvolcan.png",bg=t) zi <- volcano[ -1,-1] + volcano[ -1,-61] + volcano[-87,-1] + volcano[-87,-61] / 4 fcol[-i1,-i2] <- hcl.colors(20,"Plasma")[cut(zi,stats::quantile(zi, seq(0,1, length.out = 21)), include.lowest = TRUE)] Persp("paleta Plasma","",fcol) #----------------------- #Example 3 fill[ , i2 <- c(1,ncol(fill))] <- "green3" fill[i1 <- c(1,nrow(fill)) , ] <- "green3" fcol <- fill png(file="Jvolcan.png",bg=t) zi <- volcano[ -1,-1] + volcano[ -1,-61] + volcano[-87,-1] + volcano[-87,-61] / 4 fcol[-i1,-i2] <- terrain.colors(20)[cut(zi, stats::quantile(zi, seq(0,1, length.out = 21)), include.lowest = TRUE)] nota="Suggestions are welcome, if you like please upvote and comment" Persp("Terrain color",nota,fcol) #---------------------- cat("

PART 3: Rendering

By Jesus Eduardo CanulK

") #By Jesus Eduardo Canulk

⛰️⛰️ Volcano 3D Part 3 | R ⛰️⛰️