根据经纬度计算两点间距离

# 根据经纬度计算两点间距离
# 经度 long 纬度 lat
def GetDistance( lng1, lat1, lng2, lat2):
u '''''计算两点间球面距离 单位为m'''
EARTH_RADIUS = 6378.137 # 地球周长/2*pi 此处地球周长取40075.02km pi=3.1415929134165665
from math import asin,sin,cos,acos,radians, degrees, pow ,sqrt, hypot,pi
# 方法0
# 最简单的求平面两点间距离 误差比较大
d = hypot(lng2 - lng1,lat2 - lat1) * 40075.02 / 360 * 1000
print 'd0=' ,d
# 方法1
# ''' d=111.12cos{1/[sinΦAsinΦB+cosΦAcosΦBcos(λB—λA)]} '''
# ''' 其中A点经度、纬度分别为λA和ΦA,B点的经度、纬度分别为λB和ΦB,d为距离。'''
m = cos(radians(lat1)) * cos(radians(lat2)) * cos(radians(lng2 - lng1))
x = 1 / (sin(radians(lat1)) * sin(radians(lat2)) + m)
d = 111.12 * cos(radians( 1 / x))
print 'd1=' ,d
# 方法2
# 据说来源于 google maps 的脚本
radLat1 = radians(lat1) # a点纬度(单位是弧度)
radLat2 = radians(lat2) # b点纬度(单位是弧度
a = radLat1 - radLat2 # 两点间的纬度弧度差
b = radians(lng1) - radians(lng2) # 两点间的经度弧度差
s = 2 * asin(sqrt( pow (sin(a / 2 ), 2 ) + cos(radLat1) * cos(radLat2) * pow (sin(b / 2 ), 2 ))) # 两点间的弧度
s = s * EARTH_RADIUS
# s = round(s * 10000) / 10000 # 四舍五入保留小数点后4位
d = s * 1000
print 'd2=' ,d
# 方法3
# ''' 经纬坐标为P(x1,y1) Q(x2,y2) '''
# D=arccos[cosy1*cosy2*cos(x1-x2)+siny1*siny2]*2*PI*R/360
d = acos(cos(radians(lat1)) * cos(radians(lat2)) * cos(radians(lng1 - lng2)) + sin(radians(lat1)) * sin(radians(lat2))) * EARTH_RADIUS * 1000
print 'd3=' ,d
return d
GetDistance( 116.95400 , 39.95400 , 116.95300 , 39.95300 )
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值