python经纬度转换xy坐标公式_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:
pass
3 验证
选取深圳野生动物园(22.599578, 113.973129)为起点,深圳坪山站(22.6986848, 114.3311032)为终点,结合百度地图、谷歌地图等进行效果验证。
程序运行结果如下:
百度测距为38.3km
Google地图手动测距为39.31km
距离与角度均无问题。
到此这篇关于Python经纬度坐标转换为距离及角度的实现的文章就介绍到这了,更多相关Python经纬度坐标转换为距离及角度内容请搜索我们以前的文章或继续浏览下面的相关文章希望大家以后多多支持我们!
本文标题: Python经纬度坐标转换为距离及角度的实现
本文地址: http://www.cppcns.com/jiaoben/python/360192.html
python经纬度转换xy坐标公式_Python经纬度坐标转换为距离及角度的实现相关推荐
- python xml转换键值对_Python 提取dict转换为xml/json/table并输出
#!/usr/bin/python #-*- coding:gbk -*- #设置源文件输出格式 import sys import getopt import json import createD ...
- python坐标表示_python经纬度转xy坐标系
经纬度转换成xy坐标def millerToXY (lon, lat): """ :param lon: 经度 :param lat: 维度 :return: " ...
- 经纬度转换XY坐标,并计算距离
import pandas as pd import numpy as np from pyproj import Transformer import matplotlib.pyplot as pl ...
- 经纬度坐标转换为距离及角度(Python)
最近项目上有这样的需求,需要依据设备的经纬度坐标计算距离及角度.经验证后效果较好,并分享. 1 经纬度转换距离代码 #!/usr/bin/env python # -*- coding: utf-8 ...
- 腾讯地图api php经纬度转换地址,腾讯地图经纬度转换为百度地图经纬度
利用微信来获取用户当前的详细地址是非常方便的,但是利用微信获取到用户的经纬度转换为具体地址时,老是出现误差很大,因此,决心把腾讯地图经纬度转换城百度地图经纬度,然后再利用百度地图api来获取详细 ...
- python文件转换成jar包_Python一键转Jar包,Java调用Python新姿势!
粉丝朋友们,不知道大家看故事看腻了没(要是没腻可一定留言告诉我^_^),今天这篇文章换换口味,正经的来写写技术文.言归正传,咱们开始吧! 今天的这篇文章,聊一个轩辕君之前工作中遇到的需求:如何在Jav ...
- python oop求三角形面积公式_Python面向对象编程-OOP
Python面向对象编程-OOP 20141216 Chenxin整理 OOP的3大特性: 封装,继承,多态 一.封装 OOP目的,OOP为了代码重用 :分解代码 ,最小化代码的冗余以及对现在的代码进 ...
- python输入字符串转换为公式_python从字符串转换为boolean ?
实际上,你只需将字符串与期望接受的内容进行比较,这样你就可以这样做: s== 'True'或者检查一组值: 复制代码 sin ['true', '1', 't', 'y', 'yes', 'yeah' ...
- python数字转换成中文大写_python初学者笔记(2):阿拉伯数字转换成中文大写
题:输入一个数字,转换成中文大写的写法 可运行的程序(Python 2.7.9): 1 #-*- coding: utf-8 -*- #在python2的py文件里面写中文,必须要添加一行声明文件编码 ...
- python单位转换编程英寸厘米_Python字节单位转换实例
我就废话不多说了,直接上代码! from enum import Enum class Values(): values={'B':1} @staticmethod def getValues(): ...
最新文章
- python调用c++返回带成员指针的类指针
- .Net Discovery系列文章阅读索引--带你探索未知的.Net世界
- h5移动端如何让页面强制横屏
- 前端:jQuery笔记
- php中函数参数个数问题
- 新网站SEO要做的事情有哪些
- Kubernetes的系统架构与设计理念
- 学习笔记:Oracle的trace文件可见性
- C语言通过枚举网卡,API接口可查看man 7 netdevice--获取接口IP地址
- python输入一个序列_Python序列合并,python
- wxParse解析iframe播放视频
- 【华为机试题 HJ102】字符统计
- matlab中uigetfile的使用
- Word 关闭拼写检查 (去掉Word中拼写检查的所有红色和绿色的浪线)
- 标准差SD、相对标准偏差RSD学习和python实现
- 插入(希尔)排序时间、空间复杂度
- com.alibaba.fastjson.JSONException: syntax error, pos 1, line 1, column 2
- Android Retrofit详解(retrofit:2.3.0)
- linux ftp下载文件到本地,如何使用get指令利用FTP下载文件到本地
- VS SDK更新问题(error MSB8036: 找不到 Windows SDK 版本10.0.14393.0)