[Qgis-developer] reading GRIB data

[Even Rouault]

On jeudi 5 janvier 2017 16:59:35 CET Paolo Cavallini wrote:
> Il 05/01/2017 16:31, Even Rouault ha scritto:
> > This is a bit vague as a description of issues ;-) gdalinfo reports on
> > datasets would be a start. GRIB support is generally a pain due to the
> > diversity of the format.
> 
> Hi Even,
> sorry for being vague - I hoped I was just missing something obvious,
> and an howto was already available. Life is mostly more complicated :)
> So, symptoms:
> * the grib is loaded with a custom CRS:
>  * Generated CRS (+proj=longlat +a=6371200 +b=6371200 +no_defs)

That's not really an issue per se, is it? GRIB data comes with non standard CRS, but that's how 
the data is registered. Can be reprojected to WGS84 later if that makes life easier

> * it occupies a thin vertical stretch, from about 1.86884,-1418.31515 to
> about 18.142,67.489

That's more an issue. There's properly a lack of support for the projection method, or 
something wrong when reading the extent.

> * a series of warning is issued, of the type:
> Warning 1: Band 895 of GRIB dataset is 33x23, while the first band and
> dataset is 65x45.  Georeferencing of band 895 may be incorrect, and data
> access may be incomplete.

That's a defect of the driver I'm aware of. It should expose those as subdatasets to fit the 
model of a GDAL dataset.


-- 
Spatialys - Geospatial professional services
http://www.spatialys.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/qgis-developer/attachments/20170105/982acbd4/attachment.html>