plotGEOmapXY {GEOmap}R Documentation

Plot a GEO map

Description

High Level plot of GEO map

Usage

plotGEOmapXY(MAP, LIM = c(-180, -90, 180, 90), PROJ = list(), PMAT=NULL,
add = TRUE,  SEL=NULL , GRID = NULL, GRIDcol = 1,
MAPcol = NULL, cenlon = 0, shiftlon = 0, linelty = 1, linelwd = 1, NUMB = FALSE, ...)

Arguments

MAP Map Structure
LIM Lat-Lon limits
PROJ Projection list
PMAT Perspective matrix conversion
add logical, TRUE= add to existing plot
SEL Index vector of strokes to be used in plotting,default=NULL(use all that pass other tests)
GRID logical, TRUE=add grid lines
GRIDcol color for grid lines
MAPcol color for maps
cenlon center longitude of plot
shiftlon degrees, rotate longitude
linelty Line type
linelwd line width
NUMB logical, number the strokes on the map
... graphical parameters

Details

plotGEOmapXY includes projection of the data, though plotGEOmap does not.

Value

Graphical Side Effects

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

See Also

DOTOPOMAPI, addLLXY

Examples


PLOC=list(LON=c(137.008, 141.000), LAT=c(34.000, 36.992), x=c(137.008, 141.000), y=c(34.000, 36.992) )
  PROJ = setPROJ(type=2, LAT0=mean(PLOC$y) , LON0=mean(PLOC$x) )
      
gxy = GLOB.XY(PLOC$LAT, PLOC$LON, PROJ)
 PLAT =  pretty(PLOC$LAT)
    PLAT = c(min(PLOC$LAT),  PLAT[PLAT>min(PLOC$LAT) & PLAT<max(PLOC$LAT)],max(PLOC$LAT)) 
  PLON  = pretty(PLOC$LON)
        PLON = c(min(PLOC$LON), PLON[PLON>min(PLOC$LON) & PLON<max(PLOC$LON)], max(PLOC$LON))
     

plot(gxy$x, gxy$y,  asp=TRUE)
addLLXY(PLAT,  PLON, PROJ=PROJ, LABS=TRUE, PMAT=NULL, TICS=c(.1,.1) )



[Package GEOmap version 1.0-8 Index]