gpt4 book ai didi

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

转载 作者:行者123 更新时间:2023-11-30 22:04:12 24 4
gpt4 key购买 nike

我在 xarray 数据集中加载了 Netcdf 文件,我想制作没有闰日的每日气候学,即不包含 2 月 29 日。我正在尝试 syntaxDataset.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 天来计算每日气候,那就太好了,但我也想知道如何在需要时从特定时间步长中删除数据。

最佳答案

此处使用 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)))

关于python - 从 xarray 数据集中的闰年中删除闰日,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/53376113/

24 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com