<div class="gmail_quote"><div><div class="h5"><div class="gmail_quote"><br><div class="gmail_quote"><div>Hey</div><div> </div><blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<div>
<div><div>
</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><div><br><div>Say, use gdalwarp to preprocess the data. Please try with one channel and</div></div><div>
post what gdalinfo reports on the resulting file.</div></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><div>gdalinfo HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR</div><div><br></div></div></div><div>I get:</div><div><div><div></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["WGS 84",</div><div> DATUM["WGS_1984",</div><div> SPHEROID["WGS 84",6378137,298.257223563,</div>
<div> AUTHORITY["EPSG","7030"]],</div><div> TOWGS84[0,0,0,0,0,0,0],</div><div> AUTHORITY["EPSG","6326"]],</div><div> PRIMEM["Greenwich",0,</div>
<div> AUTHORITY["EPSG","8901"]],</div><div> UNIT["degree",0.0174532925199433,</div><div> AUTHORITY["EPSG","9108"]],</div><div> AXIS["Lat",NORTH],</div>
<div> AXIS["Long",EAST],</div><div> AUTHORITY["EPSG","4326"]]</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></div><div>If I do :</div><div><div>gdalwarp HDF5:"/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000"://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" target="_blank">http://grass.osgeo.org/wiki/MODIS</a>) Gdal Warp is not used. Instead, it'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><div></div><div><div><br></div><div> </div><blockquote class="gmail_quote" style="border-left:1px solid rgb(204, 204, 204);margin:0pt 0pt 0pt 0.8ex;padding-left:1ex">
<font color="#888888"><br>
Markus<br>
</font><div><br>
> Because there is not much inofmration regarding how to do<br>
> this. I mean, it seems that if I do r.in.gdal I loose all GCP and<br>
> coordinates system information...<br>
><br>
> Sorry for these lammer questions but I think I'm lost<br>
><br>
> Thanks<br>
><br>
> Pedro<br>
><br>
> 2010/2/17 Pedro Roma <<a href="mailto:pedroroma1982@gmail.com" target="_blank">pedroroma1982@gmail.com</a>><br>
>><br>
>> Thank you all.<br>
>> Regarding Antonio's suggestion I will see if I can do that. It seems a<br>
>> good Idea. Is it possible to have a button to browse files?<br>
>> Thanks<br>
>> Best regards<br>
>> Pedro<br>
>><br>
>> 2010/2/17 António Rocha <<a href="mailto:antonio.rocha@deimos.com.pt" target="_blank">antonio.rocha@deimos.com.pt</a>><br>
>>><br>
>>> Hello Pedro<br>
>>><br>
>>> My suggestion, in case you have a consistent dataset, would be to create<br>
>>> your own script to read import HDF5 files. My question regarding this is it<br>
>>> possible to have a "Browse files" button in a GRASS scripts?<br>
>>><br>
>>> Antonio<br>
>>><br>
>>><br>
>>> __________ Information from ESET NOD32 Antivirus, version of virus<br>
>>> signature database 4872 (20100216) __________<br>
>>><br>
>>> The message was checked by ESET NOD32 Antivirus.<br>
>>><br>
>>> <a href="http://www.eset.com" target="_blank">http://www.eset.com</a><br>
>>><br>
>>><br>
>><br>
><br>
><br>
</div><div><div></div><div>> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
><br>
</div></div></blockquote></div></div></div><br>
</div><br>
</div></div></div><br>