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
投稿猜你喜欢
1.折线图 plt.plot()常用的一些参数：颜色（color)：‘c’ 青红（cyan）&- 比如，在使用DWR的时候，如果你想传递下拉框的参数到后台的话，此时就需要先获取到下拉框的值了。 其实想要获取到下拉框的值是很简单的。 最关键
- python数据化运营数据化运营的核心是运营，所有数据工作都是围绕运营工作链条展开的，逐步强化数据对于运营工作的驱动作用。数据化运营的价值体
- 为了方便使用分类，我定义了一个分类表category，里面字段是id(自动编号) cat_name(分类名) pare
- 读写分离，顾名思义，我们可以把读和写两个操作分开，减轻数据的访问压力，解决高并发的问题。那么我们今天就Python两大框架来做这个读写分离的
- 当我们使用电脑的时候，不可避免的会碰到重复点击的任务。所以我们可以使用python制作一个简单的连点器进行使用，同时可以提升写代码的水平。第
- 本文实例讲述了JS添加或删除HTML dom元素的方法。分享给大家供大家参考，具体如下：创建新的 HTML 元素如需向 HTML DOM 添
- Vue.js是当下很火的一个JavaScript MVVM库，它是以数据驱动和组件化的思想构建的。相比于Angular.js，Vue.js提
- 本文实例讲述了PHP面向对象程序设计继承用法。分享给大家供大家参考，具体如下：<?php//继承是子类(派生类)从父类(基类，超类)继
- 这两天一直在看tensorflow中的读取数据的队列，说实话，真的是很难懂。也可能我之前没这方面的经验吧，最早我都使用的theano，什么都
- 首先呢我去安装了一个那个pytorch，然后导入一下发现：连numpy都找不到，于是我表示很生气重新安装，它说安装过了，地址是balabal
- Dump ALL MySQL Databasesmysqldump --user=XXXXXXXX --password=XXXXXXX -
- 本文实例讲述了Python实现将蓝底照片转化为白底照片功能。分享给大家供大家参考，具体如下：import cv2import numpy a
- 质数又称素数。一个大于1的自然数，除了1和它自身外，不能被其他自然数整除的数叫做质数；否则称为合数。一、判断一个数是否为素数：基于定义def
- 今天想说的是内容和容器的关系，顺便把之前设计中碰到的问题和大家一起探讨下。我们从软件的设置说起。（这里以QQ的设置举例）一个软件的设置（常称
- python程序运行中，可由程序抛出异常。异常触发：使用raise命令抛出异常，即可使用异常基类Exception，也可使用自定义异常类（继
- 首先导入包含apriori算法的mlxtend库，pip install mlxtend调用apriori进行关联规则分析，具体代码如下，其
- 今天我们来一起看看Dreamweaver MX 2004在加密FTP 传送 方面的新功能。我们一般在做
- 前文介绍了Oracle 中实现数据透视表的几种方法，今天我们来看看在 MySQL/MariaDB 中如何实现相同的功能。本文使用的示例数据可
- 今天来分享python学习的一个小例子，使用python暴力破解mysql数据库，实现方式是通过UI类库tkinter实现可视化面板效果，在
