如何提取不规则年份的数据

灰色节能君

今天想用提取厄尔尼诺年的数据，使用xarray读取了一个nc文件，时间是1870-2020年逐月的。
nc=xr.open_dataset('D:\data_english\HadISST_sst.nc')

具体时间格式是这样的：

'1870-01-16T11:59:59.505615234', '1870-02-14T23:59:59.340820312', '1870-03-16T11:59:59.340820312',.........

然后提取春季，

sst_spring =nc['sst'].sel(time=nc['time.season']=='MAM')


然后问题来了，

如果单纯是提取1951年到1979年每年春季的资料，那么

ElNino1_spring=sst_spring.sel(time=slice('1951','1979'))


可是我要提取1951年到1979年当中的厄尔尼诺年，那么我这样写

ElNino1_spring=sst_spring.sel(time=['1951','1957','1963','1965','1968','1972','1976','1977','1979'])



出现这样的错误

KeyError: "not all values found in index 'time'"

时间索引不对？ 难道说必须把时间索引写成 [ 1951-01-16T11:59:59.505615234，1957-01-16T11:59:59.505615234.........]。这也太麻烦了

能不能只输入年份，就像time=slice('1951','1979')一样简单得把不连续的的年份数据提取出来。

关键是单独输入一个年份时
ElNino1_spring=sst_spring.sel(time=‘1951’)
又不会报错
求高手指教要怎么输入索引才能把['1951','1957','1963','1965','1968','1972','1976','1977','1979']年的数据提取出来。(避免写成
1951-01-16T11:59:59.505615234  这种datetime64格式的)

Alexander

试试time.year isin[1979,1980,2010..]这种？

灰色节能君

Alexander 发表于 2020-11-7 20:25
试试time.year isin[1979,1980,2010..]这种？

哇 可以了！ 谢谢！ 之前把年份写成字符的形式一直没成功，原来问题在这

Nafisa6

灰色节能君 发表于 2020-11-8 12:00
哇 可以了！ 谢谢！ 之前把年份写成字符的形式一直没成功，原来问题在这

同学请问一下最后怎么写的？我写成time.year isin[1979,1980,2010..]这种报错 'DataArray' object has no attribute 'year'

灰色节能君

Nafisa6 发表于 2020-11-8 16:38
同学请问一下最后怎么写的？我写成time.year isin[1979,1980,2010..]这种报错 'DataArray' object has no ...

ElNino1_spring=sst_spring.sel(time=sst_spring['time.year'].isin([1951,1957]))