GEOmap (version 2.1)

lamaz.eqarea: 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

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