getGEOperim {GEOmap} | R Documentation |
Get rectangular perimeter of region defined by set of Lat-Lon
getGEOperim(lon, lat, PROJ, N)
lon |
vector of lons |
lat |
vector of lats |
PROJ |
projection structure |
N |
number of points per side |
perimeter is used for antipolygon
List:
x |
x-coordinates projected |
y |
y-coordinates projected |
Jonathan M. Lees<jonathan.lees.edu>
### target region PLOC= list(LON=c(138.3152, 139.0214), LAT=c(35.09047, 35.57324)) PLOC$x =PLOC$LON PLOC$y =PLOC$LAT #### set up projection PROJ = setPROJ(type=2, LAT0=mean(PLOC$y) , LON0=mean(PLOC$x) ) perim= getGEOperim(PLOC$LON, PLOC$LAT, PROJ, 50)