Suponho que essa seja uma pergunta básica, mas não consigo encontrar ou reconhecer a solução.
Este site retorna
Point:
X: -11705274.6374
Y: 4826473.6922
ao pesquisar com o primeiro valor da chave 000090 como exemplo. Eu acho que essa é uma referência espacial e eu meio que entendo o que é isso.
Estou procurando instruções ou exemplos de como converter isso para Latitude e longitude usando Python.
python
coordinate-system
convert
Vincent
fonte
fonte
Por padrão, o site ao qual você vinculou usa o Sistema de Referência Espacial EPSG 3857 (WGS84 Web Mercator). Encontrei esta informação aqui .
Você pode especificar outro Sistema de Referência Espacial inserindo o EPSG desejado no formulário em
Spatial Reference
ou pode converter as coordenadas retornadas com Python.Por exemplo, você pode usar as ligações GDAL Python para converter esse ponto do sistema de coordenadas projetadas (EPSG 3857) em um sistema de coordenadas geográficas (EPSG 4326).
Isso retorna para o seu ponto as coordenadas de
-105.150271116 39.7278572773
.fonte
O afalciano tem a resposta certa, mas queria incluir um uso variante do pyproj.
Ele faz necessário que você sabe a corda proj4 e é um pouco mais rápido.
fonte
p = pyproj.Proj(init='epsg:3857')
e o resultado é o mesmo.A saída não é um sistema de referência espacial / coordenada , é um par de coordenadas. Você precisa saber qual é a referência espacial para reprojetar as coordenadas.
No entanto, isso não é necessário neste caso. Basta passar uma referência espacial de saída apropriada para o serviço e ele retornará as coordenadas em Lon / Lat.
Aqui está a página com coordenadas de saída no formato Lon / Lat usando o sistema de referência geográfica geográfica WGS-84 ( EPSG 4326 ).
fonte
Tentei o código sugerido por Marcel Wilson e é mais rápido:
fonte
Encontrei este post ao procurar maneiras de fazer isso no QGIS. Conforme descrito aqui , o método usado tem a seguinte aparência:
fonte
v3.x
você precisará usarxform = QgsCoordinateTransform(crsSrc, crsDest, QgsProject.instance())
Observe que a
transform
função de acceptpyproj
tambémarrays
, o que é bastante útil quando se trata de quadros de dadosfonte
fonte