使用GDAL python读取WRF NetCDF文件中的一维变量

Question

我的问题很简单。
例如，带有wrfout文件“ out.nc”。
该文件包含Geo2D，Geo3D和1D变量。

使用Python 2.7中的GDAL包，我可以像这样轻松提取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输出变量的任意维！

Answer 1

您也可以在scipy.io中使用NetCDF阅读器：

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是一个非常受欢迎且安装广泛的软件包，同时在某些方面类似于打开文件。