Koordinatenumwandlung mit pyproj

Question

ich nicht in gis erlebt bin Koordinatenumwandlung aber gelungen, diese Seite mit: http://all-geo.org/volcan01010/2012/11/change-coordinates-with-pyproj/ zu konvertieren Shape-Datei von EPSG Koordinaten: 28992 bis EPSG: 4326 pyproj mit diesen Aussagen der Python-Modul mit:

wgs84=pyproj.Proj("+init=EPSG:4326") 
epsg28992=pyproj.Proj("+init=EPSG:28992") 
pyproj.transform(epsg28992, wgs84,x,y) 

Wenn ich diese Koordinaten in Google Maps umkehre und sie gebe, geben sie mir korrekte Standorte. Das funktioniert also gut.

Jetzt habe ich eine andere Shapefile (s) und ich schaue auf die Shapefile.prj-Datei, um zu bestimmen, welche Projektion verwendet wurde. Der ESRI WKT entspricht ESRI: 102686 was ich hier finde: http://epsg.io/102686 Da der ESRI: 102686 Code von pyproj nicht bekannt ist (gibt Fehler), muss ich die proj4-Notation verwenden, die ich von der gleichen Seite (http://epsg.io/102686) erhalten habe:

wgs84=pyproj.Proj("+init=EPSG:4326") 
esri102686=pyproj.Proj("+proj=lcc +lat_1=41.71666666666667 +lat_2=42.68333333333333 +lat_0=41 +lon_0=-71.5 +x_0=200000 +y_0=750000.0000000001 +datum=NAD83 +units=us-ft +no_defs") 
pyproj.transform(esri102686, wgs84,x,y) 

Ich bekomme z und Koordinaten verwenden diese in Google Maps: 60,275122729462495, -61,873986125999316 , die irgendwo im Ozean ist ...

Aber meine Ergebnisse sollten in Cambridge, MA in den USA sein, so mehr um: 41,00000, -71, 5000000

Was mache ich falsch?

Answer 1

gelöst, hinzugefügt preserve_units = True, wie folgt aus:

esri102686 = pyproj.Proj("+proj=lcc +lat_1=41.71666666666667 +lat_2=42.68333333333333    +lat_0=41 +lon_0=-71.5 +x_0=200000 +y_0=750000.0000000001 +datum=NAD83 +units=us-ft  +no_defs",preserve_units= True) 

Jetzt funktioniert es gut. Wenn das optionale Schlüsselwort 'preserve_units' den Wert 'True' hat, werden die Einheiten in den Kartenprojektionskoordinaten nicht zu Zählern gezwungen. Siehe here.