numpy fichier binaire

Version imprimable

Bonjour,

J'ai un fichier binaire qui a ce format :
Code:

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:

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:

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:

1 2 3 4 [52428800] 1e-6 280.0 1e-6 40.0
Merci si vous pouvez m'aiguiller :-)

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:

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:

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:

25/05/2011, 08h10
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à :)
25/05/2011, 13h54
Matthieu Brucher

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