python3-经纬度的相关计算

#from math import radians, cos, sin, asin, sqrt
from math import *
def haversine(lon1, lat1, lon2, lat2): # 经度1，纬度1，经度2，纬度2 （十进制度数）
# 将十进制度数转化为弧度
# math.degrees(x):为弧度转换为角度
# math.radians(x):为角度转换为弧度
lon1, lat1, lon2, lat2 = map(radians, [lon1, lat1, lon2, lat2])
# haversine公式
dlon = lon2 - lon1
dlat = lat2 - lat1
a = sin( dlat /2 ) **2 + cos(lat1) * cos(lat2) * sin( dlon /2 ) **2
c = 2 * asin(sqrt(a))
r = 6371 # 地球平均半径，单位为公里
return c * r*1000

if __name__=='__main__':
a=haversine(103.85012,36.03653,103.84085,36.03842)
print(a)
#注：该结果的数值存在问题

#通过经纬度，找出中心点
from math import *
def center_geolocation(geolocations):
x = 0
y = 0
z = 0
lenth = len(geolocations)
for lon, lat in geolocations:
lon = radians(float(lon))
lat = radians(float(lat))

x += cos(lat) * cos(lon)
y += cos(lat) * sin(lon)
z += sin(lat)

x = float(x / lenth)
y = float(y / lenth)
z = float(z / lenth)

return (degrees(atan2(y, x)), degrees(atan2(z, sqrt(x * x + y * y))))

if __name__ == '__main__':
locations = [[116.568627, 39.994879], [116.564791, 39.990511], [116.575012, 39.984311]]
print(center_geolocation(locations))
##结果
(116.56947685224392, 39.98990040970147)