Discussion :

[Potté75]

BonjourBonjour!

Je vous ecris pour un petit soucis concernant deux boucles et une avec subplot.
en fait je voudrais afficher 12 plots (1 par mois)
bon je peux le faire "manuellement" mais c'est pas très glamour et sa rend mon script tres lourd
je pense qu'il est possible de faire des boucles mais après plusieurs tests rien ne fonctionne!

donc j'ai un critère de base: le mois

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
 
mois=str2num(datestr(DATEVRAC,5));
id=find(mois==1);   %1=janvier, 2=fevrier.....

et je voudrais afficher un plot de ma commande, par mois.
il faut donc deux boucles: une pour avancer le mois (faire un +1 sur id) et une boucle pour afficher a chaque fois un subplot
et j'ai un peu de mal dans l'organisation...
je tente avec ca:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
i0=id;
for k = 1:12
    subplot(2, 6, k)
	for i0:(id+1)
		MA MEGA COMMANDE
	end
end

mais sans succes...
Pourrize vous me donner un coup de pouce s'il vous plait?
Merci a vous
Bonne journée

[Gooby]

Bonjour,

Je n'ai pas compris à quoi sert ta boucle sur io, peux tu expliquer?

[Potté75]

Hello!

ma boucle sur i0 sert a faire avancer de "un mois" a chaque tour de boucle: je pars de janvier et je veux un plot pour janvier, un plot pour fevrier...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
id=find(mois==1);   %1=janvier, 2=fevrier.....
i0=id;
	for i0:(id+1)
		MA MEGA COMMANDE
end

[Gooby]

dans ce cas c'est la même boucle que celle pour k?

[Potté75]

Hum..
je pensais qu'il fallait UNE boucle pour i0 (faire +1mois a chaque tour) ET UNE boucle pour k, qui m'ajoute un plot a chaque tour de boucle...

in fine je voudrais une page graphique de matlab avec 12 petites cartes.. donc un plot par mois
je ne veux pas faire apparaitre une carte pour janvier, la supprimer et faire aparaitre la carte pour fevrier...
je veux un graph avec les 12 mois... est-ce plus clair?
dans ce cas la, il ne faut pas deux boucles??

[Gooby]

Non, une boucle suffit.

avant de pouvoir te répondre, j'aimerai savoir ce que contiens ta variable mois et quelle est cette super commande?

[Potté75]

ok!
alors la variable "mois" est construite comme dans le 1er post de la discussion...
j'utilise la fonction "datestr" sur DATEVRAC qui est un vecteur de 200000 lignes
et donc ca convertie les dates de DATEVRAC en chiffres représentant le mois

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
mois(1:10)
ans =
 
     2
     4
     4
     4
     5
     5
     5
     5
     5
     7

donc la première ligne, c'est un mois de février, la seconde c'est un mois d'avril ....

ensuite je sélectionne par exemple tous les mois de janvier

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

id=find(mois==1)

et je fais tourner la "méga" commande et a la fin de cette commande j'obtiens une carte pour les mois de janvier (par exemple, ou fevrier si je mets mois==2)...
La méga commande n'a rien de spéciale, c'est juste qu'elle est longue ^^
c'est un m_pcolor sur lequel tu m'avais déjà aidé la semaine derniere Gooby!

[Gooby]

En fait je ne comprends pas bien, quelle est ta "variable d'entrée" de ta super commande? Je pense qu'une seule boucle suffit mais sans savoir ce que tu veux faire une fois que tu as repéré où se trouve le mois de Janvier (par exemple), je ne puis en être sûr.

Qu'est-ce que tu essayes de faire en faisant

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
id=find(mois==1)
i0=id;
	for i0:(id+1)

i0 contient un vecteur, de même que id (puisqu'id = i0).

[Potté75]

ca donne ca:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
mois=str2num(datestr(DATEVRAC,5));
id=find(month==1);
lon=LONGITUDE(id); %LONGITUDE est un vecteur
lat=LATITUDE(id);      %LATITUDE est un vecteur
Z=PROF(id);               %PROF est un vecteur
warning off
m_proj('stereographic','latitude',-90,'long',-30,'radius',55,'rotangle',25);
[X,Y]=m_ll2xy(lon,lat);
m_coast('patch',[0.3 0.3 0.3],'edgecolor',[0.3 0.3 0.3]);
m_grid('box','fancy','tickdir','out','xtick',7, 'yticklabel',[ ],'linestyle','none','fontsize',16,'XaxisLocation','top');
hold on
gx=min(lon):max(lon);
gy=min(lat):max(lat);
[xm,ym]=meshgrid(gx,gy);
ZI=griddata(lon,lat,Z,xm,ym);
m_pcolor(gx,gy,ZI);
shading interp
c=colorbar;load colormapbicolor; colormap(colormapbicolor);
set(c,'Position',[0.9 0.2583 0.0446 0.4167]);

Voila grosso-modo la commande. Donc a la fin j'obtiens une carte.
et je voudrais mettre tout ca dans une boucle ou i0 augmente de 1 a chaque tour pour finir a 12: de cette facon on aura une carte par mois...
et en meme temps un sublot de type

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

subplot(2,6,1);

pour afficher les 12 cartes

[Gooby]

Etant donné que tu as déjà une boucle k allant de 1 à 12, pourquoi refaire une boucle qui elle aussi parcourra 1 à 12?

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
for k = 1:12
    subplot(2, 6, k)
    idx_temp = find(mois==k)
    lon = LONGITUDE(idx_temp);
    lat = LATTITUDE(idx_temp);
    %la suite du code...
    clear idx_temp
end

[Potté75]

ROOOH si tu voyais la beauté de mes 12 cartes sous les yeux!!
Merveilleux! Merci

Oui je ne sais pas pourquoi, je réfléchis pas encore en langage matlab.. Pour moi il fallait absolument deux boucles..
mais en voyant tes lignes de codes, c'est vrai que c'est très clair
Merci une fois de plus