取一段时间内每个季节的最大降雨量（xarray）

Question

我试图找出10年内每个季节的最大降雨量。我使用netcdf数据和xarray来尝试这样做。数据包括降雨量（每3小时记录一次）、lat和lon数据。现在我有以下代码：

ds.groupby('time.season).max('time')

但是，当我这样做的时候，输出的形状是(4,145,192)，表示它在整个周期内每个季节的值都是最大的。我想每年每个季节的最大值。换句话说，输出应该具有类似(40,145,192)（每年4个值×10年）

我也试过用dataset.resample来做这个，也用time=3M作为频率，但是它没有正确地分割月份。如果必须的话，我可以改变数据集，这样它就可以从正确的位置开始，但是我希望有一种更简单的方法来考虑已经有一个函数可以正确地对它进行分组。在

谢谢，如果你需要更多的细节请告诉我！在

Answer 1

重采样将是这项工作最简单的工具。您接近时间频率，但您可能希望使用带有偏移量的季度频率：

ds.resample(time='QS-Mar').max('time')

可以按照Pandas文档中的描述进一步配置这些偏移：http://pandas.pydata.org/pandas-docs/stable/timeseries.html#offset-aliases