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

Wood, Alexander awood at agi.com
Mon Jun 30 10:30:11 PDT 2014


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?

Thanks in advance,
Alex
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20140630/3d4e7b17/attachment-0001.html>


More information about the gdal-dev mailing list