Program: MKLAMBERT
mklambert.......Lambert's equal area azimutha projection
call: [x,y]=mklambert(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
Martin Knapmeyer, 07.07.2003
Read M-File Source Code
eof.