python求球面距离
今天我们用python代码写一个求球面距离的公式,可以通过经纬度,计算球面间两点的距离。
代码如下:
def geodistance(lng1,lat1,lng2,lat2):
lng1, lat1, lng2, lat2 = map(radians, [float(lng1), float(lat1), float(lng2), float(lat2)]) # 经纬度转换成弧度
dlon=lng2-lng1
dlat=lat2-lat1
a=sin(dlat/2)**2 + cos(lat1) * cos(lat2) * sin(dlon/2)**2
distance=2*asin(sqrt(a))*6371*1000 # 地球平均半径,6371km
distance=round(distance/1000,3)
return distance
lng1,lat1,lng2,lat2分别是两对经纬度。通过上述代码可以计算这两对经纬度地点之间的距离。