Извлечение данных из файла NetCDF4 с файлом региональной маски - PullRequest
Новая
       31

Извлечение данных из файла NetCDF4 с файлом региональной маски

0 голосов
/ 26 апреля 2020

Я работаю над несколькими файлами 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 из файла, чтобы я мог провести дальнейшую проверку набора данных. Есть ли способ, который может позволить мне сделать это?

1 Ответ

0 голосов
/ 26 апреля 2020

Теперь я могу ответить ... Только не вычисляйте средние значения в l oop, а непосредственно добавляйте (маскируйте) "данные" к результатам, а затем выполняйте любую дополнительную постобработку.

...