lamaz.eqarea {GEOmap} | R Documentation |
Map Projection (Lambert-Azimuthal Equal Area) for global plots.
lamaz.eqarea(phi1, lam0, phi, lam)
phi1 |
Central Latitude, radians |
lam0 |
Central Longitude |
phi |
vector of Latitude, points for plotting, radians |
lam |
vector of Longitude, points for plotting , radians |
x |
position on the plot |
y |
position on the plot |
This is a projection routine that does not need to be set in advance.
Jonathan M. Lees<jonathan.lees@unc.edu>
setPROJ
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")