Python netcdf - среднемесячная медиана за все время ежедневных данных

Question

У меня есть файл NetCDF input.nc.Этот файл представляет почти 18 лет данных, взятых каждые 4 дня.Из этого файла я хотел бы рассчитать месячное значение медианы за весь период времени.Таким образом, выходной файл должен содержать только 12 временных шагов.

Я использую следующий код:

import xarray as xr
os.chdir(inbasedir)
data = xr.open_dataset('input.nc')
monthly_data = data.resample(freq='m', dim ='time', how = 'median')
monthly_data.to_netcdf("test.nc")

К сожалению, когда я смотрю на выходной файл, мой код принял значение медианыза каждый месяц серии все время и я получаю более 200 значений.Как я могу изменить свой код, чтобы рассчитать 12 месячных медиан за все время?

Answer 1

Вы хотите использовать метод groupby :

monthly_data = data.groupby('time.month').median()

Вот несколько хороших примеров использования xarray с данными временных рядов: http://xarray.pydata.org/en/stable/time-series.html