Извлечение Lon и Lat диапазона файла grib

Question

Я довольно новичок в формате GRIB, но сейчас мне нужно с этим справиться. С помощью wgrib мне удалось извлечь нужные поля, соответствующее время и т. Д.

Но я скучаю по координатам точек сетки или по широте и долготе диапазона сетки.

Как я могу получить это? Я бы предпочел использовать только python и wgrib, но если есть какой-либо другой простой способ, пожалуйста, дайте мне знать! То, что я нашел в Интернете, кажется, в основном касается формата wgrib2.

Answer 1

Как уже упоминалось в комментариях, вы должны использовать модуль pygrib.

Я публикую пример кода для получения массивов lat и lon из всех сообщений файла grib методом latlons. После извлечения я конвертирую данные в массивы:

import pygrib
import numpy as np

grib = pygrib.open('test.grib')
for g in grib:
    lt, ln = g.latlons() # lt - latitude, ln - longitude
    lt, ln = np.array(lt), np.array(ln)
    print (lt, ln)