// 獲取兩個經(jīng)緯度的距離 (km)
getDistance(long1,lat1,long2,lat2) {
var earth_radius = 6378.137;
var a,b,distance,sa2,sb2;
lat1 = this.rad(lat1);
lat2 = this.rad(lat2);
a = lat1-lat2;
// b = this.rad(long1-long2);
b = this.rad(long1)-this.rad(long2);
sa2 = Math.sin(a/2.0);
sb2 = Math.sin(b/2.0);
distance = 2earth_radiusMath.asin(Math.sqrt(sa2 * sa2 + Math.cos(lat1)
* Math.cos(lat2)sb2sb2));
return distance;
}
rad(x){
return x*Math.PI/180.0;
}
// 適用于近距離
GetShortDistance(lon1,lat1,lon2,lat2) {
var DEF_PI = 3.14159265359; // PI
var DEF_2PI= 6.28318530712; // 2*PI
var DEF_PI180= 0.01745329252; // PI/180.0
var DEF_R =6370693.5; // radius of earth
var ew1,ns1,ew2,ns2;
var dx,dy,dew;
var distance;
// 角度轉(zhuǎn)換為弧度
ew1 = lon1 * DEF_PI180;
ns1 = lat1 * DEF_PI180;
ew2 = lon2 * DEF_PI180;
ns2 = lat2 * DEF_PI180;
// 經(jīng)度差
dew = ew1 - ew2;
// 若跨東經(jīng)和西經(jīng)180 度，進(jìn)行調(diào)整
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); // 南北方向長度(在經(jīng)度圈上的投影長度)
// 勾股定理求斜邊長
distance = Math.sqrt(dx * dx + dy * dy);
return distance;
}
// 適用于遠(yuǎn)距離
GetLongDistance(lon1,lat1,lon2,lat2)
{
var DEF_PI = 3.14159265359; // PI
var DEF_2PI= 6.28318530712; // 2*PI
var DEF_PI180= 0.01745329252; // PI/180.0
var DEF_R =6370693.5; // radius of earth
var ew1, ns1, ew2, ns2;
var distance;
// 角度轉(zhuǎn)換為弧度
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);
// 調(diào)整到[-1..1]范圍內(nèi)蜂桶，避免溢出
if (distance > 1.0)
distance = 1.0;
else if (distance < -1.0)
distance = -1.0;
// 求大圓劣弧長度
distance = DEF_R * Math.acos(distance);
return distance;
}
