TOPOCOL function

Create Topography ColorMAP

Create Topography ColorMAP

Given an x-y-Z create a matrix of colors for plotting in persp

TOPOCOL(IZ, calcol)

Arguments

  • IZ: Matrix of values
  • calcol: Color mapping of elevations to rgb colors

Details

colors are interpolated between boundaries in the color map

Returns

Matrix of colors suitable for insertion to persp

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

See Also

persp

Examples

colk1 = 50 colk2 = 210 colk3 = 220 colk4 = 250 BWpal2 = list(z1=c(-3000, 0, 2000, 3500), r1=c(0,colk1, colk3, colk4), g1=c(0,colk1, colk3, colk4), b1=c(0,colk1, colk3, colk4), z2=c(0, 2000, 3500, 5000), r2=c(0,colk2,colk4,255), g2=c(0,colk2,colk4,255), b2=c(0,colk2,colk4,255), note=c("black, black", "grey, grey", "white, white", "white, white") ) data(volcano) MYCOLL = TOPOCOL(volcano, BWpal2) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA par(bg = "slategray") Dcol = attr( MYCOLL , "Dcol") persp(x, y, z, theta = 135, phi = 30, col = MYCOLL[1:(Dcol[1]-1), 1:(Dcol[2]-1)], scale = FALSE, ltheta = -120, shade = 0.75, border = NA, box = FALSE) calcol=settopocol() MYCOLL = TOPOCOL(volcano, calcol$calcol) Dcol = attr( MYCOLL , "Dcol") K <- 8 *volcano MYCOLL = TOPOCOL(K, calcol$calcol) persp(x, y, z, theta = 135, phi = 30, col = MYCOLL[1:(Dcol[1]-1), 1:(Dcol[2]-1)], scale = FALSE, ltheta = -120, shade = 0.75, border = NA, box = FALSE)
  • Maintainer: Jonathan M. Lees
  • License: GPL (>= 2)
  • Last published: 2024-07-09

Useful links

    Downloads (last 30 days):