Suite

Effectuer une transformation avec pyproj ?


Cela semble devoir être simple, mais j'ai du mal à obtenir la bonne transformation de mon plan d'état en lat long.

Mon code est le suivant

# projection de spatialreference.org/ref/epsg/2992/ # oregon est espg ref 2992 projection = r'+proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=399999.9999984 +y_0 =0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048 +no_defs' osp2992 = Proj(projection)

Maintenant, lorsque j'essaie de vérifier les entrées/sorties, elles s'affichent dans le coin supérieur droit

# Coordonnées d'entrée : -119.88481361812, 43.453242550659 # Coordonnées de sortie : 1475673.108308, 621492.180223 >>>osp2992(-119.88481361812, 43.453242550659) (449785.16341222316, 189430.81653184193)

C'est loin.

Je ne sais pas où je me suis trompé.

Edit: suppression de la conversion to_meter et toujours une mauvaise sortie

projection = r'+proj=lcc +lat_1=43 +lat_2=45.5 +lat_0=41.75 +lon_0=-120.5 +x_0=399999.9999984 +y_0=0 +ellps=GRS80 +datum=NAD83 +no_defs' >>>osp2992(- 119.88481361812, 43.453242550659) (550261.2548834425, 157725.95224630416)

Vos coordonnées de sortie de test sont en pieds, tandis que pyproj vous donne la réponse en mètres, ce qu'il fait par défaut

Définissez plutôt votre projection de cette façon :

osp2992 = Proj(projection, preserve_units=True) osp2992(-119.88481361812, 43.453242550659) (1475673.1083078186, 621492.18022257846)