gpt4 book ai didi

python - 用 python 平均多个 netCDF4 文件

转载 作者:行者123 更新时间:2023-12-05 04:48:52 27 4
gpt4 key购买 nike

我是 python 菜鸟中的 netCDF,所以请原谅这个菜鸟问题。

我有一个文件夹,里面装满了大约 3650 个 netCDF4 文件。十年来每天一个文件。奈尔斯被命名为 yyyymmdd.nc(例如 20100101、20100102、20100103 等)。每个 .nc 文件都包含同一区域(汤加专属经济区的一部分)在同一时间点的纬度、经度和温度。

我想做的是计算所有文件中每个纬度和经度的平均温度,即我想最终得到一个 .nc 文件,该文件具有相同的纬度和经度以及 10 年的平均温度。

我尝试过不同的东西/代码版本,通常,它们最终看起来像这样......

files = glob('*.nc')
ds = xr.open_mfdataset(files,)
mean = np.mean(ds['temp'][:, 0].values)

......这段代码会给我一个 .nc 文件中所有 .nc 文件的平均温度,而不是基于纬度和经度的十年文件的平均温度。

非常感谢所有帮助。

谢谢。

最佳答案

假设您在 linux/macOS 上工作,这可以使用我的 nctoolkit 包轻松完成(参见详细信息 here)。

以下将计算所有文件的平均值,然后绘制结果:

import nctoolkit as nc
files = glob('*.nc')
ds = nc.open_data(file)
ds.ensemble_mean()
ds.plot()

nctoolkit 默认使用 CDO 作为后端,但也可以使用 NCO,从而提高性能。所以以下可能会更快:

import nctoolkit as nc
files = glob('*.nc')
ds = nc.open_data(file)
ds.ensemble_mean(nco=True)
ds.plot()

关于python - 用 python 平均多个 netCDF4 文件,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/67868777/

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