[gdal-dev] Can't read subdatasets list from SeaWiFS HDF4

[Ivan]

Hi there,

I tried to read HDF4 files downloaded from this website:

http://oceancolor.gsfc.nasa.gov/

More specifically, the SeaWiFS L3 BIN:

http://oceandata.sci.gsfc.nasa.gov/SeaWiFS/L3BIN/2009/312/S2009312.L3b_DAY_CHL.main.bz2

But running gdalinfo on that file doesn't show any SUBDATASET:

% gdalinfo S2009312.L3b_DAY_CHL.main
Driver: HDF4/Hierarchical Data Format Release 4
Files: S2009312.L3b_DAY_CHL.main
Size is 512, 512
Coordinate System is `'
Metadata:
  Product Name=S2009312.L3b_DAY_CHL.main
  Title=SeaWiFS Level-3 Binned Data
  Sensor Name=SeaWiFS
  Data Center=NASA/GSFC SeaWiFS Data Processing Center
  Mission=SeaStar SeaWiFS
  Mission Characteristics=Nominal orbit: inclination = 98.2 (Sun-synchronous); node = 12 noon local (descending); eccentricity = <0.002; altitude = 705 km; ground speed = 6.75 km/sec
  Sensor=Sea-viewing Wide Field-of-view Sensor (SeaWiFS)
  Sensor Characteristics=Number of bands = 8; number of active bands = 8; wavelengths per band (nm) = 412, 443, 490, 510, 555, 670, 765, 865; bits per pixel = 10; instantaneous field-of-view = 1.5835 mrad; pixels per scan = 1285; scan rate = 6/sec; sample rate = 7710/sec
  Product Type=D
  Processing Version=2009.1
  Software Name=l2bin
  Software Version=2.3.3
  Orbit=65360
  Start Orbit=65340
  End Orbit=65360
  Processing Time=2009322051028000
  Processing Control=l2bin par=S2009312.L3b_DAY_CHL.param
  Input Parameters=INFILE = /data1/sdpsoper/vdc/vpu0/workbuf/l2inlist.dat|OFILE = S2009312.L3b_DAY_CHL.main|FILEUSE = S2009312.L3b_DAY_CHL.contrib|SDAY = 2009312|EDAY = 2009312|RESOLVE = 9|ROWGROUP = 2700|FLAGUSE = ATMFAIL,LAND,HILT,HISATZEN,STRAYLIGHT,CLDICE,COCCOLITH,LOWLW,CHLWARN,CHLFAIL,NAVWARN,MAXAERITER,ATMWARN,HISOLZEN,NAVFAIL,FILTER,HIGLINT|L3BPROD = chlor_a|PRODTYPE = D|PVERSION = 2009.1|SUITE = CHL|AVERAGE = standard|NIGHT = 0|VERBOSE = 0|MINOBS = 0|QUALPROD = NONE|QUAL MAX = N/A|
  Input Files=/data1/sdpsoper/vdc/vpu0/workbuf/S2009311192733.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009311193509.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009311210841.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009311224722.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312002604.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312020446.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312034327.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312052209.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312070053.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312083932.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312101814.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312115655.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312133537.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312151418.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312165300.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312183023.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312183706.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312201024.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312214906.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009312232747.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009313010819.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009313024700.L2_GAC_OC,/data1/sdpsoper/vdc/vpu0/workbuf/S2009313042542.L2_GAC_OC
  L2 Flag Names=ATMFAIL,LAND,HILT,HISATZEN,STRAYLIGHT,CLDICE,COCCOLITH,LOWLW,CHLWARN,CHLFAIL,NAVWARN,MAXAERITER,ATMWARN,HISOLZEN,NAVFAIL,FILTER,HIGLINT
  Period Start Year=2009
  Period Start Day=312
  Period End Year=2009
  Period End Day=312
  Start Time=2009311192732633
  End Time=2009313050540954
  Start Year=2009
  Start Day=311
  Start Millisec=70052633
  End Year=2009
  End Day=313
  End Millisec=18340954
  Latitude Units=degrees North
  Longitude Units=degrees East
  Northernmost Latitude=45.625
  Southernmost Latitude=-74.95834
  Easternmost Longitude=179.9583
  Westernmost Longitude=-179.9583
  Data Bins=384375
  Percent Data Bins=6.4705
  Station Name=GSFC HRPT, NASA, MD
  Station Latitude=38.9958
  Station Longitude=-76.8511
  Units=chlor_a:mg m^-3
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

Does anybody have any experience with this data provider and how to get raster images from that file?

My best regards,

Ivan