Discussion :

[Heimdall]

Bonjour,

J'ai un fichier binaire qui a ce format :

Code :

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

1
2
3
4
5
6
7
 
1 int
np float
np float
np float
np float
np float

"np float" étant un ceratin nombre de floats, ici np = 52428800.

J'ai fait le code suivant pour le lire, mais je ne sais pas pourquoi, ça ne marche pas, il semble pas bien lire les données....

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
 
 
 
#! /use/bin/env python
 
 
import numpy as np
import struct as st
 
filename = '/ou/mes/donnees/sont/fichier.dat'
nbpart   = 0
 
 
file = open(filename,"rb")
nbpart=np.fromfile(filename,dtype='int32',count=1)
 
print nbpart
 
 
x  = np.fromfile(file, np.float32, nbpart[0])
y  = np.fromfile(file, np.float32, nbpart[0])
vx = np.fromfile(file, np.float32, nbpart[0])
vy = np.fromfile(file, np.float32, nbpart[0])
vz = np.fromfile(file, np.float32, nbpart[0])
 
print x.min(), x.max()
print y.min(), y.max()
 
 
file.close()

voici ce que ça m'affiche :

Code :

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

1
2
3
4
5
 
me@here:~/code/visu/particules$ python derive_z.py
[52428800]
4.70197740329e-37 280.0
7.68104314375e-07 275.521087646

alors que ça devrait m'afficher :

Code :

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

1
2
3
4
 
[52428800]
1e-6 280.0
1e-6 40.0

Merci si vous pouvez m'aiguiller :-)

[Heimdall]

Hello,

j'ajoute deux petits codes de test qui me laissent touours dans l'incompréhension...

le code C qui génère un fichier binaire :

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
 
#include <stdio.h>
#include <stdlib.h>
 
 
 
int main(void)
{
float x[10] = {0,1,2,3,4,5,6,7,8,9};
int np;
FILE *fp;
 
np = 10;
 
fp = fopen("c2python.dat", "w");
 
fwrite (&np, sizeof(int), 1, fp);
fwrite (x, sizeof(float), 10, fp);
fwrite (x, sizeof(float), 10, fp);
fwrite (x, sizeof(float), 10, fp);
fwrite (x, sizeof(float), 10, fp);
 
fclose(fp);
 
return 0;
}

et le code python qui doit le lire :

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
 
#! /use/bin/env python
 
 
import numpy as np
import struct as st
 
filename = './c2python.dat'
nbpart   = 0
 
 
file = open(filename,"rb")
nbpart=np.fromfile(filename,dtype='int32',count=1)
 
print nbpart
 
 
x  = np.fromfile(file, np.float32, nbpart[0])
y  = np.fromfile(file, np.float32, nbpart[0])
vx = np.fromfile(file, np.float32, nbpart[0])
vy = np.fromfile(file, np.float32, nbpart[0])
 
 
print x.min(), x.max()
print y.min(), y.max()
 
 
file.close()

je pige vraiment pas pourquoi j'obtiens :

Code :

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

1
2
3
4
 
[10]
0.0 8.0
0.0 9.0

[Heimdall]

j'ai trouvé l'erreur, je la donne au cas où quelqu'un fasse la même un jour.

En regardant de près le code que j'ai donné, vous pourrez voir que j'ouvre le fichier 'filename', et qu'un objet 'file' est ainsi créé. La fonction numpy.fromfile() prend une string ou un objet fichier en argument.

Le hic étant que j'ai passé lors de mon premier appel à fromfile() la string filename., autrement dit fromfile a ouvert et fermé le fichier pour moi....

or les appels suivant donnent 'file' à fromfile(), du coup le pointeur dans ce fichier est à 0....


Voilà

[Matthieu Brucher]

Ah oui... Dans ce genre de cas, il faut pas hésiter à bien relire la doc