我有一个数据集,其中包含几年的每日 x,y 网格气象数据。我只对计算冬季数据的年度平均值感兴趣,即。也不包括夏季数据。
我认为我需要使用该resample
命令,例如,AS-OCT
将时间序列重新采样为年度频率,冬季从每年 10 月开始(它是北纬地区)。
我无法解决的是如何指定我只想使用从 10 月到 4 月/5 月的数据,而忽略 6 月、7 月和 8 月。
由于 resample 函数适用于ndarray
对象,我想出了一个相当不便携的方法来做这个总和:
def winter(x,axis):
# Only use data from 1 October to end of April (day 211)
return np.sum(x[0:211,:,:],axis=0)
win_sum = all_data.resample('AS-OCT',how=winter,dim='TIME')
但我觉得应该有一个更优雅的解决方案。有任何想法吗?