Python经纬度坐标转换为距离及角度的实现_1km是多少经纬度资源-CSDN文库

经纬度坐标转换

72 浏览量 2021-01-19 23:31:14 上传评论 1 收藏 361KB PDF 举报

资源推荐

资源详情

资源评论

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:

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余2页未读，立即下载

内容反馈

weixin_38637878

粉丝: 3
资源: 926

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈

feedback-tip