Discussion :

[Hallu33]

Bonjour,
je souhaite extraire de différentes figures que j'ai faites (en omettant de sauvegarder pour chaque figure les données) l'histogramme qui y est tracé afin de calculer l'écart type. J'ai essayé différentes méthodes, sans succès :

Code :

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

1
2
3
4
s=hgload('mafigure.fig')
h = findobj(s,'Type','hggroup')
y=get(h,'ydata');
x=get(h,'xdata');

Code :

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

1
2
3
4
s=hgload('mafigure.fig')
h = findobj(s,'Type','line')
y=get(h,'ydata');
x=get(h,'xdata');

Code :

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

1
2
3
h=findall(gca,'type','line');
xx=get(h,'xdata');
yy=get(h,'ydata');

Je précise que je me place évidemment dans le bon répertoire de la figure avec la commande 'cd' mais à chaque fois, h est vide... Est-ce le 'type' qui n'est pas bon ? Que dois-je mettre ? Merci d'avance de votre aide.

[Jerome Briot]

Quelle commande as-tu utilisée pour dessiner les histogramme ?
Sur quelle version de MATLAB travailles-tu ?

[Hallu33]

Je suis sous la R2008a, et j'ai tracé les histos avec la commande histo(y,x), y étant une série de valeurs, x étant la discrétisation voulue de ces valeurs. histo va trier les valeurs y, et chaque fois qu'une de ces valeurs est comprise entre (x(i)-x(i-1))/2 et (x(i)-x(i+1))/2 il ajoute 1 à la barre de x(i), on obtient donc ensuite un histogramme. Je veux extraire les données des figures pour ne pas avoir à tous les retracer pour avoir l'écart type de chacun.

[Jerome Briot]

la commande histo(y,x)

Tu parles plutôt de HIST, non ?

Si c'est bien le cas, essaie comme ceci :

Code :

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

1
2
3
4
5
6
7
8
s = hgload('mafigure.fig');
h = findobj(s,'Type','patch');
 
x = get(h,'xdata');
x = mean(x);
 
y = get(h,'ydata');
y = mean(y(2:3,:));

[Hallu33]

Voilà c'était le 'patch' qui me manquait merci, entretemps j'avais trouvé une autre méthode :

Code :

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

1
2
3
4
5
6
axs = get(gcf, 'Children');
pos = get(axs(1), 'Children');
get(pos);
X = get(pos, 'XData');
Y = get(pos, 'YData');
histo=Y(2,:);

Mais j'avoue que j'ai du mal avec les X et Y obtenus, ce sont des matrices de 4 lignes... Autant dans X, on comprend qu'il y ait 2 lignes différentes (l'une est la position du début de barres l'autre la position du milieu des barres de l'histogramme), mais pourquoi ces 2 lignes sont-elles dupliquées pour en donner 4 ? Pour Y, on a deux lignes de 0 et deux lignes étant l'histogramme dupliqué... Etrange. Enfin merci de ton aide en tout cas.
Sinon autre problème rencontré : hgload n'accepte pas les "-" dans le nom des fichiers... Du coup j'ai du me débrouiller avec "open" qui lui les accepte... Je pensais que ce genre d'idioties de matlab était terminé en voyant qu'il acceptait enfin les espaces dans la commande 'cd' mais on dirait que non...