<div class="gmail_quote"><div>Hey</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div class="h5">
</div></div>Not necessarily. It is also possible that GCPs are stored and gdalwarp<br>
needs to be<br>
used first to obtain a really geocoded file.<br>
<div class="im"><br>Say, use gdalwarp to preprocess the data. Please try with one channel and</div>
post what gdalinfo reports on the resulting file.</blockquote><div><br></div><div>Since GRASS documentation extensively preprocess to GEOTIFF, I will try that. Is it the best thing to do? </div><div>ok here it goes</div>
<div><br></div><div>By using</div><div><div>gdalinfo HDF5:&quot;/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000&quot;://FAPAR</div><div><br></div></div><div>I get:</div><div><div>Driver: HDF5Image/HDF5 Dataset</div><div>
Files: none associated</div><div>Size is 1701, 651</div><div>Coordinate System is:</div><div>GEOGCS[&quot;WGS 84&quot;,</div><div>    DATUM[&quot;WGS_1984&quot;,</div><div>        SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,</div>
<div>            AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],</div><div>        TOWGS84[0,0,0,0,0,0,0],</div><div>        AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],</div><div>    PRIMEM[&quot;Greenwich&quot;,0,</div>
<div>        AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],</div><div>    UNIT[&quot;degree&quot;,0.0174532925199433,</div><div>        AUTHORITY[&quot;EPSG&quot;,&quot;9108&quot;]],</div><div>    AXIS[&quot;Lat&quot;,NORTH],</div>
<div>    AXIS[&quot;Long&quot;,EAST],</div><div>    AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</div><div>Corner Coordinates:</div><div>Upper Left  (    0.0,    0.0)</div><div>Lower Left  (    0.0,  651.0)</div><div>Upper Right ( 1701.0,    0.0)</div>
<div>Lower Right ( 1701.0,  651.0)</div><div>Center      (  850.5,  325.5)</div><div>Band 1 Block=1701x1 Type=Int16, ColorInterp=Undefined</div><div>  Metadata:</div><div>    FAPAR:CLASS=Data</div><div>    FAPAR:PRODUCT=product</div>
<div>    FAPAR:PRODUCT_ID=0 </div><div>    FAPAR:N_COLS=1701 </div><div>    FAPAR:N_LINES=651 </div><div>    FAPAR:NB_BYTES=2 </div><div>    FAPAR:SCALING_FACTOR=10000 </div><div>    FAPAR:OFFSET=0 </div><div>    FAPAR:MISSING_VALUE=-10 </div>
<div>    FAPAR:UNITS=N/A</div><div>    FAPAR:CAL_SLOPE=1 </div><div>    FAPAR:CAL_OFFSET=0 </div><div><br></div><div>If I do :</div><div><div>gdalwarp HDF5:&quot;/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000&quot;://FAPAR  /mnt/GIS/testing.tif</div>
<div>I get:</div><div><div>ERROR 1: Unable to compute a transformation between pixel/line</div><div>and georeferenced coordinates for HDF5:/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000://FAPAR.</div><div>There is no affine transformation and no GCPs.</div>
<div><br></div></div></div><div>Is this helpful to solve this question?</div><div>At MODIS Wiki (<a href="http://grass.osgeo.org/wiki/MODIS">http://grass.osgeo.org/wiki/MODIS</a>) Gdal Warp is not used. Instead, it&#39;s create a intermediate file (TIF) using gdal_translate and only then is r.in.gdal is used.</div>
<div><br></div><div>Thank you Markus and Antonio and Hamish for your help</div></div><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><font color="#888888"><br>

Markus<br>
</font><div class="im"><br>
&gt; Because there is not much inofmration regarding how to do<br>
&gt; this. I mean, it seems that if I do r.in.gdal I loose all GCP and<br>
&gt; coordinates system information...<br>
&gt;<br>
&gt; Sorry for these lammer questions but I think I&#39;m lost<br>
&gt;<br>
&gt; Thanks<br>
&gt;<br>
&gt; Pedro<br>
&gt;<br>
&gt; 2010/2/17 Pedro Roma &lt;<a href="mailto:pedroroma1982@gmail.com">pedroroma1982@gmail.com</a>&gt;<br>
&gt;&gt;<br>
&gt;&gt; Thank you all.<br>
&gt;&gt; Regarding Antonio&#39;s suggestion I will see if I can do that. It seems a<br>
&gt;&gt; good Idea. Is it possible to have a button to browse files?<br>
&gt;&gt; Thanks<br>
&gt;&gt; Best regards<br>
&gt;&gt; Pedro<br>
&gt;&gt;<br>
&gt;&gt; 2010/2/17 António Rocha &lt;<a href="mailto:antonio.rocha@deimos.com.pt">antonio.rocha@deimos.com.pt</a>&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Hello Pedro<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; My suggestion, in case you have a consistent dataset, would be to create<br>
&gt;&gt;&gt; your own script to read import HDF5 files. My question regarding this is it<br>
&gt;&gt;&gt; possible to have a &quot;Browse files&quot; button in a GRASS scripts?<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; Antonio<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; __________ Information from ESET NOD32 Antivirus, version of virus<br>
&gt;&gt;&gt; signature database 4872 (20100216) __________<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; The message was checked by ESET NOD32 Antivirus.<br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt; <a href="http://www.eset.com" target="_blank">http://www.eset.com</a><br>
&gt;&gt;&gt;<br>
&gt;&gt;&gt;<br>
&gt;&gt;<br>
&gt;<br>
&gt;<br>
</div><div><div></div><div class="h5">&gt; _______________________________________________<br>
&gt; grass-user mailing list<br>
&gt; <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
&gt; <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
&gt;<br>
&gt;<br>
</div></div></blockquote></div><br>