[GRASS-user] HDF-5 parse

[Markus Neteler]

2010/2/17 Pedro Roma <pedroroma1982 at gmail.com>:
> Hello all
>
> Regarding this question, I’m using a MSG SEVIRI product called FAPAR. And
> I’m a bit lost on import data to GRASS and rectify.
>
>
>
> By making gldainfo I get:
>
> Driver: HDF5/Hierarchical Data Format Release 5
>
> Files: /mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000
>
> Size is 512, 512
>
> Coordinate System is `'
>
> Subdatasets:
>
>
> SUBDATASET_1_NAME=HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR
>
>   SUBDATASET_1_DESC=[651x1701] //FAPAR (8-bit integer)
>
>
> SUBDATASET_2_NAME=HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR_QF
>
>   SUBDATASET_2_DESC=[651x1701] //FAPAR_QF (8-bit character)
>
>
> SUBDATASET_3_NAME=HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR_err
>
>   SUBDATASET_3_DESC=[651x1701] //FAPAR_err (8-bit integer)
>
> Corner Coordinates:
>
> Upper Left  (    0.0,    0.0)
>
> Lower Left  (    0.0,  512.0)
>
> Upper Right (  512.0,    0.0)
>
> Lower Right (  512.0,  512.0)
>
> Center      (  256.0,  256.0)
>
>
>
> I get a list of all subsets. Apparentely there is no coordinate system
> right?
>
>
>
> But if I make
>
> gdalinfo HDF5:"/mnt/GIS/DATA/HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR
> (one of the SUBSETS) I get:
>
> Driver: HDF5Image/HDF5 Dataset
>
> Files: none associated
>
> Size is 1701, 651
>
> Coordinate 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
>
>
>
> This means that it’s georreferenced right?

Not necessarily. It is also possible that GCPs are stored and gdalwarp
needs to be
used first to obtain a really geocoded file.


> But if I do
>
> r.in.gdal -o --quiet input=HDF5:"/mnt/GIS/ DATA
> /HDF5_DATA_MSG_FAPAR_Euro_200702220000"://FAPAR output=out00
>
> I get the error:
>
> Over-riding projection check
>
> G_set_window(): Illegal latitude for NorthOver-riding projection check
>
> G_set_window(): Illegal latitude for North
>
>
>
> So it doesn’t’ work. I mean shouldn’t I be using gdalwarp in order to import
> and rectify image?

Say, use gdalwarp to preprocess the data. Please try with one channel and
post what gdalinfo reports on the resulting file.

Markus

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