Get Lat-Lon Perimeter
Get rectangular perimeter of region defined by set of Lat-Lon
getGEOperim(lon, lat, PROJ, N)
lon
: vector of lonslat
: vector of latsPROJ
: projection structureN
: number of points per sideperimeter is used for antipolygon
List: - x: x-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)
Useful links