Как понять формат файла netCDF - PullRequest
0 голосов
/ 14 октября 2019

Я пытаюсь создать график пространственной средней (т.е. 2x2 градуса) высоты волны и скорости ветра каждого моря для различных морей в течение 5 лет с годовым или месячным интервалом времени (т. Е. Север / юг Атлантики, Тихий океанАрабское море, Северное море, Балтийское море, Средиземное море и т. Д.). Есть ли где-нибудь руководство о том, как это сделать?

Также я пытаюсь определить, является ли скорость ветра скоростью поверхностного ветра или ветра на высоте 10 м над поверхностью моря. Я знаю, что это зависит от средств наблюдения (то есть, если это корабль буй). Есть ли объяснение от NOAA или другой организации о том, как данные форматируются. Единственный документ, который я нашел, был https://icoads.noaa.gov/e-doc/imma/R3.0-imma1_short.pdf.

Заранее спасибо.

С уважением, Джордж

1 Ответ

0 голосов
/ 15 октября 2019

Из беседы Gitter об этом вопросе самое простое решение выглядит так:

ds = xr.tutorial.open_dataset('air_temperature')
ds_proc = ds.sel(lat=slice(90, 0), lon=slice(200, 250)).groupby('time.month').mean('time').mean(['lat', 'lon']))
ds_proc['air'].plot()
plt.savefig('test.png')

Если вы хотите включить полную географическую карту, мы рекомендуем использовать Cartopy . Что касается ваших метаданных, вам придется взглянуть на атрибуты, связанные с переменной ветра в вашем файле. Я не вижу ничего конкретного в предоставленном вами файле ICOADS. Если вы используете xarray, это может быть что-то вроде:

import xarray as xr
ds = xr.open_dataset('filename.nc')
print(ds['wind_variable'].attrs)
...