[gdal-dev] Truncated Coordinate System WKT post GDAL-1.9.2

[Even Rouault]

Le lundi 30 juin 2014 19:30:11, Wood, Alexander a écrit :
> Hello,
> 
> I have recently upgraded from GDAL 1.9.2 to GDAL 1.11.0.  With this
> upgrade, I've experienced difficulties in extracting the projection
> definition from the swisstopo dataset
> (http://www.swisstopo.admin.ch/internet/swisstopo/en/home/products/height/
> swissALTI3D/swissALTI3D.html).
> 
> In GDAL-1.9.2, running gdalinfo on the lv95 GeoTIFF, 2m data returns the
> following Coordinate System: PROJCS["Hotine Oblique Mercator Azimuth
> Center",
>     GEOGCS["CH1903+ (L+T1997)",
>         DATUM["CH1903+ (L+T1997)",
>             SPHEROID["Bessel",6377397.155,299.1528131060786],
>             TOWGS84[674.374,15.056,405.346,0,0,0,0]],
>         PRIMEM["Greenwich",0],
>         UNIT["degree",0.0174532925199433]],
>     PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],
>     PARAMETER["latitude_of_center",46.95240555555557],
>     PARAMETER["longitude_of_center",7.439583333333334],
>     PARAMETER["azimuth",90.00000000000003],
>     PARAMETER["scale_factor",1],
>     PARAMETER["false_easting",2600000],
>     PARAMETER["false_northing",1200000],
>     UNIT["meters",1]]
> 
> However in GDAL-1.11.0, gdalinfo on the same file returns the following
> (truncated) WKT for the Coordinate System: PROJCS["Hotine Oblique Mercator
> Azimuth Center",
>     GEOGCS["CH1903+ (L+T1997)",
>         DATUM["CH1903+ (L+T1997)",
>             SPHEROID["Bessel",6377397.155,299.1528131060786],
>             TOWGS84[674.374,15.056,405.346,0,0,0,0]],
>         PRIMEM["Greenwich",0],
>         UNIT["degree",0.0174532925199433]],
>     UNIT["meters",1]]
> 
> Notice that PROJECTION and related PARAMETER fields are missing in the WKT.
>  I've run gdalinfo from GDAL-1.10.0 and 1.10.1 with results identical to
> GDAL-1.11.0.  I suspect that the issue lies with the fact that the
> coordinate system is defined using ESRI WKT.
> 
> Is it possible in 1.11.0 to obtain a complete Projection WKT from
> GDALGetProjectionRef() that would include the
> PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"] in order to construct
> a OGRSpatialReference from the ESRI WKT?

Alex

It was actually a regression introduced in 1.10. Now fixed in 1.10, 1.11 
branches and trunk per ticket http://trac.osgeo.org/gdal/ticket/5551

Best regards,

Even

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html