GIS measure 测距公式和C#,JavaScript的实现(转)

单位：m

Util.distanceByLnglat = function(lng1,lat1,lng2,lat2)
{
    var radLat1 = Util.Rad(lat1);
    var radLat2 = Util.Rad(lat2);
    var a = radLat1 - radLat2;
    var b = Util.Rad(lng1) - Util.Rad(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 * 6378137.0;
    s = Math.round(s * 10000) / 10000;
    return s;
}

Util.Rad = function(d)
{
    return d * Math.PI / 180.0;
}

note
Math.sqrt(x)    :X的平方根
Math.pow(x, n)    :返回x的n次方的数值
Math.cos(x)    :X的余弦函数
Math.sin(x)    :X的正弦函数
Math.asin(y)    :X的反正弦函数
Math.round(x)    :取最接近整数x的数值
6378137        :长半径

坐标系
WGS—84世界大地坐标系(Wor1d Geodetic System)。属于地心坐标系

从Google Map上弄来的根据经纬度求地球表面两点间距离的实现， 稍微改编了一下，对于我国境内空间距离计算，该实现已经够用，以米为单位。.Net2.0，C#实现。

       public static double DistanceOfTwoPoints(double lng1,double lat1,  double lng2, double lat2, GaussSphere gs)
        {          
            double radLat1 = Rad(lat1);
            double radLat2 = Rad(lat2);
            double a = radLat1 - radLat2;
            double b = Rad(lng1) - Rad(lng2);
       double 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))); 字串8
            s = s * (gs == GaussSphere.WGS84 ? 6378137.0 : (gs == GaussSphere.Xian80 ? 6378140.0 : 6378245.0));
            s = Math.Round(s * 10000) / 10000;
            return s;
        }
        private static double Rad(double d)
        {
            return d * Math.PI / 180.0;
        }

    GaussSphere 为自定义枚举类型
    /// <summary>
    /// 高斯投影中所选用的参考椭球
    /// </summary>
    public enum GaussSphere
    {
        Beijing54,
        Xian80,
        WGS84,
    }