"""
Created on Thu Jun 1 14:41:06 2017
@author: Jipon
"""
from math import *
# input Lat_A 纬度A
# input Lng_A 经度A
# input Lat_B 纬度B
# input Lng_B 经度B
# output distance 距离(km)
"""
基于googleMap中的算法得到两经纬度之间的距离
"""
def calcDistance(Lat_A, Lng_A, Lat_B, Lng_B):
EARTH_RADIUS = 6378.137;#赤道半径(单位km)
radLat1 = radians(Lat_A)
radLat2 = radians(Lat_B)
a=radLat1-radLat2
b=radians(Lng_A)-radians(Lng_B)
s = 2 * asin(sqrt(pow(sin(a/2),2)+cos(radLat1)*cos(radLat2)*pow(sin(b/2),2)));
s = s * EARTH_RADIUS;
#s = Math.round(s * 10000d) / 10000d;
return s;
Lat_A=32.060255; Lng_A=118.796877 # 南京
Lat_B=39.904211; Lng_B=116.407395 # 北京
distance=calcDistance(Lat_A,Lng_A,Lat_B,Lng_B)
print("----",distance)
