Program: MKGLOBUSMODELKIT
mkglobusmodelkit......produce sheet of paper that can be glued on a ball
call: [x,y]=mkglobusmodelkit(lon,lat,mapopt,direction);
lon: nearly arbitrary matrix of longitude values [deg]
longitude is expected to range from -180deg to 180deg
(see description of LAT for assumptions made)
lat: nearly arbitrary matrix of latitude values [deg]
(expected to be a rectangular matrix, longitude counted in the
second index!)
mapopt: map options structure as returned by MKMAPOPT
The .proj_specops field is interpreted as folows:
.proj_specops(1): number of lobes
.proj_specops(2): ball radius [cm]
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
This routine produces a plot, which with some scaling, can be glued on
a ball in order to get a globe from map data.
The routine does not only compute map projection data, it also sets some
figure properties to scale the map to a given ball's size\
for .proj_specops(1)=1, this projection becomes the Mercator-Sanson projection.
Before printing the figure, you have to set several figure properties in a way
that one unit of the x and y axes correspond to 1cm on paper.
Martin Knapmeyer, 02.07.2003
Read M-File Source Code
eof.