Python计算地图上两点经纬度间的距离

处理地图数据时，经常需要用到两个地理位置间的距离。比如A点经纬度（110.0123, 23.32435），B点经纬度（129.1344,25.5465），求AB两点之间的距离。我们可以用haversine()函数求出距离结果。Python版本的haversine()如下所示：

 

[python] view plain copy
    1. from math import radians, cos, sin, asin, sqrt  
    2.   
    3. def haversine(lon1, lat1, lon2, lat2): # 经度1，纬度1，经度2，纬度2 （十进制度数）  
    4.     """ 
    5.     Calculate the great circle distance between two points  
    6.     on the earth (specified in decimal degrees) 
    7.     """  
    8.     # 将十进制度数转化为弧度  
    9.     lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])  
    10.   
    11.     # haversine公式  
    12.     dlon = lon2 - lon1   
    13.     dlat = lat2 - lat1   
    14.     a = sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2  
    15.     c = 2 * asin(sqrt(a))   
    16.     r = 6371 # 地球平均半径，单位为公里  
    17.     return c * r * 1000 

posted on 2018-06-22 17:42  枫飞飞  阅读(1533)  评论(0)  编辑  收藏  举报