Daha çox

Yanlış dəyərləri qaytaran Pyproj çevirməsi (EPSG 3875 - EPSG 4326)


EPSG 3875-dən EPSG 4326-ya çevirməyə və bu kodu əldə etməyə çalışıram:

pyproj idxal Proj-dən çevirmək def convert_coordinates (input_format, output_format, x, y): input_format = "+ init = epsg:" + str (input_format) output_format = "+ init = epsg:" + str (output_format) inProj = Proj ( input_format) outProj = Proj (output_format) geri dönüş (inProj, outProj, x, y) x, y = 4152144, 219039 x2, y2 = convert_coordinates ("3875", "4326", x, y) print x2, y2

Bu qayıdır 112.264404642 154.392043378 bu EPSG 4326 üçün məhduddur. Lakin ArcGIS geri qayıdır 37.2993441706517, 1.96727415792103 hüdudlarda olan və daha yaxşı görünən.

Bu kodda açıq şəkildə səhv bir şey varmı?


Doğru prosedur budur

pyproj idxalından Proj, çevirmə inProj = Proj ("+ init = EPSG: 3857")) outProj = Proj ("+ init = EPSG: 4326") x, y = 4152144, 219039 çap çevirmə (inProj, outProj, x, y ) (37.29934417065166, 1.96727415792103)

Səhviniz budurconvert_coordinates ("3875", "4326", x, y): 3875 və 4326, sadə iplərdir, PyProj proyeksiyaları deyil

PyProj Koordinat Transformasiyasına Yanlış baxın