[gdal-dev] HowTo open NetCDF SUBDATASET and write to new file

[Steve Gaffigan]

Hello.  Here's an example of reading an entire array for a single subdataset.

[gaffigan at dp2 AEFF]$ gdalinfo latest.nc
Driver: netCDF/Network Common Data Format
...
  SUBDATASET_25_NAME=NETCDF:"latest.nc":TMP_GPML
  SUBDATASET_25_DESC=[49x11x78x141]
air_temperature_at_constant_altitude_above_m
ean_sea_level (8-bit integer)
...

[gaffigan at dp2 AEFF]$ python
>>> from osgeo import gdal
>>> ds = gdal.Open('NETCDF:"latest.nc":TMP_GPML')
>>> data = ds.ReadAsArray().reshape(49,11,78,141)

Using a more generic call of the form below you can get at particular
hyperslabs:

>>> data=ds.GetRasterBand(i).ReadAsArray(xoff=x0,yoff=y0,win_xsize=nx,win_ysize=ny)

where, for the above example, i=[1,...,539].

Steve


> Hi List,
>
> I'm would like to use the Python bindings to open a specific SUBDATASET in
> a
> NetCDF file and read the contents.  I read in
> http://www.gdal.org/gdal_datamodel.html, under the "SUBDATASETS domain"
> section, that I should be able to pass the "_NAME" parameter into
> GDALOpen()
> in order to do this, but I'm not exactly sure how.  Could someone pass me
> a
> couple lines of Python that show how this is done?
>
> Thanks in advance.
>
> Roger
> --
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev