Discussion :

[zoro_le_renard]

Bonjour

Toujours dans la même veine que mon précédent post, je dois lire une matrice faite sous Unix : matrice 1536*1536*2 de float

Elle contient 1536*1536 pixels et 2 variables à lire pour chaque pixel. J'ai essayé le bout de code suivant en supposant que mes variables sont stockées sur 4 octets mais sans succès

Code :

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

1
2
3
4
5
6
DIMENSION GRIDx(2359296),GRIDy(2359296)
 
INTEGER*4 GRIDx,GRIDy
	open(1,file='fich',access='direct',form='unformatted',recl=18874368)
	read(1,rec=1)GRIDx,GRIDy
	close(1)

Je suis toujours aussi débile en Fortran et j'ai essayé de m'inspirer du code de Da Zumba, je devrais peut-être déclarer du float plutôt que de l'integer. J'avoue que j'ai besoin d'aide ...

Merci

Zoro

[DaZumba]

Elle contient 1536*1536 pixels et 2 variables à lire pour chaque pixel. J'ai essayé le bout de code suivant en supposant que mes variables sont stockées sur 4 octets mais sans succès

Quel est le probleme exactement ? Message d'erreur ? Plantage ? Valeurs aberrantes?

D'abord il faut que tu saches le type de donnees contenu dans le fichier (float? integer?). Si c'est un float, il faut declarer gridx et gridy comme real. Sinon, c'est effectivement integer*4.

Ensuite, il faut faire attention: tu dois lire 1536*1536*2 variables de 4 octets chacune. Ca fait 18 Mo --- c'est un gros tableau, et il est possible qu'il depasse les capacites de la pile de ta machine. Dans ce cas, ton programme plante au demarrage en emettant une segmentation fault.

Egalement, il y a le recl. La longueur d'un enregistrement dans ton cas est 1536*1536*4, soit 9437180 (et pas le double comme tu le mets, puisque tu lis les deux grilles separement).
Pour faire plus clair, je ferais comme cela (Fortran90):

Code :

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

1
2
3
4
5
6
7
 
integer*4, dimension (1536*1536) :: gridx, gridy
 
open(unit=1, file='fich', access='direct', form='unformatted', recl=1536*1536*4)
read(1, rec=1) gridx
read(1, rec=2) gridy
close(1)

Enfin, si le fichier binaire a ete cree sous une machine Unix (certainement en big-endian) et que tu travailles sur une station a base de x86 (little-endian), alors tu vas devoir inverser les octets de chaque enregistrement.

[zoro_le_renard]

Le problème ce sont les valeurs aberrantes obtenues (en théorie coordonnées en m ou hm ou km incrémentées de +- 1000 m).

Le type de données est du float.

J'ai donc essayé ton code avec la déclaration float :

Code :

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

1
2
3
4
5
6
DIMENSION GRIDx(2359296),GRIDy(2359296)
real*4 GRIDx,GRIDy
open(1,file='MATRICE_1536_LAMBERT',access='direct',form='binary',recl=9437184)
read(1,rec=1)GRIDx
read(1,rec=2)GRIDy	
close(1)

De nouveau j'ai des valeurs aberrantes.

Concernant les endians, je pense que les données ont été générées sur une station unix. De mon côté je travaille sur un PC mais je ne sais pas comment y sont gérés les endians.

Merci pour ton aide

[DaZumba]

Concernant les endians, je pense que les données ont été générées sur une station unix. De mon côté je travaille sur un PC mais je ne sais pas comment y sont gérés les endians.

Si les valeurs sont aberrantes, alors c'est surement un probleme d'endianness. Tu peux utiliser la routine donnee en bas de page ici. Pour passer d'une endianness a l'autre (c'est du fortran 90 - je crains qu'en fortran 77 mvbits() n'existe pas).

[zoro_le_renard]

Je viens d'y jeter un coup d'oeil rapide. Ca a l'air de marcher . Confirmation + mon bout de code si tout est OK ds le courant de la semaine prochaine.

Encore une fois merci