[GRASS-user] HDF-5 parse

[Hamish]

Pedro Roma wrote:
> By usinggdalinfo
> HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR
> I get:Driver:
> HDF5Image/HDF5 Dataset
> Files: none associatedSize is 1701,
> 651Coordinate System
> is: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"]],  
>  AXIS["Lat",NORTH],
>     AXIS["Long",EAST],  
>  AUTHORITY["EPSG","4326"]]Corner
> Coordinates:Upper Left  (    0.0,  
>  0.0)Lower Left  (    0.0,
>  651.0)Upper Right ( 1701.0,    0.0)
> Lower Right ( 1701.0,  651.0)Center    
>  (  850.5,  325.5)Band 1 Block=1701x1
> Type=Int16,
> ColorInterp=Undefined  Metadata:  
>  FAPAR:CLASS=Data  
>  FAPAR:PRODUCT=product
>     FAPAR:PRODUCT_ID=0   
>  FAPAR:N_COLS=1701   
>  FAPAR:N_LINES=651   
>  FAPAR:NB_BYTES=2   
>  FAPAR:SCALING_FACTOR=10000   
>  FAPAR:OFFSET=0   
>  FAPAR:MISSING_VALUE=-10 
>     FAPAR:UNITS=N/A  
>  FAPAR:CAL_SLOPE=1   
>  FAPAR:CAL_OFFSET=0 
>
> If I do :gdalwarp
> HDF5:"/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000"://FAPAR
>  /mnt/GIS/testing.tif
> I get:ERROR 1: Unable to compute a
> transformation between pixel/lineand
> georeferenced coordinates for
> HDF5:/mnt/GIS/LSASAF/HDF5_LSASAF_MSG_FAPAR_Euro_200702220000://FAPAR.

*****
> There is no affine transformation and no GCPs.
*****

the above tells you that it is in lat/lon WGS84, but the image itself
is not georeferenced.

are there possibly additional data bands named "lat" and "long" or so?
do you know what the north/south/east/west & resolution should be?


maybe there is another HDF5 viewer (ISTR a java one from NASA?) which
will export the info if it is a gdal or libhdf problem which is stopping
you from seeing that information (if it is there at all)



> Is this helpful to solve this
> question?At MODIS Wiki (http://grass.osgeo.org/wiki/MODIS)
> 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.

gdalwarp is used if there is no north/south/east/west info beyond the
number of image pixels, but there are a series of GCPs listed.



Hamish