2023-3-13 经纬度距离计算公式

ihmhm12345

已于 2023-03-14 11:00:36 修改

阅读量389

点赞数

分类专栏：Ｃ＋＋ QT 文章标签： c++

于 2023-03-13 23:43:28 首次发布

本文链接：https://blog.csdn.net/ihmhm12345/article/details/129511104

版权

QT 同时被 2 个专栏收录

38 篇文章 6 订阅

订阅专栏

Ｃ＋＋

12 篇文章 1 订阅

订阅专栏

在网上找了好久才找到正确的经纬度距离计算公式，而且是较高精度的计算

公式原自：https://docs.gmt-china.org/latest/option/distcal/

大圆路径距离

该方法将地球近似为一个半径为R的球，地球上任意两点 A 和 B 的大圆路径距离可以用 Haversine 公式计算：

该方法是 GMT 默认使用的距离计算方法，适用于大多数情况。

$2Rsin^{-1}\sqrt{sin^2\frac{\theta_A - \theta_B}{2} + cos\theta_Acos\theta_Bsin^2\frac{\lambda_A-\lambda_B}{2}}$

$\theta_A$ 为A点的纬度
$\theta_B$ 为B点的纬度
$\lambda_A$ 为A点的经度
$\lambda_B$ 为B点的经度
$sin^{-1}$ 为反正弦
$sin^2$ 求出正弦值后再求平方
$R$ 为地球半径

以下代码来自Qt QGeoCoorditude 的代码，与上述公式吻合

static const double qgeocoordinate_EARTH_MEAN_RADIUS = 6371.0072;
double distance(double lat1, double lng1, double lat2, double lng2)
{
    //    // Haversine formula
    double dlat = qDegreesToRadians(lat1 - lat2);
    double dlon = qDegreesToRadians(lng1 - lng2);
    double haversine_dlat = sin(dlat / 2.0);
    haversine_dlat *= haversine_dlat;
    double haversine_dlon = sin(dlon / 2.0);
    haversine_dlon *= haversine_dlon;
    double y = haversine_dlat
               + cos(qDegreesToRadians(lat1)) * cos(qDegreesToRadians(lat2)) * haversine_dlon;
    double x = 2 * asin(sqrt(y));
    return qreal(x * qgeocoordinate_EARTH_MEAN_RADIUS * 1000);
    return 0.0;
}