python 三个月滑动平均判断厄尔尼诺年

感觉啥都不会

我想用三个月滑动平均大于0.5，且维持5个月以上来判断厄尔尼诺年，但是我用这个条件算出来的厄尔尼诺年与有记录的厄尔尼诺年对不上。比如1976年，它是一个强厄尔尼诺年。但是我算出来的三个月滑动平均只有最后两个值是大于0.5的，其余值都小于0.5，都不能判定为一个厄尔尼诺年。想请问大家这是为什么啊，我该怎么改呢



#调取EINO3.4区的数据
sst34=f1['sst'].loc['1960-01-01':'2022-12-01'].loc[:,5:-5,190:240]
lon=f1['lon'].loc[190:240]
lat=f1['lat'].loc[5:-5]
#time=f1.variables['time']
#print(f1)
#print(sst34)#(744,5,26)
#将数据转换成四维数组（year，month，lat，lon）对年求平均场及距平场
sst34=np.array(sst34).reshape(63,12,5,26)
sst34_1=np.mean(np.mean(sst34,2),2)#平均场 (63,12)
sst34_2=np.mean(sst34_1,0) #(12)

sst34_jp=np.zeros([63,12])
for m in range(63):
    for n in range(12):
        sst34_jp[m,n]=sst34_1[m,n]-sst34_2[n]
#找出三个月滑动平均大于0.5，且维持5个月及以上的年份
sst34_move=np.zeros([63,12])
for m in range(63):
    sst34_move[m,:]=gaussian_filter1d(sst34_jp[m,:], 3)

cookie-o-o

这可能与数据有关，可以换一套海温资料试试。其实，CPC定义的1976/1977年冬季的El Nino也不强。见https://origin.cpc.ncep.noaa.gov ... nsostuff/ONI_v5.php

感觉啥都不会

cookie-o-o 发表于 2023-12-17 08:38
这可能与数据有关，可以换一套海温资料试试。其实，CPC定义的1976/1977年冬季的El Nino也不强。见https://o ...

好的，非常感谢

puck66

没看到你挑选“三个月滑动平均大于0.5，且维持5个月及以上的年份”对应的代码