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 _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
