将 ESRI 投影坐标转换为 lat-lng答案

【问题标题】：Convert ESRI projection coordinates to lat-lng将 ESRI 投影坐标转换为 lat-lng
【发布时间】：2010-11-12 03:14:42
【问题描述】：

我在“NAD 1983 StatePlane Michigan South FIPS 2113 Feet”（又名 ESRI 102690）中有一个大型 x,y 坐标数据集。我想将它们转换为 lat-lng 点。

理论上，proj 可以处理这种情况，但文档并没有给我任何线索——它似乎描述了更复杂的情况。

我尝试过使用python interface，如下所示：

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

【问题讨论】：

标签： mapping esri proj4js proj

【解决方案1】：

如果您在 TXT、CSV 或 XLS 文件中有坐标，您可以按 CTRL+C 并将它们插入到http://cs2cs.mygeodata.eu，您可以在其中设置适当的输入和所需的输出坐标系。可以插入数千个各种格式的坐标...

【讨论】：

【解决方案2】：

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

例如，使用 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}

所以在你的情况下是这样的：

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/

【讨论】：

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