Mapping 将ESRI投影坐标转换为lat lng

我有一个很大的x，y坐标数据集，在“NAD 1983密歇根州南部FIPS 2113英尺”（又名）中。我想把它们转换成lat lng点

从理论上讲，这是为处理而构建的，但文档并没有给我任何线索——它似乎描述了更复杂的案例

我试过用一个，像这样：

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

13304147.06410000000 288651.94040000000
...

但这给出了极不正确的输出

有一个Javascript库，但我要处理大约50000个点，所以这似乎不合适

---

对我有用的是：

我创建了一个名为

ptest

的文件，每一对都在自己的行上，x和y坐标用空格分隔，如下所示：

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

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

---

如果您只需要重新投影，并且可以对文本文件进行一些数据挖掘，请使用您喜欢的任何内容作为参考

例如，使用Proj4并将其存储在shell/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}

---

如果您在TXT、CSV或XLS文件中有坐标，您可以执行CTRL+C并将它们插入到可以设置适当输入和所需输出坐标系的位置。可以以各种格式插入数千个坐标…

我仍然对cs2cs的用法感到困惑。我如何输入我的两个样本点？如何格式化一个包含多个点的文件？如何将该文件交给cs2cs？文档似乎建议它是一个带有值标签的文件，每对在自己的行上，最后一行是“EOF”。我在某处有一个工作示例，但我记得它比人们想象的要简单，在所有参数之后，我使用-E将输入坐标获取到输出。第一个参数是源文件。所有的东西都是标准输出的，所以只需使用cs2cs+++多个参数++++源文件>输出文件谢谢，这比看起来容易。谢谢