Suite

Extraction de l'emplacement géographique à partir du grib NOAA NARR

Extraction de l'emplacement géographique à partir du grib NOAA NARR


J'ai téléchargé un fichier grib NOAA NARR dont je souhaite extraire la composante éolienne ainsi que l'emplacement géographique. Je suis capable d'extraire la composante du vent du grib, mais je n'arrive pas à trouver d'indication qu'il existe une valeur de latitude/longitude spécifique pour chaque point de données. En fin de compte, je souhaite extraire les données et les charger dans une base de données relationnelle avec heure, emplacement, wind_data. Le fichier grib que j'utilise se trouve ici.

http://nomads.ncdc.noaa.gov/data/narr/201410/20141001/narr-a_221_20141001_0000_000.grb


Vous pouvez traduire la première (ou toute autre) bande du raster grib au format XYZ en utilisant :

gdal_translate -b 1 -of "XYZ" narr-a.gsb narr-a.xyz

Cela vous donne les coordonnées dans la projection lcc utilisée.

Ensuite, vous pouvez charger ce fichier sous forme de texte délimité dans QGIS, en utilisantchamp_1comme X etchamp_2comme coordonnée Y, et le CRS personnalisé suivant (dérivé de ce que gdalinfo rapporte sur le grb) :

+proj=lcc +lat_1=50 +lat_2=50 +lat_0=50 +lon_0=-107 +x_0=0 +y_0=0 +a=6371200 +b=6371200 +units=m +no_defs

Cette couche de points peut facilement être reprojetée sur WGS84 et les coordonnées WGS84 ajoutées à la table attributaire.

Alternativement à QGIS, la reprojection devrait également être possible en utilisant GDAL pur.


Voir la vidéo: pre-code hollywood