Python：计算经纬度距离

噶噶~

已于 2022-07-22 10:38:43 修改

阅读量5.3k

点赞数 3

分类专栏：机器学习--数据处理统计学文章标签： python

于 2022-07-22 10:18:57 首次发布

本文链接：https://blog.csdn.net/weixin_42587745/article/details/125926903

版权

机器学习--数据处理同时被 2 个专栏收录

48 篇文章 4 订阅

订阅专栏

统计学

19 篇文章 4 订阅

订阅专栏

这篇博客介绍了两种使用Python计算经纬度之间距离的方法。第一种利用Haversine公式，将经纬度转换为弧度并进行计算；第二种则借助Scipy库的cdist函数，通过欧氏距离计算。两种方法的结果以米为单位，适用于地理定位和距离估算。

摘要由CSDN通过智能技术生成

Python计算经纬度的方法：
这里介绍两种方式，其最终结果都是以【米】作为单位，虽然结果值稍微有误差，但误差不大。

1. List item

方式1：

维基百科公式（要求的是公式中的d）：
在这里插入图片描述
Python代码如下（返回结果是单位是：米）：

from math import radians,sin,cos,asin,sqrt
def haversine_dis(lon1, lat1, lon2, lat2):
    #将十进制转为弧度
    lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])
    
    #haversine公式
    d_lon = lon2 - lon1
    d_lat = lat2 - lat1
    aa = sin(d_lat/2)**2 + cos(lat1)*cos(lat2)*sin(d_lon/2)**2
    c = 2 * asin(sqrt(aa))
    r = 6371 # 地球半径，千米
    return c*r*1000

方式2：

from scipy.spatial.distance import cdist
lonlat1=[[102.8083989,24.8962294]]    # 经度、维度
lonlat2=[[102.471938,24.912491]]   # 经度、维度
distance=cdist(lonlat1,lonlat2,metric='euclidean')
distance=list(list(list(distance))[0]*100000)[0]