python实现两个经纬度点之间的距离和方位角的方法
作者：bboyzqh 发布时间：2022-03-15 02:41:27
标签：python,经纬度,距离,方位角
最近做有关GPS轨迹上有关的东西，花费心思较多，对两个常用的函数总结一下，求距离和求方位角，比较精确，欢迎交流！
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
来源：https://blog.csdn.net/zhuqiuhui/article/details/53180395
0
投稿猜你喜欢
前言在启动 Django 项目时，Django 默认监听的端口号为 8000，设置的默认 IP 地址为 127.0.0.1 。如果需要修改默- 楼主在做公司项目的时候遇到url重定向的问题，因此上网简单查找，作出如下结果由于使用的是语言是python所以以下是python的简单解决方
- 当在函数中定义默认值时，值初始化只会进行一次，就是执行到def methodname时执行。看下面代码:from datetime impo
- 此次遇到的是一个函数使用不熟练造成的问题，但有了分析工具后可以很快定位到问题（此处推荐一个非常棒的抓包工具fiddler） 正文如
- 就目前互联网上大小网站而言，大部分都是采用ASP+ACCESS/SQL Server或者PHP+MySQL来编写；事实上，ASP和MySQL
- Socket是什么Socket是一个抽象层，连接了应用层与传输层，是一组接口。在设计模式中，Socket其实就是一个门面模式，它把复杂的TC
- 前言人类都是视觉动物，不管是男生还是女生看到漂亮的小姐姐、小哥哥就想截图保存下来。可是截图会对画质会产生损耗，截取的画面不规整，像素不高等问
- 介绍还记得你在小学时学习如何加减数字吗？现在，你也可以对图像做同样的事情!输入图像可以进行算术运算，例如加法、减法和按位运算（AND、OR、
- 本文通过将同一个数据集在三种不同的简便项窗口部件中显示。三个窗口的数据得到实时的同步,数据和视图分离。当添加或删除数据行，三个不同的视图均保
- 无聊,写了一个图片后加载的东东,感谢asfman帮忙<!DOCTYPE html PUBLIC "-//W3C//DTD X
- 1，判断图像清晰度，明暗，原理，Laplacian算法。偏暗的图片，二阶导数小，区域变化小；偏亮的图片，二阶导数大，区域变化快。import
- 学习目的 接触验证控件 昨天介绍了SQL SERVER插入数据，但是我们省略了验证输入这一步。以前的做法是用语句逐个判断输入的正确性，如是否
- 由于工作需要,所以前一阵子将IE升级到了8.0,结果今天发现出现一个问题,eWebEditor的在线编辑器不好用了,仔细想想,肯定是IE8搞
- 本文实例讲述了Python使用time模块实现指定时间触发器。分享给大家供大家参考，具体如下：其实很简单，指定某个时间让脚本处理一个事件，比
- 首先说明，伪造访问来路不是什么光明正大的事情，目的就是为了欺骗服务器。原本以为给 XMLHTTP 对象增加一个 Referer 的heade
- 比较好奇python对于多进程中copy on write机制的实际使用情况。目前从实验结果来看，python 使用multiprocess
- 已知有一个XML文件（bookstore.xml）如下：<?xml version="1.0" encoding=
- 本文较为深入的探究了php中in_array函数用法。分享给大家供大家参考。具体如下：今天突然想到php中的in_array函数有个其怪的用
- 这篇文章主要介绍了基于Python批量生成指定尺寸缩略图代码实例,文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值
- 【原文地址】 Recipe: Deploying a SQL Database to a Remote Hosting Environmen
