plotGEOmapXY {GEOmap}R Documentation

Plot a projected 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, MAPstyle = NULL, border = NA, 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 override color for maps
MAPstyle override plotting style for maps
border color, add border to polygons, NA=no border
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, plotGEOmap does not. MAPcol and MAPstyle can be used to override the colors and style in the map-list. These are applied to all the strokes.

Value

Graphical Side Effects

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

See Also

DOTOPOMAPI, addLLXY, plotGEOmap

Examples


data(japmap)

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, ann=FALSE , axes=FALSE)

plotGEOmapXY(japmap, LIM=c(PLOC$LON[1], PLOC$LAT[1],PLOC$LON[2],
     PLOC$LAT[2]) , PROJ=PROJ, add=TRUE )

addLLXY(PLAT,  PLON, PROJ=PROJ, LABS=TRUE, PMAT=NULL, TICS=c(.1,.1) )


[Package GEOmap version 1.4-0 Index]