GEOmap (version 2.1)

getGEOperim: 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

Value

• List:
• x
  x-coordinates projected
• y
  y-coordinates projected

Details

perimeter is used for antipolygon

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)