extracción de la localización geográfica a partir del grib NOAA NARR

Question

He descargado un archivo grib NOAA NARR del que quiero extraer la componente de viento así como la localización geográfica. Puedo extraer la componente del viento del grib pero no encuentro ninguna indicación de que haya un valor específico de latitud/longitud para cada punto de datos. En última instancia, quiero extraer los datos y cargarlos en una base de datos relacional con tiempo, ubicación, datos_de_viento. El archivo grib que estoy utilizando se puede encontrar aquí.

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

Answer 1

Puedes traducir la primera banda (o cualquier otra) del ráster grib al formato XYZ utilizando:

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

Esto le da las coordenadas en la proyección lcc utilizada.

A continuación, puede cargar este archivo como texto delimitado en QGIS, utilizando field_1 como X y field_2 como coordenada Y, y el siguiente CRS personalizado (derivado de lo que informa gdalinfo en el 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

Esta capa de puntos puede ser fácilmente reproyectada a WGS84, y las coordenadas WGS84 añadidas a la tabla de atributos.

Como alternativa a QGIS, la reproyección también debería ser posible utilizando GDAL puro.