通过经纬度计算地球两点间的距离

设两点A、B的经、纬度分别为(jA,wA)(jB,wB)，则半径为R的球面上两点间的最短距离(大圆弧)为：

弧AB=R*arccos[sin(wA)sin(wB)+cos(wA)cos(wB)*cos(jA-jB)] .

// 赤道半径

        private const double EARTH_RADIUS = 6378137.0;   

        // 平均半径
        //private const double EARTH_RADIUS = 6371004.0;

        // 把角度变换为弧度

        private double rad(double d)
        {
            return d * Math.PI / 180.0;
        }

        ///

        /// 通过经纬度计算地球两点间的距离
        ///

        /// A点纬度
        /// A点经度
        /// B点纬度
        /// B点经度
        /// 返回两点间的距离，单位公里
        public double GetDistance(double lat1, double lng1, double lat2, double lng2)
        {
            double radLat1 = rad(lat1);
            double radLat2 = rad(lat2);
            double b = rad(lng1) - rad(lng2);
            double d = Math.Acos(Math.Cos(radLat1) * Math.Cos(radLat2) * Math.Cos(b) + Math.Sin(radLat1) * Math.Sin(radLat2));
            return d * EARTH_RADIUS; ;
        }

地图基础知识
地图应用接口中使用的经纬度单位是NTU，您需要把其它单位的经纬度值转化过来，以下为常用的转换关系：
基本转换：
NTU = 度*100000
度 = NTU/100000
例如：
经度 = 116.21345°
纬度 = 39.445875°
NTU经度：116.21345 * 100000 = 11621345（NTU）
NTU纬度：39.445874*100000=3944587(NTU)

度分转换：
将度分单位数据转换为度单位数据
度=度+分/60
例如：
经度 = 116°20.12’
纬度 = 39°12.34’
经度 = 116 + 20.12 / 60 = 116.33533°
纬度 = 39 + 12.34 / 60 = 39.20567°
NTU经度 = 116.33533 * 100000 = 11633533（NTU）
NTU纬度 = 39.20567 * 100000 = 3920567（NTU）

度分秒转换：
将度分秒单位数据转换为度单位数据
度 = 度 + 分 / 60 + 秒 / 60 / 60
例如：
经度 = 116°20’43”
纬度 = 39°12’37”
经度 = 116 + 20 / 60 + 43 / 60 / 60 = 116.34528°
纬度 = 39 + 12 / 60 + 37 / 60 / 60 = 39.21028°
NTU经度 = 116.34528 * 100000 = 11634528（NTU）
NTU纬度 = 39.21028 * 100000 = 3921028（NTU）
其格式的经纬值先转换度，再采用上面的运算。
与经纬度有关的基础知识：
NTU 相当于十万分之一度。
如何计算经度方向距离
比如: 点A的纬度为3995400，点B的纬度为3995300，则这两个点在沿着纬线的方向相差100米。
点A的经度为11695400，点B的经度为11695300，则这两个点在沿着经线的方向相差大约为77米,
该值的计算方法可以用近似公式：
经度方向距离 = 经度差 * cos(纬度值) = 100 * cos(39) = 77米。

如何计算两点间距离
比如：点A的经度为11695400，纬度为3995400。点B的经度为11695300，纬度为3995300。
公式:两点间距离 = [ (A点经度 - B点经度)^2 + (A点纬度 - B点纬度)^2 ] ^ (1/2) = [ (11695400 - 11695300)^2
+ (3995400 - 3995300)^2 ] ^(1/2) =(10000+10000) ^ (1/2) =141米.

公式说明:加法 + ，减法 - ，乘法 *，除法 /，幂运算 ^(1/2)表示平方根 ^2表示平方。
适用范围：未加密经纬度/加密经纬度