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

[Pierre Leman]

Hello,


I'm trying to read a WAFS GRIB2 file using gdal in order to produce Geotiff
image displayable in OsgEarth (file WAFS_blended_2016090700f06.grib2 from
http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/gfs.2016090700/).
When i do, i get an error.

The gdalinfo command give the following informations :
gdalinfo --debug on WAFS_blended_2016090700f06.grib2
GRIB: This was only designed for templates 0, 1, 2, 5, 8, 9, 10, 11, 12,
20, 30

ERROR 4: WAFS_blended_2016090700f06.grib2 is a grib file, but no raster
dataset was successfully identified.
gdalinfo failed - unable to open 'WAFS_blended_2016090700f06.grib2'.


A bug seems to be already opened for this error :
https://trac.osgeo.org/gdal/ticket/5768


Using the wgrib2 tool, i can extract the 1-3 band that are accepted by gdal
resulting in a success for these bands :
wgrib2 ..\WAFS_blended_2016090700f06.grib2
1:0:d=2016090700:ICAHT:cumulonimbus base:6 hour fcst:
2:68039:d=2016090700:ICAHT:cumulonimbus top:6 hour fcst:
3:146518:d=2016090700:CBHE:entire atmosphere:6 hour fcst:
4:167577:d=2016090700:CTP:700 mb:6 hour fcst:spatial ave
5:204299:d=2016090700:CTP:700 mb:6 hour fcst:spatial max
6:251461:d=2016090700:CTP:600 mb:6 hour fcst:spatial ave
7:288183:d=2016090700:CTP:600 mb:6 hour fcst:spatial max
8:340565:d=2016090700:CTP:500 mb:6 hour fcst:spatial ave
9:366847:d=2016090700:CTP:500 mb:6 hour fcst:spatial max
10:403569:d=2016090700:CTP:400 mb:6 hour fcst:spatial ave
11:424631:d=2016090700:CTP:400 mb:6 hour fcst:spatial max
12:461353:d=2016090700:CTP:300 mb:6 hour fcst:spatial ave
13:482415:d=2016090700:CTP:300 mb:6 hour fcst:spatial max
14:513917:d=2016090700:CAT:400 mb:6 hour fcst:spatial ave
15:555859:d=2016090700:CAT:400 mb:6 hour fcst:spatial max
16:603021:d=2016090700:CAT:350 mb:6 hour fcst:spatial ave
17:644963:d=2016090700:CAT:350 mb:6 hour fcst:spatial max
18:692125:d=2016090700:CAT:300 mb:6 hour fcst:spatial ave
19:734067:d=2016090700:CAT:300 mb:6 hour fcst:spatial max
20:786449:d=2016090700:CAT:250 mb:6 hour fcst:spatial ave
21:828391:d=2016090700:CAT:250 mb:6 hour fcst:spatial max
22:875553:d=2016090700:CAT:200 mb:6 hour fcst:spatial ave
23:912275:d=2016090700:CAT:200 mb:6 hour fcst:spatial max
24:959437:d=2016090700:CAT:150 mb:6 hour fcst:spatial ave
25:996159:d=2016090700:CAT:150 mb:6 hour fcst:spatial max
26:1043321:d=2016090700:ICIP:800 mb:6 hour fcst:spatial ave
27:1080043:d=2016090700:ICIP:800 mb:6 hour fcst:spatial max
28:1116765:d=2016090700:ICIP:700 mb:6 hour fcst:spatial ave
29:1153487:d=2016090700:ICIP:700 mb:6 hour fcst:spatial max
30:1190209:d=2016090700:ICIP:600 mb:6 hour fcst:spatial ave
31:1226931:d=2016090700:ICIP:600 mb:6 hour fcst:spatial max
32:1263653:d=2016090700:ICIP:500 mb:6 hour fcst:spatial ave
33:1300375:d=2016090700:ICIP:500 mb:6 hour fcst:spatial max
34:1337097:d=2016090700:ICIP:400 mb:6 hour fcst:spatial ave
35:1373819:d=2016090700:ICIP:400 mb:6 hour fcst:spatial max
36:1410541:d=2016090700:ICIP:300 mb:6 hour fcst:spatial ave
37:1426383:d=2016090700:ICIP:300 mb:6 hour fcst:spatial max

wgrib2 ../WAFS_blended_2016090700f06.grib2 -for 1:3 -grib
../grib2extractblock.grib2


Is there a workaround for the 4-37 bands ?
Thanks for any help you can give
Pierre
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20160914/5192a3fe/attachment-0001.html>