[GRASS-user] HDF-5 parse

[Pedro Roma]

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?

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