Чтение одномерных переменных в файлах WRF NetCDF с помощью GDAL python

У меня простой вопрос.
Например, с файлом wrfout "out.nc".
Файл содержит переменные Geo2D, Geo3D и 1D.

Используя пакет GDAL в Python 2.7, я могу легко извлечь переменные Geo2D следующим образом:

## T2 is 2-d variable means temperature 2 m above the ground
temp = gdal.Open('NETCDF:"'+"out.nc"+'":T2')

Но когда я хочу использовать этот код для извлечения массива 1d, это не удалось.

## Time is 1-d array represent the timeseries throught the simulation period
time = gdal.Open('NETCDF:"'+"out.nc"+'":Time')

Ничего не произошло! Хотелось бы, чтобы кто-нибудь дал совет, как легко читать выходные переменные WRF любого измерения!


python netcdf matplotlib-basemap gdal
person Han Zhengzu    schedule 24.03.2016    source источник
comment
Вы пытались использовать python-netcdf4 вместо gdal?   -  person kakk11    schedule 27.03.2016

Ответы (1)


arrow_upward
3
arrow_downward

Вы также можете использовать программу чтения NetCDF в scipy.io:

import scipy.io.netcdf as nc

# Open a netcdf file object and assign the data values to a variable
time = nc.netcdf_file('out.nc', 'r').variables['Time'][:]

Его преимущество заключается в том, что scipy является очень популярным и широко устанавливаемым пакетом, хотя в некоторых отношениях он работает аналогично открытию файлов.

person Luis Ortiz    schedule 03.05.2016
comment
Спасибо за ваш ответ. Я постараюсь! - person Han Zhengzu; 04.05.2016