<!doctype html public "-//W3C//DTD W3 HTML//EN">
<html><head><style type="text/css"><!--
blockquote, dl, ul, ol, li { padding-top: 0 ; padding-bottom: 0 }
--></style><title>Re: [Gdal-dev] GDAL 1.2.1 Released - HDF
error</title></head><body>
<div>Until GDAL 1.2.0 I have HDF4 support and it works.</div>
<div>Now I've compiled new GDAL with new HDF-EOS support and I have
this error:</div>
<div><br></div>
<div>When I use gdalinfo for HDF file the output is similar:</div>
<div><br></div>
<div>OUTPUT FROM GDAL 1.2.0</div>
<div><br></div>
<blockquote><font face="Helvetica" size="-1" color="#000000">[lor:~]
lor% gdalinfo
/Volumes/ASTER_files/AST_L1B_003_01102002101254_05222002150244.hdf</font
></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Driver:
HDF4/Hierarchical Data Format Release 4</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Size is
512, 512</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Coordinate System is `'</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Metadata:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
HDFEOSVersion=HDFEOS_V2.6</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
IDOFASTERGDSDATAGRANULE=ASTL1B
0201101012540205220094</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
QAPERCENTINTERPOLATEDDATA=0.000000</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
REPROCESSINGACTUAL=not reprocessed</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PGEVERSION=04.00R06</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PROCESSINGLEVELID=1B</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
MAPPROJECTIONNAME=Universal Transverse Mercator</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Subdatasets:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_1_NAME=HDF4_SDS:ASTER_L1B:"/Volumes/PB_APPLI/LAVORI/M<span
></span
>APPE/ASTER_FO_gen02/AST_L1B_003_01102002101254_05222002150244.hdf<span
></span>":2</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_1_DESC=[4200x4980] ImageData1 (8-bit unsigned
integer)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_2_NAME=HDF4_SDS:ASTER_L1B:"/Volumes/PB_APPLI/LAVORI/M<span
></span
>APPE/ASTER_FO_gen02/AST_L1B_003_01102002101254_05222002150244.hdf<span
></span>":3</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_2_DESC=[4200x4980] ImageData2 (8-bit unsigned
integer)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Corner
Coordinates:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Left ( 0.0,
0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Left ( 0.0, 512.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Right ( 512.0, 0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Right ( 512.0, 512.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Center ( 256.0,
256.0)</font></blockquote>
<div><br></div>
<div><br></div>
<div>OUTPUT FROM GDAL 1.2.1</div>
<div><br></div>
<blockquote><font face="Helvetica" size="-1" color="#000000">[lor:~]
lor% gdalinfo
/Volumes/ASTER_files/AST_L1B_003_01102002101254_05222002150244.hdf</font
></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Driver:
HDF4/Hierarchical Data Format Release 4</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Size is
512, 512</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Coordinate System is `'</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Metadata:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
HDFEOSVersion=HDFEOS_V2.6</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
IDOFASTERGDSDATAGRANULE=ASTL1B
0201101012540205220094</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
QAPERCENTINTERPOLATEDDATA=0.000000</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
REPROCESSINGACTUAL=not reprocessed</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PGEVERSION=04.00R06</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PROCESSINGLEVELID=1B</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
MAPPROJECTIONNAME=Universal Transverse Mercator</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Subdatasets:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_1_NAME=HDF4_EOS:EOS_SWATH:"/Volumes/PB_APPLI/LAVORI/M<span
></span
>APPE/ASTER_FO_gen02/AST_L1B_003_01102002101254_05222002150244.hdf<span
></span>":VNIR_Swath:ImageData1</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_1_DESC=[4200x4980] ImageData1 (8-bit unsigned
integer)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_2_NAME=HDF4_EOS:EOS_SWATH:"/Volumes/PB_APPLI/LAVORI/M<span
></span
>APPE/ASTER_FO_gen02/AST_L1B_003_01102002101254_05222002150244.hdf<span
></span>":VNIR_Swath:ImageData2</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SUBDATASET_2_DESC=[4200x4980] ImageData2 (8-bit unsigned
integer)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Corner
Coordinates:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Left ( 0.0,
0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Left ( 0.0, 512.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Right ( 512.0, 0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Right ( 512.0, 512.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Center ( 256.0,
256.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<div><font face="Helvetica" size="-1"
color="#000000"><br></font></div>
<div>THE DIFFERENCES ARE AT THE END.</div>
<div>Now I have a band description.</div>
<div>But the error is after when I select one band.</div>
<div><font face="Helvetica" size="-1"
color="#000000"><br></font></div>
<div><br></div>
<div>OUTPUT FROM GDAL 1.2.0</div>
<div><font face="Helvetica" size="-1"
color="#000000"><br></font></div>
<blockquote><font face="Helvetica" size="-1" color="#000000">[lor:~]
lor% gdalinfo
HDF4_SDS:ASTER_L1B:"/Volumes/ASTER_files/AST_L1B_003_01102002101<span
></span>254_05222002150244.hdf":2</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Driver:
HDF4Image/HDF4 Dataset</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Size is
4980, 4200</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Coordinate System is:</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">PROJCS["UTM Zone 33, Northern
Hemisphere",</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"> GEOGCS["WGS
84",</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
DATUM["WGS_1984",</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
>
SPHEROID["WGS 84",6378137,298.257223563,</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
> <span
></span>
AUTHORITY["EPSG","7030"]],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
>
TOWGS84[0,0,0,0,0,0,0],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
>
AUTHORITY["EPSG","6326"]],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PRIMEM["Greenwich",0,</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
>
AUTHORITY["EPSG","8901"]],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
UNIT["degree",0.0174532925199433,</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"
>
AUTHORITY["EPSG","9108"]],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
AXIS["Lat",NORTH],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
AXIS["Long",EAST],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
AUTHORITY["EPSG","4326"]],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PROJECTION["Transverse_Mercator"],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PARAMETER["latitude_of_origin",0],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PARAMETER["central_meridian",15],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PARAMETER["scale_factor",0.9996],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PARAMETER["false_easting",500000],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
PARAMETER["false_northing",0],</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
UNIT["Meter",1]]</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">GeoTransform =</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
273342.0580804553, 14.62007381662245,
-3.354636379426534</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
4820040.394499515, -3.354618932788324,
-14.6200685003371</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">GCP
Projection = GEOGCS["WGS
84",DATUM["WGS_1984",SPHEROID["WGS
84",6378137,298.257223563,AUTHORITY["EPSG","7030<span
></span
>"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","632<span
></span
>6"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG",<span
></span
>"8901"]],UNIT["degree",0.0174532925199433,AUTHOR<span
></span
>ITY["EPSG","9108"]],AXIS["Lat",NORTH],<span
></span
>AXIS["Long",EAST],AUTHORITY["EPSG","4326<span
></span>"]]</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">GCP[
0]: Id=, Info=</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
(0.5,0.5) -> (12.1965,43.4988,0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">GCP[
1]: Id=, Info=</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
(498.5,0.5) -> (12.2872,43.486,0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">GCP[
2]: Id=, Info=</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">
(996.5,0.5) -> (12.3778,43.4731,0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PROJECTIONPARAMETERS14=6378137.000000, 6356752.300000, 0.999600,
0.000000, 0.261799, 0.000000, 500000.000000, 0.000000, 0.000000,
0.000000, 0.000000, 0.000000, 0.000000</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
UTMZONECODE14=33</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
INCL14=0.005225</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
INSTRUMENTSHORTNAME=ASTER</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
WESTBOUNDINGCOORDINATE=12.049274</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
NORTHBOUNDINGCOORDINATE=43.498844</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
EASTBOUNDINGCOORDINATE=13.101212</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
SOUTHBOUNDINGCOORDINATE=42.811432</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
TIMEOFDAY=101254443000Z</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">.....</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
PROCESSINGLEVELID=1B</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">
MAPPROJECTIONNAME=Universal Transverse Mercator</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Corner
Coordinates:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Left ( 273342.058, 4820040.394) ( 12d11'47.31"E,
43d29'56.12"N)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Left ( 259252.585, 4758636.107) ( 12d 2'57.12"E,
42d56'32.39"N)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Right ( 346150.026, 4803334.392) ( 13d 6'4.10"E,
43d22'1.61"N)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Right ( 332060.553, 4741930.105) ( 12d56'45.26"E,
42d48'41.44"N)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Center (
302701.305, 4780985.250) ( 12d34'23.91"E, 43d
9'21.21"N)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Band 1
Block=4980x1 Type=Byte, ColorInterp=Gray</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<div><br></div>
<div>OUTPUT FROM GDAL 1.2.1</div>
<div><font face="Helvetica" size="-1"
color="#000000"><br></font></div>
<blockquote><font face="Helvetica" size="-1" color="#000000">[lor:~]
lor% gdalinfo
HDF4_EOS:EOS_SWATH:"/Volumes/ASTER_files/AST_L1B_003_01102002101<span
></span>254_05222002150244.hdf":VNIR_Swath:ImageData1</font></blockquote
>
<blockquote><font face="Helvetica" size="-1" color="#000000">Driver:
HDF4Image/HDF4 Dataset</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Size is
4980, 4200</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Coordinate System is `'</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Corner
Coordinates:</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Left ( 0.0,
0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Left ( 0.0, 4200.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Upper
Right ( 4980.0, 0.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Lower
Right ( 4980.0, 4200.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000">Center ( 2490.0,
2100.0)</font></blockquote>
<blockquote><font face="Helvetica" size="-1" color="#000000">Band 1
Block=4980x1 Type=Byte, ColorInterp=Gray</font></blockquote>
<blockquote><font face="Helvetica" size="-1"
color="#000000"><br></font></blockquote>
<div><font face="Helvetica" size="-1"
color="#000000"><br></font></div>
<div>Now I've lost the projection and now I can not translate with
gdalwarp for importing in grass. In Grass 5.7 there is a module
r.in.aster from Michael Burton based on this method and it's very easy
method for importing ASTER files in GRASS 5.7.</div>
<div><br></div>
<div>What is the solution ?</div>
<div><br></div>
<x-sigsep><pre>--
</pre></x-sigsep>
<div><font face="Verdana"
color="#000000"
>____________________________________________________________________<span
></span>____</font></div>
<div><font face="Courier New"
color="#000000">|| Lorenzo
Moretti e-mail:
lorenzo.moretti@bologna.enea.it </font></div
>
<div><font face="Courier New" color="#000000">||/|/| ENEA
prot/idr Web:
http://wwwamb.bologna.enea.it/ </font></div>
<div><font face="Courier New" color="#000000">|| |
via Don Fiammelli, 2 FTP:
ftp://ftpamb.bologna.enea.it/ (ris.)</font></div>
<div><font face="Courier New" color="#000000">~~~~~~ 40128
BOLOGNA - ITALY Ph: +39-0516098086 Fax:
+39-0516098131</font></div>
<div><font face="Courier New"
color="#000000"
>____________________________________________________________________<span
></span>____</font></div>
</body>
</html>