网格数据的时间轴上的滑动窗口阈值计数

时间:2019-02-09 06:58:55

标签: python numpy python-xarray

我有时间lat和lon(324,72,144)的每月网格化降水数据,并且我想计算在60个月的滑动窗口中超过某个阈值的次数。假设1-60的阈值是5,则5以下的月份数和2-61、3-63等。

我已经创建了一个for循环来执行此操作,但它似乎无法正常工作,因为它在运行后仅显示空列表,我很乐意就如何处理此问题提出建议。

a = []
for i in range(len(lon)):
    for j in range(len(lat)):
        aa = AS[:,j,i].where(AS[:,j,i]<-50).rolling(time=60).count()

我希望通过频率计数将aa变为aa.shape = 324,72,144。 任何有关创建此倾斜窗口循环的更好方法的建议将不胜感激。

0 个答案:

没有答案