<div dir="ltr"><p class="gmail-MsoNormal">Hi there,</p><p class="gmail-MsoNormal"> </p>

<p class="gmail-MsoNormal">I was trying to re-project a EPSG:4326 coordinate into EPSG:3857 using GRASS, but I found GDAL gave me a different result. </p>

<p class="gmail-MsoNormal">GRASS (7.2.2):<span></span></p>

<p class="gmail-MsoNormal">Re-project this " -111.648251 40.246882" lon-lat pair
into a existing location with EPSG:3857 projection<span></span></p>

<p class="gmail-MsoNormal">m.proj -i input=in4326.txt<span></span></p>

<p class="gmail-MsoNormal">-12428626.45 4874271.67<span></span></p>

<p class="gmail-MsoNormal"><span> </span></p>

<p class="gmail-MsoNormal">GDAL (2.2):<span></span></p>

<p class="gmail-MsoNormal">gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857<span></span></p>

<p class="gmail-MsoNormal">-111.648251 40.246882<span></span></p>

<p class="gmail-MsoNormal">-12428626.4492796 4901883.62422139 0<span></span></p>

<p class="gmail-MsoNormal"><span> </span></p>

<p class="gmail-MsoNormal">Results on X axis are almost the same, but the offset along Y is pretty big.<span></span></p>

<p class="gmail-MsoNormal">Any suggestion would be appreciated.<span></span></p>

<p class="gmail-MsoNormal"><span> </span></p>

<p class="gmail-MsoNormal">Thanks & Happy New Year!<span></span></p>

<p class="gmail-MsoNormal"><br></p><p class="gmail-MsoNormal">Drew<span></span></p><p class="gmail-MsoNormal">Brigham Young University</p><p class="gmail-MsoNormal">Provo, UT USA</p></div>