Zhiyu (Drew) Li wrote > 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.
have a look into a related ticket: https://trac.osgeo.org/grass/ticket/2229 in my experience EPSG:3857 isn't the best projection for analysis. ----- best regards Helmut -- Sent from: http://osgeo-org.1560.x6.nabble.com/Grass-Users-f3884509.html _______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
