C# 经纬度计算

#region 计算经纬度距离
public const double EARTH_RADIUS = 6377.830;//地球半径(千米)

/// <summary>
/// 角度换算成弧度
/// </summary>
/// <param name="d">角度</param>
/// <returns>弧度</returns>
public static double rad(double d)
{
return d * Math.PI / 180.0;
}
/// <summary>
/// 计算距离
/// </summary>
/// <param name="lat1">纬度1</param>
/// <param name="lng1">经度1</param>
/// <param name="lat2">纬度2</param>
/// <param name="lng2">经度2</param>
/// <returns>距离（千米）</returns>
public static double GetDistance(double lat1, double lng1, double lat2, double lng2)
{
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)));
s = s * EARTH_RADIUS;
s = Math.Round(s * 10000) / 10000;
return s;
}
#endregion