1

如果我有两个 GPS 位置，例如 51.507222、-0.1275 和 48.856667、2.350833，我可以使用什么公式来计算两者之间的距离？我听说过很多关于半正弦公式的信息，但找不到任何关于它的信息，或者如何将它应用到 C 中。

我已经编写了以下代码，但是，它非常不准确。有人知道为什么吗？我想不通。问题来自函数本身，但我不知道它是什么。

float calcDistance(float A, float B, float C, float D) 
{
    float dLat;
    float dLon;
    dLat = (C - A);
    dLon = (D - B);
    dLat /= 57.29577951;
    dLon /= 57.29577951; 
    float v_a;
    float v_c;
    float distance;

    v_a = sin(dLat/2) * sin(dLat/2) + cos(A) * cos(C) * sin(dLon/2) * sin(dLon/2);
    v_c = 2 * atan2(sqrt(v_a),sqrt(1-v_a));
    distance = r * v_c;
    return distance;
}
4

4 回答 4

4

你的代码是绝对正确的。

我会将参数A、和重命名为B、和。 CDlat1lon1lat2lon2

它以公里为单位返回距离。您需要定义r为 6371，大致为地球的半径。

于 2011-05-04T11:36:37.907 回答
2

您正在寻找的是大圆距。

于 2011-04-29T23:56:19.320 回答
1

这个页面有一个 Javascript 块，我相信你可以很容易地在 C 中适应它：

var R = 6371e3; // metres
var φ1 = lat1.toRadians();
var φ2 = lat2.toRadians();
var Δφ = (lat2-lat1).toRadians();
var Δλ = (lon2-lon1).toRadians();

var a = Math.sin(Δφ/2) * Math.sin(Δφ/2) +
        Math.cos(φ1) * Math.cos(φ2) *
        Math.sin(Δλ/2) * Math.sin(Δλ/2);
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));

var d = R * c;
于 2011-04-30T00:08:04.833 回答
1

也许为时已晚，但这就是计算不准确的原因：-> A 和 C（lat1 和 lat2）也必须从度数转换为弧度

dLat /= 57.295779513082325;
dLon /= 57.295779513082325;
lat1 /= 57.295779513082325;
lat2 /= 57.295779513082325;

问候， 欧米尔

于 2017-01-16T22:19:56.280 回答