新手求教各位师兄师姐，对300个nc文件求均值

z15980268996

因为是小白，所以正在一步一步学习，不知道大家有脚本吗？能否发一个？拜托~
300个nc文件  0.25x0.25分辨率
求某变量各格点的均值。

andrewsoong

CDO适合你

苏拉苏拉

grads批量描述

demonisgood

这个类似问题在我处理海冰数据的时候遇到过，建议使用python下面的OS库简单编写一下就可以了（当然你的文件名必须有规律）
举个例子：
import numpy as np
import netCDF4 as nc
import os
import Ngl
ice_thickness_jan=np.empty((896,608,31))
file_all=os.listdir(r'/mnt/c/Users/59799/Desktop/smos/01')
#file_all=os.listdir(r'C:\Users\59799\Desktop\smos\01')
t=0
for i in file_all:
    domin=os.path.abspath(r'/mnt/c/Users/59799/Desktop/smos/01')
    #domin=os.path.abspath(r'C:\Users\59799\Desktop\smos\01')
    file_path=os.path.join(domin,i)
    data=nc.Dataset(file_path)
    ice_thickness_jan[:,:,t]=np.array(data.variables['sea_ice_thickness'])[0,:,:]
    lat_grid=np.array(data.variables['latitude'])
    lon_grid=np.array(data.variables['longitude'])
    t=t+1

Eremite

用python的glob批量读取
import glob
import xarray as xr
import warnings
warnings.filterwarnings("ignore")

hgt_li=glob.glob(r'f:/data2/hgt.*.*.nc') #文件名不同的部分用*代替
hgt_data=xr.open_mfdataset(hgt_li,concat_dim='time',chunks={"time": 20}) #用dask多线程读取,连接时间维度
hgt_data.z.loc[:,:,:].mean().values#这是多年平均所有格点平均，最后只剩一个值

南气斑织逍遥

装个cdo吧，方便的很

z15980268996

demonisgood 发表于 2021-1-8 22:36
这个类似问题在我处理海冰数据的时候遇到过，建议使用python下面的OS库简单编写一下就可以了（当然你的文件 ...

谢谢