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;
}