Преобразование «GLNX/Y» в долготу, широту

Question

У меня есть координаты данных, состоящие из двух 6-значных чисел (например, 300 000 250 000), которые я хочу преобразовать в long., Lat. координаты. Мне сказали, что данные координат (называемые GLNX, GLNY) исходят из системы координат штата Мичиган, номер EPSG ESRI: 102121 Однако, когда я передаю этот номер 102121 в gdal.ImportFromEPSG, он жалуется, что он не знает, , Два вопроса:

1. Как создать SpatialReference для ESRI: 102121
2. Могу ли я передать мой 6-значный номер пары непосредственно перепроецировать, или мне нужно «подогнать» его, например, масштаб некоторых мощность 10, или конвертировать из футов в метры, или что?

Answer 1

Я не знаю, как вы используете GDAL API, но с GDAL 2.0 через Python, это работает для меня:

from osgeo import osr 
osr.UseExceptions() 
sr = osr.SpatialReference() 
sr.ImportFromEPSG(102121) # returns 0 for success, which I get 

Но я подозреваю, что это не работает, как описано в вашем вопрос. Таким образом, вы можете импортировать из proj.4 кода вместо, который вы можете получить от http://epsg.io/102121 или добавить .proj4 расширения к исходному коду:

import urllib2 
srid = 102121 
response = urllib2.urlopen('http://epsg.io/%d.proj4' % (srid,)) 
sr.ImportFromProj4(response.read()) # returns 0 for success 
print(sr.ExportToPrettyWkt()) # shows that it is understood 

Код proj.4, в конечном счете используется libproj сделать фактическую проекцию , а не WKT.