最近做一個項目挪蹭，需要計算計算地球表面兩點間距離超燃，然后就用JS實現(xiàn)了一下郭怪。
翻看了下資料計算地球表面兩點間的距離大概有兩種辦法。
第一種是默認(rèn)地球是一個光滑的球面宽档，然后計算任意兩點間的距離，這個距離叫做大圓距離(The Great Circle Distance)并级。
公式如下：
使用JS來實現(xiàn)為：
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;
}
這個公式在大多數(shù)情況下比較正確巴柿，只有在處理球面上的相對點的時候，會出現(xiàn)問題死遭，有一個修正的公式，因為沒有需要凯旋，就沒有找出來呀潭，可以在wiki上查到。
當(dāng)然至非，我們都知道钠署，地球其實并不是一個真正的圓球體，而是橢球荒椭，所以有了下面的公式：
/**
* 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));
}
這個公式計算出的結(jié)果要比第一個好一些谐鼎，當(dāng)然，最后結(jié)果的精度實際上還取決于傳入的坐標(biāo)的精度趣惠。
