Suite

Transformation Pyproj renvoyant des valeurs incorrectes (EPSG 3875 à EPSG 4326)


J'essaie de convertir EPSG 3875 en EPSG 4326 et j'ai ce code :

de pyproj import Proj, transform 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) return transform(inProj,outProj,x,y) x, y = 4152144, 219039 x2,y2 = convert_coordinates("3875", "4326", x, y) print x2, y2

Cela revient 112.264404642 154.392043378 qui est hors limites pour EPSG 4326. Mais ArcGIS renvoie 37.2993441706517, 1.96727415792103 qui est dans les limites et a l'air mieux.

Y a-t-il quelque chose qui ne va pas dans ce code ?


La bonne procédure est

de pyproj import Proj, transform inProj = Proj("+init=EPSG:3857")) outProj = Proj("+init=EPSG:4326") x, y = 4152144, 219039 print transform(inProj,outProj,x,y ) (37.29934417065166, 1.96727415792103)

Votre erreur estconvertir_coordonnées("3875", "4326", x, y): 3875 et 4326 sont des chaînes simples et non des projections PyProj

Regardez la transformation de coordonnées pyProj incorrecte