我有一个NetCDF文件input.nc
。该文件代表每4天抽样近18年的数据。从这个文件,我想计算所有时间段的月中值。因此输出文件应该只包含12个步骤。
我使用以下代码:
import xarray as xr
os.chdir(inbasedir)
data = xr.open_dataset('input.nc')
monthly_data = data.resample(freq='m', dim ='time', how = 'median')
monthly_data.to_netcdf("test.nc")
不幸的是,当我查看输出文件时,我的代码已经完成了所有时间系列的每个月的中值,最终我得到了超过200个值。如何更改我的代码,以便计算所有时间段内的12个月中位数?
您想使用groupby方法:
monthly_data = data.groupby('time.month').median()
有一些很好的例子,说明如何在这里使用xarray和时间序列数据:qazxsw poi