source: trunk/toolbox/plot_image.m @ 23

function [Image2d]=plot_image(data,Valid_Pixels,Val_Pix_Nuages,Val_Pix_Terre,CLim)
%[Image2d]=plot_image(data,Valid_Pixels,Val_Pix_Nuages,Val_Pix_Terre,CLim)
%Affiche la donnée data sur l'image de la mer Méditerranée avec le masque
%de Terre et de Nuage.
%Renvoie la donnée dans une matrice à 2 dimensions de taille (512,1024).
%Les données manquantes (nuages/terre) sont mises à la valeur
%nan;
%CLim représente l'échelle de couleur Clim=[val_min val_max]

Image2d=nan*ones(512,1024);
Image2d(Valid_Pixels)=data;
Cdirect=nan*ones(512,1024);

%CLim=[0 0.05];

%cmap=[data;manquantes;land;cloud];
N=64;
cmap=[jet(N);[0.8 0.8 0.8]; [0.8 0.6 0.4] ; [1 1 1]];

data(data>=CLim(2))=CLim(2);
data(data<=CLim(1))=CLim(1);

X=floor((N-1)*(data-CLim(1))/(CLim(2)-CLim(1)))+1;

Cdirect(Valid_Pixels)=X;
Cdirect(Val_Pix_Nuages)=N+3;
Cdirect(Val_Pix_Terre)=N+2;
Cdirect(isnan(Cdirect))=N+1;
clf
image(linspace(-9.5,41.5,1024),linspace(49,29,512),Cdirect);
colormap(cmap);
axis image;
axis xy;
xlabel('longitude')
ylabel('latitude')
axis image

%Colorbar
pos=[0.2 0.1 0.6 0.04];
lev=linspace(CLim(1),CLim(2),10);
lev_pas=2;
clip_lev=CLim;
hpal=cmap;
orien=0;
[h0]=colorbartype1(pos,lev,lev_pas,clip_lev,hpal,orien);