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,
    }

Shenzhen Universiade - Shenzhen 2011 Summer Universiade