R: Land and Sea Colors

LandSeaCol {GEOmap} | R Documentation

## Land and Sea Colors

### Description

Color pixels with two palettes, one for land the other for sea.

### Usage

    LandSeaCol(IZ, coastmap, PROJ, calcol = NULL)

### Arguments

`IZ` | list of x, y, z suitable for plotting with image or contour.
`coastmap` | coastal map from GEOmap
`PROJ` | projection list
`calcol` | color map for the land

### Details

The program uses closed polygons in the map list to separate the pixels into land versus sea. Sea is colored with a palette of blues, land is colored according to topographic color scheme extracted from palettes similar to GMT palettes.

All map and pixel coordinates are projected with the same projection parameters. calculations are done in XY coordinates.

### Value

`Cmat` | Matrix of colors for each pixel
`UZ` | Under water
`AZ` | Above Sea Level

### Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

### See Also

settopocol, TOPOCOL

### Examples

    ## Not run:

    Lat.range = c(-10, 30)
    Lon.range = c(65, 117)
    ######

    ########  load up the important libraries
     library(RFOC)

     library(geomapdata)

     data(coastmap)

         data(ETOPO5)
     PLOC=list(LON=Lon.range,LAT=Lat.range,lon=Lon.range,lat=Lat.range,
                     x=Lon.range, y=Lat.range )

    #####   set up topography colors
    COLS = settopocol()

    ####  set the projection ##   utm
     PROJ = setPROJ(type=2, LAT0=mean(PLOC$y) , LON0=mean(PLOC$x) )
     NK = 300

       ###   extract topography from the etopo5 data base in geomapdata
         JMAT = GEOTOPO(ETOPO5, PLOC, PROJ, COLS$calcol,nx=NK, ny=NK )
    #####   select relevant earthquakes

    IZ = list(x=JMAT$xo, y=JMAT$yo, z=JMAT$IZ$z)

    CMAT = LandSeaCol(IZ, coastmap, PROJ, calcol=NULL)

    Mollist =CMAT$Cmat
    dMol = attr(Mollist, "Dcol")

         ####  Under water
    UZ = CMAT$UZ
        #####  above water
    AZ = CMAT$AZ
        ####   blues for underwater:
       blues = shade.col(100, acol=as.vector(col2rgb("darkblue")/255),
            bcol= as.vector(col2rgb("paleturquoise")/255))

    plot(x=range(IZ$x), y=range(IZ$y),
             type='n', asp=1, axes=FALSE, ann=FALSE)

    image(x=IZ$x, y=IZ$y, z=(UZ), col=blues, add=TRUE)

    image(x=IZ$x, y=IZ$y, z=(AZ), col=terrain.colors(100) , add=TRUE)

    plotGEOmapXY(coastmap,
           LIM = c(Lon.range[1],Lat.range[1] ,Lon.range[2] ,Lat.range[2]),
           PROJ =PROJ,MAPstyle =2,MAPcol ="black" ,   add=TRUE  )

    ## End(Not run)

[Package GEOmap version 2.4-0 Index]