如何将全年每天3小时的数据导入xarray

Question

我试图将每小时3次的netCDF数据导入到xarray中，并找到它的每日平均值。我可以简单地选择某一天并得到每日平均值，但如何在一整年内完成这个过程呢？如何做到这一点

import xarray as xr

fname='/home/atmosphere/GRIDSAT-B1.2010.12.01.*.v02r01.nc'         
ds=xr.open_mfdataset(fname))
varlist=list(ds.variables)
imr=ds.sel(lat=slice(0,30),lon=slice(25,45)) #subsetting dataset
bt=imr['irwin_cdr']                          #making into a data array  
dailymean1=imrbt.mean('time')

在上面的代码中，我可以得到2010年12月1日的每日平均值，但我如何才能找到全年的每日平均值？谢谢你的帮助

Answer 1

我想您可能想使用一个名为^{}的方法：

daily_mean = ds.resample(time='D').mean()

它的一个主要用例是将时间序列下采样到一个规则的粗略频率（例如，从三小时数据到每日平均数据）。确保ds包含您希望重新采样的所有数据（您可能需要更改glob字符串）