# -*- coding:utf-8 -*-
# /usr/bin/python
import warnings

warnings.filterwarnings("ignore")

from math import radians, cos, sin, asin, sqrt
方法1
def geodistance(lat1,lng1,lat2,lng2):

''''数学转换'''

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
方法2
from geopy.distance import geodesic
def fromGeodesic(lat1,lng1,lat2,lng2):

'''调用包'''

distance=geodesic((lat1,lng1), (lat2,lng2))

return distance
方法3
def simpleDistance(lat1,lng1,lat2,lng2):

dx = lng1 - lng2 # 经度差值

dy = lat1 - lat2 # 纬度差值

b = (lat1 + lat2) / 2.0 # 平均纬度

Lx = radians(dx) * 6367.0 * cos(radians(b))# 东西距离

Ly = 6367.0 * radians(dy) #南北距离

distance = sqrt(Lx * Lx + Ly * Ly)

return distance
print(geodistance(39.941,116.45,39.94,116.451))

print(fromGeodesic(39.941,116.45,39.94,116.451))

print(simpleDistance(39.941,116.45,39.94,116.451))