Я предполагаю, что это основной вопрос, но я не могу найти или распознать решение.
Этот сайт возвращается
Point:
X: -11705274.6374
Y: 4826473.6922
при поиске с первым значением ключа 000090 в качестве примера. Я предполагаю, что это пространственная привязка, и я понимаю, что это такое.
Я ищу инструкции или примеры того, как преобразовать это в широту и долготу, используя Python.
python
coordinate-system
convert
Винсент
источник
источник
По умолчанию сайт, на который вы ссылаетесь, использует систему пространственной привязки EPSG 3857 (WGS84 Web Mercator). Я нашел эту информацию здесь .
Вы можете либо указать другую систему пространственной привязки, введя нужный EPSG в форму,
Spatial Reference
либо вы можете преобразовать возвращенные координаты с помощью Python.Например, вы можете использовать привязки GDAL Python для преобразования этой точки из проецированной системы координат (EPSG 3857) в географическую систему координат (EPSG 4326).
Это возвращает для вашей точки координаты
-105.150271116 39.7278572773
.источник
У afalciano правильный ответ, но он хотел включить вариант использования pyproj.
Это действительно требует , чтобы Вы знаете строку Proj4 и чуть - чуть быстрее.
источник
p = pyproj.Proj(init='epsg:3857')
и результат будет таким же.Выходные данные не являются пространственной / координатной системой координат , это пара координат. Вы должны знать, что такое пространственная привязка, чтобы перепроектировать координаты.
Тем не менее, это не требуется в этом случае. Просто передайте соответствующую выходную пространственную привязку службе, и она вернет координаты в Lon / Lat.
Вот страница с выходными координатами в формате Lon / Lat с использованием системы географической пространственной привязки WGS-84 ( EPSG 4326 ).
источник
Попробовал код, предложенный Марселем Уилсоном, и он быстрее:
источник
Я нашел этот пост, когда искал способы сделать это в QGIS. Как описано здесь , используемый метод выглядит так:
источник
v3.x
вы используете его, вам нужно будет его использоватьxform = QgsCoordinateTransform(crsSrc, crsDest, QgsProject.instance())
Обратите внимание, что
transform
функцияpyproj
принимает такжеarrays
, что весьма полезно, когда речь идет о фреймах данныхисточник
источник