Convertir proyección ESRI coordenadas latitud y la longitud

Question

Tengo un gran conjunto de datos de coordenadas x, y en "NAD 1983 stateplane Michigan Sur FIPS 2113 Pies" (también conocido como ESRI 102690 ). Me gustaría convertirlos en puntos lat-lng.

En teoría, esto es algo proj está construida para manejar, pero la documentación no ha dado yo ni idea -. que parece describir casos mucho más complicados

He intentado usar un pitón interfaz , así:

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

Pero eso da salida tremendamente incorrecto.

Hay una biblioteca de Javascript, pero tengo ~ 50.000 puntos de manejar, por lo que no parece apropiado.

---

Lo que funcionó para mí:

I creado un archivo llamado ptest con cada par en su propia línea, coordenadas x e y separadas por un espacio, así:

13304147.06410000000 288651.94040000000
...

A continuación, di de comer a ese archivo en el comando para el abastecimiento de los resultados a un archivo de salida:

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

Si sólo necesita reproyectar y puede hacer algo de la minería de datos en sus archivos de texto usar lo que quiera y usar http://spatialreference.org/ref/esri/102690/ como referencia.

Por ejemplo utilizar el Proj4 y almacenarla en un archivo de cáscara / cmd y llamar a cabo el archivo de entrada con proj4 (/ Linux versión de Windows está disponible) no hay problema con el tamaño del conjunto de datos.

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}

por lo que en su caso algo como esto:

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

Si tiene coordenadas en TXT, CSV o XLS, puede hacerlo CTRL + C e insertarlos a http: // cs2cs.mygeodata.eu donde se puede establecer de entrada apropiada y la salida deseada del sistema de coordenadas. Es posible insertar miles de coordenadas en varios formatos ...