Converter ESRI coordenadas de projecção de latitude e longitude

Question

Eu tenho um grande conjunto de dados de coordenadas x, y em "NAD 1983 stateplane Michigan Sul FIPS 2113 Pés" (aka ESRI 102690 ). Eu gostaria de convertê-los em pontos de latitude e longitude.

Em teoria, isso é algo proj é construído para lidar, mas a documentação não tenha dado -me uma pista -. parece descrever casos muito mais complicado

Eu tentei usar um python interface de , assim:

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

Mas isso dá saída descontroladamente incorreta.

Há uma biblioteca Javascript, mas eu tenho ~ 50.000 pontos ao cabo, de modo que não parece apropriado.

---

O que funcionou para mim:

Eu criei um arquivo chamado ptest com cada par em sua própria linha, coordenadas xey separados por um espaço, assim:

13304147.06410000000 288651.94040000000
...

Então eu alimentei esse arquivo para o comando canalizado e os resultados para um arquivo de saída:

$>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 você só precisa reprojetar e pode fazer alguma mineração de dados em seus arquivos de texto usar o que você gosta e usa http://spatialreference.org/ref/esri/102690/ como referência.

Por exemplo utilização do Proj4 e armazená-lo em um arquivo de shell / cmd e chamar o seu arquivo de entrada com proj4 (/ versão Linux Windows disponível) nenhum problema com o tamanho do seu conjunto de dados.

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}

Assim, no seu caso, algo como isto:

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 você tem coordenadas em TXT, CSV ou arquivo XLS, você pode fazer CTRL + C e inseri-los para http: // cs2cs.mygeodata.eu onde você pode definir entrada apropriado e saída desejada sistema de coordenadas. É possível inserir milhares de coordenadas em vários formatos ...