14

我尝试了以下方法，输入：纬度/经度数据，然后我将计算它周围的一个框，比如说 50 m，所以东/北值是 +/- 50 m。

现在我将它重新转换为纬度/经度并使用脚本：

http://robotics.ai.uiuc.edu/~hyoon24/LatLongUTMconversion.py我得到的结果是不可能的，lon 之前是 7 左右，之后是 2 左右。

zone, easting, northing = LLtoUTM(23, location.get_lat(), location.get_lon()) 

topUTM = northing + error
bottomUTM = northing - error
leftUTM = easting - error
rightUTM = easting + error
left, top = UTMtoLL(23, leftUTM, topUTM, zone)

是我的代码中的错误，还是脚本有缺陷？

所以我尝试使用 pyproj，只是 lat/lon 到 utm 到 lat/lon 看看会发生什么

>>> p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
>>> p
<pyproj.Proj object at 0x7ff9b8487dd0>
>>> x,y = p(47.9941214, 7.8509671)
>>> print x,y
5159550.36822 1114087.43925
>>> print p(x,y,inverse=True)
(47.971558538495991, 7.8546573140162605)

在这里，它并不像上面的脚本那么遥远，但它似乎仍然非常不正确，以至于无法使用它。怎么来的？我该怎么做才能获得更准确的结果？

编辑：

我运行了 test() 并且所有测试都通过了。

在 epsg 文件中没有这样的东西。我发现的最接近的是：

<32632> +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs <>

没有 tmerc。另外我需要什么将 towgs84 作为参数传递？上面那些？

4

4 回答 4

52

我上周为 Python 创建了一个小型 UTM 转换库，并将其上传到 Python 包索引： http: //pypi.python.org/pypi/utm

我将它与使用 pyproj 进行了比较，它更快、更准确。给定您的样本数据，结果如下：

>>> import utm

>>> u = utm.from_latlon(47.9941214, 7.8509671)
>>> print u
(414278, 5316285, 32, 'T')

>>> print utm.to_latlon(*u)
(47.994157948891505, 7.850963967574302)

更新：下面的理查兹回答描述了这个问题的真正解决方案。

于 2013-01-05T20:23:06.297 回答
29

错误在您的代码中。

首先，其他答案之一中列出的 PyProj 问题是真实的。您应该检查您的 epsg 文件并确保它包含该行

<2392> +proj=tmerc +lat_0=0 +lon_0=24 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=intl +towgs84=-90.7,-106.1,-119.2,4.09,0.218,-1.05,1.37 +units=m +no_defs no_defs <>

注意towgs84参数。

您对 PyProj 的问题源于误用了投影命令。

如果我们取 47.9941214N, 7.8509671E 并转换为 UTM，我们会得到 Zone 32, 414278 Easting, 5316286 Northing。

您执行以下 PyProj 操作：

p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
>>> x,y = p(47.9941214, 7.8509671)
>>> print x,y
5159550.36822 1114087.43925
>>> print p(x,y,inverse=True)
(47.971558538495991, 7.8546573140162605)

但是，如果我们查阅 PyProj文档，我们会看到以下内容：

使用参数 lon 调用 Proj 类实例，lat 会将 lon/lat（以度为单位）转换为 x/y 本机地图投影坐标（以米为单位）。

让我们再次尝试运行 OP 的 PyProj 操作，但切换 lon/lat 参数的顺序：

p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
>>> x,y = p(7.8509671, 47.9941214)
>>> print x,y
414278.16731 5316285.59492
>>> print p(x,y,inverse=True)
(7.850967099999812, 47.994121399999784)

操作完全反转（几乎）！

要回答您问题的第一部分，如果您查看http://robotics.ai.uiuc.edu/~hyoon24/LatLongUTMconversion.py的定义UTMtoLL，您会发现以下内容：

UTMtoLL(ReferenceEllipsoid, northing, easting, zone)

然而，您使用UTMtoLL(23, leftUTM, topUTM, zone)where leftUTM 是 Easting 而 topUTM 是 Northing。

因此，对于您的第一个脚本和 PyProj，您使用了错误的参数顺序。

这是一个很好的提醒，在暗示别人的错误之前，一定要反复（或三次）检查你的工作。也就是说，Python 的文档并不是最好的，而 PyProj 的文档在这种情况下充其量是神秘的。这个命令的一个很好的基于网络的解释并附有它的使用示例可能会阻止你的焦虑。

于 2013-09-04T18:10:30.803 回答
4

我对pyproj没有问题，试试下面的代码

from pyproj import Proj

Lat = 52.063098675
Lon = -114.132980348 #Calgary

ZoneNo = "11" #Manually input, or calcuated from Lat Lon
myProj = Proj("+proj=utm +zone="+ZoneNo+",\
+north +ellps=WGS84 +datum=WGS84 +units=m +no_defs") #north for north hemisphere
UTMx, UTMy = myProj(Lon, Lat)

########################################

#UTM ==> Lat Lon:
ZoneNo = "11" #Manually input or from other sources
myProj = Proj("+proj=utm +zone="+\
ZoneNo+", +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
Lon2, Lat2 = myProj(UTMx, UTMy,inverse=True)

print Lat2
print Lon2
于 2017-05-14T22:01:36.773 回答
2

您对 pyProj 的问题听起来就像这里描述的那样：

http://code.google.com/p/pyproj/issues/detail?id=3

已解决：

solved! in epsg file there must be

<2392> +proj=tmerc +lat_0=0 +lon_0=24 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=intl +towgs84=-90.7,-106.1,-119.2,4.09,0.218,-1.05,1.37 +units=m +no_defs no_defs <>

note the towgs84 parameter!

如果您想继续使用 pyproj，请检查该线程。

另外，test()模块的功能是否有效？您是否尝试过test目录中随附的任何脚本？

于 2011-07-21T15:28:54.070 回答