getGEOperim {GEOmap}

Get Lat-Lon Perimeter

Description

Get rectangular perimeter of region defined by set of Lat-Lon

Usage

getGEOperim(lon, lat, PROJ, N)

Arguments

lon

vector of lons

lat

vector of lats

PROJ

projection structure

N

number of points per side

Details

perimeter is used for antipolygon

Value

List:

x

x-coordinates projected

y

y-coordinates projected

Author(s)

Jonathan M. Lees<jonathan.lees.edu>

Examples

###  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)