我有时间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。 任何有关创建此倾斜窗口循环的更好方法的建议将不胜感激。