转换ESRI突起坐标，以经纬度-eng

Question

我有一个大数据集的x、y坐标"NAD1983年StatePlane密歇根南智2113脚"(aka ESRI102690).我想把他们纬度lng点。

在理论上，这是什么东西 proj 建立来处理，但该文件没有给我一个线索--这似乎描述的复杂得多的情况。

我已经尝试采用一个 蟒蛇口, 是这样的：

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

但是，让疯狂不正确的输出。

有一个Javascript图书馆,但是我有-50,000点处理，因此，似乎并不合适。

---

什么工作对我来说：

我创建了一个文件叫 ptest 每对对其自己的路线、x和y坐标分隔的空间，就像这样：

13304147.06410000000 288651.94040000000
...

然后我送这个文件成的命令和管道的结果为一个输出的文件:

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

如果你仅仅需要重新投影和可以做一些数据挖掘在你的文本文件使用的任何你喜欢的和使用 http://spatialreference.org/ref/esri/102690/ 作为参考。

例如使用Proj4并将其存储在一个壳/cmd文件和喊出你的输入文件proj4(linux/windows版)没有任何问题的大小的数据集。

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}

所以你的情况是这样的：

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

如果你有坐标在TXT，CSV或XLS文件，你可以做的CTRL+C和插入他们 http://cs2cs.mygeodata.eu 在这里你可以设置适当的输入和输出所需的坐标系统。它可以插入成千上万的坐标在各种不同的格式...