Forgot to copy the list. ---------- Forwarded message ---------- From: Scott Sinclair <scott.sinclair...@gmail.com> Date: 2009/6/30 Subject: Re: [gdal-dev] can't open grib2 files To: Alberto Pettazzi <alberto.petta...@meteogalicia.es>
>2009/6/29 Alberto Pettazzi <alberto.petta...@meteogalicia.es>: > I just have executed the command you suggested me, > > gdal-config --formats > > > and that's what I got: > > gxf gtiff hfa aigrid aaigrid ceos ceos2 iso8211 xpm sdts raw dted mem jdem > envisat elas fit vrt usgsdem l1b nitf bmp pcidsk airsar rs2 ilwis rmf > leveller sgi srtmhgt idrisi gsg ingr ers jaxapalsar dimap gff cosar pds adrg > coasp tsx terragen blx msgn wcs wms grib bsb jpeg2000 hdf5 hdf4 gif jpeg png > netcdf pcraster rik > > So, I guess that I have already built the jpeg2000 support. Any Suggestion? When you run gdalinfo on the input file you'll see that the GRIB file has 3 bands. The last two are quality flags at a different spatial resolution to the rainfall estimates (it seems that's what causes the GRIB driver to choke). See http://www.eumetsat.int/groups/ops/documents/document/pdf_metprod_mpe_quality_info.pdf Using the -b flag to gdal_translate works for me: $ gdalinfo --version GDAL 1.6.1, released 2009/05/11 $ gdal_translate -b 1 MPE_20090629_0900_M9_00.grb MPE_20090629_0900_M9_00.tif $ gdalinfo MPE_20090629_0900_M9_00.tif Driver: GTiff/GeoTIFF Files: MPE_20090629_0900_M9_00.tif Size is 3712, 3712 Coordinate System is: LOCAL_CS["unnamed", UNIT["metre",1, AUTHORITY["EPSG","9001"]]] Origin = (-5568748.275999999605119,5568748.275999999605119) Pixel Size = (3000.403165948275728,-3000.403165948275728) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left (-5568748.276, 5568748.276) Lower Left (-5568748.276,-5568748.276) Upper Right ( 5568748.276, 5568748.276) Lower Right ( 5568748.276,-5568748.276) Center ( 0.0000000, -0.0000000) Band 1 Block=3712x1 Type=Float64, ColorInterp=Gray Metadata: GRIB_UNIT=[kg/(m^2*s)] GRIB_COMMENT=Instantaneous rain rate [kg/(m^2*s)] GRIB_SHORT_NAME=0 undefined GRIB_REF_TIME=1246266000 sec UTC GRIB_VALID_TIME=1246266000 sec UTC GRIB_FORECAST_SECONDS=0 sec GRIB_PDS_PDTN=30 GRIB_PDS_TEMPLATE_NUMBERS=1 1 8 7 1 1 77 0 56 207 255 255 255 255 255 Cheers, Scott _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev