plotUTM {GEOmap} | R Documentation |
Plot UTM
plotUTM(proj, LIM, shiftlon = 0)
proj |
projection |
LIM |
Limit vector |
shiftlon |
rotation around z axiz, default=0 |
Graphical Side Effect
Jonathan M. Lees<jonathan.lees@unc.edu>
GLOB.XY
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)