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
Afficher la suite