Discussion :

[niepoc]

Bonjour

dans un programme, je lis une 18 fichiers .mat en boucle (6 par année) où j'effectue un calcul

je souhaiterais enregistrer les fichiers correspondants à une année sur une figure avec un subplot

voilà le code de la lecture en boucle des fichiers .mat

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
clear all
close all
 
%% FICHIERS A LIRE
rep = 'G:\AeRoNeT\comptage des évènements poussières et histo\histogramme';
ext = '*.mat';
 
chemin = fullfile(rep,ext);
 
listAERONET = dir(chemin);
Nom=numel(listAERONET);
p=0;
figure('Name','Tous les sites année 2005','NumberTitle','off') 
%figure('Name','Tous les sites année 2006','NumberTitle','off') 
%figure('Name','Tous les sites année 2007','NumberTitle','off') 
%% LECTURE DU FICHIER ET INITIALISATION DE LA MATRICE
for f=1:Nom
    z=f;
    fichier =listAERONET(f).name
 
%récupération des données
%fichier='05Agoufou.mat';
load(fichier)
    %nbjouraot
    %AOTfort
    %AOTintense
    %moyjanjuin
    %moyjuidec
    %moyfevmai
    M;
    joursintenses;
    joursforts;
    AOTfort;
    AOTintense;

ensuite vient un calcul, je vous epargne le code, le résultat est 'AOT'

puis l'enregistrement sur la figure avec ce code
code

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
20
21
22
23
24
if str2double(fichier(1:2))==05   %permet d'identifier les fichiers de 2005 et de les enregistrer sur la figure 'Tous les sites année 2005'
p=p+1;   
subplot (3,2,p)
bar(AOT,'group')
title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')
%AJOUTER LES MOIS
set(gca,'xlim',[0 13],'xtick',1:12)
Xt = 1:12;
% Add the months as tick labels.
months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];
ax = axis; % Current axis limits
axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual
Yl = ax(3:4); % Y-axis limits
% Place the text labels
t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));
set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...
'Rotation',45);
% Remove the default labels
set(gca,'XTickLabel','')
h = legend('Evènements forts','Evènements intenses',1);
set(h,'Interpreter','none')
 
end
end

le deuxième "end" correspond à la fin de ma boucle for du début

je n'arrive à créer dans mon programme une commande pour demander d'enregistrer sur la figure 'Tous les sites de année 2006' les fichiers de 2006 et idem pour l'année 2007
actuellement il enregistre mes résultat sur la derniere appelée, ici 'Tous les sites année 2005'

merci pour l'aide en espérant d'avoir été clair.

[christophe_halgand]

Salut,

Tu veux enregistrer tes 3 figures ou les données te permettant de faire ses figures ?

Je suppose que si tu arrives à le faire pour le dernier indice dans ta boucle, c'est que tu as un problème qui est souvent comme celui-ci :

Code :

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

1
2
3
4
for t=1:5
vect=t;
vect(t)=t;
end

Regardes si tu n'aurais pas fait ce type d'erreur...

[niepoc]

bonjour merci pour la réponse

l'idée c'est de pouvoir enregistrer les données de 2005 sur la figure correspondant à l'année 2005, les données de 2006 sur la figure de 2006 et celles de 2007 sur la figure de 2007

avec un code comme celui ci

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
if str2double(fichier(1:2))==05   %permet d'identifier les fichiers de 2005
% enregistrer sur la figure 'Tous les sites année 2005'
%je ne sais pas comment dire d'enregistrer les données sur la figure 'Tous les sites année 2005
p=p+1;   
subplot (3,2,p)
bar(AOT,'group')
title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')
%AJOUTER LES MOIS
set(gca,'xlim',[0 13],'xtick',1:12)
Xt = 1:12;
% Add the months as tick labels.
months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];
ax = axis; % Current axis limits
axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual
Yl = ax(3:4); % Y-axis limits
% Place the text labels
t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));
set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...
'Rotation',45);
% Remove the default labels
set(gca,'XTickLabel','')
h = legend('Evènements forts','Evènements intenses',1);
set(h,'Interpreter','none')

if str2double(fichier(1:2))==06   %permet d'identifier les fichiers de 2006
% enregistrer sur la figure 'Tous les sites année 2006'
%[COLOR="Red"]je ne sais pas comment dire d'enregistrer les données sur la figure 'Tous les sites année 2006'/COLOR]
p=p+1;   
subplot (3,2,p)
bar(AOT,'group')
title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')
%AJOUTER LES MOIS
set(gca,'xlim',[0 13],'xtick',1:12)
Xt = 1:12;
% Add the months as tick labels.
months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];
ax = axis; % Current axis limits
axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual
Yl = ax(3:4); % Y-axis limits
% Place the text labels
t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));
set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...
'Rotation',45);
% Remove the default labels
set(gca,'XTickLabel','')
h = legend('Evènements forts','Evènements intenses',1);
set(h,'Interpreter','none')

%....puis 2007
end

merci pour l'aide

bonne journée

[Caro-Line]

Je dirais :
Là :

Code :

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

figure('Name','Tous les sites année 2005','NumberTitle','off')

Mettre le handle en sortie :

Code :

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

hFug2005=figure('Name','Tous les sites année 2005','NumberTitle','off')

Et dans ton if :

Code :

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

1
2
3
4
if...%2005
   figure(hFig2005);%force La figure 2005 a être la figure courante
   %....
end

Normalement cela fonctionne.
Ce n'est pas forcément le plus propre, mais c'est plus rapide.
La plus propre c'est pour chaque commande qui "écrité dans la figure de lui passer à un moment ou un autre le handle de la figure.
Ex avec subplot :

Code :

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

subplot(3,2,p,'Parent',hFig2005)

[HS]Un peu de vocabulaire :
On ne dit pas "enregistrer" sur la figure.
"enregistrer" (pour moi en tout cas) sous-entend une notion de sauvegarde dans un contenant fichier ou variable.
On dit plutôt "tracer"[/HS]

[niepoc]

merci pour l'aide

ok pour le vocabulaire, désolé

l'idée est bien de forcer les figures pour qu'elles soient les figures courantes

malheureusement ça ne trace pas mes résultats avec les deux tests

code pour appeler les figures

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
20
21
22
23
24
25
26
27
28
%% FICHIERS A LIRE
 
rep = 'G:\AeRoNeT\comptage des évènements poussières et histo\histogramme';
 
ext = '*.mat';
 
chemin = fullfile(rep,ext);
 
listAERONET = dir(chemin);
 
Nom=numel(listAERONET);
 
p5=0;
p6=0;
p7=0;
 
hFig2005=figure('Name','Tous les sites année 2005','NumberTitle','off');
hFig2006=figure('Name','Tous les sites année 2006','NumberTitle','off'); 
hFig2007=figure('Name','Tous les sites année 2007','NumberTitle','off'); 
 
%% LECTURE DU FICHIER ET INITIALISATION DE LA MATRICE
 
for f=1:Nom
    z=f;
    fichier =listAERONET(f).name
 
%récupération des données
load(fichier)

premier test

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
if str2double(fichier(1:2))==05  
 
figure(hFig2005)
 
p5=p5+1;   
 
subplot (3,2,p5)%,'Parent',hFig2005)
 
bar(AOT,'group')
 
title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')
 
%AJOUTER LES MOIS
 
set(gca,'xlim',[0 13],'xtick',1:12)
 
Xt = 1:12;
 
% Add the months as tick labels.
 
months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];
 
ax = axis; % Current axis limits
 
axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual
 
Yl = ax(3:4); % Y-axis limits
 
% Place the text labels
 
t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));
 
set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...
 
'Rotation',45);
 
% Remove the default labels
 
set(gca,'XTickLabel','')
 
h = legend('Evènements forts','Evènements intenses',1);
 
set(h,'Interpreter','none')
 
end

deuxième test:

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
if str2double(fichier(1:2))==05  
 
%figure(hFig2005)
 
p5=p5+1;   
 
subplot (3,2,p5,'Parent',hFig2005)
 
bar(AOT,'group')
 
title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')
 
%AJOUTER LES MOIS
 
set(gca,'xlim',[0 13],'xtick',1:12)
 
Xt = 1:12;
 
% Add the months as tick labels.
 
months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];
 
ax = axis; % Current axis limits
 
axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual
 
Yl = ax(3:4); % Y-axis limits
 
% Place the text labels
 
t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));
 
set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...
 
'Rotation',45);
 
% Remove the default labels
 
set(gca,'XTickLabel','')
 
h = legend('Evènements forts','Evènements intenses',1);
 
set(h,'Interpreter','none')
 
end

dans les deux cas les fenêtres des figures s'ouvrent mais rien n'est tracé

merci pour l'aide

[niepoc]

j'ai trouvé

j'avais fait une petite erreur d'étourderie dans mon programme

voilà le code en entier, j'ai souligné e rouge les parties qui nous intéressent

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
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
clear all

close all

%% FICHIERS A LIRE

rep = 'Y:\Lexar\AeRoNeT\comptage des évènements poussières et histo\histogramme';

ext = '*.mat';

chemin = fullfile(rep,ext);

listAERONET = dir(chemin);

Nom=numel(listAERONET);

p5=1; %place dans le subplot p5 pour l'année 2005

p6=1;

p7=1;

Fig2005=figure('Name','Tous les sites année 2005','NumberTitle','off');

Fig2006=figure('Name','Tous les sites année 2006','NumberTitle','off'); 

Fig2007=figure('Name','Tous les sites année 2007','NumberTitle','off'); 

%% LECTURE DU FICHIER 

for f=1:Nom

    z=f;

    fichier=listAERONET(f).name

    

%récupération des données

load(fichier)

    M;

    joursintenses;

    joursforts;

    AOTfort;

    AOTintense;

    

AOTjan=0;AOTfev=0;AOTmar=0;AOTavr=0;AOTmai=0;AOTjuin=0;

AOTjui=0;AOTaou=0;AOTsep=0;AOToct=0;AOTnov=0;AOTdec=0;

    

for n=1:length(joursforts)

    if joursforts(n,2)==1

        AOTjan=AOTjan+1;

    elseif joursforts(n,2)==2

        AOTfev=AOTfev+1;

    elseif joursforts(n,2)==3

        AOTmar=AOTmar+1;

    elseif joursforts(n,2)==4

        AOTavr=AOTavr+1;

    elseif joursforts(n,2)==5

        AOTmai=AOTmai+1;

    elseif joursforts(n,2)==6

        AOTjuin=AOTjuin+1;

    elseif joursforts(n,2)==7

        AOTjui=AOTjui+1;

    elseif joursforts(n,2)==8

        AOTaou=AOTaou+1;

    elseif joursforts(n,2)==9

        AOTsep=AOTsep+1;

    elseif joursforts(n,2)==10

        AOToct=AOToct+1;

    elseif joursforts(n,2)==11

        AOTnov=AOTnov+1;

    elseif joursforts(n,2)==12

        AOTdec=AOTdec+1;

    end

end

    AOTf=[AOTjan;AOTfev;AOTmar;AOTavr;AOTmai;AOTjuin;AOTjui;AOTaou;AOTsep;AOToct;AOTnov;AOTdec];

    

AOTjan=0;AOTfev=0;AOTmar=0;AOTavr=0;AOTmai=0;AOTjuin=0;

AOTjui=0;AOTaou=0;AOTsep=0;AOToct=0;AOTnov=0;AOTdec=0;

    

for n=1:length(joursintenses)

    if joursintenses(n,2)==1

        AOTjan=AOTjan+1;

    elseif joursintenses(n,2)==2

        AOTfev=AOTfev+1;

    elseif joursintenses(n,2)==3

        AOTmar=AOTmar+1;

    elseif joursintenses(n,2)==4

        AOTavr=AOTavr+1;

    elseif joursintenses(n,2)==5

        AOTmai=AOTmai+1;

    elseif joursintenses(n,2)==6

        AOTjuin=AOTjuin+1;

    elseif joursintenses(n,2)==7

        AOTjui=AOTjui+1;

    elseif joursintenses(n,2)==8

        AOTaou=AOTaou+1;

    elseif joursintenses(n,2)==9

        AOTsep=AOTsep+1;

    elseif joursintenses(n,2)==10

        AOToct=AOToct+1;

    elseif joursintenses(n,2)==11

        AOTnov=AOTnov+1;

    elseif joursintenses(n,2)==12

        AOTdec=AOTdec+1;

    end

end

    

AOTi=[AOTjan;AOTfev;AOTmar;AOTavr;AOTmai;AOTjuin;AOTjui;AOTaou;AOTsep;AOToct;AOTnov;AOTdec];



AOT=[AOTf AOTi];



%Affichage des résultats

if str2double(fichier(1:2))==05  

figure(Fig2005)

subplot (3,2,p5)%,'Parent',hFig2005)

bar(AOT,'group')

title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')

%AJOUTER LES MOIS

set(gca,'xlim',[0 13],'xtick',1:12)

Xt = 1:12;

% Add the months as tick labels.

months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];

ax = axis; % Current axis limits

axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual

Yl = ax(3:4); % Y-axis limits

% Place the text labels

t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));

set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...

'Rotation',45);

% Remove the default labels

set(gca,'XTickLabel','')

h = legend('Ev.forts','Ev.intenses',1);

set(h,'Interpreter','none')
p5=p5+1;
end



if str2double(fichier(1:2))==06  

figure(Fig2006)
   
subplot (3,2,p6)%,'Parent',Fig2006)

bar(AOT,'group')

title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')

%AJOUTER LES MOIS

set(gca,'xlim',[0 13],'xtick',1:12)

Xt = 1:12;

% Add the months as tick labels.

months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];

ax = axis; % Current axis limits

axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual

Yl = ax(3:4); % Y-axis limits

% Place the text labels

t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));

set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...

'Rotation',45);

% Remove the default labels

set(gca,'XTickLabel','')

h = legend('Ev.forts','Ev.intenses',1);

set(h,'Interpreter','none')
p6=p6+1;  
end



if str2double(fichier(1:2))==07  

figure(Fig2007)
  

subplot (3,2,p7)%,'Parent',Fig2007)

bar(AOT,'group')

title([fichier(3:end-4) ' 20' fichier(1:2)],'fontsize',14,'fontweight','b','color','b')

%AJOUTER LES MOIS

set(gca,'xlim',[0 13],'xtick',1:12)

Xt = 1:12;

% Add the months as tick labels.

months = ['Jan ';'Fev ';'Mars';'Avr ';'Mai ';'Juin';'Juil';'Août';'Sep ';'Oct ';'Nov ';'Dec '];

ax = axis; % Current axis limits

axis(axis); % Set the axis limit modes (e.g. XLimMode) to manual

Yl = ax(3:4); % Y-axis limits

% Place the text labels

t = text(Xt,Yl(1)*ones(1,length(Xt)),months(1:12,:));

set(t,'HorizontalAlignment','right','VerticalAlignment','top', ...

'Rotation',45);

% Remove the default labels

set(gca,'XTickLabel','')

h = legend('Ev.forts','Ev.intenses',1);

set(h,'Interpreter','none')
p7=p7+1;  
end

end

merci à caro pour son aide précieuse