¿Cómo convertir coordenadas proyectadas a lat/lon utilizando Python?

Question

Supongo que es una pregunta básica pero no consigo encontrar o reconocer la solución.

Este sitio devuelve

Point:
X: -11705274.6374
Y: 4826473.6922

al buscar con el primer valor de la clave 000090 como ejemplo. Supongo que se trata de una referencia espacial y más o menos entiendo lo que es.

Estoy buscando instrucciones o ejemplos de cómo convertir esto en Latitud y Longitud usando Python.

Answer 1

La forma más sencilla de transformar coordenadas en Python es pyproj es decir, el Interfaz de Python para la biblioteca PROJ.4 . De hecho:

from pyproj import Proj, transform

inProj = Proj(init='epsg:3857')
outProj = Proj(init='epsg:4326')
x1,y1 = -11705274.6374,4826473.6922
x2,y2 = transform(inProj,outProj,x1,y1)
print x2,y2

devuelve -105.150271116 39.7278572773

---

EDITADO a partir del comentario de Marc:

pyproj 2.4 da un FutureWarning sobre el deprecated Proj inicialización con el init= sintaxis. La sintaxis actualizada es idéntica pero sin el init= . Así:

inProj = Proj('epsg:3857')
outProj = Proj('epsg:4326')

Answer 2

Por defecto, el sitio al que ha enlazado utiliza el Sistema de Referencia Espacial EPSG 3857 (WGS84 Web Mercator). He encontrado esta información aquí .

Puede especificar otro Sistema de Referencia Espacial introduciendo el EPSG deseado en el formulario bajo Spatial Reference o puede convertir las coordenadas devueltas con Python.

Por ejemplo, puede utilizar el Enlaces GDAL Python para convertir este punto del sistema de coordenadas proyectadas (EPSG 3857) a un sistema de coordenadas geográficas (EPSG 4326).

import ogr, osr

pointX = -11705274.6374 
pointY = 4826473.6922

# Spatial Reference System
inputEPSG = 3857
outputEPSG = 4326

# create a geometry from coordinates
point = ogr.Geometry(ogr.wkbPoint)
point.AddPoint(pointX, pointY)

# create coordinate transformation
inSpatialRef = osr.SpatialReference()
inSpatialRef.ImportFromEPSG(inputEPSG)

outSpatialRef = osr.SpatialReference()
outSpatialRef.ImportFromEPSG(outputEPSG)

coordTransform = osr.CoordinateTransformation(inSpatialRef, outSpatialRef)

# transform point
point.Transform(coordTransform)

# print point in EPSG 4326
print point.GetX(), point.GetY()

Esto devuelve para su punto las coordenadas de -105.150271116 39.7278572773 .

Answer 3

afalciano tiene la respuesta correcta pero quería incluir una variante de uso de pyproj.

Es hace requiere que conozcas la cadena proj4 y es un poco más rápido.

import pyproj
p = pyproj.Proj("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
lon, lat = p(x, y, inverse=True)
print lat, lon

Answer 4

La salida no es un sistema de referencia espacial/de coordenadas es un par de coordenadas. Es necesario saber cuál es la referencia espacial para reproyectar las coordenadas.

Sin embargo, eso no es necesario en este caso. Basta con pasar una referencia espacial de salida adecuada al servicio y éste devolverá las coordenadas en Lon/Lat.

Aquí está el página con coordenadas de salida en formato Lon/Lat utilizando el sistema de referencia espacial geográfico WGS-84 ( EPSG 4326 ).

Answer 5

He probado el código sugerido por Marcel Wilson y es más rápido:

from pyproj import Proj, transform
import time
import pyproj

# Test 1 - 0.0006158 s
start=time.time()
inProj = Proj(init='epsg:3857')
outProj = Proj(init='epsg:4326')
x1,y1 = -11705274.6374,4826473.6922
x2,y2 = transform(inProj,outProj,x1,y1)
end=time.time()
print(y2,x2)
print('%.7f' % (end-start))

# Test 2 - 0.0000517 s --- factor 11,9
start=time.time()
x,y = -11705274.6374,4826473.6922
p = pyproj.Proj("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")
lon, lat = p(x, y, inverse=True)
end=time.time()
print(lat, lon)
print('%.7f' % (end-start))
-----------------

39.72785727727918 -105.15027111593008
0.0006158
39.72785727727918 -105.15027111593008
0.0000517