<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 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@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;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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">Hi Even,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Sorry about my lack of clarity. GetEPSGGeogCS didn’t error for me either, but it returned -1 instead of a valid EPSG code (4326) as it used to. I have generally used WKT2 when possible, but this GDAL feature seems to internally use WKT1_GDAL
which led me to ask. Does this mean that GetEPSGGeogCS will no longer return the expected result due to this side effect? I was under the impression that I could represent this CRS previously even if I could not use it for transformations, but it sounds
like that isn’t correct.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Jonathan<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in">
<p class="MsoNormal"><b>From:</b> Even Rouault <even.rouault@spatialys.com> <br>
<b>Sent:</b> Monday, December 28, 2020 11:48 AM<br>
<b>To:</b> proj@lists.osgeo.org<br>
<b>Cc:</b> Jonathan Gale <jgale@mathworks.com><br>
<b>Subject:</b> Re: [PROJ] WKT1 Vertical Perspective CRS regression<o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Jonathan<br>
<br>
This is a side effect of <a href="https://github.com/OSGeo/PROJ/pull/1971">
https://github.com/OSGeo/PROJ/pull/1971</a><br>
<br>
The past mapping of ESRI:54049 to the below WKT was incorrect. GDAL 2.4 (which <br>
is what WKT1_GDAL means) has never recognized the nsper projection. If you <br>
tried the below WKT with it, it would fail.<br>
<br>
You have to use WKT2 with that projection.<br>
<br>
I tried quickly with GDAL (master and 3.2 branch) to reproduce the issue by <br>
doing in Python<br>
<br>
from osgeo import osr<br>
sr = osr.SpatialReference()<br>
sr.SetFromUserInput('ESRI:54049')<br>
sr.AutoIdentifyEPSG() (which calls GetEPSGGeogCS())<br>
<br>
but I couldn't reproduce any error when doing so.<br>
<br>
Even<br>
<br>
> Hi PROJ list,<br>
> <br>
> I'm not sure which list to report this since my initial issue was expressed<br>
> through GDAL, but since the issue is reproducible in PROJ alone, I'll start<br>
> here. I've encountered a regression upgrading from PROJ 6.2.1 to 7.2.0<br>
> when using GDAL's GetEPSGGeogCS with the Vertical Perspective CRS expressed<br>
> by ESRI: 54089. Based on an initial investigation, I think the issue is in<br>
> WKT1:GDAL generation.<br>
> <br>
> Using PROJ 6.2.1, I can do the following:<br>
> <br>
> projinfo ESRI:54049 -o WKT1:GDAL<br>
> WKT1_GDAL:<br>
> PROJCS["World_Vertical_Perspective",<br>
> GEOGCS["WGS 84",<br>
> DATUM["WGS_1984",<br>
> SPHEROID["WGS 84",6378137,298.257223563,<br>
> AUTHORITY["EPSG","7030"]],<br>
> AUTHORITY["EPSG","6326"]],<br>
> PRIMEM["Greenwich",0],<br>
> UNIT["Degree",0.0174532925199433]],<br>
> PROJECTION["Vertical_Near_Side_Perspective"],<br>
> PARAMETER["False_Easting",0],<br>
> PARAMETER["False_Northing",0],<br>
> PARAMETER["Longitude_Of_Center",0],<br>
> PARAMETER["Latitude_Of_Center",0],<br>
> PARAMETER["Height",35800000],<br>
> UNIT["metre",1,<br>
> AUTHORITY["EPSG","9001"]],<br>
> AXIS["Easting",EAST],<br>
> AXIS["Northing",NORTH],<br>
> AUTHORITY["ESRI","54049"]]<br>
> <br>
> But in PROJ 7.2.0, I get the following:<br>
> <br>
> projinfo ESRI:54049 -o WKT1:GDAL<br>
> WKT1:GDAL string:<br>
> Error when exporting to WKT1:GDAL: Unsupported conversion method: Vertical<br>
> Perspective<br>
> <br>
> <br>
> I'm aware of release notes for GDAL that impacted Vertical Perspective (and<br>
> I haven't found any corresponding in PROJ), but I suspect the GDAL issue is<br>
> rooted in PROJ too.<br>
> <br>
> <br>
> Thanks,<br>
> Jonathan<br>
<br>
<br>
-- <br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com">http://www.spatialys.com</a><o:p></o:p></p>
</div>
</body>
</html>