JS-根据经纬度计算地球上任意两点之间的距离

   最近做一个项目，需要计算计算地球表面两点间距离，然后就用JS实现了一下。

翻看了下资料计算地球表面两点间的距离大概有两种办法。
第一种是默认地球是一个光滑的球面，然后计算任意两点间的距离，这个距离叫做大圆距离(The Great Circle Distance)。
公式如下：

使用JS来实现为：
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;

}

这个公式在大多数情况下比较正确，只有在处理球面上的相对点的时候，会出现问题，有一个修正的公式，因为没有需要，就没有找出来，可以在wiki上查到。
当然，我们都知道，地球其实并不是一个真正的圆球体，而是椭球，所以有了下面的公式：
/**

 * 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));

}

这个公式计算出的结果要比第一个好一些，当然，最后结果的精度实际上还取决于传入的坐标的精度。