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 <pedroroma1...@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.ro...@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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user