计算两个地理点的距离（以公里为单位）c#

我想计算两个地理点的距离。这些点以经度和纬度给出。

坐标是：

点 1：36.578581，-118.291994

点 2：36.23998，-116.83171

这是一个比较结果的网站：

http://www.movable-type.co.uk/scripts/latlong.html http://www.movable-type.co.uk/scripts/latlong.html

这是我从此链接中使用的代码：计算谷歌地图V3中两点之间的距离 https://stackoverflow.com/questions/1502590/calculate-distance-between-two-points-in-google-maps-v3

    const double PIx = Math.PI;
    const double RADIO = 6378.16;

    /// <summary>
    /// Convert degrees to Radians
    /// </summary>
    /// <param name="x">Degrees</param>
    /// <returns>The equivalent in radians</returns>
    public static double Radians(double x)
    {
        return x * PIx / 180;
    }

    /// <summary>
    /// Calculate the distance between two places.
    /// </summary>
    /// <param name="lon1"></param>
    /// <param name="lat1"></param>
    /// <param name="lon2"></param>
    /// <param name="lat2"></param>
    /// <returns></returns>
    public static double DistanceBetweenPlaces(double lon1, double lat1, double lon2, double lat2)
    {
        double R = 6371; // km
        double dLat = Radians(lat2 - lat1);
        double dLon = Radians(lon2 - lon1);
        lat1 = Radians(lat1);
        lat2 = Radians(lat2);

        double a = Math.Sin(dLat / 2) * Math.Sin(dLat / 2) + Math.Sin(dLon / 2) * Math.Sin(dLon / 2) * Math.Cos(lat1) * Math.Cos(lat2);
        double c = 2 * Math.Atan2(Math.Sqrt(a), Math.Sqrt(1 - a));
        double d = R * c;

        return d;
    }


Console.WriteLine(DistanceAlgorithm.DistanceBetweenPlaces(36.578581, -118.291994, 36.23998, -116.83171));

问题是我得到两个不同的结果。

我的结果：163,307 公里

网站结果：136公里

有什么建议么？？？

torti

您的公式几乎是正确的，但是您必须交换经度和纬度的参数

Console.WriteLine(DistanceAlgorithm.DistanceBetweenPlaces(-118.291994, 36.578581, -116.83171, 36.23998)); // = 136 km

我正在使用简化的公式：

// cos(d) = sin(φА)·sin(φB) + cos(φА)·cos(φB)·cos(λА − λB),
//  where φА, φB are latitudes and λА, λB are longitudes
// Distance = d * R
public static double DistanceBetweenPlaces(double lon1, double lat1, double lon2, double lat2)
{
    double R = 6371; // km

    double sLat1 = Math.Sin(Radians(lat1));
    double sLat2 = Math.Sin(Radians(lat2));
    double cLat1 = Math.Cos(Radians(lat1));
    double cLat2 = Math.Cos(Radians(lat2));
    double cLon = Math.Cos(Radians(lon1) - Radians(lon2));

    double cosD = sLat1*sLat2 + cLat1*cLat2*cLon;

    double d = Math.Acos(cosD);

    double dist = R * d;

    return dist;
}

Testing:

（赤道距离）：经度 0、100；纬度 = 0,0;地点之间的距离(0, 0, 100, 0) = 11119.5 公里

（北极距离）：经度 0、100；纬度 = 90,90；地点之间的距离(0, 90, 100, 90) = 0 公里

经度：-118.291994、-116.83171；纬度：36.578581, 36.23998 = 135.6 公里

经度：36.578581、36.23998；纬度：-118.291994，-116.83171 = 163.2 公里

此致

P.S. At web site http://www.movable-type.co.uk/scripts/latlong.html您用于结果比较，对于每个点，第一个文本框是纬度，第二个文本框是经度