<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span></span></div><div><br></div>  <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> ----- Mensaje reenviado -----<br>  <b><span style="font-weight:bold;">De:</span></b> David Ortiz &lt;dvdov@yahoo.com&gt;<br> <b><span style="font-weight: bold;">Para:</span></b> "chaitanya.ch@gmail.com" &lt;chaitanya.ch@gmail.com&gt; <br> <b><span style="font-weight: bold;">Enviado:</span></b> Miércoles, febrero 8, 2012 12:12 A.M.<br> <b><span style="font-weight: bold;">Asunto:</span></b> Re: [gdal-dev] HDF-EOS AIRS CO2 images without georeference using gdal_translate<br> </font> </div> <br><div id="yiv1004822053"><div><div style="color: rgb(0, 0, 0);
 background-color: rgb(255, 255, 255); font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "><div><span>Hi Chaitanya,</span></div><div><span><br></span></div><div><span>I have some questions:</span></div><div><span><br></span></div><div>1. &nbsp;Is there some way to get the residuals or errors from georeferencing when i use gdalwarp with any image?</div><div><br></div><div>2. &nbsp;Currently , i have a particular Geotiff image without spatial reference but i know it has the same features of AIRS radiances product (AIRIBRAD product &nbsp; &nbsp; &nbsp;from&nbsp;http://mirador.gsfc.nasa.gov/cgi-bin/mirador/homepageAlt.pl?keyword=AIRIBRAD for 2009.01.18. date) like&nbsp;coordinate system,resolution etc, but &nbsp; &nbsp; &nbsp;with other values for each pixel. I want set geolocation information to Geotiff from radiances file so,</div><div>&nbsp;&nbsp; &nbsp; Is posible associate &nbsp;geolocation arrays from hdf file
 (latitude,longitude subdatasets) with pixels of Geotiff image without using "-gcp" or "--optfile"options?</div><div>&nbsp;&nbsp; &nbsp; may be extracting arrays to GeoTiff format first?</div><div><br></div><div>3. &nbsp;I tried the command:</div><div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; gdal_translate -a_srs "EPSG:4326"&nbsp;--optfile GCPs.txt inputfile.tif outputfile.tif, &nbsp; then</div><div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; gdalwarp outputfile.tif outputfile_1.if</div><div>&nbsp;&nbsp; &nbsp; where GCPs.txt file has 12200 gcps. I get:<br></div><div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;Input file size is 90, 135</div><div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0ERROR 1: TIFFFetchNormalTag:Incorrect count for "GeoTiePoints"; tag ignored</div><div>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...10...20...30...40...50...60...70...80...90...100 - done.</div><div>&nbsp;&nbsp; &nbsp; &nbsp;it worked with 10800
 gcps.</div><div>&nbsp;&nbsp; &nbsp; &nbsp;Is there a limit number of gcps &nbsp;??</div><div><br></div><div>4. Image with 10800 gcps does not have the same size of radiances image provided by HEG tool (Geotiff) when both are loaded into GIS software.</div><div>&nbsp;&nbsp; &nbsp;Do you know why ?</div><div><br></div><div><br></div><div>Below i describe some lines of AIRS radiances and particular geotiff images:</div><div><br></div><div>gdalinfo HDF4_EOS:EOS_SWATH:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":L1B_AIRS_Science:radiances<br></div><div><div>SWATH_TITLE=
 =L1B_AIRS_Science</div><div>...</div><div>...</div><div>Geolocation:</div><div>&nbsp;&nbsp;LINE_OFFSET=0</div><div>&nbsp;&nbsp;LINE_STEP=1</div><div>&nbsp;&nbsp;PIXEL_OFFSET=0</div><div>&nbsp;&nbsp;PIXEL_STEP=1</div><div>&nbsp;&nbsp;SRS=</div><div>&nbsp;&nbsp;X_BAND=1</div><div>&nbsp;&nbsp;X_DATASET=HDF4_EOS:EOS_SWATH_GEOL:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":L1B_AIRS_Science:Longitude</div><div>&nbsp;&nbsp;Y_BAND=1</div><div>&nbsp;&nbsp;Y_DATASET=HDF4_EOS:EOS_SWATH_GEOL:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":L1B_AIRS_Science:Latitude</div><div>Corner Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div><div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; 90.0)</div><div>Upper Right ( 2378.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( 2378.0, &nbsp; 90.0)</div><div>Center &nbsp; &nbsp; &nbsp;( 1189.0, &nbsp;
 45.0)</div><div>&nbsp;--- Here i note the dimensions are wrong because radiances is a dataset with 135*90 &nbsp;pixels and 2378 bands</div><div><br></div><div>gdalinfo&nbsp;<br></div><div><div>HDF4_SDS:UNKNOWN:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":3 (radiances)</div><div>...</div><div>...</div><div>Geolocation:</div><div>&nbsp;&nbsp;LINE_OFFSET=0</div><div>&nbsp;&nbsp;LINE_STEP=1</div><div>&nbsp;&nbsp;PIXEL_OFFSET=0</div><div>&nbsp;&nbsp;PIXEL_STEP=1</div><div>&nbsp;&nbsp;SRS=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"]]</div><div>&nbsp;&nbsp;X_BAND=1</div><div>&nbsp;&nbsp;X_DATASET=HDF4_SDS:UNKNOWN:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":1</div><div>&nbsp;&nbsp;Y_BAND=1</div><div>&nbsp;&nbsp;Y_DATASET=HDF4_SDS:UNKNOWN:"/home/david/Descargas/2009.01.18/HDFs/AIRS.2009.01.18.236.L1B.AIRS_Rad.v5.0.0.0.G09020142509.hdf":0</div><div>Corner Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div><div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; 90.0)</div><div>Upper Right ( 2378.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( 2378.0, &nbsp; 90.0)</div><div>Center &nbsp; &nbsp; &nbsp;( 1189.0, &nbsp; 45.0)</div><div><br></div><div><br></div><div><div># gdalinfo
 /home/david/Descargas/2009.01.18/Sin_referencia/Airs_referencia/GCPs/CO2_particularGTiff_2009.01.18.tif</div><div>Driver: GTiff/GeoTIFF</div><div>Files: /home/david/Descargas/2009.01.18/Sin_referencia/uwairs_Airs_referencia/GCPs/CO2_particularGTiff_2009.01.18.tif</div><div>Size is 90, 135</div><div>Coordinate System is `'</div><div>Image Structure Metadata:</div><div>&nbsp;&nbsp;INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div><div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp;135.0)</div><div>Upper Right ( &nbsp; 90.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( &nbsp; 90.0, &nbsp;135.0)</div><div>Center &nbsp; &nbsp; &nbsp;( &nbsp; 45.0, &nbsp; 67.5)</div><div>Band 1 Block=90x22 Type=Float32, ColorInterp=Gray</div><div><br></div><div><br></div><div>Thanks in advance</div><div>Best Regards</div><div><br></div><div><br></div><div>David A. Ortiz
 V</div><div><br></div></div></div></div><div><br></div><div><br></div><div><br></div>  <div style="font-size: 12pt; font-family: times, serif; "> <div style="font-size: 12pt; font-family: times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1">  <b><span style="font-weight:bold;">De:</span></b> David Ortiz &lt;dvdov@yahoo.com&gt;<br> <b><span style="font-weight:bold;">Para:</span></b> Chaitanya kumar CH &lt;chaitanya.ch@gmail.com&gt; <br><b><span style="font-weight:bold;">CC:</span></b> "gdal-dev@lists.osgeo.org" &lt;gdal-dev@lists.osgeo.org&gt; <br> <b><span style="font-weight:bold;">Enviado:</span></b> Sábado, enero 28, 2012 10:59 P.M.<br> <b><span style="font-weight:bold;">Asunto:</span></b> Re: [gdal-dev] HDF-EOS AIRS CO2 images without georeference using gdal_translate<br> </font> </div> <br><div id="yiv1004822053"><div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-size: 12pt; font-family: times,
 serif; "><div><span>Chaitanya,</span></div><div><span><br></span></div><div><span>i &nbsp;used gdalwarp with several values for <i>tolerance</i> and <i>minimum_gcps </i><span>respectively</span>&nbsp;:</span></div><div><span>-refine_gcps 1 10, &nbsp;1 35, 1 41, 200 10, 200 35, 200 41, 329 10, 329 35, 329 41, 9000 10, 9000 41&nbsp;<br></span></div><div><span>and loaded the images in Gis software with layer of points whose coordinates are from geolocation fields (Latitude, Longitude) in hdf file in order to compare</span></div><div>pixels against points.</div><div><br></div><div>when parameters are 1 41, 200 41, 329 41,9000 41, 9000 10 , most of points fall within its corresponding pixel. for other values, the images have a displacement towards west. Then,images with 1 10, 200 10, 329 10 values for parameters&nbsp;<span>&nbsp;are overlapped
 as well as images with 1 35, 200 35, 329 35, and 1 41, 200 41, 329 41, 9000 10, 9000 41. &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;</span>&nbsp;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;</div><div><br></div><div><span><br></span></div><div><span><br></span></div><div><span><br></span></div><div><span><br></span></div><div><span><br></span></div><div><br></div>  <div style="font-size: 12pt; font-family: times, serif; "> <div style="font-size: 12pt; font-family: times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1">  <b><span style="font-weight:bold;">De:</span></b> Chaitanya kumar CH &lt;chaitanya.ch@gmail.com&gt;<br> <b><span style="font-weight:bold;">Para:</span></b> David Ortiz &lt;dvdov@yahoo.com&gt; <br><b><span style="font-weight:bold;">CC:</span></b> "gdal-dev@lists.osgeo.org" &lt;gdal-dev@lists.osgeo.org&gt; <br> <b><span style="font-weight:bold;">Enviado:</span></b> Viernes, enero 27, 2012 11:09 P.M.<br>
 <b><span style="font-weight:bold;">Asunto:</span></b> Re: [gdal-dev] HDF-EOS AIRS CO2 images without georeference using gdal_translate<br> </font> </div> <br><div id="yiv1004822053">David,<br><br>gdalwarp should stop eliminating outliers when it can't find any more. But to be sure you can pass the minimum number of GCPs. Also, the argument <i>tolerance</i> is the tolerance in GCP's SRS units. Since the outliers are indicated by (-9999,-9999) and our GCPs can go between -180 and 180, anything between 360 and 9000 can be given as tolerance.<br>
<br>BTW, there is no need for -s_srs because we are using GCPs.<br><br>Do let us know how it goes.<br><br><div class="yiv1004822053gmail_quote">On Sat, Jan 28, 2012 at 6:57 AM, David Ortiz <span dir="ltr">&lt;<a rel="nofollow" ymailto="mailto:dvdov@yahoo.com" target="_blank" href="mailto:dvdov@yahoo.com">dvdov@yahoo.com</a>&gt;</span> wrote:<br>
<blockquote class="yiv1004822053gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="font-size: 12pt; font-family: times, serif; "><div style="font-size:12pt;">
<span><span style="font-size: 13px; font-family: Arial; ">Chaitanya,</span><br></span></div><div><font face="Arial"><br></font></div><div><font face="Arial">This image have 329 GCPs and 41 of them are not outliers , can i put&nbsp;
<span style="font-size:16px;">-refine_gcps 329 41 </span></font><span style="font-size: 13px; font-family: Arial; ">in the command ?</span></div><div><span style="font-size: 13px; font-family: Arial; "><br>
</span></div><div><span style="font-size: 13px; font-family: Arial; "><br></span></div><div><span style="font-size: 13px; font-family: Arial; "><br></span></div><div><span style="font-size: 13px; font-family: Arial; ">--</span></div><div>
<span style="font-size: 13px; font-family: Arial; ">David A. Ortiz V</span></div><div style="font-size:12pt;"><br></div><div style="font-size:12pt;">
<br></div><div style="font-size:12pt;"> <div style="font-size:12pt;">
 <div dir="ltr"> <font face="Arial"><div class="yiv1004822053im"> <hr size="1">  <b><span style="font-weight:bold;">De:</span></b> Chaitanya kumar CH &lt;<a rel="nofollow" ymailto="mailto:chaitanya.ch@gmail.com" target="_blank" href="mailto:chaitanya.ch@gmail.com">chaitanya.ch@gmail.com</a>&gt;<br>
 <b><span style="font-weight:bold;">Para:</span></b> David Ortiz &lt;<a rel="nofollow" ymailto="mailto:dvdov@yahoo.com" target="_blank" href="mailto:dvdov@yahoo.com">dvdov@yahoo.com</a>&gt; <br><b><span style="font-weight:bold;">CC:</span></b> "<a rel="nofollow" ymailto="mailto:gdal-dev@lists.osgeo.org" target="_blank" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>" &lt;<a rel="nofollow" ymailto="mailto:gdal-dev@lists.osgeo.org" target="_blank" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a>&gt; <br>
 </div><b><span style="font-weight:bold;">Enviado:</span></b> Jueves, enero 26, 2012 10:49 P.M.<div class="yiv1004822053im"><br> <b><span style="font-weight:bold;">Asunto:</span></b> Re: [gdal-dev] HDF-EOS AIRS CO2 images without georeference using
 gdal_translate<br> </div></font> </div><div><div class="yiv1004822053h5"> <br><div>David,<br><br>I checked the images I created now. The GCPs were not considered. So, I made sure they were used with gdalwarp.<br><span style="">gdalwarp -s_srs "EPSG:4326" -t_srs "EPSG:4326" -refine_gcps 200 10 -of GTiff HDF4_EOS:EOS_SWATH:"AIRS.2008.03.12.228.L2.CO2_Sup.v5.4.11.0.CO2.T09254115637.hdf":CO2:CO2ret&nbsp; CO2_sup_Gtiff_EOS_warp.tif</span><br>


I used the -refine_gcps. It will help in excluding the dummy gcps 
pointing to (-9999,-9999,0). In this option, the number of minimum_gcps 
had to be given. I chose the number 10 based on the number of usable 
GCPs available. It should not be more than the number of usable GCPs. If I did not exclude the dummy gcps, I got black image.<br>
<br>Also note that the output image dimensions changed.<br><br><div>On Fri, Jan 27, 2012 at 7:26 AM, David Ortiz <span dir="ltr">&lt;<a rel="nofollow" ymailto="mailto:dvdov@yahoo.com" target="_blank" href="mailto:dvdov@yahoo.com">dvdov@yahoo.com</a>&gt;</span> wrote:<br>

<blockquote style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="font-size:12pt;"><div style="font-size:12pt;">
<span>Hi Chaitanya</span></div><div style="font-size:12pt;"><span><br></span></div><div style="font-size:12pt;">
<span style="font-size:12pt;"><br></span></div><div style="font-size:12pt;"><span style="font-size:12pt;">Applying&nbsp;</span><span style="font-size:12pt;">gdalinfo -stats to input and output images i get coherent values &nbsp;for statistics in both .img &nbsp;and GTiff &nbsp;images. The problem was in data display in GIS since QGIS show a Gtiff image with all pixel values equal to zero when i use "identify spatial objects" tool and ArcGIS
 &nbsp;doesn't display it correctly , just two lines and no data, but&nbsp;</span><font size="3">gvSIG &nbsp;allows see the&nbsp;</font>appropriate<font size="3">&nbsp;values of CO2 amounts.</font></div><div style="font-size:12pt;">
<br></div><div style="font-size:12pt;">However, neither the .img nor GTiff have geographical spatial reference in WGS84 datum coordinates. Even opening &nbsp;subdatasets directly&nbsp;from hdf file&nbsp;<span style="font-size:12pt;">with ArcGis , the images are missing georeference. As it is seen below, hdf image is between -133.57, -157.94 longitudes and 22.23, 45.88 latitudes. I mentioned HDF_SDS:UNKNOWN datasets on previous post because it have some info about coordinate reference system:</span></div>

<div><div style="font-size:12pt;">SRS=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"]]<br>

</div></div><div style="font-size:12pt;"><span style="font-size:12pt;">Actually , when i load Geotiff generated with HEGTool&nbsp;</span><span style="font-size:12pt;">(with other AIRS products)&nbsp;</span><span style="font-size:12pt;">in any gis software , automatically have wgs84 Datum coordinates</span><br>

</div><div style="font-size:12pt;"><span style="font-size:12pt;"><br></span></div><div style="font-size:12pt;">
<span style="font-size:12pt;">Here the results of gdalinfo command:</span></div><div style="font-size:12pt;"><span style="font-size:12pt;">1. Input image</span><br>
</div><div><div><br></div><div>gdalinfo -stats HDF4_EOS:EOS_SWATH:"/home/david/Descargas/pruebas/AIRS.2008.03.12.228.L2.CO2_Sup.v5.4.11.0.CO2.T09254115637.hdf":CO2:CO2ret&nbsp;</div><div>Driver: HDF4Image/HDF4 Dataset</div>

<div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...</div><div><div>Corner Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div><div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; 22.0)</div><div>Upper Right ( &nbsp; 15.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( &nbsp; 15.0, &nbsp; 22.0)</div>

<div>Center &nbsp; &nbsp; &nbsp;( &nbsp; &nbsp;7.5, &nbsp; 11.0)</div><div>Band 1 Block=15x22 Type=Float32, ColorInterp=Gray</div></div><div>&nbsp; Min=379.197 Max=393.697&nbsp;</div><div>&nbsp; Minimum=379.197, Maximum=393.697, Mean=386.463,
 StdDev=3.071</div><div>&nbsp; NoData Value=-9999</div><div>&nbsp; Metadata:</div><div>&nbsp; &nbsp; STATISTICS_MAXIMUM=393.6969909668</div><div>&nbsp; &nbsp; STATISTICS_MEAN=386.46261742001</div><div>&nbsp; &nbsp; STATISTICS_MINIMUM=379.1969909668</div><div>&nbsp; &nbsp; STATISTICS_STDDEV=3.0709672735453</div>

<div><br></div><div>2. &nbsp;with&nbsp;Geotiff&nbsp;output &nbsp;image :&nbsp;</div><div><div>gdalinfo -stats &nbsp;/home/david/Descargas/pruebas/EOS/CO2_sup_Gtiff_EOS_traslate.tif</div><div>Driver: GTiff/GeoTIFF</div><div>Files: /home/david/Descargas/pruebas/EOS/CO2_sup_Gtiff_EOS_traslate.tif</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp;/home/david/Descargas/pruebas/EOS/CO2_sup_Gtiff_EOS_traslate.tif.aux.xml</div><div><div>Size is 15, 22</div><div>Coordinate System is `'</div></div><div>GCP Projection =&nbsp;</div><div>GCP[ &nbsp;0]: Id=1, Info=</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (0.5,0.5) -&gt; (-9999,-9999,0)</div><div>GCP[ &nbsp;1]: Id=2,
 Info=</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (1.5,0.5) -&gt; (-9999,-9999,0)</div><div>GCP[ &nbsp;2]: Id=3, Info=</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (2.5,0.5) -&gt; (-9999,-9999,0)</div><div>GCP[ &nbsp;3]: Id=4, Info=</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (3.5,0.5) -&gt; (-145.229995727539,23.3500003814697,0)</div>

<div>GCP[ &nbsp;4]: Id=5, Info=</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; (4.5,0.5) -&gt; (-9999,-9999,0)</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...&nbsp;</div><div><div>GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= HORIZONTALSPATIALDOMAINCONTAINER</div><div>&nbsp; &nbsp; &nbsp; GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= BOUNDINGRECTANGLE</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = EASTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp;
 &nbsp; &nbsp; &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= -133.537399291992</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = EASTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = WESTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= -157.941314697266</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = WESTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = SOUTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
 &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 22.2317676544189</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = SOUTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = NORTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 45.889762878418</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = NORTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; END_GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= BOUNDINGRECTANGLE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...&nbsp;</div></div>

<div>&nbsp; INTERLEAVE=BAND</div><div><div>Corner
 Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div><div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; 22.0)</div><div>Upper Right ( &nbsp; 15.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( &nbsp; 15.0, &nbsp; 22.0)</div><div>Center &nbsp; &nbsp; &nbsp;( &nbsp; &nbsp;7.5, &nbsp; 11.0)</div><div>

Band 1 Block=15x22 Type=Float32, ColorInterp=Gray</div></div><div>&nbsp; Minimum=379.197, Maximum=393.697, Mean=386.463, StdDev=3.071</div><div>&nbsp; NoData Value=-9999</div><div>&nbsp; Metadata:</div><div>&nbsp; &nbsp; STATISTICS_MAXIMUM=393.6969909668</div>

<div>&nbsp; &nbsp; STATISTICS_MEAN=386.46261742001</div><div>&nbsp; &nbsp; STATISTICS_MINIMUM=379.1969909668</div><div>&nbsp; &nbsp; STATISTICS_STDDEV=3.0709672735453</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;&nbsp;</div><div>3. .img image :</div><div><br></div><div>gdalinfo -stats /home/david/Descargas/pruebas/EOS/CO2_sup_img_EOS.img</div>

<div>Driver: HFA/Erdas Imagine
 Images (.img)</div><div>Files: /home/david/Descargas/pruebas/EOS/CO2_sup_img_EOS.img</div><div>&nbsp; &nbsp; &nbsp; &nbsp;/home/david/Descargas/pruebas/EOS/CO2_sup_img_EOS.img.aux.xml</div><div><div>Size is 15, 22</div><div>Coordinate System is `'</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...&nbsp;</div></div><div><div>GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= SPATIALDOMAINCONTAINER</div><div>&nbsp; &nbsp; GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= HORIZONTALSPATIALDOMAINCONTAINER</div><div>&nbsp; &nbsp; &nbsp; GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= BOUNDINGRECTANGLE</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = EASTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp;
 &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= -133.537399291992</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = EASTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = WESTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= -157.941314697266</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = WESTBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = SOUTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;
 &nbsp; &nbsp; &nbsp;= 22.2317676544189</div><div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = SOUTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = NORTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; NUM_VAL &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 1</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; VALUE &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= 45.889762878418</div>

<div>&nbsp; &nbsp; &nbsp; &nbsp; END_OBJECT &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = NORTHBOUNDINGCOORDINATE</div><div>&nbsp; &nbsp; &nbsp; END_GROUP &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;= BOUNDINGRECTANGLE</div><div>&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;...&nbsp;</div></div><div><div>Corner Coordinates:</div><div>Upper Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp; &nbsp;0.0)</div>

<div>Lower Left &nbsp;( &nbsp; &nbsp;0.0, &nbsp;
 22.0)</div><div>Upper Right ( &nbsp; 15.0, &nbsp; &nbsp;0.0)</div><div>Lower Right ( &nbsp; 15.0, &nbsp; 22.0)</div><div>Center &nbsp; &nbsp; &nbsp;( &nbsp; &nbsp;7.5, &nbsp; 11.0)</div></div><div>Band 1 Block=64x64 Type=Float32, ColorInterp=Undefined</div><div>&nbsp; Description = Layer_1</div>

<div>&nbsp; Minimum=379.197, Maximum=393.697, Mean=386.463, StdDev=3.071</div><div>&nbsp; NoData Value=-9999</div><div>&nbsp; Metadata:</div><div>&nbsp; &nbsp; LAYER_TYPE=athematic</div><div>&nbsp; &nbsp; STATISTICS_MAXIMUM=393.6969909668</div><div>&nbsp; &nbsp; STATISTICS_MEAN=386.46261742001</div>

<div>&nbsp; &nbsp; STATISTICS_MINIMUM=379.1969909668</div><div>&nbsp; &nbsp; STATISTICS_STDDEV=3.0709672735453</div><div><br></div><div><br></div><div><br></div><div>What is wrong?. I tried with gdal 1.8 and 1.9 with hdf4 support.</div><div>

<br></div><div><br></div><div>Best regards.</div></div></div></div><span><font color="#888888"><div style="font-size:12pt;">--</div><div style="font-size:12pt;">
<br></div><div style="font-size:12pt;">David Ortiz V.</div></font></span><div style="font-size:12pt;">
<br></div></div></div></blockquote></div><br>
</div> </div></div></div> </div>  </div></div></blockquote></div><br clear="all"><br>-- <br>Best regards,<br>Chaitanya kumar CH.<br><br>+91-9494447584<br>17.2416N 80.1426E<br>
</div><br><br> </div> </div>  </div></div></div><br>_______________________________________________<br>gdal-dev mailing list<br><a rel="nofollow" ymailto="mailto:gdal-dev@lists.osgeo.org" target="_blank" href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br><a rel="nofollow" target="_blank" href="http://lists.osgeo.org/mailman/listinfo/gdal-dev">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br><br> </div> </div>  </div></div></div><br><br> </div> </div>  </div></body></html>