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.