Lecture de fichiers GRIB avec pyGrib

Version imprimable

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:

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

bonjour

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

Code:

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

Bonjour papajoker
Merci d'avoir répondu

Citation:

Envoyé par papajoker

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:

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:

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:

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:

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...

20/03/2023, 19h18
wiztricks

Citation:

Envoyé par bmayer

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
21/03/2023, 11h01
bmayer

Citation:

Envoyé par wiztricks

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.
21/03/2023, 11h20
wiztricks

Citation:

Envoyé par bmayer

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