Tôi đang viết một tập lệnh python đọc nhiều tệp XML chứa tọa độ x và y và kết hợp tất cả chúng thành một tệp csv. Vĩ độ và Kinh độ là các trường bắt buộc trong csv, nhưng tôi gặp khó khăn khi chuyển đổi tọa độ x, y trong Máy bay Bang Bắc Ohio usFt thành 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)
Cả hai phương pháp trên đều trả về cùng một kết quả, tuy nhiên thời gian dài này là một nơi nào đó trong Vịnh Hudson. Khi tôi vẽ tọa độ trong ArcMap, độ dài chính xác là: -81.142311,41.688205.
* Lưu ý tất cả các độ dài lat được cung cấp dài, lat vì đây là thứ tự Proj sử dụng
Có ai biết tại sao tôi lại nhận được tọa độ sai từ Proj.4 và pyproj không?