GEOmap (version 2.1)

plotUTM: Plot UTM

Description

Plot UTM

Usage

plotUTM(proj, LIM, shiftlon = 0)

Arguments

proj

projection

LIM

Limit vector

shiftlon

rotation around z axiz, default=0

Value

• Graphical Side Effect

See Also

GLOB.XY

Examples

data(USAmap)


proj = setPROJ(type=3, LAT0=33.75, LON0= fmod(-79., 360) ,LAT1=34.333333,  LAT2=36.166667, LATS=NULL, LONS=NULL, DLAT=NULL, DLON=NULL,FE=0,FN=0)



ALOC=list(lon=c(274.5,288), lat=c(31, 38), LON=c(274.5, 288), LAT=c(31, 38), shiftlon=0)

plotGEOmapXY(USAmap, LIM=c(ALOC$LON[1], ALOC$lat[1], ALOC$LON[2], ALOC$lat[2]) , PROJ=proj, add=FALSE, shiftlon=0)


plotUTM(proj, c(ALOC$LON[1], ALOC$lat[1], ALOC$LON[2], ALOC$lat[2]))

##############  larger scale

p = plotusa()

plotUTM(p$PROJ, LIM=p$LIM)