获取两点（经纬度表示）间距离【转】

转自：http://www.cnblogs.com/3echo/archive/2010/04/28/1722982.html

经常有朋友或同事征询有关如何获取两点（经纬度表示）间距离的问题，现将解决方法写出来，以供自己和朋友查询、学习。

解决此问题的关键是要理解空间几何模型，在理解空间几何模型的基础上再利用数学公式求取两点之间的值。其解决步骤如下：

1、设两点分别为P1、P2，如果其值是用度分秒形式表示，则需将其转换成十进制度的形式，如P1点纬度为23度30分，则其纬度值转换成十进制度的形式为23.5度。如果值为十进制度的形式，则直接进入第二步。

2、分别将两点的经度、纬度值转换成弧度制形式，如P1纬度为23.5度，转换成弧度制则为：23.5*PI / 180。分别用 P1latInRad、P1LongInRad、P2latInRad、P2LongInRad表示。

3、分别求取两点间的纬度差（dlat）与经度差（dlon）；

4、求取两点间的正弦与余弦值，公式如下：

A=sin²(dlat/2) + cos(P1LatInRad)*cos(P2LatInRad)*Sin²(dlon/2)       (1)

5、求取两点的正切值，公式如下：

C=2*Math.Atan2(Math.Sqrt(A), Math.Sqrt(1-A))                     (2)

6、返回两点间的距离：公式如下：

D=EarthRadiusKm * C                                          (3)

 

附C#源代码：

 

 1       /// <summary>
 2       /// 获取两点(经纬度表示)间的距离
 3       /// </summary>
 4       /// <param name="p1Lat">第一点纬度值</param>
 5       /// <param name="p1Lng">第一点经度值</param>
 6       /// <param name="p2Lat">第二点纬度值</param>
 7       /// <param name="p2Lng">第二点经度值</param>
 8       /// <returns>返回两点间距离</returns>
 9       public double GetDistance(double p1Lat, double p1Lng, double p2Lat, double p2Lng)
10       {
11          double dLat1InRad = p1Lat * (Math.PI / 180);
12          double dLong1InRad = p1Lng * (Math.PI / 180);
13          double dLat2InRad = p2Lat * (Math.PI / 180);
14          double dLong2InRad = p2Lng * (Math.PI / 180);
15          double dLongitude = dLong2InRad - dLong1InRad;
16          double dLatitude = dLat2InRad - dLat1InRad;
17          double a = Math.Pow(Math.Sin(dLatitude / 2), 2) + Math.Cos(dLat1InRad) * Math.Cos(dLat2InRad) * Math.Pow(Math.Sin(dLongitude / 2), 2);
18          double c = 2 * Math.Atan2(Math.Sqrt(a), Math.Sqrt(1 - a));
19          double dDistance = EarthRadiusKm * c;
20          return dDistance;
21       }
22 
23       /// <summary>
24       /// Radius of the Earth
25       /// </summary>
26       public double EarthRadiusKm = 6378.137; // WGS-84