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.