python – lat / lon to utm to lat / lon是非常有缺陷的,怎么样?
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了python – lat / lon to utm to lat / lon是非常有缺陷的,怎么样?,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含1563字,纯文字阅读大概需要3分钟。
内容图文
![python – lat / lon to utm to lat / lon是非常有缺陷的,怎么样?](/upload/InfoBanner/zyjiaocheng/708/229ed7e641f44f598d06c342f5860acb.jpg)
我试过以下,
输入: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 – lat / lon to utm to lat / lon是非常有缺陷的,怎么样?全部内容,希望文章能够帮你解决python – lat / lon to utm to lat / lon是非常有缺陷的,怎么样?所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。