快捷導(dǎo)航

根據(jù)經(jīng)緯度計(jì)算地球上兩點(diǎn)之間的距離js實(shí)現(xiàn)代碼

更新時(shí)間：2013年03月05日 14:13:45 作者：

最近用到了根據(jù)經(jīng)緯度計(jì)算地球表面兩點(diǎn)間距離的公式，然后就用JS實(shí)現(xiàn)了一下，根據(jù)地球是一個(gè)光滑的球面以及并不是一個(gè)真正的圓球體，而是橢球延伸出兩種方法，感興趣的你可以參考下

利用JS實(shí)現(xiàn)的根據(jù)經(jīng)緯度計(jì)算地球上兩點(diǎn)之間的距離
最近用到了根據(jù)經(jīng)緯度計(jì)算地球表面兩點(diǎn)間距離的公式，然后就用JS實(shí)現(xiàn)了一下。
計(jì)算地球表面兩點(diǎn)間的距離大概有兩種辦法。

第一種是默認(rèn)地球是一個(gè)光滑的球面，然后計(jì)算任意兩點(diǎn)間的距離，這個(gè)距離叫做大圓距離(The Great Circle Distance)。
公式如下：
使用JS來(lái)實(shí)現(xiàn)為：

復(fù)制代碼代碼如下:

 
var EARTH_RADIUS = 6378137.0; //單位M 
var PI = Math.PI; 

function getRad(d){ 
return d*PI/180.0; 
} 

/** 
* caculate the great circle distance 
* @param {Object} lat1 
* @param {Object} lng1 
* @param {Object} lat2 
* @param {Object} lng2 
*/ 
function getGreatCircleDistance(lat1,lng1,lat2,lng2){ 
var radLat1 = getRad(lat1); 
var radLat2 = getRad(lat2); 

var a = radLat1 - radLat2; 
var b = getRad(lng1) - getRad(lng2); 

var s = 2*Math.asin(Math.sqrt(Math.pow(Math.sin(a/2),2) + Math.cos(radLat1)*Math.cos(radLat2)*Math.pow(Math.sin(b/2),2))); 
s = s*EARTH_RADIUS; 
s = Math.round(s*10000)/10000.0; 

return s; 
} 

這個(gè)公式在大多數(shù)情況下比較正確，只有在處理球面上的相對(duì)點(diǎn)的時(shí)候，會(huì)出現(xiàn)問(wèn)題，有一個(gè)修正的公式，因?yàn)闆](méi)有需要，就沒(méi)有找出來(lái)，可以在wiki上查到。
當(dāng)然，我們都知道，地球其實(shí)并不是一個(gè)真正的圓球體，而是橢球，所以有了下面的公式：

復(fù)制代碼代碼如下:

 
/** 
* approx distance between two points on earth ellipsoid 
* @param {Object} lat1 
* @param {Object} lng1 
* @param {Object} lat2 
* @param {Object} lng2 
*/ 
function getFlatternDistance(lat1,lng1,lat2,lng2){ 
var f = getRad((lat1 + lat2)/2); 
var g = getRad((lat1 - lat2)/2); 
var l = getRad((lng1 - lng2)/2); 

var sg = Math.sin(g); 
var sl = Math.sin(l); 
var sf = Math.sin(f); 

var s,c,w,r,d,h1,h2; 
var a = EARTH_RADIUS; 
var fl = 1/298.257; 

sg = sg*sg; 
sl = sl*sl; 
sf = sf*sf; 

s = sg*(1-sl) + (1-sf)*sl; 
c = (1-sg)*(1-sl) + sf*sl; 

w = Math.atan(Math.sqrt(s/c)); 
r = Math.sqrt(s*c)/w; 
d = 2*w*a; 
h1 = (3*r -1)/2/c; 
h2 = (3*r +1)/2/s; 

return d*(1 + fl*(h1*sf*(1-sg) - h2*(1-sf)*sg)); 
} 