xarray怎样从逐月数据文件挑选部分年的冬季数据做逐年季节平均？

Nafisa6

使用xarray读取了一个nc文件，数据是逐月的，尝试了选取特定年份的冬季数据，目前只能按下面的方法选，有点太繁琐，求助有没有简单点的方法？
times = xr.DataArray(
pd.to_datetime(['1966-12-01','1967-01-01','1967-02-01',……,'1999-12-01','2000-01-01','2000-02-01']), dims="new_time")
air_DJF=(air_DJF.sel(time=times))

然后要求逐年冬季平均怎么求呢？用groupby做只能求总的平均
[size=13.0667px]air.groupby(air.time.dt.season).mean())

Nafisa6

Nafisa6 发表于 2020-11-8 22:19
找了很久按照先重采样再做季节平均来算，可以出来结果，但是跟直接取某一年的冬季结果不一样，有没有大神帮 ...

找到原因了，改成air_season=(air.resample(time = '3MS').mean())就可以了

Nafisa6

尽头的尽头 发表于 2020-11-8 20:50
季节平均以后对时间做平均就是四个季节了.mean('time')

请问可以说的具体点怎么做么，我想要的结果是1966冬、1967冬……2000冬这样的

尽头的尽头

季节平均以后对时间做平均就是四个季节了.mean('time')

Nafisa6

找了很久按照先重采样再做季节平均来算，可以出来结果，但是跟直接取某一年的冬季结果不一样，有没有大神帮忙看看是什么问题？

1. air=air.loc['1960-12-01':'2019-02-28']
   
2. air_season=(air.resample(time = '3M').mean())
   
3. air_DJF=air_season.loc[air_season.time.dt.month==12,:,:]
   
4. air_DJF_mean=air_DJF.mean(dim=['time'])
   
5. 
   
6. air_DJF=air_DJF.sel(time=air_DJF['time.year'].isin([1960]))
   

复制代码

冬青老先生

Nafisa6 发表于 2020-11-8 22:04
请问可以说的具体点怎么做么，我想要的结果是1966冬、1967冬……2000冬这样的

噗，不好意思，刚刚点错点到反对了

ahtoffff

Nafisa6 发表于 2020-11-8 22:19
找了很久按照先重采样再做季节平均来算，可以出来结果，但是跟直接取某一年的冬季结果不一样，有没有大神帮 ...

请问一下，我需要的是79年冬到18年冬天，将您loc后的那两个时间改成我这个时间，但却求出来的都是nan
filename = r'E:\pyclass\mslp.mon.mean.nc'
dataset = xr.open_dataset(filename)
mslp = dataset['mslp']
lon = dataset['lon']
lat = dataset['lat']
slp = mslp.loc['1979-01-01':'2019-02-28']
slp_season = (slp.resample(time='3MS').mean())
slp_DJF = slp_season.loc[slp_season.time.dt.month == 12, :, :]
slp_DJF_mean = slp_DJF.mean(dim=['time'])
这个slp_DJP都是nan