如何使用 python 從 netcdf 數據集中找到氣候月平均值？

Question

我有一個 netcdf 數據集（每月時間x緯度x經度）。 月數是 216，所以總共 18 年。 第一個索引代表 1 月數據，第 2 個索引代表 2 月數據……第 12 個索引代表 12 月數據，第 13 個索引代表 1 月數據，第 14 個索引代表 2 月數據，依此類推。 現在我想找到所有 1 月數據的平均值，所有 2 月數據的平均值......所有 12 月數據的平均值。 因此，在執行此操作后，數組的大小將減少到 12xlatxlon（以前是 216xlatxlon）。 我怎樣才能做到這一點？

Answer 1

如果您正在使用 linux 或 macOS，您可以使用我的 nctoolkit package ( Z5E056C500A1C4B6A7110B50D807BADE .

import nctoolkit as nc

ds = nc.open_data("infile.nc")
# get the mean using month as the averaging period
ds.tmean("month")
# plot to visualize
ds.plot()

Answer 2

您可以使用xarray / pandas中非常方便的resample()選項。

import xarray as xr 

ds = xr.open_dataset('file.nc')
ds_monthly_mean = ds.resample(time="1MS").mean()

print(ds.time.shape)
(5619,)

print(ds_monthly_mean.time.shape)
(12,)

如果您不知道要使用什么重采樣選項，請查看偏移別名。 在定義了重新采樣數據的維度之后，只需 select 您的重采樣方法（平均值、總和等....）