World two scales

## =====================================================================
## World map with two scales; uses dataset Hypsometry from OceanView
## implemented by Karline Soetaert
## =====================================================================
library(OceanView)

DD <- Hypsometry$z; DD[DD > 0] <- NA
HH <- Hypsometry$z; HH[HH < 0] <- NA
ocean.col <- ramp.col (c("lightblue", "darkblue"))
land.col  <- ramp.col (c(gray(0.7), gray(0.2)))

image2D(-DD, col = ocean.col, rasterImage = TRUE,  
  axes = FALSE, xlab = "", ylab = "",  shade = 0.05, 
  zlim = c(0, 8000),
  colkey = list(length = 0.3, shift = 0.15, width = 0.5,
  cex.axis = 0.6))
image2D(HH, add = TRUE, NAcol = "transparent", 
  col =  land.col, rasterImage = TRUE, shade = 0.05,  
  colkey = list(length = 0.3, shift = -0.15, width = 0.5,
  cex.axis = 0.6))

Go back