Konwertowanie rzutowanych współrzędnych na długość / długość przy użyciu Pythona?

50

Zakładam, że jest to podstawowe pytanie, ale nie mogę znaleźć ani rozpoznać rozwiązania.

Ta strona powraca

Point:
X: -11705274.6374
Y: 4826473.6922

podczas wyszukiwania z pierwszą wartością klucza 000090 jako przykładem. Myślę, że to odniesienie przestrzenne i trochę rozumiem, co to jest.

Szukam instrukcji lub przykładów, jak przekonwertować to na szerokość i długość geograficzną za pomocą Pythona.

Vincent
źródło

Odpowiedzi:

100

Najprostszym sposobem transformacji współrzędnych w Pythonie jest pyproj , tj. Interfejs Pythona do biblioteki PROJ.4 . W rzeczywistości:

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

zwroty -105.150271116 39.7278572773

Antonio Falciano
źródło
4
Tak. Pyproj do końca.
sgillies,
To działa dla mnie
lenhhoxung
36

Domyślnie witryna, z którą się łączysz, korzysta z systemu odniesienia przestrzennego EPSG 3857 (WGS84 Web Mercator). Znalazłem tę informację tutaj .

Możesz albo określić inny System Odniesienia Przestrzennego, wpisując żądany EPSG w formularzu pod, Spatial Referencealbo możesz przekonwertować zwrócone współrzędne za pomocą Pythona.

Na przykład możesz użyć powiązań GDAL Python do konwersji tego punktu z rzutowanego układu współrzędnych (EPSG 3857) na geograficzny układ współrzędnych (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()

Zwraca to dla twojego punktu współrzędne -105.150271116 39.7278572773.

ustroetz
źródło
6

Wynik nie jest przestrzennym / współrzędnym układem odniesienia , to para współrzędnych. Musisz wiedzieć, czym jest odniesienie przestrzenne, aby ponownie rzutować współrzędne.

W tym przypadku nie jest to jednak wymagane. Wystarczy przekazać odpowiednie odniesienie przestrzenne wyjściowe do usługi, a zwróci współrzędne w Lon / Lat.

Oto strona ze współrzędnymi wyjściowymi w formacie Lon / Lat przy użyciu geograficznego systemu odniesienia przestrzennego WGS-84 ( EPSG 4326 ).

użytkownik2856
źródło
6

afalciano ma poprawną odpowiedź, ale chciał uwzględnić wariant użycia pyproj.

To nie wymaga wiesz ciąg proj4 i jest odrobinę szybciej.

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
Marcel Wilson
źródło
2
Nie potrzebujesz ciągu proj4, zastąp go drugim wierszem, p = pyproj.Proj(init='epsg:3857')a wynik będzie taki sam.
alphabetasoup
1
Rezultat jest taki sam, ale ostatnio sprawdziłem, że to było trochę szybciej.
Marcel Wilson,
3

Wypróbowałem kod sugerowany przez Marcela Wilsona i jest on szybszy:

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
użytkownik1174460
źródło
1

Znalazłem ten post, gdy szukam sposobów na zrobienie tego w QGIS. Jak opisano tutaj , zastosowana metoda wygląda tak:

def convertProjection(self,x,y,from_crs,to_crs):
    crsSrc = QgsCoordinateReferenceSystem(from_crs)
    crsDest = QgsCoordinateReferenceSystem(to_crs)
    xform = QgsCoordinateTransform(crsSrc, crsDest)
    pt = xform.transform(QgsPoint(x,y))
    return pt.x, pt.y

# Remove the "EPSG:" part
from_crs = 3857
to_crs = 4326
x = -11705274.6374    
y = 4826473.6922
lon, lat = self.convertProjection(x,y,from_crs, to_crs)
Toivo Säwén
źródło
2
Uwaga, w QGIS 3 nastąpiła przełomowa zmiana API, więc jeśli używasz v3.x, musisz użyćxform = QgsCoordinateTransform(crsSrc, crsDest, QgsProject.instance())
Jonny
0

Należy pamiętać, że transformfunkcja pyprojakceptuje również arrays, co jest dość przydatne, jeśli chodzi o ramki danych

import pandas as pd
from pyproj import Proj, transform

df = pd.DataFrame({'x': [-11705274.6374]*100, 
                   'y': [4826473.6922]*100})
inProj, outProj = Proj(init='epsg:3857'), Proj(init='epsg:4326')
df['x2'], df['y2'] = transform(inProj, outProj, df['x'].tolist(), df['y'].tolist())
J. Doe
źródło
-1
import cv2
import numpy as np
def onMouse(event, x, y, flags, param):
    if event == cv2.EVENT_LBUTTONDOWN:
       # draw circle here (etc...)
       print('x = %f, y = %f'%((x*2/100),(y*2/100)))
a=float(input("enter length of original dimension in cm"))
b=float(input("enter width of original dimension in cm"))
print("origional image coordinates")
im=cv2.imread('mask1.jpg',0)
re_im=cv2.resize(im,(round(a*100/2),round(b*100/2)))
cv2.imshow('image',re_im)
cv2.setMouseCallback('image', onMouse)
cv2.waitKey(0)
cv2.destroyAllWindows()
Subhash Verma
źródło