# TOPOCOL: Create Topography ColorMAP In GEOmap: Topographic and Geologic Mapping

## Description

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

## Usage

 `1` ```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

## Value

Matrix of colors suitable for insertion to persp

## Author(s)

Jonathan M. Lees<jonathan.lees.edu>

 ``` 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44``` ```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) ```