<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
Even,</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
Thank you for your prompt and detailed reply. You've given us a lot to digest, but 2 things seem to be apparent: 1) use EPSG:XXXX notation to specify proj parameters, and 2) be glad we are not attempting any survey-precision applications with those 2 map projections.
In fact, the data sets we analyze (which are distributed in those projections) are Arctic-wide, coarse pixel-resolution rasters. Notwithstanding a desire to be as spatially accurate as possible, perhaps the imprecision of transforming coordinates between these
2 projections is inconsequential with respect to the broad scales across which we are quantifying spatial information. </div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
Thanks again.</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
Be well.</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
David</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Arial, Helvetica, sans-serif; font-size: 10pt; color: rgb(0, 0, 0);">
<br>
</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Even Rouault <even.rouault@spatialys.com><br>
<b>Sent:</b> Wednesday, April 22, 2020 11:24<br>
<b>To:</b> Douglas, David C <ddouglas@usgs.gov><br>
<b>Cc:</b> proj@lists.osgeo.org <proj@lists.osgeo.org>; Fischbach, Anthony S <afischbach@usgs.gov>; r.hijmans@gmail.com <r.hijmans@gmail.com>; Durner, George M <gdurner@usgs.gov><br>
<b>Subject:</b> Re: [EXTERNAL] Re: [PROJ] package 'geosphere' -- transforming spherical-based geographic coordinates to WGS84</font>
<div> </div>
</div>
<style type="text/css">
<!--
p, li
{white-space:pre-wrap}
-->
</style>
<div style="font-family:'monospace'; font-size:9pt; font-weight:400; font-style:normal">
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
David,</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
That could be the matter of days of discussions. I'll try to make it short :-)</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
- the +proj= syntax you mention for EPSG definitions has been deprecated since PROJ 6, since it was a rather imperfect way of capturing datum expressions. With GDAL 3 & PROJ >= 6, we encourage strongly to use directly EPSG:XXXX syntax instead of deprecated
PROJ.4 strings</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
- ... but for the code you mention, the datum behind EPSG:3408 is datum EPSG:6053 "Not specified (based on International 1924 Authalic Sphere)" and the one for EPSG:3411 is EPSG:6054 "Not specified (based on Hughes 1980 ellipsoid)". So even with PROJ 6, there's
no hope to get any significant datum transformation</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
- EPSG:3408 has been recently deprecated by EPSG in favor of EPSG:6931, and EPSG:3411 in favor of EPSG:3413. Both use the WGS84 datum, and are thus not equivalent to their non-deprecated code.</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
- With PROJ.4 syntax, there's no toXYZ qualifier, only towgs84. WGS84 in PROJ.4 era was used as the universal pivot when transforming between 2 datums (sometimes it was abused a bit to transform directly between a CRS A and B when A or B was sufficiently close
to GRS80/WGS84 for the approximation to be valid). Starting with PROJ 6, one can do direct transformations, but with your rather unspecified datums, I've no idea what would be the right way of transforming them...</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Even</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
-- </p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
Spatialys - Geospatial professional services</p>
<p style="margin-top:0px; margin-bottom:0px; margin-left:0px; margin-right:0px; text-indent:0px">
http://www.spatialys.com</p>
</div>
</body>
</html>