[gdal-dev] NetCDF rotated grid support

[Rousseau Lambert2, Louis-Philippe (EC)]

Hi all,

I have a NetCDF file with a rotated grid that I try to manipulate with GDAL (master build GDAL 3.1.0dev-c6bff27467, released 2020/04/14) and open in QGIS (nightly build) and it looks like the rotated grid is not recognized or supported:

ERROR 1: Only OGC WKT GEOGCS and PROJCS Projections supported for writing to NetCDF.  GEOGCRS["unnamed",BASEGEOGCRS["unknown",DATUM["unknown",ELLIPSOID["unknown",6371229,0,LENGTHUNIT["metre",1,ID["EPSG",9001]]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433],ID["EPSG",8901]]],DERIVINGCONVERSION["unknown",METHOD["PROJ ob_tran o_proj=longlat"],PARAMETER["o_lon_p",0,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],PARAMETER["o_lat_p",31.758312,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],PARAMETER["lon_0",-92.402969,ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]],CS[ellipsoidal,2],AXIS["longitude",east,ORDER[1],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]],AXIS["latitude",north,ORDER[2],ANGLEUNIT["degree",0.0174532925199433,ID["EPSG",9122]]]] not supported.

I also tried converting a valid GRIB2 file with a rotated grid to a NetCDF and got the same error. Note that the GRIB2 file is correctly read and opened in QGIS and with GDAL.

Here is the proj4 definition of the GRIB2 projection: "+proj=ob_tran +o_proj=longlat +o_lon_p=-0 +o_lat_p=31.758312 +lon_0=-92.402969 +R=6371229 +no_defs +type=crs"

Is this a known issue (or expected), or should the rotated grid be supported with in NetCDF ?

Thanks

LP
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20200414/fc5b6349/attachment.html>