netCDF: files show up in wrong place in QGIS

[rduivenvoorde]

Loading netcdf files from our national weather institute KNMI, for example this one:

​ftp://data.knmi.nl/download/KNW-NetCDF-3D/1.0/noversion/2013/11/14/

ends up in QGIS in the totally wrong place.

Gdalinfo ends with: Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 512.0) Upper Right ( 512.0, 0.0) Lower Right ( 512.0, 512.0) Center ( 256.0, 256.0)

QGIS debug info shows me:

Warning 6: NetCDF driver import of LCC-1SP with standard_parallel1 != latitude_of_projection_origin (which forces a computation of scale_factor) is experimental (bug #3324)

Using Panoply (java program from nasa ​https://www.giss.nasa.gov/tools/panoply/download/) the data is ok (see screenshot), and one of the arrays contains:

int Lambert_Conformal;

:grid_mapping_name = "lambert_conformal_conic"; :standard_parallel = "54.0 54.0"; :longitude_of_central_meridian = "2.0"; :latitude_of_projection_origin = "60"; :proj = "+proj=lcc +a=6371220.0 +b=6371220.0 +lon_0=2.0 +lat_2=54.0 +lat_1=54.0 +lat_0=60";

I'm totally not enough into this knowledge to understand ​https://trac.osgeo.org/gdal/ticket/3324 but hope it is possible to make this work in QGIS :-)

[rduivenvoorde]

QGIS showing grid on wrong place

[rduivenvoorde]

Panoply showing the right place

[Even Rouault]

In 40423:

netCDF: be more tolerant on the formatting of standard parallel (space separated instead of {x,y,...} syntax), and accept up to 2/1000 error on spacing to consider a regular grid, to be able to read files provided by the national weather institute of Netherlands (KNMI) (fixes #7086)

[Even Rouault]

In 40424:

netCDF: be more tolerant on the formatting of standard parallel (space separated instead of {x,y,...} syntax), and accept up to 2/1000 error on spacing to consider a regular grid, to be able to read files provided by the national weather institute of Netherlands (KNMI) (fixes #7086)