Aide sur matlab

Fermé
debutant_matlab - 8 avril 2010 à 14:46
j'ai ecris ce programme pour lire un fichier et tracer l'image mais ça me affiche des image bizzarre.
j'ai fourni tout mes effort pour resoudre le probleme mais jusqu'a présent rien n'a changer.je voulais joindre le fichier mais je ne vois la ou on peut joindre un fichier.je suis prés pour discuiter sur email



close all;
clear all;

clc;
Nb_jour=31;
Nb_annee=12;

Nb_lignes = 241;
Nb_colonnes = 201;

Nb_lignes1 = 548;
Nb_colonnes1 = 456;

Nb_pixels = Nb_lignes*Nb_colonnes;
Nb_pixels_48 = Nb_lignes1*Nb_colonnes1;
sst=zeros(Nb_pixels,12);


mois={'janv', 'fev', 'mars', 'avril', 'mai', 'juin', 'juillet', 'aout', 'sept', 'oct','nov', 'dec'};
disp('Programme de lecture des données SST');

for j=1:12
%adapter le chemin d'acces
s = ['fid=load([''C:\MATLAB6p1\work\data_perou_sst\SST\2003\sst_' mois{j} '_2003.txt''], ''r'' );'];
eval(s);

sst_old=fid(2:end,2:end);
lon_old=fid(1,2:end);
lat_old=fid(2:end,1);
sst_old(sst_old==-1.E+07)=NaN;
lon_old=-lon_old;

%%%%%%%%%%%%%%%%%%%%%%%
%transformation de la matrice en 201x241
xmin = min(lon_old);
xmax = max(lon_old);
ymin = min(lat_old);
ymax = max(lat_old);

xi = -90:0.1:-70;
yi = -20:0.1:4;

xi = linspace(xmin, xmax, 241);
yi = linspace(ymin, ymax, 201);
zi = griddata(lon_old, lat_old, sst_old, xi', yi);

%sauvegarde des matrices
%adapter le chemin de sauvegarde
eval(['sst_' mois{j} '=zi;'])
eval(['save C:\MATLAB6p1\work\data_perou_sst\SST\2003\sst_' mois{j} ' ' 'sst_' mois{j} ';'])

sst(:,j)=zi(:);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%images de l'annee 2002

valeur_titre={'janvier','fevrier','mars','avril','mai','juin','juillet','aout','septembre','octobre','novembre','decembre'};

for j=1:12

figure(1);
subplot(3,4,j);
zi=reshape(sst(:,j),201,241);
CLim=[min(zi(:)) max(zi(:))];
xi = -90:0.1:-70;
yi = -20:0.1:4;

imagesc(xi,yi,zi,CLim);
titre=['SST 'valeur_titre{j} ' 2003'];
set(gca,'xtick',[-90:4:-70],'xticklabel',{'-90', '-86', '-82', '-78', '-74', '-70'}, ...
'ytick',[-20:4:4],'yticklabel',{'4','0', '-4', '-8', '-12', '-16', '-20' })
set(get(gca,'title'),'string',titre);
h=colorbar;
localcolor=jet(255);
MAP=[0.7 0.7 0.7;localcolor];
colormap(MAP)
axis image
set(get(h,'ylabel'),'String','SST(°C)');
xlabel('longitude');
ylabel('latitude');
%grid on;

end