Я пишу сценарій python, який читає кілька XML-файлів, що містять координати x і y, і об'єднує їх у один файл csv. Широта та Довгота - обов'язкові поля в csv, але у мене виникають труднощі при перетворенні координат x, y в штаті Огайо, штат Північна Держава usFt, до WGS84.
>>> p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters
>>> p(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)
>>> p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code
>>> p1(739400.91,2339327.3,inverse=True)
(-80.138057868777224, 60.278230707978487)
Обидва вищезазначені методи дають однаковий результат, однак ця давній час знаходиться десь у Гудзоновій бухті. Коли я будувати координати в ArcMap, правильна довжина ширини: -81.142311,41.688205.
* Зверніть увагу на те, що всі тривалість Lat надається довга, lat, оскільки такий порядок використовує Proj
Хтось знає, чому я отримав би неправильні координати від Proj.4 та pyproj?