Get Lat-Lon Perimeter

Description

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

Usage

1
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
###  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)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.