首先是公式,百度即可知道,需要注意的点是python中不能直接使用三角函数,需要import math库,并且在三角函数中输入的变量并不是角度,而是弧度,由于经纬度从某种意义上来说是角度,所以要换算成弧度再进行计算
C = math.sin(Mlat_A*pi/180) * math.sin(Mlat_B*pi/180) *
math.cos((Mlng_A - Mlng_B)*pi/180) +math.cos(Mlat_A*pi/180) * math.cos(Mlat_B*pi/180)
import math
lng_A = 113.37935836
lat_A = 31.71785761
lng_B = 114.3933937
lat_B = 30.63068048
R = 6371.004
pi = 3.141592654
#经度:113.379