Fwd: [gdal-dev] can't open grib2 files

[Scott Sinclair]

Forgot to copy the list.

---------- Forwarded message ----------
From: Scott Sinclair <scott.sinclair.za at gmail.com>
Date: 2009/6/30
Subject: Re: [gdal-dev] can't open grib2 files
To: Alberto Pettazzi <alberto.pettazzi at meteogalicia.es>


>2009/6/29 Alberto Pettazzi <alberto.pettazzi at 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