2016-03-24 10 views
0

Meine Frage ist einfach.
Mit einer Wrfout-Datei "out.nc" zum Beispiel.
Die Datei enthält Geo2D-, Geo3D- und 1D-Variablen.Lesen von 1-D-Variablen in WRF NetCDF-Dateien mit GDAL Python

Mit GDAL Paket in Python 2.7, kann ich die Geo2D Variablen leicht wie folgt extrahieren:

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

Aber wenn ich diesen Code verwenden möchten 1d Array zu extrahieren, es ist fehlgeschlagen.

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

Nichts ist passiert! Wünschen Sie ein Angebot, um eine beliebige Dimension der WRF-Ausgangsvariablen leicht zu lesen!

+1

Haben Sie versucht, Python-Netcdf4 anstelle von Gdal? – kakk11

Antwort

3

Sie können auch die NetCDF Leser in scipy.io verwenden:

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'][:] 

Es den Vorteil scipy hat ein sehr beliebtes und weit installierte Paket ist, während ähnliche Arbeitsdateien in mancher Hinsicht zu öffnen.

+0

Danke für Ihre Antwort. Ich werde es versuchen! –