Lecture de fichier avec architecture particulière

Version imprimable

Bonjour,
Je souhaite lire un fichier de résultats issu d'un logiciel de calcul seulement l'architecture est assez particulière:

<list><x row="631" col="1" list="-94.847176|
-93.62446|-92.40155|-91.17844|-89.955124|-88.731606|
-87.50789|-86.283966|-85.05985|-83.83553|-82.61101|
-81.38629|-80.16137|-78.93625|-77.71092|-76.4854|
-75.259674|-74.03375|-72.80763|-71.58131|-70.35478|
-69.12806|-67.90113|-66.67401|-65.446686|-64.21916|
-62.991436|-61.76351|-60.53539|-59.307064|-58.078545|
-56.84982|-55.6209|-54.391777|-53.16246|
607.81995|609.0398|610.2598|611.4801|612.7005|
613.9212|615.142|616.3631|617.58435|618.80585|
620.02747|621.2493|622.47144|623.69366|624.91614"/><y row="631" col="1" list="1.8219452|
1.8576648|1.8948541|1.9335837|1.9739263|2.0159614|
2.0597715|2.1054447|2.1530752|2.2027586|2.2546024|
2.3087146|2.3652134|2.4242232|2.485872|2.5503006|
2.6176577|2.688096|2.7617812|2.8388884|2.9196048|
3.0041225|3.0926516|3.1854088|3.282627|3.384555|
1.6867067|1.7169546|1.74841|1.7811378|1.815192"/></list></data></plot><plot id="25002" label1="hioqshihqsihioqottype="xy" sn="temp" xaxtype="engine-cyclic" xcycle="4-stroke" xlabel="sdvsv="61" ylabel="Temperature [UNIT]" yunit="32" zunit="0">
<data datatype="xy" sdg2" name="cylàiàvq" typesdff="32" zunit="0">
<list><x row="631" col="1" list="-94.847176|
-93.62446|-92.40155|-91.17844|-89.955124|-88.731606|
-87.50789|-86.283966|-85.05985|-83.83553|-82.61101|
<list><x row="631" col="1" list="-94.847176|
-93.62446|-92.40155|-91.17844|-89.955124|-88.731606|
-87.50789|-86.283966|-85.05985|-83.83553|-82.61101|
-81.38629|-80.16137|-78.93625|-77.71092|-76.4854|
-75.259674|-74.03375|-72.80763|-71.58131|-70.35478|
-69.12806|-67.90113|-66.67401|-65.446686|-64.21916|
-62.991436|-61.76351|-60.53539|-59.307064|-58.078545|
-56.84982|-55.6209|-54.391777|-53.16246|
607.81995|609.0398|610.2598|611.4801|612.7005|
613.9212|615.142|616.3631|617.58435|618.80585|
620.02747|621.2493|622.47144|623.69366|624.91614"/><y row="631" col="1" list="1.8219452|
1.8576648|1.8948541|1.9335837|1.9739263|2.0159614|
2.0597715|2.1054447|2.1530752|2.2027586|2.2546024|
2.3087146|2.3652134|2.4242232|2.485872|2.5503006|
2.6176577|2.688096|2.7617812|2.8388884|2.9196048|
3.0041225|3.0926516|3.1854088|3.282627|3.384555|
1.6867067|1.7169546|1.74841|1.7811378|1.815192"/></list></data></plot><plot id="25002" label1="hioqshihqsihioqottype="xy" sn="temp" xaxtype="engine-cyclic" xcycle="4-stroke" xlabel="sdvsv="61" ylabel="Temperature [UNIT]" yunit="32" zunit="0">
<data datatype="xy" sdg2" name="cylàiàvq" typesdff="32" zunit="0">
<list><x row="631" col="1" list="-94.847176|
-93.62446|-92.40155|-91.17844|-89.955124|-88.731606|
-87.50789|-86.283966|-85.05985|-83.83553|-82.61101|

Le format est quasi identique au précédent. Ce qui est récurent c'est que les nombres que je souhaite extraire sont toujours après le string list=".
Sinon il y très peu de récurrence. Je ne sais donc pas trop comment m'y prendre.

Merci d'avance de votre aide.

05/03/2014, 12h18
Jerome Briot

Ceci ressemble à un format XML.

Tu peux utiliser la fonction de lecture xmlread (voir la fonction parseXML données dans la documentation) pour lire les fichiers puis extraire les valeurs numériques des champs XML lus

[Edit] En y regardant de plus près, il y a des soucis de balisage pour la lecture XML
Pourrais-tu nous dire quel logiciel génère ces fichiers ?
Pourrais-tu nous transmettre un de ces fichiers (dans une archive zip ou rat) ?
05/03/2014, 13h58
odjapo

1 pièce(s) jointe(s)

Merci pour votre réponse, j'ai mis un fichier en pièce jointe.

Voici une idée :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
xDoc = xmlread('exemple.gu');
 
allx = xDoc.getElementsByTagName('x');
ally = xDoc.getElementsByTagName('y');
 
for n = 0:allx.getLength-1
 
    x_txt{n+1} = char(allx.item(n).getAttribute('list'));
    y_txt{n+1} = char(ally.item(n).getAttribute('list'));
 
end
 
x_txt = strrep(x_txt, '|', ' ');
y_txt = strrep(y_txt, '|', ' ');
 
x = cellfun(@str2num, x_txt, 'UniformOutput', false);
y = cellfun(@str2num, y_txt, 'UniformOutput', false);

Les valeurs numériques sont au final stockées dans deux tableaux de cellules

Par curiosité, quel logiciel génère ces fichiers ?

05/03/2014, 18h59
odjapo

Merci cela marche plutôt bien.
Le logiciel est GT-Power un logiciel de simulation numérique.

Est-il également possible de sortir de la même façon les "string" en gras ci-dessous (L259 et L513) pour chaque série de donnée ?:

L259
1.6867067|1.7169546|1.74841|1.7811378|1.815192"/></list></data></plot><plot id="25002" label1="EngCylinder part cylinder-01" name="Temperature" plottype="xy" sn="temp" xaxtype="engine-cyclic" xcycle="4-stroke" xlabel="Crank Angle [UNIT]" xunit="61" ylabel="Temperature [UNIT]" yunit="32" zunit="0">

L260
<data datatype="xy" dbname="" filetype="2" name="cylinder-01" type="exp" xunit="61" yunit="32" zunit="0">

L513
411.51013|413.35733|415.2496|417.1884|419.1745"/></list></data></plot><plot id="25003" label1="EngCylinder part cylinder-01" name="Volume" plottype="xy" sn="vol" xaxtype="engine-cyclic" xcycle="4-stroke" xlabel="Crank Angle [UNIT]" xunit="61" ylabel="Volume [UNIT]" yunit="17" zunit="0">
L514
<data datatype="xy" dbname="" filetype="2" name="cylinder-01" type="exp" xunit="61" yunit="17" zunit="0">

Merci de votre aide :lol:
05/03/2014, 19h05
Jerome Briot

Il suffit de reprendre exactement la même idée en changeant simplement les arguments d'entrée de getElementsByTagName et getAttribute

Essaie et tiens nous au courant ;)

J'ai essayé cela:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
xDoc = xmlread('exemple.gu');
 
allx = xDoc.getElementsByTagName('x');
ally = xDoc.getElementsByTagName('y');
allname = xDoc.getElementsByTagName('name');
 
for n = 0:allx.getLength-1
 
    x_txt{n+1} = char(allx.item(n).getAttribute('list'));
    y_txt{n+1} = char(ally.item(n).getAttribute('list'));
    name_txt{n+1} = char(ally.item(n).getAttribute('name'));
 
end

Mais cela ne fonctionne pas :(

C'est parce que tu n'as pas le même nombre d'entité x et plot
Code:

1 2 3 4 5 6 7 xDoc = xmlread('exemple.gu'); allx = xDoc.getElementsByTagName('x'); allx.getLength allplot = xDoc.getElementsByTagName('plot'); allplot.getLength
donne :
Code:

1 2 3 4 5 6 7 8 9 ans = 39 ans = 35
Il te faut 2 boucles for-end indépendantes

J'ai donc essayé ça:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
xDoc = xmlread('exemple.gu');
 
allx = xDoc.getElementsByTagName('x');
ally = xDoc.getElementsByTagName('y');
allname = xDoc.getElementsByTagName('name');
 
for n = 0:allx.getLength-1
 
    x_txt{n+1} = char(allx.item(n).getAttribute('list'));
    y_txt{n+1} = char(ally.item(n).getAttribute('list'));
 
end
for n = 0:allname.getLength-1
    name_txt{n+1} = char(allname.item(n).getAttribute('name'));
end

Mais cela ne fonctionne toujours pas :(

05/03/2014, 19h43
Invité
Bonsoir,

'name' est un attribut des éléments 'plot', il faut donc récupérer d'abord ces derniers éléments, puis leur attribut 'name'.
Comme ceci :
Code:

1 2 3 4 5 allplot = xDoc.getElementsByTagName('plot'); for n = 0:allplot.getLength-1 name_txt{n+1} = char(allplot.item(n).getAttribute('name')); end

J'ai trouvé:

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
xDoc = xmlread('exemple.gu');
 
allx = xDoc.getElementsByTagName('x');
ally = xDoc.getElementsByTagName('y');
allplot = xDoc.getElementsByTagName('plot');
 
for n = 0:allx.getLength-1
 
    x_txt{n+1} = char(allx.item(n).getAttribute('list'));
    y_txt{n+1} = char(ally.item(n).getAttribute('list'));
 
end
for n = 0:allplot.getLength-1
    name_txt{n+1} = char(allplot.item(n).getAttribute('name'));
end

Merci :lol:

05/03/2014, 20h12
odjapo

Encore une question. Je cherche maintenant à lier les vecteurs x et y avec leurs "plot" respectif. En effet il arrive que des graphiques comportent 2 ou 3 matrices mais je ne sais pas comment savoir quels graphiques comportent plusieurs matrices. (C'est l'explication du fait qu'il y ai 35 éléments "plot" et 39 éléments "x" et "y").

Je pense que vous l'avez compris le but est de retracer les graphiques enregistrés.

Merci encore de votre aide.

Ayant la hiérarchie suivante :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
<plot ...>
    <data ...>
        <list ...>
            <x ...>
            <y ...>
        </list>
    </data>
    <data ...>
        <list ...>
            <x ...>
            <y ...>
        </list>
    </data>
    ...
</plot>
<plot ...>
    ...
</plot

On peut faire comme suit :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
xDoc = xmlread('exemple.gu');
 
allplot = xDoc.getElementsByTagName('plot');
 
for n = 0:allplot.getLength-1
    plotElement = allplot.item(n);
    name = char(plotElement.getAttribute('name'));
 
    allDataOfElement = plotElement.getElementsByTagName('data');
    for m = 0:allDataOfElement.getLength-1
        dataElement = allDataOfElement.item(m);
        xElement = dataElement.getElementsByTagName('x').item(0);
        yElement = dataElement.getElementsByTagName('y').item(0);
 
        x = char(xElement.getAttribute('list'));
        y = char(yElement.getAttribute('list'));
 
    end
end

Les plots étant eux-mêmes englobés par des éléments <group> ... </group>, tu peux procéder de la même façon.

05/03/2014, 20h57
odjapo

Merci beaucoup pour ce code seulement, je ne peux plus utiliser la fonction:

Code:

x = cellfun(@str2num, x_txt, 'UniformOutput', false);

Car je dois extraire les "double" présent dans un "char".

str2num ne convient pas ici, les valeurs sont séparées par des '|', or cette fonction ne prend que des espaces ou virgules.
Tu peux par contre utiliser sscanf ou textscan :
Code:

1 2 3 4 5 6 7 8 9 10 >> sscanf('516.91156|518.1161|519.32086|520.52576|521.7309', '%f|') ans = 516.9116 518.1161 519.3209 520.5258 521.7309
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 >> textscan('516.91156|518.1161|519.32086|520.52576|521.7309|', '%n', 'delimiter', '|') ans = [5x1 double] >> ans{1} ans = 516.9116 518.1161 519.3209 520.5258 521.7309

05/03/2014, 23h45
odjapo

Merci beaucoup pour toute votre aide.:P

Bonjour, j'ai de nouveau deux questions pour mon code:

Code:

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
 
clear all
close all
[FileName,PathName] = uigetfile('*.gu','Sélectionnez un fichier GT-Power', 'MultiSelect', 'Off');
xDoc = xmlread(fullfile(PathName,FileName));
 
allplot = xDoc.getElementsByTagName('plot');
alldata = xDoc.getElementsByTagName('data');
 
temp=0;
for n = 0:allplot.getLength-1
    plotElement = allplot.item(n);
    name = char(plotElement.getAttribute('name'));
    xlegend = char(plotElement.getAttribute('xlabel'));
    ylegend = char(plotElement.getAttribute('ylabel'));
    group = char(plotElement.getAttribute('label1'));
    allDataOfElement = plotElement.getElementsByTagName('data');
    figure('Name',name,'NumberTitle','off') 
    for m = 0:allDataOfElement.getLength-1
        data = allDataOfElement.item(m);
        xElement = data.getElementsByTagName('x').item(0);
        yElement = data.getElementsByTagName('y').item(0);
 
        x = char(xElement.getAttribute('list'));
        y = char(yElement.getAttribute('list'));
 
        x2 =sscanf(x, '%f|');
        y2 =sscanf(y, '%f|');
 
        hold all
        p(m+1) = plot(x2,y2);
        datname = char(alldata.item(n+temp).getAttribute('name'));
        leg{m+1}=datname;
        if (m~=0)
            datname = char(alldata.item(n+m+temp).getAttribute('name'));
            leg{m+1}=sprintf(datname,m+1);
        end
        hold on
    end
    temp=temp+m;
    title(name)
    xlabel(xlegend);
    ylabel(ylegend);
    legend(p,leg)
    clearvars leg
    grid on
end

1) J'ai un message d'erreur pour mes legendes:

Warning: Ignoring extra legend entries.
> In legend at 286
In dep_gt at 44

2) Jaimerais ranger mes figures par groupes de figures qui sont fonction de la variable group dans mon code mais je ne sais pas comment faire.

Merci d'avance de votre aide. ;)

06/03/2014, 12h57
Jerome Briot

Citation:

Envoyé par odjapo

1) J'ai un message d'erreur pour mes legendes

Il faut aussi nettoyer la variable p :

Code:

clearvars leg p

Citation:

Envoyé par odjapo

2) Jaimerais ranger mes figures par groupes de figures qui sont fonction de la variable group

C'est à dire ?

Veux-tu une fenêtre par groupe avec un seul objet Axes et toutes les courbes du groupe superposées ?
Ou bien une fenêtre par groupe avec autant d'objets Axes (via subplot) que de courbes du groupe ?
06/03/2014, 13h16
odjapo

Pour les groupes je pensais à une fenêtre graphique par groupe et avoir la possibilité (en cochant une case ou en selectionnant dans un menu déroulant), d'afficher la série xy de mon choix.
Exemple afficher la série "Temperature" du groupe "EngCylinder part cylinder-01".