pyproj - Comment décider les points initiaux de toute zone UTM

Question

Je suis en train de convertir X, système Y à la latitude et la longitude, et il semble que beaucoup de gens recommandent pyproj

Dans leur exemple de code, ils utilisent p1 = Proj(init='epsg:26915') et p2 = Proj(init='epsg:26715')

Mais c'est pour UTM Zone15 Voici mon échantillon x = 492893,53, y = 5458938,38 et mon exemple de code pour UTM Zone10N:

p1 = Proj(init='epsg:26910') 
p2 = Proj(init='epsg:26710') 
x,y = 492893.53, 5458938.38 
x2,y2 = transform(p1, p2, x, y) 

Le x2 de sortie, y2 ne sont pas la latitude droite et la longitude

J'ai essayé de vérifier la valeur epsg pour UTM Zone 10N, il y a 3 versions en ligne: version 1, version 2, version 3 J'ai essayé tous, aucun d'entre eux travaillent. Que veut vraiment dire p1, p2 et comment puis-je trouver les valeurs de p1, p2 pour toute zone UTM?

Answer 1

Si vous convertissez en lat-long, vous avez besoin de la projection pour votre système de coordonnées lat-long. Dans la plupart des cas est WGS84 Lat-temps utilisé par GPS, qui est le code EPSG 4326.

import pyproj 
ll = pyproj.Proj(init="epsg:4326") 

alors vous avez besoin d'une projection de vos points de données. Dans ce cas, vous dites son UTM10, mais vous ne dites pas laquelle des variantes, mais il n'y a aucun moyen de dire à partir des coordonnées. Tout d'abord permet de essayer la version WGS84:

>>> utm10 = pyproj.Proj(init="epsg:32610") 
>>> pyproj.transform(utm10, ll, 492893.53, 5458938.38) 
(-123.0977143330712, 49.28315708900033) 

Cela ressemble à un possible lat long (en fait (-123, 49) est long de latitude) pour la région valide de UTM10. Donc, je ne suis pas sûr de savoir ce que «aucun d'eux ne fonctionne» dans votre question.

Les deux autres codes de projection donné des résultats proches ou identiques. Lequel vous utilisez ne peut pas être déterminées à partir des seuls chiffres - vous devez demander celui qui a fourni les chiffres que la projection qu'ils sont dans

Vous ne pouvez pas déterminer même quelle zone UTM points sont donné votre coordonnées x et y.. Ils pourraient être n'importe quelle zone de la planète, puisque les nombres x et y sont dupliqués pour chaque zone dans la bande de longitude de chaque zone.