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
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
