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

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

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

大圆路径距离

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

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

d = 2 R s i n − 1 s i n 2 θ A − θ B 2 + c o s θ A c o s θ B s i n 2 λ A − λ B 2 d = 2Rsin^{-1}\sqrt{sin^2\frac{\theta_A - \theta_B}{2} + cos\theta_Acos\theta_Bsin^2\frac{\lambda_A-\lambda_B}{2}} d=2Rsin1sin22θAθB+cosθAcosθBsin22λAλB

θ A \theta_A θA 为A点的纬度
θ B \theta_B θB 为B点的纬度
λ A \lambda_A λA 为A点的经度
λ B \lambda_B λB 为B点的经度
s i n − 1 sin^{-1} sin1 为反正弦
s i n 2 sin^2 sin2 求出正弦值后再求平方
R R 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;
}

你可能感兴趣的:(C++,QT,c++)