如何在python中计算网格数据集的平均值

Question

您好我有.nc格式的网格化数据集。我尝试使用下面给出的代码读取.nc文件中不同的纬度和经度组合，但我知道这不是理想的方法。

import netCDF4
from netCDF4 import Dataset
f= Dataset('data.nc')

data1 =  f.variables['asa'][:,3,23]
data2 =  f.variables['asa'][:,3,24]
data3 =  f.variables['asa'][:,3,25]
data4 =  f.variables['asa'][:,4,23]
data5 =  f.variables['asa'][:,4,24]
data6 =  f.variables['asa'][:,4,25]

最后

data = np.vstack((data1,data2,data3,data4...))

Data = np.average(data, axis = 0)    

有没有简单的方法可以做到这一点而不是读取每一个网格并计算平均值？ 您的帮助将受到高度赞赏。

Answer 1

我不知道nc格式化，但我希望你能够做到

import netCDF4
from netCDF4 import Dataset

f= Dataset('data.nc')
Data = np.average(f.variables['asa'][:,3:5,23:26], axis = 0)