[Qgis-developer] Problems with gdal and hdf files

Agustin Lobo alobolistas at gmail.com
Wed May 2 02:41:34 PDT 2012


Hi!

The gdal info of a multiband satellite image in hdf is reported by
gdal as "Size is 512, 512"
in contradiction with the size of the subdatasets. QGIS and R
(raster+rgdal) report an error at opening the file,
while envi and HDFLook can open it with no problem. I include the
gdalinfo output and the
info reported by HDFLook. Any help appreciated.
Thanks
Agus

$ gdalinfo --version
GDAL 1.8.0, released 2011/01/12

$ gdalinfo lndsr.L5169054_05420101216.hdf
Driver: HDF4/Hierarchical Data Format Release 4
Files: lndsr.L5169054_05420101216.hdf
Size is 512, 512
Coordinate System is `'
Metadata:
  DataProvider=UMD
  Satellite=LANDSAT_5
  Instrument=TM
  AcquisitionDate=2010-12-16T00:00:00.000000Z
  Level1ProductionDate=2012-02-22T00:00:00.000000Z
  SolarZenith=42.02009
  SolarAzimuth=141.4236
  WRS_System=2
  WRS_Path=169
  WRS_Row=54
  NumberOfBands=6
  BandNumbers=1, 2, 3, 4, 5, 7
  ShortName=L5TSR
  LocalGranuleID=L5TSR.a2010350.w2p169r054.020.2012117105207.hdf
  ProductionDate=2012-04-26T10:52:07Z
  PGEVersion=1.0.9
  ProcessVersion=1.0.4
  WestBoundingCoordinate=36.6199876177544
  EastBoundingCoordinate=38.780676500528
  NorthBoundingCoordinate=9.61826703329522
  SouthBoundingCoordinate=7.72285510894074
  OrientationAngle=0
  PixelSize=30
  HDFEOSVersion=HDFEOS_V2.4
  Cloud Mask Algo Version=CMReflectanceBasedv1.0
Subdatasets:
  SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band1
  SUBDATASET_1_DESC=[6961x7901] band1 Grid (16-bit integer)
  SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band2
  SUBDATASET_2_DESC=[6961x7901] band2 Grid (16-bit integer)
  SUBDATASET_3_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band3
  SUBDATASET_3_DESC=[6961x7901] band3 Grid (16-bit integer)
  SUBDATASET_4_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band4
  SUBDATASET_4_DESC=[6961x7901] band4 Grid (16-bit integer)
  SUBDATASET_5_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band5
  SUBDATASET_5_DESC=[6961x7901] band5 Grid (16-bit integer)
  SUBDATASET_6_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band7
  SUBDATASET_6_DESC=[6961x7901] band7 Grid (16-bit integer)
  SUBDATASET_7_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:atmos_opacity
  SUBDATASET_7_DESC=[6961x7901] atmos_opacity Grid (16-bit integer)
  SUBDATASET_8_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:lndsr_QA
  SUBDATASET_8_DESC=[6961x7901] lndsr_QA Grid (16-bit integer)
  SUBDATASET_9_NAME=HDF4_EOS:EOS_GRID:"lndsr.L5169054_05420101216.hdf":Grid:band6
  SUBDATASET_9_DESC=[6961x7901] band6 Grid (16-bit integer)
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)


HDFLooK:

Global attributes: 26
DataProvider: UMDLANDSAT_5TM2010-12-16T00:00:00.000000Z2012-02-22T00:00:00.000000Z	SolarZenith:	42.0201
	SolarAzimuth:	141.424
2	WRS_Path:	169
	WRS_Row:	54
	NumberOfBands:	6
	BandNumbers:	1,  2,  3,  4,  5,  7
L5TSRL5TSR.a2010350.w2p169r054.020.2012117105207.hdf2012-04-26T10:52:07Z1.0.91.0.4	WestBoundingCoordinate:	36.62
	EastBoundingCoordinate:	38.7807
	NorthBoundingCoordinate:	9.61827
	SouthBoundingCoordinate:	7.72286
	OrientationAngle:	0
	PixelSize:	30
HDFEOS_V2.4-----------------------------------------------------------------------------------
GROUP=SwathStructure
END_GROUP=SwathStructure
GROUP=GridStructure
	GROUP=GRID_1
		GridName="Grid"
		XDim=7901
		YDim=6961
		UpperLeftPointMtrs=(238785.000000,1063215.000000)
		LowerRightMtrs=(475815.000000,854385.000000)
		Projection=GCTP_UTM
		ZoneCode=37
		SphereCode=12
		GridOrigin=HDFE_GD_UL
		GROUP=Dimension
		END_GROUP=Dimension
		GROUP=DataField
			OBJECT=DataField_0
				DataFieldName="band1"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_0
			OBJECT=DataField_1
				DataFieldName="band2"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_1
			OBJECT=DataField_2
				DataFieldName="band3"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_2
			OBJECT=DataField_3
				DataFieldName="band4"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_3
			OBJECT=DataField_4
				DataFieldName="band5"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_4
			OBJECT=DataField_5
				DataFieldName="band7"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_5
			OBJECT=DataField_6
				DataFieldName="atmos_opacity"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_6
			OBJECT=DataField_7
				DataFieldName="lndsr_QA"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_7
			OBJECT=DataField_9
				DataFieldName="band6"
				DataType=DFNT_INT16
				DimList=("YDim","XDim")
			END_OBJECT=DataField_9
		END_GROUP=DataField
		GROUP=MergedFields
		END_GROUP=MergedFields
	END_GROUP=GRID_1
END_GROUP=GridStructure
GROUP=PointStructure
END_GROUP=PointStructure
END
-----------------------------------------------------------------------------------
CMReflectanceBasedv1.0


More information about the Qgis-developer mailing list