[gdal-dev] can't handle grib file

[Even Rouault]

> So it looks like the information is actually there, but the gdal driver
> can't handle it. Any hint?

Margherita,

There's a possibility that your data file is of a type that isn't understood by 
the underlying library that GDAL uses to decode GRIB files. Or perhaps GDAL 
isn't using it correctly with that particular file. Difficult to know.

GDAL uses the g2clib library (v1.0.4), that has evolved ( 
http://www.nco.ncep.noaa.gov/pmb/docs/grib2/download/g2clib.changes ) since it 
was integrated in GDAL. So perhaps a newer version of it better supports your 
data file. (but merging a newer version to GDAL might be non trivial, since a 
lot of fixes were done in GDAL copy (to address cross-platform portability, 
etc..)).

> 
> Thanks,
> 
> Margherita
> 
> [1] http://www.ecmwf.int/products/data/software/grib_api.html
> 
> On Thu, Aug 23, 2012 at 8:53 AM, Margherita Di Leo 
<diregola at gmail.com>wrote:
> > Hi all,
> > 
> > I have a problem handling some grib files, I get a warning message and I
> > can't understand how to deal with it..
> > 
> > running gdalinfo, it says:
> > 
> > Un-handled possible ensemble section center 78 subcenter 255
> > [... several times, cut]
> > Un-handled possible ensemble section center 78 subcenter 255
> > Driver: GRIB/GRIdded Binary (.grb)
> > Files:
> > /home/leomarg/Desktop/Fire_event_test_case_Valencia/grib/dwd_grib1_ispra_
> > gmi202_ 2012063000
> > Size is 201, 161
> > Coordinate System is:
> > GEOGCS["Coordinate System imported from GRIB file",
> > 
> >     DATUM["unknown",
> >     
> >         SPHEROID["Sphere",6367470,0]],
> >     
> >     PRIMEM["Greenwich",0],
> >     UNIT["degree",0.0174532925199433]]
> > 
> > Origin = (-10.000000000000000,74.960000000000008)
> > Pixel Size = (0.000000000000000,-0.000000000000000)
> > Corner Coordinates:
> > Upper Left  ( -10.0000000,  74.9600000) ( 10d 0' 0.00"W, 74d57'36.00"N)
> > Lower Left  ( -10.0000000,  74.9600000) ( 10d 0' 0.00"W, 74d57'36.00"N)
> > Upper Right ( -10.0000000,  74.9600000) ( 10d 0' 0.00"W, 74d57'36.00"N)
> > Lower Right ( -10.0000000,  74.9600000) ( 10d 0' 0.00"W, 74d57'36.00"N)
> > Center      ( -10.0000000,  74.9600000) ( 10d 0' 0.00"W, 74d57'36.00"N)
> > Band 1 Block=201x1 Type=Float64, ColorInterp=Undefined
> > 
> >   Description = 10[m] HTGL (Specified height level above ground)
> >   
> >   Metadata:
> >     GRIB_COMMENT=u-component (zonal) of wind [m/s]
> >     GRIB_ELEMENT=U
> >     GRIB_FORECAST_SECONDS=0 sec
> >     GRIB_REF_TIME=  1341014400 sec UTC
> >     GRIB_SHORT_NAME=10-HTGL
> >     GRIB_UNIT=[m/s]
> >     GRIB_VALID_TIME=  1341014400 sec UTC
> > 
> > Band 2 Block=201x1 Type=Float64, ColorInterp=Undefined
> > 
> >   Description = 10[m] HTGL (Specified height level above ground)
> >   
> >   Metadata:
> >     GRIB_COMMENT=v-component (merdional) of wind [m/s]
> >     GRIB_ELEMENT=V
> >     GRIB_FORECAST_SECONDS=0 sec
> >     GRIB_REF_TIME=  1341014400 sec UTC
> >     GRIB_SHORT_NAME=10-HTGL
> >     GRIB_UNIT=[m/s]
> >     GRIB_VALID_TIME=  1341014400 sec UTC
> > 
> > [...]
> > 
> > Furthermore, when I try with gdal_translate to convert one of the bands
> > in tif format, the size of the pixels created is 0x0.
> > 
> > Thanks in advance for any hint.
> > 
> > --
> > Dr. Margherita Di Leo