[gdal-dev] WCS driver, GetSubdatasets (Python) returns empty list

[Jachym Cepicky]

Hi,

using Python, gdal 1.8.0, Ubuntu

I need to access some Temporal information about available WCS layer

http://www.gdal.org/frmt_wcs.html says at the end, that

"Each time position available for the coverage will be treated as a subdataset."

I have following service description file (with service, I need to deal with):

<WCS_GDAL>
  <ServiceURL>http://briseide01.ingr.briseide.eu:8080/geoserver/ows?</ServiceURL>
  <CoverageName>20101012T210222_pig</CoverageName>
  <Timeout>5000</Timeout>
</WCS_GDAL>

The DescribeCoveridge request seems to be passed correctly, the file
is "extended" with aditional informations, among others it has

      <temporalDomain>
        <timePosition>2010-10-17T00:00:00.000Z</timePosition>
        <timePosition>2010-10-17T01:00:00.000Z</timePosition>
        <timePosition>2010-10-17T02:00:00.000Z</timePosition>
        <timePosition>2010-10-17T03:00:00.000Z</timePosition>
        <timePosition>2010-10-17T04:00:00.000Z</timePosition>
      </temporalDomain>

in it

So, I would assume, there is Temporal dimension available in the WCS dataset.

In my Python code, I approach:

from osgeo import gdal
ds = gdal.Open("the_wcs_def_file.xml")
ds.GetMetadata()
# returns empty list
ds.GetSubDatasets()
# returns empty list

So, my question would be: how to access the temporal information from
WCS dataset, using Python, according to WCS format description
(http://www.gdal.org/frmt_wcs.html)?

Thanks

Jachym

P.S. Congrats, Frank! Good luck in your new job.
-- 
Jachym Cepicky
e-mail: jachym.cepicky gmail com
URL: http://les-ejk.cz
GPG: http://les-ejk.cz/pgp/JachymCepicky.pgp