plotm.m

主元分析法的相关算法

%PLOTM  Trace un champ de donnee geographique
% [CS,H]=PLOTM(C,[ZONE]) 
%   Trace le champ C(LONGITUDE,LATITUDE) sur une carte 
%   du monde en projection Mercator par defaut.
%   Le parametre optionnel ZONE peut prendre les valeurs
%     ZONE = 1 : Projection Stereographique centree sur le pole sud
%     ZONE = 2 : Projection Mercator de l'hemisphere sud
%     ZONE = 3 : Projection Mercator du globe
%     ZONE = 4 : Pas de projection (imagesc)
%     ZONE = 5 : Projection Mercator de l'hemisphere sud complet
%     ZONE = 6 : Projection Mercator du globe complet
%     ZONE = 7 : Projection Stereographique centree sur le pole nord
%     ZONE = 8 : Idem 2 mais sans label sur les axes
%     ZONE = 9 : Projection Mercator de l'hemisphere sud aux limites du Fmask
%   Returns contour matrix CS as described in CONTOURC and a vector H
%   of handles to PATCH objects (for use by CLABEL).
%
%   La definition des grilles (longitude,latitude) est celle du
%   modele par defaut, mais elle peut etre modifiee en passant en
%   global les nouvelles variables nommees lon et lat.
%
% NECESSITE LE PACKAGE M_MAP
%
%
% (nov. 2003 -> [e-mail])

function varargout = plotm(C,varargin)

if (nargin>2)|(nargin<1)
     help plotm.m
     error('plotm.m : Wrong number or bad parameter(s)')
     return
elseif (nargin==2)
   arg = varargin{:};
   if (arg>14)|(arg<1)
     help plotm.m
     error('plotm.m : Wrong number or bad parameter(s)')
     return
   end
end %if

%--------------------------------------------------
% Variabes
%--------------------------------------------------
[lx,ly]=size(C);
% Si lon et lat ne sont pas definie comme dans le modele
% on prend en global n'importe quel definition
if ((lx~=64)|(ly~=32)) 
   global lon lat
else
   lon= linspace(0.,360.,64);
   lat= [85.76 80.27 74.74 69.21 63.68 58.14 52.61 47.07 41.53 ...
      36.00 30.46 24.92 19.38 13.84 8.31 2.77 -2.77 -8.31 ...
      -13.84 -19.38 -24.92 -30.46 -36.00 -41.53 -47.07 -52.61 ...
      -58.14 -63.68 -69.21 -74.74 -80.27 -85.76];
end

% Le nombre de contour est par defaut a 10, mais on peut le
% modifier en declarant nbcont en global
global nbcont
if isempty(nbcont)==1
   nbcont=10;
end

%-----------------------------------------------
% Projection
%-----------------------------------------------
if (nargin==2)
   arg = varargin{:};
   switch arg(1)
     case 1 % Projection stereo 
       m_proj('stereo','lon',0,'lat',-90,'rad',71);
     case 2 % Projection mercator (HS)
       m_proj('miller','lon',[0 360],'lat',[-80 -20]);
       m_proj('miller','lon',[0 360],'lat',[-60 -20]);
     case 3 % Projection mercator (Tout)  -> Defaut
       m_proj('miller','lon',[0 360],'lat',[-80 80]);
     case 5 % Projection mercator (HS-complet)
       m_proj('miller','lon',[0 360],'lat',[-85 -5]);
     case 6 % Projection mercator (Tout-complet)
       m_proj('miller','lon',[0 360],'lat',[-85 85]);
     case 7 % Projection stereo 
       m_proj('stereo','lon',0,'lat',90,'rad',71);
     case 8 % Projection mercator (HS) sans Label sur les axes
       m_proj('miller','lon',[0 360],'lat',[-80 -20]);
       m_proj('miller','lon',[0 360],'lat',[-75 -30]);
     case 9 % Projection mercator (HS) sans Label sur les axes aux limtes du FMASK
       m_proj('miller','lon',[0 360],'lat',[-69 -25]);
     case 10 % Projection stereo recentree sur le pole sud
       m_proj('stereo','lon',0,'lat',-90,'rad',51);
     case 11 % Indien
       m_proj('lambert','long',[11 125],'lat',[-60 -30]);
     case 12 % Pacifique
       m_proj('lambert','long',[120 290],'lat',[-60 -30]);
     case 13 % Atlantique
       m_proj('lambert','long',[-70 20],'lat',[-60 -30]);
     case 14 % Centr