var
DEF_PI
=
Math
.
PI
;
// 3.14159265359
var DEF_2PI = 2 * Math . PI ; // 6.28318530712
var DEF_PI180 = Math . PI / 180.0 ; // 0.01745329252
var DEF_R = 6370996.81 ; // radius of earth
// 利用弧夹角方式计算地图上两点之间的距离,和百度地图计算结果不一致
function GetShortDistance ( lon1 , lat1 , lon2 , lat2 ) {
var ew1 , ns1 , ew2 , ns2 ;
var dx , dy , dew ;
var distance ;
// 角度转换为弧度
ew1 = lon1 * DEF_PI180 ;
ns1 = lat1 * DEF_PI180 ;
ew2 = lon2 * DEF_PI180 ;
ns2 = lat2 * DEF_PI180 ;
// 经度差
dew = ew1 - ew2 ;
// 若跨东经和西经180 度,进行调整
if ( dew > DEF_PI )
dew = DEF_2PI - dew ;
else if ( dew < - DEF_PI )
dew = DEF_2PI + dew ;
dx = DEF_R * Math . cos ( ns1 ) * dew ; // 东西方向长度(在纬度圈上的投影长度)
dy = DEF_R * ( ns1 - ns2 ); // 南北方向长度(在经度圈上的投影长度)
// 勾股定理求斜边长
distance = Math . sqrt ( dx * dx + dy * dy );
return distance ;
}
// 利用弧夹角方式计算地图上两点之间的距离,与百度方式一致
function GetLongDistance ( lon1 , lat1 , lon2 , lat2 ) {
var ew1 , ns1 , ew2 , ns2 ;
var distance ;
// 角度转换为弧度
ew1 = lon1 * DEF_PI180 ;
ns1 = lat1 * DEF_PI180 ;
ew2 = lon2 * DEF_PI180 ;
ns2 = lat2 * DEF_PI180 ;
// 求大圆劣弧与球心所夹的角(弧度)
distance = Math . sin ( ns1 ) * Math . sin ( ns2 ) + Math . cos ( ns1 ) * Math . cos ( ns2 ) * Math . cos ( ew1 - ew2 );
// 调整到[-1..1]范围内,避免溢出
if ( distance > 1.0 )
distance = 1.0 ;
else if ( distance < - 1.0 )
distance = - 1.0 ;
// 求大圆劣弧长度
distance = DEF_R * Math . acos ( distance );
return distance ;
}