[gdal-dev] Difficulties reading a ODIM HDF5 file
Sjur Kolberg
Sjur.A.Kolberg at sintef.no
Thu Mar 14 05:48:07 PDT 2013
Hello again!
Thanks for your help. I downloaded the pre-built GDAL binaries from Tamas Szekeres' www.gisinternals.com and managed to read the HDF5 file.
Almost.
And actually better with GDAL 1.7.3 than with GDAL 1.9.2, both downloaded prebuilt from www.gisinternals.com. 1.7.3 can read the projection info, 1.9.2 cannot.
They also seem to handle metadata (and subdataset?) differently; 1.9.2 seem to attribute the first 7 metadata items to a dataset1, whereas 1.7.3 does not. 1.9.2 also finds the "Conventions" attribute.
None of the two detects the what:nodata attribute.
I can see that the corner coordinates are wrong, I'll discuss that with the provider.
But does anyone see what causes 1.9.2 to miss the refsystem?
Grateful for any clue. The attached file may be different from the one in the previous mail.
Sjur :-)
GDALINFO output from 1.9.2:, command window spawn from SDKShell.bat:
Setting environment for using the GDAL and MapServer tools.
Hiding the OCI plugin library.
C:\Devel_VS2010\GDAL-1.9.2>gdalinfo norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000
039Z.hdf
Driver: HDF5Image/HDF5 Dataset
Files: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Size is 480, 480
Coordinate System is `'
Metadata:
Conventions=ODIM_H5/V2_1
dataset1_what_gain=0.5
dataset1_what_nodata=255
dataset1_what_offset=-32
dataset1_what_prodpar=0.5
dataset1_what_product=PPI
dataset1_what_quantity=DBZH
dataset1_what_undetect=0
how_endepochs=1180656039
how_software=RAINBOW
how_startepochs=1180656009
how_system=GEMA
what_date=20070601
what_object=IMAGE
what_source=NOD:norsa
what_source_plc=PLC:Rissa
what_source_wmoid=WMO:00
what_time=000039
what_version=H5rad 2.1
where_height=620
where_lat=63.69
where_lat_sphere=63.536783
where_LL_lat=61.303028
where_LL_lon=5.70656
where_lon=10.204
where_lon_sphere=10.204
where_projdef=+proj=aeqd +lat_0=63.536783 +lon_0=10.204000 +R=6371000 +datum=W
GS84
where_UR_lat=65.59814
where_UR_lon=15.411344
where_xscale=1000
where_xsize=480
where_yscale=1000
where_ysize=480
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 480.0)
Upper Right ( 480.0, 0.0)
Lower Right ( 480.0, 480.0)
Center ( 240.0, 240.0)
Band 1 Block=480x480 Type=Byte, ColorInterp=Undefined
Metadata:
dataset1_data1_data_CLASS=IMAGE
dataset1_data1_data_IMAGE_VERSION=1.2
C:\Devel_VS2010\GDAL-1.9.2>
GDALINFO output from 1.7.3:, command window spawn from SDKShell.bat:
Setting environment for using the GDAL and MapServer tools.
Hiding the OCI plugin library.
C:\Devel_VS2010\GDAL-1.7.3>gdalinfo norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000
039Z.hdf
Driver: HDF5Image/HDF5 Dataset
Files: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Size is 480, 480
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9108"]],
AUTHORITY["EPSG","4326"]]
Metadata:
what:product=PPI
what:prodpar=0.5
what:quantity=DBZH
what:gain=0.5
what:offset=-32
what:nodata=255
what:undetect=0
how:system=GEMA
how:software=RAINBOW
how:startepochs=1180656009
how:endepochs=1180656039
what:object=IMAGE
what:version=H5rad 2.1
what:date=20070601
what:time=000039
what:source=NOD:norsa
what:source_wmoid=WMO:00
what:source_plc=PLC:Rissa
where:LL_lat=61.303
where:LL_lon=5.70656
where:UR_lat=65.5981
where:UR_lon=15.4113
where:projdef=+proj=aeqd +lat_0=63.536783 +lon_0=10.204000 +R=6371000 +datum=W
GS84
where:lat=63.69
where:lon=10.204
where:lat_sphere=63.5368
where:lon_sphere=10.204
where:height=620
where:xsize=480
where:ysize=480
where:xscale=1000
where:yscale=1000
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 480.0)
Upper Right ( 480.0, 0.0)
Lower Right ( 480.0, 480.0)
Center ( 240.0, 240.0)
Band 1 Block=480x480 Type=Byte, ColorInterp=Undefined
Metadata:
data:CLASS=IMAGE
data:IMAGE_VERSION=1.2
C:\Devel_VS2010\GDAL-1.7.3>
GDALINFO output from 1.7.3 on Linux at met.no, with h5ls: Version 1.8.7. This is equal to the result I get with GDAL 1.7.3 on Windows.
christofferae at pc2752:~/src/prorad-clean/processing/prorad-process$ gdalinfo norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Driver: HDF5Image/HDF5 Dataset
Files: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Size is 480, 480
Coordinate System is:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9108"]],
AUTHORITY["EPSG","4326"]]
Metadata:
what:product=PPI
what:prodpar=0.5
what:quantity=DBZH
what:gain=0.5
what:offset=-32
what:nodata=255
what:undetect=0
how:system=GEMA
how:software=RAINBOW
how:startepochs=1180656009
how:endepochs=1180656039
what:object=IMAGE
what:version=H5rad 2.1
what:date=20070601
what:time=000039
what:source=NOD:norsa
what:source_wmoid=WMO:00
what:source_plc=PLC:Rissa
where:LL_lat=61.303
where:LL_lon=5.70656
where:UR_lat=65.5981
where:UR_lon=15.4113
where:projdef=+proj=aeqd +lat_0=63.536783 +lon_0=10.204000 +R=6371000 +datum=WGS84
where:lat=63.69
where:lon=10.204
where:lat_sphere=63.5368
where:lon_sphere=10.204
where:height=620
where:xsize=480
where:ysize=480
where:xscale=1000
where:yscale=1000
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 480.0)
Upper Right ( 480.0, 0.0)
Lower Right ( 480.0, 480.0)
Center ( 240.0, 240.0)
Band 1 Block=480x480 Type=Byte, ColorInterp=Undefined
Metadata:
data:CLASS=IMAGE
data:IMAGE_VERSION=1.2
christofferae at pc2752:~/src/prorad-clean/processing/prorad-process$ h5ls --version
h5ls: Version 1.8.7
christofferae at pc2752:~/src/prorad-clean/processing/prorad-process$ gdalinfo --version GDAL 1.7.3, released 2010/11/10
> -----Original Message-----
> From: gdal-dev-bounces at lists.osgeo.org [mailto:gdal-dev-
> bounces at lists.osgeo.org] On Behalf Of Frank Warmerdam
> Sent: 13. mars 2013 16:46
> To: Rutger
> Cc: gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] Difficulties reading a ODIM HDF5 file
>
> Sjur / Rutger,
>
> It also seems to work fine for me with trunk on linux. I suspect the problem
> is some sort of mixup at build time between versions of the HDF include files
> or something like that. Very hard to debug such a thing remotely.
>
> Best regards,
> Frank
>
> On Wed, Mar 13, 2013 at 12:51 AM, Rutger <kassies at gmail.com> wrote:
> > Hey,
> >
> > For what its worth; it works fine on Windows for me, using verion:
> > GDAL 1.9.2, released 2012/10/08
> >
> > Gdalinfo output looks equal to your Linux output, and gdal_translate
> > is able to convert the layer to a GeoTiff for example.
> >
> > Im using the pre-build cmd utilities from
> > http://www.gisinternals.com/sdk/ <http://www.gisinternals.com/sdk/> .
> > Maybe you can peak over there to see if there is any difference in building.
> >
> > Regards,
> > Rutger
> >
> >
> >
> > --
> > View this message in context:
> > http://osgeo-org.1560.n6.nabble.com/gdal-dev-Difficulties-reading-a-OD
> > IM-HDF5-file-tp5039864p5040058.html
> > Sent from the GDAL - Dev mailing list archive at Nabble.com.
> > _______________________________________________
> > gdal-dev mailing list
> > gdal-dev at lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
> --
> ---------------------------------------+--------------------------------
> ---------------------------------------+------
> I set the clouds in motion - turn up | Frank Warmerdam, warmerdam at pobox.com
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush | Geospatial Software Developer
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
-------------- next part --------------
A non-text attachment was scrubbed...
Name: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Type: application/x-hdf
Size: 15868 bytes
Desc: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20130314/19d79748/attachment-0001.hdf>
More information about the gdal-dev
mailing list