lamaz.eqarea {GEOmap}R Documentation

Lambert-Azimuthal Equal Area

Description

Map Projection (Lambert-Azimuthal Equal Area) for global plots.

Usage

lamaz.eqarea(phi1, lam0, phi, lam)

Arguments

phi1 Central Latitude, radians
lam0 Central Longitude
phi vector of Latitude, points for plotting, radians
lam vector of Longitude, points for plotting , radians

Value

x position on the plot
y position on the plot

Note

This is a projection routine that does not need to be set in advance.

Author(s)

Jonathan M. Lees<jonathan.lees@unc.edu>

See Also

setPROJ

Examples

data(coastmap)
#########  coastmap is a GEOmap list
DEGRAD = pi/180

phicen  = -90*DEGRAD
lamcen  = 0*DEGRAD

i = 7
j1 = coastmap$STROKES$index[i]+1
j2 = j1+ coastmap$STROKES$num[i]-1
lat = coastmap$POINTS$lat[j1:j2]*DEGRAD
lon = coastmap$POINTS$lon[j1:j2]*DEGRAD

 xy = lamaz.eqarea(phicen, lamcen,lat, lon)

plot(xy, asp=1, type='n')

 polygon(xy, col=grey(.8))

title("Antarctica")




[Package GEOmap version 1.4-0 Index]