<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hello,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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 (<a href="http://www.swisstopo.admin.ch/internet/swisstopo/en/home/products/height/swissALTI3D/swissALTI3D.html">http://www.swisstopo.admin.ch/internet/swisstopo/en/home/products/height/swissALTI3D/swissALTI3D.html</a>).<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">In GDAL-1.9.2, running gdalinfo on the lv95 GeoTIFF, 2m data returns the following Coordinate System:<o:p></o:p></p>
<p class="MsoNormal">PROJCS["Hotine Oblique Mercator Azimuth Center",<o:p></o:p></p>
<p class="MsoNormal"> GEOGCS["CH1903+ (L+T1997)",<o:p></o:p></p>
<p class="MsoNormal"> DATUM["CH1903+ (L+T1997)",<o:p></o:p></p>
<p class="MsoNormal"> SPHEROID["Bessel",6377397.155,299.1528131060786],<o:p></o:p></p>
<p class="MsoNormal"> TOWGS84[674.374,15.056,405.346,0,0,0,0]],<o:p></o:p></p>
<p class="MsoNormal"> PRIMEM["Greenwich",0],<o:p></o:p></p>
<p class="MsoNormal"> UNIT["degree",0.0174532925199433]],<o:p></o:p></p>
<p class="MsoNormal"> PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["latitude_of_center",46.95240555555557],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["longitude_of_center",7.439583333333334],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["azimuth",90.00000000000003],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["scale_factor",1],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["false_easting",2600000],<o:p></o:p></p>
<p class="MsoNormal"> PARAMETER["false_northing",1200000],<o:p></o:p></p>
<p class="MsoNormal"> UNIT["meters",1]]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">However in GDAL-1.11.0, gdalinfo on the same file returns the following (truncated) WKT for the Coordinate System:<o:p></o:p></p>
<p class="MsoNormal">PROJCS["Hotine Oblique Mercator Azimuth Center",<o:p></o:p></p>
<p class="MsoNormal"> GEOGCS["CH1903+ (L+T1997)",<o:p></o:p></p>
<p class="MsoNormal"> DATUM["CH1903+ (L+T1997)",<o:p></o:p></p>
<p class="MsoNormal"> SPHEROID["Bessel",6377397.155,299.1528131060786],<o:p></o:p></p>
<p class="MsoNormal"> TOWGS84[674.374,15.056,405.346,0,0,0,0]],<o:p></o:p></p>
<p class="MsoNormal"> PRIMEM["Greenwich",0],<o:p></o:p></p>
<p class="MsoNormal"> UNIT["degree",0.0174532925199433]],<o:p></o:p></p>
<p class="MsoNormal"> UNIT["meters",1]]<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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. <o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">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?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks in advance,<o:p></o:p></p>
<p class="MsoNormal">Alex<o:p></o:p></p>
</div>
</body>
</html>