#1 Wed 22 September 2010 08:48
- Mropop
- Juste Inscrit !
- Date d'inscription: 22 Sep 2010
- Messages: 5
Gdal Grib
Bonjour,
Avant tout excusez moi pour ma minuscule connaissance en l'outil gdal et donc surement aux bêtises que je vais pouvoir dire sur mon premier sujet... ^^
En fait j'ai une question surement assez simple pour un connaisseur, mais alors pour moi... j'y passe déjà 3h sur google pour essayer d'avoir une petite explication.
Bien sur j'y ai vu quelques pistes mais rien de fameux...
Je vais vous détailler mon soucis, j'ai installer gdal comme il faut sur mon Linux, et une fois près à travailler. PAF comprends pas, je dois simplement lire un fichier GRIB récupérer sur un site de prévisions météo. Mais maintenant j'arrive a la super question, y'a t il un moyen de récupérer par exemple une infos concernant le vent ou la houle par rapport à des coordonnées latitude longitude sans faire 500 initialisation de trucs dont je ne me servirais pas?
Voila, si jamais vous avez besoin de détails, demandez moi.
Excusez moi encore pour mon débutantisme. 
Merci beaucoup.
Hors ligne
#2 Wed 22 September 2010 10:14
- Mropop
- Juste Inscrit !
- Date d'inscription: 22 Sep 2010
- Messages: 5
Re: Gdal Grib
Euh, au cas ou ma question soit incompréhensible dites le moi ^^
Merci.
Dernière modification par Mropop (Wed 22 September 2010 11:14)
Hors ligne
#3 Wed 22 September 2010 11:54
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: Gdal Grib
C'est tout à fait faisable, en supposant que vent et houle sont des données des bandes du fichier GRIB en question. Mais même s'il n'y a pas besoin de "500 initialisations", il va falloir te plonger un peu dans l'API GDAL et mettre en oeuvre quelques capacités de codage/scriptage.
Je vais supposer que tu connais un peu le Python...
Des généralités sur les concepts utilisés par GDAL : http://gdal.org/gdal_datamodel.html
Et un tutorial les mettant en pratique (ex C/C++/Python) : http://gdal.org/gdal_tutorial.html
Le point d'entrée général pour trouver des infos sur l'API Python pour GDAL : http://trac.osgeo.org/gdal/wiki/GdalOgrInPython
Concernant ton besoin particulier, voici quelques raccourcis :
* extraction d'une valeur à partir des coordonnées exprimées dans le SRS d'origine : http://www.gis.usu.edu/~chrisg/python/2 … lides4.pdf
* Si ton fichier GRIB est projeté et non géographique, il faudra de plus reprojeter préalablement tes coordonnées (long, lat) dans le SRS du fichier GRIB. Des infos intéressantes : http://www.gis.usu.edu/~chrisg/python/2 … lides2.pdf à partir du slide 26
Sinon le fichier d'exemple http://svn.osgeo.org/gdal/trunk/gdal/sw … latlong.py contient une partie de dont tu as besoin, sauf qu'il faut "inverser" ce qu'il fait. En particulier l'ordre des SRS à la création de l'objet CoordinateTransformation et manipuler l'inverse de la matrice de géotransformation (obtenue via gdal.InvGeoTransform(mat_src, mat_inv_dst)) afin de convertir de coordonnées projetées vers coordonnées image (x,y).
A l'occasion je verrais pour contribuer un petit script d'exemple car ce besoin(extraction d'une valeur à partir de coordonnées (long, lat)) doit être effectivement assez courant.
Dernière modification par rouault (Wed 22 September 2010 11:56)
Hors ligne
#4 Wed 22 September 2010 19:52
- rouault
- Participant assidu
- Date d'inscription: 26 Apr 2009
- Messages: 168
Re: Gdal Grib
Voilà un petit script qui extrait les valeurs des bandes d'un fichier à la coordonnée spécifiée : http://svn.osgeo.org/gdal/trunk/gdal/sw … t_coord.py
Hors ligne
#5 Mon 07 March 2011 17:07
- Mropop
- Juste Inscrit !
- Date d'inscription: 22 Sep 2010
- Messages: 5
Re: Gdal Grib
Bonjour,
6 mois plus tard je découvre ta réponse et je te remercie 
J'ai reussi à récuperer des informations météo sur mon grib, mais j'aimerai savoir comment récuperer l'infos....
Y a 500 couches sur mon grib mais j'arrive pas a trouver celle que je veux...
Par exemple sur ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/ je lis gfs.2011030700/gfs.t00z.pgrb2f00 et je cherche la force et la direction du vent... mais j'ai aucune idée de ou trouver ses valeurs parmis les couches... je cherche juste la vitesse du vent en noeud ou km/h et sa direction en degrés....
Merci beaucoup de l'aide
Hors ligne
