Python根据经纬度求两点间距离

痛定思痛。 2022-08-20 14:13 460阅读 0赞

Python程序根据经纬度求两点之间距离的两种方法

from math import *
    
    '''计算两点之间直线距离'''
    def get_distance1(lon_a, lat_a, lon_b, lat_b):
        radlat1 = radians(lat_a)
        radlat2 = radians(lat_b)
        a = radlat1 - radlat2
        b = radians(lon_a) - radians(lon_b)
        s = 2 * asin(sqrt(pow(sin(a/2),2) + cos(radlat1) * cos(radlat2)*pow(sin(b/2),2)))
        earth_radius = 6378137
        s = s * earth_radius
        return s
    
    '''计算两点之间直线距离'''
    def get_distance2(lon_a, lat_a, lon_b, lat_b):
        if abs(lon_a - lon_b) < 0.000001 and abs(lat_a - lat_b) < 0.000001:
            return 0
        re = 6378140  # 赤道半径 (m)
        rp = 6356755  # 极半径 (m)
        oblateness = (re - rp) / re  # 地球扁率
        rad_lat_a = radians(lat_a)
        rad_lon_a = radians(lon_a)
        rad_lat_b = radians(lat_b)
        rad_lon_b = radians(lon_b)
        atan_a = atan(rp / re * tan(rad_lat_a))
        atan_b = atan(rp / re * tan(rad_lat_b))
        tmp = acos(sin(atan_a) * sin(atan_b) + cos(atan_a) * cos(atan_b) * cos(rad_lon_a - rad_lon_b))
        if tmp == 0:
            return 0
        c1 = (sin(tmp) - tmp) * (sin(atan_a) + sin(atan_b)) ** 2 / cos(tmp / 2) ** 2
        c2 = (sin(tmp) + tmp) * (sin(atan_a) - sin(atan_b)) ** 2 / sin(tmp / 2) ** 2
        dr = oblateness / 8 * (c1 - c2)
        distance = re * (tmp + dr)
        return distance
    
    def main():
        lon_a = 116.320866
        lat_a = 39.814506
        lon_b = 116.301322
        lat_b = 39.87874
        print(get_distance1(lon_a, lat_a, lon_b, lat_b))
        print(get_distance2(lon_a, lat_a, lon_b, lat_b))
    
    if __name__ == '__main__':
        main()

发表评论取消回复

表情：

评论列表（有 0 条评论，460人围观）

还没有评论，来说两句吧...

相关阅读

相关获取地理经纬度计算两点间距离

/ 距离计算 export function distance(la1, lo1, la2, lo2) \{ var La1 = la1 \* Math.PI /...

今天药忘吃喽~/ 2024年04月19日 15:08/ 0 赞/ 81 阅读

相关 mysql根据经纬度求两点距离

1.两点距离(1.4142135623730951) select st_distance(point(0,0),point(1,1)); sele...

淡淡的烟草味﹌/ 2024年04月17日 05:24/ 0 赞/ 87 阅读

相关根据两点经纬度计算距离

经纬度简介这些经纬线是怎样定出来的呢？地球是在不停地绕地轴旋转（地轴是一根通过地球南北两极和地球中心的假想线），在地球中腰画一个与地轴垂直的大圆圈，使圈上的每一点

╰+攻爆jí腚メ/ 2024年02月17日 17:09/ 0 赞/ 16 阅读

相关根据经纬度坐标计算两点的距离

现在如果有个业务需求，就是要做一个根据定位坐标实现计算距离并找到附近的店铺的推荐功能，已知数据库中会存储店铺的经纬度坐标。那么，这个需求已经很明确了，在要求不高的情况下，只要计

桃扇骨/ 2023年05月31日 06:54/ 0 赞/ 12 阅读

相关根据两点经纬度计算距离和角度

1.distance ![复制代码][48304ba5e6f9fe08f3fa1abda7d326ab.png] /

秒速五厘米/ 2023年01月13日 05:41/ 0 赞/ 266 阅读

相关根据两点经纬度计算距离

经纬度简介这些经纬线是怎样定出来的呢？地球是在不停地绕地轴旋转（地轴是一根通过地球南北两极和地球中心的假想线），在地球中腰画一个与地轴垂直的大圆圈，使圈上的每一点都

超、凢脫俗/ 2022年10月01日 02:47/ 0 赞/ 325 阅读

相关 Python根据经纬度求两点间距离

Python程序根据经纬度求两点之间距离的两种方法 from math import '''计算两点之间直线距离''' def get_d

痛定思痛。/ 2022年08月20日 14:13/ 0 赞/ 461 阅读

相关 js根据经纬度计算两点距离

js版-胡老师 google.maps.LatLng.prototype.distanceFrom = function(latlng) \{ var la

傷城~/ 2022年07月12日 09:43/ 0 赞/ 505 阅读

相关根据两点经纬度计算距离

\[b\]\[color=green\]\[size=medium\]这些经纬线是怎样定出来的呢？地球是在不停地绕地轴旋转（地轴是一根通过地球南北两极和地球中心的假想线）

r囧r小猫/ 2022年04月13日 20:10/ 0 赞/ 545 阅读

相关计算两个经纬度点间的距离

计算两个经纬度点间的距离纬度线投射在图上看似水平的平行线，但实际上是不同半径的圆。有相同特定纬度的所有位置都在同一个纬线上。赤道的纬度为0°，将行星平分为南半球和北

你的名字/ 2021年10月25日 13:30/ 0 赞/ 947 阅读