Konvertieren ESRI Projektionskoordinaten zu lat-lng

Question

Ich habe eine große Datenmenge von x, y-Koordinaten in "NAD 1983 StatePlane Michigan Süd FIPS 2113 Feet" (auch bekannt als ESRI 102690 ). Ich möchte sie zu Lat-lng Punkte konvertieren.

In der Theorie ist dies etwas, proj zu handhaben gebaut, aber die Dokumentation nicht gegeben hat ich ein Anhaltspunkt -. es scheint viel komplizierte Fälle zu beschreiben

Ich habe versucht mit einer Python-Schnittstelle , etwa so:

from pyproj import Proj
p = Proj(init='esri:102690')
sx = 13304147.06410000000 #sample points
sy = 288651.94040000000
x2, y2 = p(sx, sy, inverse=True)

Aber das gibt wild falsche Ausgabe.

Es gibt eine Javascript-Bibliothek, aber ich habe ~ 50.000 Punkte zu handhaben, so dass nicht angemessen scheint.

---

Was für mich gearbeitet:

habe ich eine Datei mit jedem Paar genannt ptest in einer eigenen Zeile, x und y Koordinaten, die durch einen Raum getrennt sind, etwa so:

13304147.06410000000 288651.94040000000
...

Dann fütterte ich die Datei in dem Befehl und geleitet, die Ergebnisse in einer Ausgabedatei:

$>cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 
+lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 
+datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to 
+proj=latlon ptest > out.txt

Answer 1

Wenn Sie nur neu projizieren müssen und können einige Data-Mining auf Textdateien tun verwenden, was auch immer Sie mögen und verwenden http://spatialreference.org/ref/esri/102690/ als Referenz.

Zum Beispiel die Proj4 verwenden und speichern sie in einer Shell / cmd-Datei und rufen Sie Ihre Eingabedatei mit proj4 (Linux / Windows-Version verfügbar) kein Problem mit der Größe Ihres Datensatzes.

cs2cs +proj=latlong +datum=NAD83 +to +proj=utm +zone=10  +datum=NAD27 -r <<EOF
cs2cs -f %.16f +proj=utm +zone=20N +to +proj=latlon - | awk '{print $1 " " $2}

so in Ihrem Fall so etwas wie folgt aus:

cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 +lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to +proj=latlon

http://trac.osgeo.org/proj/wiki/man_cs2cs

http://trac.osgeo.org/proj/

Answer 2

Wenn Sie Koordinaten in TXT, CSV oder XLS-Datei haben, können Sie STRG + C tun, und legen Sie sie auf http: // cs2cs.mygeodata.eu wo man entsprechende Eingangs- und gewünschtes Ausgangskoordinatensystem festlegen. Es ist möglich, Tausende von Koordinaten in verschiedenen Formaten einfügen ...