python坐标表示_已知经纬度坐标求两点间距离,用python表示
已知经纬度坐标求两点间距离,用python表示
已知地球上任意两点(lon1, lat1),( lon2, lat2)的经纬度坐标,求两点间的距离用haversine表示:
1.首先先将经纬度坐标的角度化成弧度(mlon1,mlat1),(mlon2,mlat2)
2.利用公式:
其中a表示两点维度的差值,即 a = mlat1 - mlat2
其中a表示两点经度的差值,即 b = mlon1 - mlon2
其中r表示地球的半径
用python表示为
def geo_distance(lon1, lat1, lon2, lat2): #地理距离
"""
Calculate the great circle distance between two points
on the earth (specified in decimal degrees)
"""
lon1, lat1, lon2, lat2 = map(radians, map(float, [lon1, lat1, lon2, lat2])) #根据提供的函数对指定序列做映,radians:将角度转换为弧度。
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
c = 2 * asin(sqrt(a))
r = 6371
return c * r
print(geo_distance(111.5,36.08,121.47,31.23)) #我老家到上海的距离
结果
1067.9651712083848 #单位 km
其实python自带haversine库,可以直接调用,这样就可以不记公式啦!
from haversine import haversine
linfen = (111.5,36.08)
shanghai = (121.47,31.23)
dis = haversine(linfen,shanghai)
print(dis)
结果
1133.5471931316163 #单位 km
可以看出还是有差距的,那那个准确呢? _
用那种方法知道了吧!嘿嘿
原文链接:https://blog.csdn.net/weixin_40264772/article/details/104952973
python坐标表示_已知经纬度坐标求两点间距离,用python表示相关推荐
- 已知经纬度坐标求两点间距离,用python表示
已知经纬度坐标求两点间距离,用python表示 已知地球上任意两点(lon1, lat1),( lon2, lat2)的经纬度坐标,求两点间的距离用haversine表示: 1.首先先将经纬度坐标的角 ...
- ArcGis中已知经纬度坐标转任意坐标系(以点和80坐标系为讲解模板)
文章首发微信公众号:ArcGisdada 一般经纬度都是以度分秒的格式显示的,如图31-1所示. 图31-1: 我们转至WGS-84坐标需要进行十进制的转化,也就是把度分秒换算为以度为单位的一串数字. ...
- wgs84 投影坐标转经纬度 代码_ArcGis中已知经纬度坐标转任意坐标系(以点和80坐标系为讲解模板)...
一般经纬度都是以度分秒的格式显示的,如图31-1所示. 图31-1: 我们转至WGS-84坐标需要进行十进制的转化,也就是把度分秒换算为以度为单位的一串数字.操作步骤:打开excel>选中代表经 ...
- lisp批量生成轴线_已知大量坐标和直径,请问如何以用autocad 的lisp 以坐标为球心批量画球...
展开全部 lisp里用62616964757a686964616fe59b9ee7ad9431333363383432command调用CAD球体命令:sphere (command "sp ...
- python 方差齐性检验_已知F和sig值_方差齐性 sig_如何判断方差齐性_怎么判断方差齐不齐...
已知F和sig值,怎么判断方差齐性 用spss做两个独立样本的t检验,首先得到Levene的结果是 F=0.063,Sig=0.806 Q1:据说只用Sig(p)值判断,如果大于0.05,就是方差齐, ...
- 已知三角形三边长怎么求面积_已知三角形三边求面积的公式——海伦公式
海伦公式又译作希伦公式.海龙公式.希罗公式.海伦-秦九韶公式,传说是古代的叙拉古国王希伦(Heron,也称海龙)二世发现的公式,是一个利用三角形的三条边长直接求三角形面积的公式.下面我们利用初中的知识 ...
- c++ 圆上任意点坐标计算_已知圆上任意三点坐标如何编程来计算这个圆的圆心和半径...
在过去的一周时间里一直在计算师兄交代的任务 ,啥都没学,就一直在干这活.其中呢,我想的一个算法中就用到了本文中的这个小知识(主要就是线性代数和高中的一些知识结合起来).师兄说做的这个东西的内容要保密, ...
- 已知法向量 求投影_已知一点坐标和平面方程,求此点在平面上的投影点的坐标。...
展开全部 平面Ax+By+cZ+D=0的法向量为n=(A,B,C),M点在平面上的投影M'的坐标e69da5e887aa3231313335323631343130323136353331333431 ...
- c++ 圆上任意点坐标计算_已知圆弧的一些信息,求得圆弧上中间点的坐标的方法(C++语言描述)...
//从圆弧一些信息得到圆弧上一个特殊的点(计算的是二维图形的情况)Point3d DBOPERATION::GetArcTangencyPoint(Point3d pStartPoint, Point ...
最新文章
- 地铁框架保护的原理_地铁屏蔽门是如何保证通讯的稳定?
- 紧急求助!配置SMTP插件出错,SMTP connect() failed
- java-多线程 String 常量池缓存功能带来的问题
- 笔记-项目质量管理-质量保证和质量控制(信管网)
- 作者:刘松灵(1988-),男,就职于中国人民银行征信中心数据部
- Android 应用中十大导航设计错误
- SVN代码更新群发邮件通知
- python漏洞扫描器爬虫_Python系列之——漏洞平台厂商列表爬虫
- 网页转圈打不开服务器,打不开网页怎么回事
- 京东广告受众基础属性预估
- STM32CubeIDE配置使用
- EasyUI TreeGrid各行调整上下位置
- 入手python时踩过的坑——缩进篇
- 计算机三级网络技术考试上机,计算机等级考试三级网络技术上机答案
- 漏洞复现篇——利用XSS漏洞实现键盘记录
- satisfy with用法
- VCO电路中的电源设计
- 数据不平衡问题解决方法——欠采样
- java 基础知识点
- 软件设计师学习笔记-知识产权