Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

Annonce

Suite à un problème technique intervenu entre le 22 et le 23 mars, nous avons du procéder dans la soirée du 25 mars, à la restauration de la base de données du 24 mars (matinée).

En clair, nous avons perdu vos contributions et inscriptions du dimanche 24 et du lundi 25 mars.
Nous vous prions de nous excuser.

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

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: 166

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: 166

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

Hors ligne

 

Pied de page des forums

Powered by FluxBB