[英]Extract data from NetCDF4 file with regional mask file
我正在處理多個 NetCDF4 文件,並希望從中提取monthly_rain
值。 這是我的代碼:
import numpy
import netCDF4
with netCDF4.Dataset('abc.nc', 'r') as mask_dataset:
mask_data = mask_dataset.variables['mask'][:]
results = []
for year in range(2010, 2019):
with netCDF4.Dataset('{:d}.monthly_rain.nc'.format(year), 'r') as dataset:
data = dataset.variables['monthly_rain'][:]
data.mask = mask_data.mask
average = numpy.mean(data)
results.append(average)
print(results)
上述代碼的結果是:
[92.82600598372804, 67.01124693612452, 54.30168356893234, 39.58771623758809, 45.30353874205824, 39.017626997972684, 50.94861235658874, 44.55133832249074, 41.7971056907917]
這就是我想要的結果。
但是,我想從文件中提取所有monthly_rain
值,以便對數據集進行進一步檢查。 有什么方法可以讓我這樣做嗎?
現在我可以回答了……只是不要在循環中計算平均值,而是直接對結果進行 append(屏蔽)“數據”,然后進行任何額外的后處理。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.