标签：python latitude-longitude utm

我试过以下,
输入：lat / lon数据
然后我会计算一个围绕它的方框,比方说50米,所以/ – 东方/北方值为50米.

现在我将它重新转换为lat / lon并使用脚本：

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来指示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作为参数？以上的？

解决方法:

我上周为Python创建了一个小的UTM转换库,并将其上传到Python Package Index：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)

更新：下面的Richards answer描述了这个问题的真正解决方案.

标签：python,latitude-longitude,utm
来源： https://codeday.me/bug/20190928/1828943.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。