更新时间：2020年11月01日 10:44:43 作者：_Seven°
这篇文章主要介绍了Python经纬度坐标转换为距离及角度的实现，文中通过示例代码介绍的非常详细，对大家的学习或者工作具有一定的参考学习价值，需要的朋友们下面随着小编来一起学习学习吧
最近项目上有这样的需求，需要依据设备的经纬度坐标计算距离及角度。经验证后效果较好，并分享。
1 经纬度转换距离代码
#!/usr/bin/env python
# -*- coding: utf-8 -*-
__author__ = 'Seven'
import math
# 计算距离
def getDistance(latA, lonA, latB, lonB):
ra = 6378140 # 赤道半径
rb = 6356755 # 极半径
flatten = (ra - rb) / ra # Partial rate of the earth
# change angle to radians
radLatA = math.radians(latA)
radLonA = math.radians(lonA)
radLatB = math.radians(latB)
radLonB = math.radians(lonB)
pA = math.atan(rb / ra * math.tan(radLatA))
pB = math.atan(rb / ra * math.tan(radLatB))
x = math.acos(math.sin(pA) * math.sin(pB) + math.cos(pA) * math.cos(pB) * math.cos(radLonA - radLonB))
c1 = (math.sin(x) - x) * (math.sin(pA) + math.sin(pB)) ** 2 / math.cos(x / 2) ** 2
c2 = (math.sin(x) + x) * (math.sin(pA) - math.sin(pB)) ** 2 / math.sin(x / 2) ** 2
dr = flatten / 8 * (c1 - c2)
distance = ra * (x + dr)
distance = round(distance / 1000, 4)
return f'{distance}km'
2 经纬度转化角度代码
#!/usr/bin/env python
# -*- coding: utf-8 -*-
__author__ = 'Seven'
import math
# 计算角度
def getDegree(latA, lonA, latB, lonB):
radLatA = math.radians(latA)
radLonA = math.radians(lonA)
radLatB = math.radians(latB)
radLonB = math.radians(lonB)
dLon = radLonB - radLonA
y = math.sin(dLon) * math.cos(radLatB)
x = math.cos(radLatA) * math.sin(radLatB) - math.sin(radLatA) * math.cos(radLatB) * math.cos(dLon)
brng = math.degrees(math.atan2(y, x))
brng = round((brng + 360) % 360, 4)
brng = int(brng)
if (brng == 0.0) or ((brng == 360.0)):
return '正北方向'
elif brng == 90.0:
return '正东方向'
elif brng == 180.0:
return '正南方向'
elif brng == 270.0:
return '正西方向'
elif 0 < brng < 90:
return f'北偏东{brng}'
elif 90 < brng < 180:
return f'东偏南{brng - 90}'
elif 180 < brng < 270:
return f'西偏南{270 - brng}'
elif 270 < brng < 360:
return f'北偏西{brng - 270}'
else:
选取深圳野生动物园(22.599578, 113.973129)为起点，深圳坪山站(22.6986848, 114.3311032)为终点，结合百度地图、谷歌地图等进行效果验证。
程序运行结果如下：
百度测距为38.3km
Google地图手动测距为39.31km
距离与角度均无问题。
到此这篇关于Python经纬度坐标转换为距离及角度的实现的文章就介绍到这了,更多相关Python经纬度坐标转换为距离及角度内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！
您可能感兴趣的文章:
-
-
python如何调用现有的matlab函数
这篇文章主要介绍了python如何调用现有的matlab函数问题，具有很好的参考价值，希望对大家有所帮助。如有错误或未考虑完全的地方，望不吝赐教
2022-12-12
-
-
-
python利用rsa库做公钥解密的方法教程
RSA是一种公钥密码算法，RSA的密文是对代码明文的数字的 E 次方求mod N 的结果。下面这篇文章主要给大家介绍了关于python利用rsa库做公钥解密的方法教程，文中通过示例代码介绍的非常详细，需要的朋友可以参考下。
2017-12-12
-
-
浅谈django 重载str 方法
这篇文章主要介绍了浅谈django 重载str 方法，具有很好的参考价值，希望对大家有所帮助。一起跟随小编过来看看吧
2020-05-05
-
-
简单的Python人脸识别系统
这篇文章主要介绍了Python人脸识别系统的实现，文中讲解非常详细，代码帮助大家更好的理解和学习，感兴趣的朋友可以了解下
2020-07-07
-
-
