python实现两个经纬度点之间的距离和方位角
内容导读
互联网集市收集整理的这篇技术教程文章主要介绍了python实现两个经纬度点之间的距离和方位角,小编现在分享给大家,供广大互联网技能从业者学习和参考。文章包含1346字,纯文字阅读大概需要2分钟。
内容图文
![python实现两个经纬度点之间的距离和方位角](/upload/InfoBanner/zyjiaocheng/1143/028f3ed0892d4e95879ca3a4ae09a39c.jpg)
from:http://blog.csdn.net/zhuqiuhui/article/details/53180395
1. 求两个经纬点的方位角,P0(latA, lonA), P1(latB, lonB)(很多博客写的不是很好,这里总结一下)
def getDegree(latA, lonA, latB, lonB): """ Args: point p1(latA, lonA) point p2(latB, lonB) Returns: bearing between the two GPS points, default: the basis of heading direction is north """ radLatA = radians(latA) radLonA = radians(lonA) radLatB = radians(latB) radLonB = radians(lonB) dLon = radLonB - radLonA y = sin(dLon) * cos(radLatB) x = cos(radLatA) * sin(radLatB) - sin(radLatA) * cos(radLatB) * cos(dLon) brng = degrees(atan2(y, x)) brng = (brng + 360) % 360 return brng
2. 求两个经纬点的距离函数:P0(latA, lonA), P1(latB, lonB)
def getDistance(latA, lonA, latB, lonB): ra = 6378140 # radius of equator: meter rb = 6356755 # radius of polar: meter flatten = (ra - rb) / ra # Partial rate of the earth # change angle to radians radLatA = radians(latA) radLonA = radians(lonA) radLatB = radians(latB) radLonB = radians(lonB) pA = atan(rb / ra * tan(radLatA)) pB = atan(rb / ra * tan(radLatB)) x = acos(sin(pA) * sin(pB) + cos(pA) * cos(pB) * cos(radLonA - radLonB)) c1 = (sin(x) - x) * (sin(pA) + sin(pB))**2 / cos(x / 2)**2 c2 = (sin(x) + x) * (sin(pA) - sin(pB))**2 / sin(x / 2)**2 dr = flatten / 8 * (c1 - c2) distance = ra * (x + dr) return distance
原文:http://www.cnblogs.com/lindaxin/p/8057999.html
内容总结
以上是互联网集市为您收集整理的python实现两个经纬度点之间的距离和方位角全部内容,希望文章能够帮你解决python实现两个经纬度点之间的距离和方位角所遇到的程序开发问题。 如果觉得互联网集市技术教程内容还不错,欢迎将互联网集市网站推荐给程序员好友。
内容备注
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 gblab@vip.qq.com 举报,一经查实,本站将立刻删除。
内容手机端
扫描二维码推送至手机访问。