Discussion :

[bmayer]

Bonjour

J'essaie de lire des fichiers grib à l'aide de la bibliothèque "pygrib"
Autant avec certaines source de grib, j'accède aux données, de façon élémentaire, autant pour d'autres sources, c'est plus problématique.
Lecture OK : ceux que je charge avec xygrib Les infos de vent et de température sont accessibles, lisibles.
Lecture KO : ceux disponible auprès de weather4d (http://grib.weather4d.com/AromeHD/) par exemple. Autant un outil de visualisation ouvre et interprète correctement le fichier (vent et temperature sont visibles) autant le parcours par la boucle de code ne révèle pas les informations attendues (vent et température)
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
import pygrib, sys
#grbs = pygrib.open("/home/bernard/MTO_gribs/20230308_213001_GFS_P25_.grb2")
#grbs = pygrib.open("/home/bernard/MTO_gribs/20230313_195809_ARPEGE_P50_.grb2")
#grbs = pygrib.open("/home/bernard/MTO_gribs/AromeHD_Atlantique_lastest.grb")
grbs = pygrib.open("/home/bernard/MTO_gribs/Brest_4km_WRF_WAM_230313-00.grb")
grbs.seek(0) # grbs.rewind()
 
for grb in grbs :
    #print("--- grb -----------------------------------------------------")
    #print(repr(grb))
    print("msg number :", grb.messagenumber)
    #print("packing :", grb.packingType)
    print("parameterName :", grb.parameterName)
    print("validDate :", grb.validDate)
   print(grb)
 
grbs.close()
sys.exit()

Quelqu'un pourrait-il éclairer ma lanterne ?
Merci d'avance

[papajoker]

bonjour

Tu peux afficher toutes les datas dans ta boucle (je suppose qu'alors, tu trouveras ce que tu désires)

Code :

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

1
2
3
4
5
6
7
8
for i, grb in enumerate(grbs) :
    print("-"*30)
    for key in grb.keys():
        try:
            print(f"{key:30}", grb[key])
        except:
            print(f"{key:30}", "!!!!!!!!! ERREUR")
    if i > 10: break

[bmayer]

Bonjour papajoker
Merci d'avoir répondu

bonjour
Tu peux afficher toutes les datas dans ta boucle (je suppose qu'alors, tu trouveras ce que tu désires)

C'est ce que je fais jusqu'à maintenant.
Et cela ne fonctionne seulement avec certaines sources de fichiers, pas avec d'autres. De façon plus simple, je fais ce que tu indiques, itérer l'objet instancié par la lecture du fichier.
Entendre "fonctionne" comme fonctionnement technique, le fichier s'ouvre, son contenu est lu, parcouru, chaque itération retourne quelque chose.
Le souci est dans la nature inatendue de ce quelque chose pour certains gribs, certaines sources de données.

Source xygrib, modèle GFS (20230308_213001_GFS_P25_.grb2)

Code :

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

1
2
3
4
5
6
319:2 metre temperature:K (instant):regular_ll:heightAboveGround:level 2 m:fcst time 237 hrs:from 202303081200
320:Wind speed (gust):m s**-1 (instant):regular_ll:surface:level 0:fcst time 237 hrs:from 202303081200
321:10 metre U wind component:m s**-1 (instant):regular_ll:heightAboveGround:level 10 m:fcst time 240 hrs:from 202303081200
322:10 metre V wind component:m s**-1 (instant):regular_ll:heightAboveGround:level 10 m:fcst time 240 hrs:from 202303081200
323:2 metre temperature:K (instant):regular_ll:heightAboveGround:level 2 m:fcst time 240 hrs:from 202303081200
324:Wind speed (gust):m s**-1 (instant):regular_ll:surface:level 0:fcst time 240 hrs:from 202303081200

source xygrib, modèle Argege (20230313_195809_ARPEGE_P50_.grb2)

Code :

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

1
2
3
4
5
6
126:10 metre V wind component:m s**-1 (instant):regular_ll:heightAboveGround:level 10 m:fcst time 93 hrs:from 202303131200
127:2 metre temperature:K (instant):regular_ll:heightAboveGround:level 2 m:fcst time 93 hrs:from 202303131200
128:Wind speed (gust):m s**-1 (max):regular_ll:heightAboveGround:level 10 m:fcst time 95-96 hrs (max):from 202303131200
129:10 metre U wind component:m s**-1 (instant):regular_ll:heightAboveGround:level 10 m:fcst time 96 hrs:from 202303131200
130:10 metre V wind component:m s**-1 (instant):regular_ll:heightAboveGround:level 10 m:fcst time 96 hrs:from 202303131200
131:2 metre temperature:K (instant):regular_ll:heightAboveGround:level 2 m:fcst time 96 hrs:from 202303131200

C'est là que ça se gâte :

source weather4D, modèle Arome (AromeHD_Atlantique_lastest.grb)

Code :

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

1
2
3
4
5
82:34:34 (instant):regular_ll:heightAboveGround:level 10:fcst time 40 hrs:from 202303131200
83:33:33 (instant):regular_ll:heightAboveGround:level 10:fcst time 41 hrs:from 202303131200
84:34:34 (instant):regular_ll:heightAboveGround:level 10:fcst time 41 hrs:from 202303131200
85:33:33 (instant):regular_ll:heightAboveGround:level 10:fcst time 42 hrs:from 202303131200
86:34:34 (instant):regular_ll:heightAboveGround:level 10:fcst time 42 hrs:from 202303131200

Là, c'est pire encore, des infos sont comme attendues (ligne 972, 973), d'autres non (ligne 970, 971)

source opensiron, modèle ? (Brest_4km_WRF_WAM_230313-00.grb)

Code :

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

1
2
3
4
970:180:180 (instant):regular_ll:surface:level 0:fcst time 54 hrs:from 202303121800
971:61:61 (accum):regular_ll:surface:level 0:fcst time 53-54 hrs (accum):from 202303121800
972:Signific.height,combined wind waves+swell:m (instant):regular_ll:surface:level 0:fcst time 48 hrs:from 202303130000
973:Mean direction of wind waves:Degree true (instant):regular_ll:surface:level 0:fcst time 48 hrs:from 202303130000

Et, pourtant, les outils de visualisation arrivent très bien à s'en sortir...

[wiztricks]

Et, pourtant, les outils de visualisation arrivent très bien à s'en sortir...

Les fichiers ne sont pas corrompus mais leur contenu semble vous étonner. Est-ce qu'il est "licite" dans le format GRIB? Sinon que devrait afficher pygrib? Si oui, le format GRIB est ce qu'il est à vous de faire avec.

- W

[bmayer]

Les fichiers ne sont pas corrompus mais leur contenu semble vous étonner. Est-ce qu'il est "licite" dans le format GRIB? Sinon que devrait afficher pygrib? Si oui, le format GRIB est ce qu'il est à vous de faire avec.

Je ne gémis pas sur un contenu défectueux.
Je constate une différence à la lecture, ne doutant pas de leur intégrité (sinon pourquoi aurais-je écris que d'autres outils les lisent ?), je cherche à comprendre où je merde, afin d'y remédier.

[wiztricks]

Je ne gémis pas sur un contenu défectueux.
Je constate une différence à la lecture, ne doutant pas de leur intégrité (sinon pourquoi aurais-je écris que d'autres outils les lisent ?), je cherche à comprendre où je merde, afin d'y remédier.

J'essaie de vous dire que si ça marche, c'est peut être votre compréhension du format grib qui est à améliorer.
Ce en quoi la programmation Python ne va pas beaucoup aider.

- W