[GRASS-user] Incorrect re-projection from EPSG:4326 to 3857

[Zhiyu (Drew) Li]

Hi there,



I was trying to re-project a EPSG:4326 coordinate into EPSG:3857 using
GRASS, but I found GDAL gave me a different result.

GRASS (7.2.2):

Re-project this " -111.648251 40.246882" lon-lat pair into a existing
location with EPSG:3857 projection

m.proj -i input=in4326.txt

-12428626.45 4874271.67



GDAL (2.2):

gdaltransform -s_srs EPSG:4326 -t_srs EPSG:3857

-111.648251 40.246882

-12428626.4492796 4901883.62422139 0



Results on X axis are almost the same, but the offset along Y is pretty big.

Any suggestion would be appreciated.



Thanks & Happy New Year!


Drew

Brigham Young University

Provo, UT USA
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20180103/bb8b9ed0/attachment.html>