hello Hamish<div><br></div><div>With some support from Ricardo I realized that there are also 2 other HDF5 files, one with latitude of each pixel and the other with longitude. So, now I have a latitude dataset, a longitude and values. How can I "merge" this and import to GRASS (georreferenced and bla bla)? (THat I wasn't expecting to have).</div>
<div><br></div><div>By the way, thanks about gdalwarp enlightment.</div><div><br></div><div>Pedro</div><div><br><div class="gmail_quote">On Thu, Feb 18, 2010 at 12:17 PM, Hamish <span dir="ltr"><<a href="mailto:hamish_b@yahoo.com">hamish_b@yahoo.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Pedro Roma wrote:<br>
> By usinggdalinfo<br>
<div class="im">> HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR<br>
> I get:Driver:<br>
> HDF5Image/HDF5 Dataset<br>
> Files: none associatedSize is 1701,<br>
</div>> 651Coordinate System<br>
<div><div></div><div class="h5">> is:GEOGCS["WGS 84",<br>
> DATUM["WGS_1984",<br>
> SPHEROID["WGS 84",6378137,298.257223563,<br>
><br>
> AUTHORITY["EPSG","7030"]],<br>
> TOWGS84[0,0,0,0,0,0,0],<br>
> AUTHORITY["EPSG","6326"]],<br>
> PRIMEM["Greenwich",0,<br>
><br>
> AUTHORITY["EPSG","8901"]],<br>
> UNIT["degree",0.0174532925199433,<br>
><br>
> AUTHORITY["EPSG","9108"]],<br>
> AXIS["Lat",NORTH],<br>
> AXIS["Long",EAST],<br>
> AUTHORITY["EPSG","4326"]]Corner<br>
> Coordinates:Upper Left ( 0.0,<br>
> 0.0)Lower Left ( 0.0,<br>
> 651.0)Upper Right ( 1701.0, 0.0)<br>
> Lower Right ( 1701.0, 651.0)Center<br>
> ( 850.5, 325.5)Band 1 Block=1701x1<br>
> Type=Int16,<br>
> ColorInterp=Undefined Metadata:<br>
> FAPAR:CLASS=Data<br>
> FAPAR:PRODUCT=product<br>
> FAPAR:PRODUCT_ID=0<br>
> FAPAR:N_COLS=1701<br>
> FAPAR:N_LINES=651<br>
> FAPAR:NB_BYTES=2<br>
> FAPAR:SCALING_FACTOR=10000<br>
> FAPAR:OFFSET=0<br>
> FAPAR:MISSING_VALUE=-10<br>
> FAPAR:UNITS=N/A<br>
> FAPAR:CAL_SLOPE=1<br>
> FAPAR:CAL_OFFSET=0<br>
><br>
> If I do :gdalwarp<br>
> HDF5:"/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000"://FAPAR<br>
> /mnt/GIS/testing.tif<br>
> I get:ERROR 1: Unable to compute a<br>
> transformation between pixel/lineand<br>
> georeferenced coordinates for<br>
> HDF5:/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000://FAPAR.<br>
<br>
*****<br>
> There is no affine transformation and no GCPs.<br>
</div></div>*****<br>
<br>
the above tells you that it is in lat/lon WGS84, but the image itself<br>
is not georeferenced.<br>
<br>
are there possibly additional data bands named "lat" and "long" or so?<br>
do you know what the north/south/east/west & resolution should be?<br>
<br>
<br>
maybe there is another HDF5 viewer (ISTR a java one from NASA?) which<br>
will export the info if it is a gdal or libhdf problem which is stopping<br>
you from seeing that information (if it is there at all)<br>
<div class="im"><br>
<br>
<br>
> Is this helpful to solve this<br>
> question?At MODIS Wiki (<a href="http://grass.osgeo.org/wiki/MODIS" target="_blank">http://grass.osgeo.org/wiki/MODIS</a>)<br>
> Gdal Warp is not used. Instead, it's create a<br>
> intermediate file (TIF) using gdal_translate and only then<br>
> is r.in.gdal is used.<br>
<br>
</div>gdalwarp is used if there is no north/south/east/west info beyond the<br>
number of image pixels, but there are a series of GCPs listed.<br>
<font color="#888888"><br>
<br>
<br>
Hamish<br>
<br>
<br>
<br>
<br>
</font></blockquote></div><br></div>