xarray - пересчет ежедневных входных данных на полугодовые данные с високосным годом

Question

У меня есть ежедневные входные данные в xarray за листовой год (например, 1972), и я хотел бы пересчитать их в полугодовые данные.

import numpy as np
import pandas as pd
import xarray as xr

time = pd.date_range('1972-01-01', freq='D', periods=366)    # a leap year
ds = xr.Dataset({'foo': ('time', np.arange(366)), 'time': time})

ds выход:

<xarray.Dataset>
Dimensions:  (time: 366)
Coordinates:
  * time     (time) datetime64[ns] 1972-01-01 1972-01-02 ... 1972-12-31
Data variables:
    foo      (time) int32 0 1 2 3 4 5 6 7 8 ... 358 359 360 361 362 363 364 365

Вот функция повторной выборки, которую я использовал для получения суммы foo каждые полгода (6 месяцев):

res = ds['foo'].resample(time='6MS', closed='left').sum('time')

Res Res:

<xarray.DataArray 'foo' (time: 3)>
array([16471., 50324.,    nan])
Coordinates:
  * time     (time) datetime64[ns] 1972-01-01 1972-07-01 1973-01-01

Я ожидал, что координаты времени будут иметь 1972-01-01 и 1972-07-01. Это ошибка для функции повторной выборки?

Answer 1

Это было идентифицировано как известная проблема и может быть исправлено обновлением панд.

pip3.exe install pandas --upgrade