[gdal-dev] GDAL fails to read projection with: Unsupported conversion method: Polar_Stereographic

[Even Rouault]

Louis-Phillippe,

this sounds like a bug. Please file a ticket at 
https://github.com/OSGeo/gdal/issues/new with all the below details

Even

Le 28/07/2021 à 17:08, Rousseau Lambert,Louis-Philippe (ECCC) a écrit :
>
> Hi,
>
> I'm trying to read a projection in proj4 string from a NetCDF and gdal 
> / proj fails to recognize the /*conversion method: Polar_Stereographic*/.
>
> The issue appears with any files in: 
> https://dd.weather.gc.ca/model_riops/netcdf/forecast/polar_stereographic/2d/06/003/
>
> example gdalinfo request:
>
> /$ gdalinfo -proj4 input.nc//
> //Driver: netCDF/Network Common Data Format//
> //...//
> //*ERROR 1: PROJ: proj_as_proj_string: Unsupported conversion method: 
> Polar_Stereographic*//*
> *//*PROJ.4 string is:*//*
> *//*''*//
> //Origin = (-2500.000000000000000,8047500.000000000000000)//
> //Pixel Size = (5000.000000000000000,-5000.000000000000000)//
> //...//
> //  SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 
> 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AXIS["Latitude",NORTH],AXIS["Longitude",EAST],AUTHORITY["EPSG","4326"]]//
> //  X_BAND=1//
> //  X_DATASET=NETCDF:"dd_riops.nc":longitude//
> //  Y_BAND=1//
> //  Y_DATASET=NETCDF:"dd_riops.nc":latitude//
> //Corner Coordinates://
> //Upper Left  (   -2500.000, 8047500.000)//
> //Lower Left  (   -2500.000,   -2500.000)//
> //Upper Right ( 8847500.000, 8047500.000)//
> //Lower Right ( 8847500.000,   -2500.000)//
> //Center      ( 4422500.000, 4022500.000)//
> //.../
>
> The gdal version I'm getting the error with: GDAL 3.1.3, and PROJ 
> version: Rel. 7.2.0.
>
> I know this command previsously worked fine, using gdal GDAL 2.4.4 and 
> proj Rel. 4.8.0 and I was able to get the PROJ4 string: '+proj=stere 
> +lat_0=90 +lat_ts=59.99993768695209 +lon_0=-100 +k=0.93301243 
> +x_0=4245000 +y_0=5295000 +a=6371229 +b=6371229 +units=m +no_defs '
>
> Also, if I try to reproject the file in EPSG:4326 with a gdalwarp 
> (GDAL 3.1.3) command, it fails, because it does not read the source 
> srs, but if I specify it using the -s_srs flag it works fine, example:
>
> |gdalwarp -of NetCDF -s_srs "+proj=stere +lat_0=90 
> +lat_ts=59.99993768695209 +lon_0=-100 +x_0=4245000 +y_0=5295000 
> +R=6371229 +units=m +no_defs" -t_srs EPSG:4326 input.nc output.nc|
>
> Is this a bug in gdal / proj, or a parameter in the projection 
> definition that is not supported anymore ?
>
> Thanks
>
> LP
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/gdal-dev

-- 
http://www.spatialys.com
My software is free, but my time generally not.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20210728/438e1ca4/attachment-0001.html>