[Gdal-dev] HDF4 format read in - command line ?
Hannes Isaak Reuter
hannes.reuter at jrc.it
Tue May 8 09:29:10 EDT 2007
Today i got an HDF4 format dataset to help on and somehow I'm to
stupid to tell my colleagues to tell how to read it in..
Usually I was reading it using the gdalinfo,
identifying the
SUBDATASET_1_NAME
and than use that for the gdalwarp process..
with that HDF file now i have no idea how to call the substring. and
after playing around for 3 hours with different combinations i
deccide to call for help..
If anybody can tell me the file/band string from the gdalinfo command
given below, i'm quite happy ..
U:\share\Luis>gdalinfo GEO2004001.hdf
Driver: HDF4Image/HDF4 Dataset
Size is 11277, 3086
Coordinate System is:
GEOGCS["Unknown datum based upon the Clarke 1866 ellipsoid",
DATUM["Not specified (based on Clarke 1866 spheroid)",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
- Ignored:
AUTHORITY["EPSG","7008"]]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]
Origin = (-29.238043888889045,70.000000000000000)
Pixel Size = (0.012964385918241,-0.012964386026500)
Metadata:
HDFEOSVersion=HDFEOS_V2.6
OldStructMetadata.0=GROUP=SwathStructure
END_GROUP=SwathStructure
GROUP=GridStructure
GROUP=GRID_1
GridName="MOD_Grid_monthly_1km_VI"
XDim=6000
YDim=4800
UpperLeftPointMtrs=(-1111950.519667,7783653.637667)
LowerRightMtrs=(4447802.078667,3335851.559000)
Projection=GCTP_SNSOID
ProjParams=(6371007.181000,0,0,0,0,0,0,0,0,0,0,0,0)
SphereCode=-1
GridOrigin=HDFE_GD_UL
GROUP=Dimension
END_GROUP=Dimension
GROUP=DataField
OBJECT=DataField_1
DataFieldName="1 km monthly NDVI"
DataType=DFNT_INT16
DimList=("YDim","XDim")
END_OBJECT=DataField_1
END_GROUP=DataField
GROUP=MergedFields
END_GROUP=MergedFields
END_GROUP=GRID_1
END_GROUP=GridStructure
GROUP=PointStructure
END_GROUP=PointStructure
END
OldArchiveMetadata.0=
GROUP = ARCHIVEDMETADATA
GROUPTYPE = MASTERGROUP
OBJECT = CHARACTERISTICBINANGULARSIZE
NUM_VAL = 1
VALUE = 30.000000
END_OBJECT = CHARACTERISTICBINANGULARSIZE
OBJECT = CHARACTERISTICBINSIZE
NUM_VAL = 1
VALUE = 926.625433
END_OBJECT = CHARACTERISTICBINSIZE
OBJECT = DATACOLUMNS
NUM_VAL = 1
VALUE = 1200
END_OBJECT = DATACOLUMNS
OBJECT = DATAROWS
NUM_VAL = 1
VALUE = 1200
END_OBJECT = DATAROWS
OBJECT = GLOBALGRIDCOLUMNS
NUM_VAL = 1
VALUE = 43200
END_OBJECT = GLOBALGRIDCOLUMNS
OBJECT = GLOBALGRIDROWS
NUM_VAL = 1
VALUE = 21600
END_OBJECT = GLOBALGRIDROWS
OBJECT = NUMBEROFDAYS
NUM_VAL = 1
VALUE = 48
END_OBJECT = NUMBEROFDAYS
OBJECT = DAYSOFYEAR
NUM_VAL = 48
VALUE = (353, 354, 355, 356, 357, 358, 359, 360,
361, 362, 36
3, 364, 365, 366, 1, 2, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,
14, 15, 16, 1
7, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32)
END_OBJECT = DAYSOFYEAR
OBJECT = NUMBEROF16DAYPERIODS
NUM_VAL = 1
VALUE = 3
END_OBJECT = NUMBEROF16DAYPERIODS
OBJECT = SPSOPARAMETERS
NUM_VAL = 4
VALUE = ("2749", "4334", "2749a", "4334a")
END_OBJECT = SPSOPARAMETERS
OBJECT = GEOANYABNORMAL
NUM_VAL = 1
VALUE = "False"
END_OBJECT = GEOANYABNORMAL
OBJECT = GEOESTMAXRMSERROR
NUM_VAL = 1
VALUE = 50.000000
END_OBJECT = GEOESTMAXRMSERROR
OBJECT = INSTRUMENTNAME
NUM_VAL = 1
VALUE = "Moderate-Resolution Imaging
SpectroRadiometer"
END_OBJECT = INSTRUMENTNAME
OBJECT = PROCESSINGCENTER
NUM_VAL = 1
VALUE = "MODAPS"
END_OBJECT = PROCESSINGCENTER
OBJECT = INPUTFILENAME
NUM_VAL = 3
VALUE =
("MOD13A2.A2004353.h17v02.004.2005008091125.hdf", "MO
D13A2.A2005001.h17v02.004.2005027190904.hdf",
"MOD13A2.A2005017.h17v02.004.20050
39040318.hdf")
END_OBJECT = INPUTFILENAME
OBJECT = PERCENTLAND
NUM_VAL = 1
VALUE = 0
END_OBJECT = PERCENTLAND
OBJECT = SEAPROCESSED
NUM_VAL = 1
VALUE = "Yes"
END_OBJECT = SEAPROCESSED
OBJECT = PROCESSINGENVIRONMENT
NUM_VAL = 1
VALUE = "IRIX64 mtvs1 6.5 10070055 IP35"
END_OBJECT = PROCESSINGENVIRONMENT
OBJECT = DESCRREVISION
NUM_VAL = 1
VALUE = "4.1"
END_OBJECT = DESCRREVISION
OBJECT = LONGNAME
NUM_VAL = 1
VALUE = "MODIS/Terra Vegetation Indices Monthly
L3 Global 1km
SIN Grid"
END_OBJECT = LONGNAME
GROUP = BOUNDINGRECTANGLE
OBJECT = NORTHBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 70.000000
END_OBJECT = NORTHBOUNDINGCOORDINATE
OBJECT = SOUTHBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 30.000000
END_OBJECT = SOUTHBOUNDINGCOORDINATE
OBJECT = EASTBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 62.24191
END_OBJECT = EASTBOUNDINGCOORDINATE
OBJECT = WESTBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = -29.238044
END_OBJECT = WESTBOUNDINGCOORDINATE
END_GROUP = BOUNDINGRECTANGLE
OBJECT = QAPERCENTPOORQ1KMMONTHNDVI
NUM_VAL = 16
VALUE = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 100)
END_OBJECT = QAPERCENTPOORQ1KMMONTHNDVI
OBJECT = QAPERCENTPOORQ1KMMONTHEVI
NUM_VAL = 16
VALUE = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 100)
END_OBJECT = QAPERCENTPOORQ1KMMONTHEVI
OBJECT = COMPOSITEMONTHNAME
NUM_VAL = 1
VALUE = "January"
END_OBJECT = COMPOSITEMONTHNAME
END_GROUP = ARCHIVEDMETADATA
END
Corner Coordinates:
Upper Left ( -29.2380439, 70.0000000) ( 29d14'16.96"W, 70d 0'0.00"N)
Lower Left ( -29.2380439, 29.9919047) ( 29d14'16.96"W, 29d59'30.86"N)
Upper Right ( 116.9613361, 70.0000000) (116d57'40.81"E, 70d 0'0.00"N)
Lower Right ( 116.9613361, 29.9919047) (116d57'40.81"E, 29d59'30.86"N)
Center ( 43.8616461, 49.9959524) ( 43d51'41.93"E, 49d59'45.43"N)
Band 1 Block=11277x1 Type=Int16, ColorInterp=Gray
U:\share\Luis>
- Done.
Received: from mx2.jrc.it (mx2.jrc.it [139.191.1.110])
by duke.maptools.org (8.12.11/8.12.10) with ESMTP id l48DQMvH008039
for <gdal-dev-request at lists.maptools.org>;
Tue, 8 May 2007 09:26:35 -0400
Received: from irelay.jrc.it (irelay.jrc.it [139.191.254.63])
by mx2.jrc.it (LMC5614Amx2) with ESMTP id l48DQEua006460
for <gdal-dev-request at lists.maptools.org>;
Tue, 8 May 2007 15:26:14 +0200 (CEST)
Received: from cheetah-1.jrc.it (cheetah-1.jrc.it [139.191.243.17])
by irelay.jrc.it (LMC5614Ainternal) with ESMTP id l48D8rq8007957
for <gdal-dev-request at lists.maptools.org>;
Tue, 8 May 2007 15:26:13 +0200 (MEST)
Received: from D01RI0302803.jrc.it (139.191.16.37) by cheetah-1.jrc.it
(7.1.006) (authenticated as hannes.reuter)
id 45C0D425001A35E6 for gdal-dev-request at lists.maptools.org;
Tue, 8 May 2007 15:25:48 +0200
Message-ID: <45C0D425001A35E6 at cheetah-1.jrc.it> (added by postmaster at jrc.it)
X-Mailer: QUALCOMM Windows Eudora Version 7.1.0.9
Date: Tue, 08 May 2007 15:26:07 +0200
To: gdal-dev-request at lists.maptools.org
From: Hannes Isaak Reuter <hannes.reuter at jrc.it>
Subject: HDF4 format read in - command line ?
Mime-Version: 1.0
Content-Type: text/plain; charset="us-ascii"; format=flowed
X-Spam-Status: No, score=-101.2 required=99.0 tests=ALL_TRUSTED,AWL,
UPPERCASE_50_75,USER_IN_WHITELIST autolearn=disabled version=3.1.1
X-Spam-Checker-Version: SpamAssassin 3.1.1 (2006-03-10) on mx2.jrc.it
Dear colleagues..
Today i got an HDF4 format dataset to help on and somehow I'm to
stupid to tell my colleagues to tell how to read it in..
Usually I was reading it using the gdalinfo,
identifying the
SUBDATASET_1_NAME
and than use that for the gdalwarp process..
with that HDF file now i have no idea how to call the substring. and
after playing around for 3 hours with different combinations i
deccide to call for help..
If anybody can tell me the file/band string from the gdalinfo command
given below, i'm quite happy ..
U:\share\Luis>gdalinfo GEO2004001.hdf
Driver: HDF4Image/HDF4 Dataset
Size is 11277, 3086
Coordinate System is:
GEOGCS["Unknown datum based upon the Clarke 1866 ellipsoid",
DATUM["Not specified (based on Clarke 1866 spheroid)",
SPHEROID["Clarke 1866",6378206.4,294.9786982139006,
AUTHORITY["EPSG","7008"]]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]]
Origin = (-29.238043888889045,70.000000000000000)
Pixel Size = (0.012964385918241,-0.012964386026500)
Metadata:
HDFEOSVersion=HDFEOS_V2.6
OldStructMetadata.0=GROUP=SwathStructure
END_GROUP=SwathStructure
GROUP=GridStructure
GROUP=GRID_1
GridName="MOD_Grid_monthly_1km_VI"
XDim=6000
YDim=4800
UpperLeftPointMtrs=(-1111950.519667,7783653.637667)
LowerRightMtrs=(4447802.078667,3335851.559000)
Projection=GCTP_SNSOID
ProjParams=(6371007.181000,0,0,0,0,0,0,0,0,0,0,0,0)
SphereCode=-1
GridOrigin=HDFE_GD_UL
GROUP=Dimension
END_GROUP=Dimension
GROUP=DataField
OBJECT=DataField_1
DataFieldName="1 km monthly NDVI"
DataType=DFNT_INT16
DimList=("YDim","XDim")
END_OBJECT=DataField_1
END_GROUP=DataField
GROUP=MergedFields
END_GROUP=MergedFields
END_GROUP=GRID_1
END_GROUP=GridStructure
GROUP=PointStructure
END_GROUP=PointStructure
END
OldArchiveMetadata.0=
GROUP = ARCHIVEDMETADATA
GROUPTYPE = MASTERGROUP
OBJECT = CHARACTERISTICBINANGULARSIZE
NUM_VAL = 1
VALUE = 30.000000
END_OBJECT = CHARACTERISTICBINANGULARSIZE
OBJECT = CHARACTERISTICBINSIZE
NUM_VAL = 1
VALUE = 926.625433
END_OBJECT = CHARACTERISTICBINSIZE
OBJECT = DATACOLUMNS
NUM_VAL = 1
VALUE = 1200
END_OBJECT = DATACOLUMNS
OBJECT = DATAROWS
NUM_VAL = 1
VALUE = 1200
END_OBJECT = DATAROWS
OBJECT = GLOBALGRIDCOLUMNS
NUM_VAL = 1
VALUE = 43200
END_OBJECT = GLOBALGRIDCOLUMNS
OBJECT = GLOBALGRIDROWS
NUM_VAL = 1
VALUE = 21600
END_OBJECT = GLOBALGRIDROWS
OBJECT = NUMBEROFDAYS
NUM_VAL = 1
VALUE = 48
END_OBJECT = NUMBEROFDAYS
OBJECT = DAYSOFYEAR
NUM_VAL = 48
VALUE = (353, 354, 355, 356, 357, 358, 359, 360,
361, 362, 36
3, 364, 365, 366, 1, 2, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13,
14, 15, 16, 1
7, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32)
END_OBJECT = DAYSOFYEAR
OBJECT = NUMBEROF16DAYPERIODS
NUM_VAL = 1
VALUE = 3
END_OBJECT = NUMBEROF16DAYPERIODS
OBJECT = SPSOPARAMETERS
NUM_VAL = 4
VALUE = ("2749", "4334", "2749a", "4334a")
END_OBJECT = SPSOPARAMETERS
OBJECT = GEOANYABNORMAL
NUM_VAL = 1
VALUE = "False"
END_OBJECT = GEOANYABNORMAL
OBJECT = GEOESTMAXRMSERROR
NUM_VAL = 1
VALUE = 50.000000
END_OBJECT = GEOESTMAXRMSERROR
OBJECT = INSTRUMENTNAME
NUM_VAL = 1
VALUE = "Moderate-Resolution Imaging SpectroRadiometer"
END_OBJECT = INSTRUMENTNAME
OBJECT = PROCESSINGCENTER
NUM_VAL = 1
VALUE = "MODAPS"
END_OBJECT = PROCESSINGCENTER
OBJECT = INPUTFILENAME
NUM_VAL = 3
VALUE =
("MOD13A2.A2004353.h17v02.004.2005008091125.hdf", "MO
D13A2.A2005001.h17v02.004.2005027190904.hdf",
"MOD13A2.A2005017.h17v02.004.20050
39040318.hdf")
END_OBJECT = INPUTFILENAME
OBJECT = PERCENTLAND
NUM_VAL = 1
VALUE = 0
END_OBJECT = PERCENTLAND
OBJECT = SEAPROCESSED
NUM_VAL = 1
VALUE = "Yes"
END_OBJECT = SEAPROCESSED
OBJECT = PROCESSINGENVIRONMENT
NUM_VAL = 1
VALUE = "IRIX64 mtvs1 6.5 10070055 IP35"
END_OBJECT = PROCESSINGENVIRONMENT
OBJECT = DESCRREVISION
NUM_VAL = 1
VALUE = "4.1"
END_OBJECT = DESCRREVISION
OBJECT = LONGNAME
NUM_VAL = 1
VALUE = "MODIS/Terra Vegetation Indices Monthly
L3 Global 1km
SIN Grid"
END_OBJECT = LONGNAME
GROUP = BOUNDINGRECTANGLE
OBJECT = NORTHBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 70.000000
END_OBJECT = NORTHBOUNDINGCOORDINATE
OBJECT = SOUTHBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 30.000000
END_OBJECT = SOUTHBOUNDINGCOORDINATE
OBJECT = EASTBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = 62.24191
END_OBJECT = EASTBOUNDINGCOORDINATE
OBJECT = WESTBOUNDINGCOORDINATE
NUM_VAL = 1
VALUE = -29.238044
END_OBJECT = WESTBOUNDINGCOORDINATE
END_GROUP = BOUNDINGRECTANGLE
OBJECT = QAPERCENTPOORQ1KMMONTHNDVI
NUM_VAL = 16
VALUE = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100)
END_OBJECT = QAPERCENTPOORQ1KMMONTHNDVI
OBJECT = QAPERCENTPOORQ1KMMONTHEVI
NUM_VAL = 16
VALUE = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 100)
END_OBJECT = QAPERCENTPOORQ1KMMONTHEVI
OBJECT = COMPOSITEMONTHNAME
NUM_VAL = 1
VALUE = "January"
END_OBJECT = COMPOSITEMONTHNAME
END_GROUP = ARCHIVEDMETADATA
END
Corner Coordinates:
Upper Left ( -29.2380439, 70.0000000) ( 29d14'16.96"W, 70d 0'0.00"N)
Lower Left ( -29.2380439, 29.9919047) ( 29d14'16.96"W, 29d59'30.86"N)
Upper Right ( 116.9613361, 70.0000000) (116d57'40.81"E, 70d 0'0.00"N)
Lower Right ( 116.9613361, 29.9919047) (116d57'40.81"E, 29d59'30.86"N)
Center ( 43.8616461, 49.9959524) ( 43d51'41.93"E, 49d59'45.43"N)
Band 1 Block=11277x1 Type=Int16, ColorInterp=Gray
U:\share\Luis>
More information about the Gdal-dev
mailing list