Convertire proiezione ESRI coordinate lat-lng

Question

Ho un grande insieme di dati di coordinate x, y in "NAD 1983 StatePlane Michigan Sud FIPS 2113 Feet" (aka ESRI 102690 ). Mi piacerebbe convertirli in punti lat-LNG.

In teoria, questo è qualcosa proj è costruito per gestire, ma la documentazione non ha me un indizio -. sembra descrivere i casi più complicati

Ho provato con un pitone interfaccia , in questo modo:

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

Ma che dà uscita selvaggiamente errato.

C'è una libreria JavaScript, ma ho ~ 50.000 punti per gestire, in modo che non sembra appropriato.

---

Che cosa ha funzionato per me:

Ho creato un file chiamato ptest con ciascuna coppia su una riga, coordinate xey separati da uno spazio, in questo modo:

13304147.06410000000 288651.94040000000
...

Poi ho dato da mangiare quel file nel comando e convogliato i risultati a un file di output:

$>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

Se avete solo bisogno di riproiettare e può fare un po 'di data mining sui vostri file di testo utilizzare quello che ti piace e utilizzare http://spatialreference.org/ref/esri/102690/ come riferimento.

Ad esempio utilizzare il Proj4 e memorizzare in un file di shell / cmd e chiamare il vostro file di input con Proj4 (Linux / Windows versione disponibile) non è un problema con la dimensione set di dati.

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}

Quindi nel tuo caso qualcosa di simile:

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

Se si dispone di coordinate in TXT, CSV o XLS, si può fare CTRL + C e inserirli a http: // cs2cs.mygeodata.eu dove è possibile impostare ingresso appropriato e output desiderato sistema di coordinate. E 'possibile inserire migliaia di coordinate in vari formati ...