从Python中的每日netcdf文件计算每月平均值

Question

您好，我有一个包含每日数据的netcdf文件。文件的形状为（5844，89，89），即16年数据。我试图从每日数据中获取每月平均值。我正在寻找类似的resample功能在熊猫数据框中运行。无论如何在python中做到这一点。据我所知，使用cdo和nco进行计算非常容易，但是我正在使用python。

我用来读取netcdf文件的示例代码是：

import netCDF4
from netCDF4 import Dataset
fh = Dataset(ncfile, mode='r')
time = fh.variables['time'][:]
lon = fh.variables['longitude'][:]
lat = fh.variables['latitude'][:]
data = fh.variables['t2m'][:]
data.shape

Answer 1

@詹姆曼谢谢你的建议xarray.resample。它比我想象的简单，我的问题的答案是：

import xarray as xr
ds = xr.open_dataset(nc_file)
monthly_data = ds.resample(freq = 'm', dim = 'time', how = 'mean')