Last active
August 22, 2024 14:22
-
-
Save springmeyer/871897 to your computer and use it in GitHub Desktop.
convert from long/lat to google mercator (or EPSG:4326 to EPSG:900913)
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
// See https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames for more details. | |
var degrees2meters = function(lon,lat) { | |
var x = lon * 20037508.34 / 180; | |
var y = Math.log(Math.tan((90 + lat) * Math.PI / 360)) / (Math.PI / 180); | |
y = y * 20037508.34 / 180; | |
return [x, y] | |
} | |
x= -77.035974 | |
y = 38.898717 | |
console.log(degrees2meters(x,y)) | |
// should result in: -8575605.398444, 4707174.018280 |
Thank you so much! It was really helpful.
console.log(degrees2meters(-180,-90)) // [-20037508.34, -Infinity]
Can you please help with converting meters to lat long?
The C++ macro posted above works well for converting meters to lat/long degrees: #define DEGREE_TO_METER_REVERSE(Y) (atan(pow(M_E, ((Y)/111319.490778)*M_PI/180.0))*360.0/M_PI-90.0)
(Thanks @jiangfeng79!)
Should be trivial to adapt to Python.
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
@MightyTedKim
Here's my general workflow for transforming a CRS in R.
With some GIS file (shp, gpkg, etc):
With a CSV: