Discussion :

[nicolasf85]

Bonjour,

j'ai créé une boucle de calcul:
6 fichiers XYZ présents dans un dossier sont chargés (un par un), 6 calculs sont fait sur la base de ces 6 fichiers, 1 calcul pour chaque fichier chargé.
Et donc 6 résultats différents.

Ma boucle s'exécute donc 6 fois.
Maintenant je voudrais que les résultats (résultats statistiques MIN, MAX, MOY) soient concaténés dans une seule et même matrice.
Malheureusement je n'y parviens pas.

2 choses que j'essaie donc de coder:

a) A chaque boucle je voudrais créer une matrice (1 ligne x 4 colonnes) contenant 4 colonnes contenant : 1.le nom du fichier XYZ importé au début de la boucle, 2.MIN des résultats, 3.MAX, 4.MOY)

Sur le code ci-dessous, j'essaie de faire une matrice "data.result_line" qui regroupe les variables : nom du fichier XYZ actuellement chargé dans la boucle, l'écart type, la moyenne, le min, le max.
Apparemment cela ne fonctionne pas, le résultat est une matrice 1x1 qui contient uniquement la première variable : nom du fichier XYZ actuellement chargé.


b) Je voudrais qu'à la fin de chaque boucle, une matrice récupère cette matrice (1x4). Mais je veux pas que les matrices résultant des 6 boucles s'écrasent les unes après les autres (et ne reste plus que la dernière à la fin).
Je veux que les 6 lignes de résultats se suivent dans une seule et même matrice.
Au final j'aurais une matrice de 6 lignes et 4 colonnes (1 ligne supplémentaire par boucle effectuée).


Merci d'avance pour votre aide.
voici le 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
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
 
 
 function import_plan(objet,event)
data=guidata(gcbf);
 
plan_load = uigetfile('*','Importer 3 pts du plan');
 
fid = fopen(plan_load,'r');
plan_import = textscan(fid,'%*s %f %f %f %f %f %f %f %f %f','delimiter',',','headerLines',0); 
fclose(fid);
 
for i=1:size(plan_import,2)
    data.plan(:,i) = plan_import{1,i}(:,1);
end
 
guidata(gcbf,data)
 
 
 
 
 
 function import_nuage(objet,event)
data=guidata(gcbf);
 
 
cd C:\MATLAB\points
str = dir('*.xyz');
 
%boucle qui permet de traiter chaque fichier un après l'autre
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
for k = 1:length(str)
    str_actuel = str(k).name;
    str_eval = ['load ', str_actuel];
    eval(str_eval);
 
 
 
 
%Calculs
%%%%%%%%%
 
P1 = [data.plan(1,1),data.plan(1,2),data.plan(1,3)]
P2 = [data.plan(1,4),data.plan(1,5),data.plan(1,6)]
P3 = [data.plan(1,7),data.plan(1,8),data.plan(1,9)]
 
H = load(str_actuel);
 
    %   1)Calculer le vecteur normal:
       N = cross(P1-P2, P1-P3);
    %   2)Extraire les valeurs a, b et c
       a=N(1,1) ; b=N(1,2) ; c=N(1,3);
    %   3)Calculer la valeur d
       d = -(a*P1(1,1)+b*P1(1,2)+c*P1(1,3));
    %   4)Calculer la distance/ecart entre le point H et le plan P [P1,P2,P3]
       data.e = ((a*H(:,1)+b*H(:,2)+c*H(:,3)+d))/sqrt(a^2+b^2+c^2);
    %   5)Calcul des éléments STATISTIQUES
       data.maximum = max(data.e);
       data.minimum = min(data.e);
       data.moyenne = mean(data.e);
       data.ecart_type = std(data.e);
 
 
 
% Création d'une matrice qui synthétise les résultats
%%%%%%%%%
 
       data.result_line = [ str_actuel data.moyenne data.ecart_type data.maximum data.minimum];
 
 
 
end;
 
guidata(gcbf,data);

[ol9245]

Code :

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

1
2
3
4
5
6
Result_Final = [] ;
for i=1:6
    result_etape = [i,  i*i, 2*i, 0]  ; 
    Result_Final = [ Result_Final ; result_etape] ;
end 
disp(Result_Final)

[Jerome Briot]

Code :

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

1
2
3
4
5
6
7
8
 
%boucle qui permet de traiter chaque fichier un après l'autre
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
for k = 1:length(str)
    str_actuel = str(k).name;
    str_eval = ['load ', str_actuel];
    eval(str_eval);

Pas besoin d'EVAL =>

Code :

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

1
2
3
4
5
6
 
%boucle qui permet de traiter chaque fichier un après l'autre
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
for k = 1:length(str)
    load(str(k).name);

Code :

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

P1 = [data.plan(1,1),data.plan(1,2),data.plan(1,3)]

s'écrit directement :

Code :

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

P1 = data.plan(1,1:3);

[Jerome Briot]

Code :

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

1
2
3
 
cd C:\MATLAB\points
str = dir('*.xyz');

Pas besoin de se déplacer dans le dossier :

Code :

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

1
2
 
str = dir('C:\MATLAB\points\*.xyz');

ou encore :

Code :

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

1
2
pathname = 'C:\MATLAB\points';
str = dir(fullfile(pathname,'*.xyz'));

[nicolasf85]

pourquoi est-ce que je n'arrive pas à créer une matrice formée de texte et de chiffres:

Code :

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

1
2
 
result_line = [ nuage_name moyenne ecart_type maximum minimum ]

sachant que
nuage_name = nuage_de_points.xyz
moyenne = 1234
ecart_type = 3345
maximum = 85443
minimum = 34599

?? impossible de créer une matrice alpha numérique ?

[Jerome Briot]

pourquoi est-ce que je n'arrive pas à créer une matrice formée de texte et de chiffres:

Parce que c'est impossible sous MATLAB.

Au mieux tu peux utiliser un tableau de cellules ou une structure

=> :

• Les tableaux de cellules (cell array)
• Les structures

[nicolasf85]

est-ce que une structure c'est exploitable ensuite ?

j'envisage de
- exporter la matrice complète de la structure dans un fichier ascii, éventuellement excel.
- générer des graphiques, sans perdre le lien entre le nom du nuage de points et ses résultats statistiques propres (min, max, etc.). Car sinon c'est n'importe quoi et je sais plus à quoi se reporte telle ou telle courbe.

[Jerome Briot]

est-ce que une structure c'est exploitable ensuite ?

Non non ça sert à rien... c'est juste la pour faire joli (à moins que ce ne soit pour justifier le prix de MATLAB)

Plus sérieusement, depuis le début tu utilises la structure data (et même str), tu sais donc que ce type de variable est très utile

- exporter la matrice complète de la structure dans un fichier ascii, éventuellement excel.

Pourquoi passer par Excel ?
Quel est le problème avec MATLAB ?

- générer des graphiques, sans perdre le lien entre le nom du nuage de points et ses
résultats statistiques propres (min, max, etc.).

Il suffit d'avoir par exemple une structure res un champs name et un autre data :

Code :

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

1
2
3
4
for n = 1:...
   res(n).name = nuage_name;
   res(n).data = [moyenne ecart_type maximum minimum];
end

Enfin... c'est l'idée

[nicolasf85]

MERCI pour ces astuces !!!

Code :

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

1
2
pathname = 'C:\MATLAB\points';
str = dir(fullfile(pathname,'*.xyz'));

J'aimerais que ce soit l'utilisateur qui spécifie le pathname (le chemin où se trouvent les fichiers).
J'ai créé un TextEdit
maintenant je ne parviens pas à rapatrier le contenu du TextEdit pour en faire une variable que je nome "pathname".

voici le 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
25
26
function chemin(objet,event)
data=guidata(gcbf);
 
%quelque chose à mettre ici ???
guidata(gcbf,data);
 
 
 
 
function import_nuage(objet,event)
data=guidata(gcbf);
 
 
%Récupérer le nom du répertoire
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
pathname = get(chemin,'string');
 
str = dir(fullfile(pathname,'*.xyz'))
 
 
%%%%%cd C:\MATLAB\points
%%%%%str = dir('*.xyz');
 
 
guidata(gcbf,data);

ERREUR:

Code :

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

1
2
3
4
Error in ==> importation_sansGUIDE_version3>import_nuage at 118
pathname = get(data.chemin,'string');
 
??? Error while evaluating uicontrol Callback

[Jerome Briot]

maintenant je ne parviens pas à rapatrier le contenu du TextEdit pour en faire une variable que je nome "pathname".

Si tu as stocké l'identifiant de l'objet Edit dans la structure data, il suffit de faire :

Code :

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

chemin = get(data.<...>,'string');

Et sinon, tu peux aussi utiliser UIGETDIR...

[nicolasf85]

OK !

j'ai ce code:

Code :

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

1
2
   res(k,1).nuage_nom = str_actuel;
   res(k,1).nuage_res = [moyenne ecart_type maximum minimum];

je souhaiterais exporter le contenu de cette structure dans un fichier Texte ASCII avec séparateur (tab ou virgule), l'aspet du fichier texte doit être formé ainsi: [nuage_nom] [nuage_res]
ce qui donne :

Code :

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

1
2
3
4
5
 
nuage1   ecart_type1    maximum1    minimum1
nuage2   ecart_type2    maximum2    minimum2
nuage3   ecart_type3    maximum3    minimum3
...

J'ai essayé ça:

Code :

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

save('resultats.txt', '-struct', 'res')

Erreur en retour:

Code :

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

1
2
??? Error using ==> save
The argument to -STRUCT must be the name of a scalar structure variable.

Et si j'essaie:

Code :

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

res_noms = res.nuage_nom

ou encore

Code :

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

res_noms = res(:,1).nuage_nom

il conserve uniquement le premier nom de la liste

pourtant quand je questionne via la commande window j'obtient ça:

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
>> res(:,1).nuage_nom
 
ans =
pil05_bb_hz0.xyz
 
ans =
pil05_br_hz0.xyz
 
ans =
pil05_b_hz0.xyz
 
ans =
pil05_bm_hz0.xyz
 
ans =
pil05_lr_hz0.xyz
 
ans =
pil05_pvc_hz0.xyz
 
ans =
pil05_sim_hz0.xyz

une idée ?

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
fid = fopen('resultats.txt','wt');
   for k = 1:...
       fprintf(fid,'%s,%.4f,%.4f,%.4f,%.4f\n',res(k).nuage_nom, res(k).nuage_res);
   end
fclose(fid);

[nicolasf85]

le fichier resultat.txt est créé, mais il est vide.

que veut dire le wt dans

Code :

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

fid = fopen('resultats.txt','wt')

?

en fait pour bien comprendre ton code,

Code :

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

fid = fopen('resultats.txt','wt');

il créer un fichier txt

Code :

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

1
2
3
   for k = 1:...
       fprintf(fid,'%s,%.4f,%.4f,%.4f,%.4f\n',res(k).nuage_nom, res(k).nuage_res);
   end

inscrit les résultats les uns après les autres ?

Code :

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

fclose(fid);

finalise le fichier txt ?

J'ai déjà une boucle de calculs qui génère à chaque tour un résultat sauvegardé dans une structure. Maintenant tout ce que j'ai besoin de faire c'est de récupérer ces données stockées dans la structure et les mettre dans un fichier txt.

[nicolasf85]

on y est presque...

Code :

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

1
2
3
fid = fopen('resultats.txt','wt');
fprintf(fid, '%s %8.6f %8.6f %8.6f\n', res.nuage_nom , res.nuage_res);
fclose(fid);

le résultat est un fichier texte avec toutes les données dedant, mais pas structurées comme voulu.

en fait je me retrouve avec les premières lignes contenant les valeurs de res.nuage_nom, avec en plus de ça des valeurs à droit de chaque nom (je sais pas du tout ce qu'elles représentent).

et les lignes suivantes correspond au contenu de res.nuage_res.

je ne parviens pas à les mettre l'un à coté de l'autre.

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
fid = fopen('resultats.txt','wt');
   for k = 1:numel(res)
       fprintf(fid,'%s,%.4f,%.4f,%.4f,%.4f\n',res(k).nuage_nom, res(k).nuage_res);
   end
fclose(fid);

Sinon, il est préférable d'utiliser une structure de tableaux plutôt qu'un tableau de structures :

Plutôt que :

Code :

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

1
2
3
4
5
res(1).nuage_nom = 'AAA';
res(2).nuage_nom = 'BB';
 
res(1).nuage_res = [0 1 2 3];
res(2).nuage_res = [4 5 6 7];

tu fais :

Code :

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

1
2
3
4
5
res.nuage_nom{1} = 'AAA';
res.nuage_nom{2} = 'BB';
 
res.nuage_res(1,1:4) = [0 1 2 3];
res.nuage_res(2,1:4) = [4 5 6 7];

ce qui donne pour l'écriture :

Code :

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

1
2
3
4
5
fid = fopen('resultats.txt','wt');
   for k = 1:numel(res.nuage_nom)
       fprintf(fid,'%s,%.4f,%.4f,%.4f,%.4f\n',res.nuage_nom{k}, res.nuage_res(k,:));
   end
fclose(fid);

[Jerome Briot]

que veut dire le wt dans

Code :

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

fid = fopen('resultats.txt','wt')

?

Par défaut, un fichier est de type binaire.

Si tu veux créer un fichier text (qui est un fichier binaire particulier), il faut spécifier l'argument t dans FOPEN pour que les commandes suivantes fonctionnent normalement.

Pour le reste, tu as tout compris

[nicolasf85]

Merci beaucoup !!! c'est excellent !