Question

我试图在10年内找到每个季节(DJF, MAM, JJA, SON)的最大降雨量值。我正在使用netcdf数据和xarray来尝试这样做。数据包括降雨量（每3小时记录一次），纬度和经度数据。现在我有以下代码：

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

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

我已尝试使用DataSet.resample尝试使用time=3M作为频率，但是它并没有正确地拆分月份。如果必须的话，我可以改变数据集，所以它从正确的位置开始，但我希望有一种更简单的方法，考虑到它已经是一个正确分组的功能。

谢谢，如果您需要更多细节，请告诉我们！

Answer 1

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

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