<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;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;}
@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="#0563C1" vlink="#954F72" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal">Hello GDAL users.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have noticed that from GDAL 1.10 – GDAL 3.1.0 that gdaltransform is producing different results for NAD27 GCS to WGS84 GCS transformations. The differences in longitude in the output is more significant. Having tested the differences
in feet between the points output from each version I can see a difference of ~5 feet primarily due to the change in the longitude value. <o:p></o:p></p>
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal">Would anyone look at the following results and share your thoughts about why the result of this transformation is different in a newer version? Could this be related to the concept of WGS 84 realization being different in newer versions
of GDAL? I am also linking GDAL 3.4.2 with PROJ 8.1.1 at this time. These results were produced with GDAL 3.4.2, but the same results occur with GDAL 3.1.0 as well. GDAL 3.4.2 is what I am currently trying to upgrade to. The results from the C++ programming
API are similar in difference. I just used gdaltransform for the purposes of this message because it was a more convenient way of producing something that others might be able to reproduce.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>GDAL 1.10.0 results<o:p></o:p></b></p>
<p class="MsoNormal">./gdaltransform -s_srs NAD27 -t_srs WGS84<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">-117.258911 32.804191<o:p></o:p></p>
<p class="MsoNormal">-117.2597<span style="color:red">8144845</span> 32.80423<span style="color:red">93819884</span> 1.32247805595398e-07<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><b>GDAL 3.4.2 results<o:p></o:p></b></p>
<p class="MsoNormal">./gdaltransform -s_srs NAD27 -t_srs WGS84<o:p></o:p></p>
<p class="MsoNormal">-117.258911 32.804191<o:p></o:p></p>
<p class="MsoNormal">-117.2597<span style="color:red">64476223</span> 32.80423<span style="color:red">88409449</span> 0<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I have read that the support for specifying epochs related to realizations is not available until GDAL 3.8.0.
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Here are some references I’ve read so far.<o:p></o:p></p>
<p class="MsoNormal"><a href="https://ascelibrary.org/doi/10.1061/%28ASCE%29SU.1943-5428.0000389">Transforming between WGS84 Realizations | Journal of Surveying Engineering | Vol 148, No 2 (ascelibrary.org)</a><o:p></o:p></p>
<p class="MsoNormal"><a href="https://gdal.org/development/rfc/rfc81_coordinate_epoch.html">RFC 81: Support for coordinate epochs in geospatial formats — GDAL documentation</a><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks,<o:p></o:p></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Tahoma",sans-serif;color:black;mso-ligatures:none">Shawn Fox<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span style="font-size:10.0pt;font-family:"Tahoma",sans-serif;color:black;mso-ligatures:none">San Diego, CA</span><span style="font-size:8.0pt;mso-ligatures:none">
</span><span style="mso-ligatures:none"><o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>