Program: MKHAMMER
mkhammer.......Hammer projection (equal area)
call: [x,y]=mkhammer(lon,lat,mapopt,direction);
lon: arbitrary matrix of longitude values [deg]
lat: arbitrary matrix of latitude values [deg]
mapopt: map options structure as returned by MKMAPOPT
direction: projection direction:
'forward': transform lon/lat to x/y
'inverse': transform x/y to lon/lat
result: x: matrix x coordinates
y: matrix of y coordinates
X and Y have same size and shape as LON and LAT (values with equal index correspond)
but are projected and scaled according to the settings in MAPOPT
The Hammer projection is an equal area projection.
It expects longitudes from -180 to 180, so better define mapopt.proj_range to be
[-180 180 -90 90]!
The Hammer projection is intended for whole-planet maps.
Martin Knapmeyer, 02.07.2003
Read M-File Source Code
eof.