[gdal-dev] Problem reading WAFS GRIB2 file with gdal

[Pierre Leman]

Hi Even,


Thanks for handling this problem so fast it's perfect.
Any information about the GDAL 2.2 release date ?

Pierre

2016-09-14 16:16 GMT+02:00 Even Rouault <even.rouault at spatialys.com>:

> Pierre,
>
> > A bug seems to be already opened for this error :
> > https://trac.osgeo.org/gdal/ticket/5768
>
> I had a look and just added support for that kind of products in GDAL
> trunk.
> The issue was that Template 4.15 (
> http://www.nco.ncep.noaa.gov/pmb/docs/grib2/grib2_temp4-15.shtml ) used
> for
> bands 4-37 wasn't supported before.
>
> I've also refreshed the GRIB tables used to get parameter names from latest
> degrib v2.0.3
>
> Now you get:
>
> $ gdalinfo WAFS_blended_2016090700f06.grib2
> Driver: GRIB/GRIdded Binary (.grb)
> Files: WAFS_blended_2016090700f06.grib2
> Size is 288, 145
> Coordinate System is:
> GEOGCS["Coordinate System imported from GRIB file",
>     DATUM["unknown",
>         SPHEROID["Sphere",6371229,0]],
>     PRIMEM["Greenwich",0],
>     UNIT["degree",0.0174532925199433]]
> Origin = (-0.625000000000000,90.625000000000000)
> Pixel Size = (1.250000000000000,-1.250000000000000)
> Corner Coordinates:
> Upper Left  (  -0.6250000,  90.6250000) (  0d37'30.00"W, 90d37'30.00"N)
> Lower Left  (  -0.6250000, -90.6250000) (  0d37'30.00"W, 90d37'30.00"S)
> Upper Right (     359.375,      90.625) (359d22'30.00"E, 90d37'30.00"N)
> Lower Right (     359.375,     -90.625) (359d22'30.00"E, 90d37'30.00"S)
> Center      ( 179.3750000,   0.0000000) (179d22'30.00"E,  0d 0' 0.01"N)
> Band 1 Block=288x1 Type=Float64, ColorInterp=Undefined
>   Description = 0[-] RESERVED(11) (Reserved)
>   Metadata:
>     GRIB_COMMENT=ICAO Standard Atmosphere Reference Height [m]
>     GRIB_ELEMENT=ICAHT
>     GRIB_FORECAST_SECONDS=21600 sec
>     GRIB_PDS_PDTN=0
>     GRIB_PDS_TEMPLATE_NUMBERS=3 3 2 0 96 0 0 0 1 0 0 0 6 11 0 0 0 0 0 255
> 0 0
> 0 0 0
>     GRIB_REF_TIME=  1473206400 sec UTC
>     GRIB_SHORT_NAME=0-RESERVED(11)
>     GRIB_UNIT=[m]
>     GRIB_VALID_TIME=  1473228000 sec UTC
> [...snip...]
> Band 4 Block=288x1 Type=Float64, ColorInterp=Undefined
>   Description = 70000[Pa] ISBL="Isobaric surface"
>   Metadata:
>     GRIB_COMMENT=In-Cloud Turbulence [%]
>     GRIB_ELEMENT=CTP
>     GRIB_FORECAST_SECONDS=21600 sec
>     GRIB_PDS_PDTN=15
>     GRIB_PDS_TEMPLATE_NUMBERS=19 21 2 0 96 0 0 0 1 0 0 0 6 100 0 0 1 17 112
> 255 0 0 0 0 0 0 3 1
>     GRIB_REF_TIME=  1473206400 sec UTC
>     GRIB_SHORT_NAME=70000-ISBL
>     GRIB_UNIT=[%]
>     GRIB_VALID_TIME=  1473228000 sec UTC
> [... snip ... ]
> Band 37 Block=288x1 Type=Float64, ColorInterp=Undefined
>   Description = 30000[Pa] ISBL="Isobaric surface"
>   Metadata:
>     GRIB_COMMENT=Icing [%]
>     GRIB_ELEMENT=ICIP
>     GRIB_FORECAST_SECONDS=21600 sec
>     GRIB_PDS_PDTN=15
>     GRIB_PDS_TEMPLATE_NUMBERS=19 20 2 0 96 0 0 0 1 0 0 0 6 100 0 0 0 117 48
> 255 0 0 0 0 0 2 3 1
>     GRIB_REF_TIME=  1473206400 sec UTC
>     GRIB_SHORT_NAME=30000-ISBL
>     GRIB_UNIT=[%]
>     GRIB_VALID_TIME=  1473228000 sec UTC
>
>
> Even
>
> --
> Spatialys - Geospatial professional services
> http://www.spatialys.com
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20160914/43e627d7/attachment.html>