[GRASS-user] HDF-5 parse

[Pedro Roma]

hello Hamish

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).

By the way, thanks about gdalwarp enlightment.

Pedro

On Thu, Feb 18, 2010 at 12:17 PM, Hamish <hamish_b at yahoo.com> wrote:

> 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
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20100218/d8847855/attachment.html