从 xarray 数据集中的闰年中删除闰日

Question

我在 xarray 数据集中加载了 Netcdf 文件，我想制作没有闰日的每日气候学，即不包括 2 月 29 日。我正在尝试 syntax 的 Dataset.drop 方法对我来说不是那么直观。这是数据集

print(ds)
>><xarray.Dataset>
Dimensions:        (lat: 1, lev: 1, lon: 720, time: 27133)
Coordinates:
* lon            (lon) float32 -180.0 -179.5 -179.0 ... 178.5 179.0 179.5
* lev            (lev) float32 1.0
* time           (time) datetime64[ns] 2000-01-02T18:00:00 ... 2018-07-30
Dimensions without coordinates: lat
Data variables:
Var1              (time, lev, lon) float32 ...
Var2              (time, lat, lon) float64 ...
Var3              (time, lat, lon) float64 ...

我试过了

ds_N_R.drop(['Var1', 'Var2', 'Var3'], time='2000-02-29')
>>TypeError: drop() got an unexpected keyword argument 'time'
##another approach
ds_N_R.sel(time='2000-02-29').drop(['Var1', 'Var2', 'Var3'])
## gives not the result I intended
<xarray.Dataset>
Dimensions:  (lev: 1, lon: 720, time: 4)
Coordinates:
* lon      (lon) float32 -180.0 -179.5 -179.0 -178.5 ... 178.5 179.0 179.5
* lev      (lev) float32 1.0
* time     (time) datetime64[ns] 2000-02-29 ... 2000-02-29T18:00:00
Data variables:
*empty*

我该如何继续？很高兴知道是否有一种直接的方法可以计算一年中仅 365 天的每日气候，但我也想知道如何在需要时从特定时间步长中删除数据。

Answer 1

此处使用 drop() 的正确方法是： ds_N_R.drop([np.datetime64('2000-02-29')], dim='time')

但我认为这实际上可以通过索引操作更干净地完成，例如， ds_N_R.sel(time=~((ds_N_R.time.dt.month == 2) & (ds_N_R.time.dt.day == 29)))

从 xarray 数据集中的闰年中删除闰日

Removing leap day from leap years in xarray dataset

python

python-xarray