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@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@pc2752:~/src/prorad-clean/processing/prorad-process$ h5ls
--version
h5ls: Version 1.8.7
christofferae@pc2752:~/src/prorad-clean/processing/prorad-process$ gdalinfo
--version GDAL 1.7.3, released 2010/11/10
> -----Original Message-----
> From: [email protected] [mailto:gdal-dev-
> [email protected]] On Behalf Of Frank Warmerdam
> Sent: 13. mars 2013 16:46
> To: Rutger
> Cc: [email protected]
> 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 <[email protected]> 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
> > [email protected]
> > http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
> --
> ---------------------------------------+--------------------------------
> ---------------------------------------+------
> I set the clouds in motion - turn up | Frank Warmerdam, [email protected]
> light and sound - activate the windows | http://pobox.com/~warmerdam
> and watch the world go round - Rush | Geospatial Software Developer
> _______________________________________________
> gdal-dev mailing list
> [email protected]
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
Description: norsa.ss.ppi-00.5-dbz.aeqd-1000.20070601T000039Z.hdf
_______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
