根据两个gps经纬度计算与正北的角度 C语言

在C语言中，要计算两个GPS经纬度之间的角度（通常表示为方向角或方位角），你可以使用Haversine公式来估算两点间的经度差对结果的影响。这里是一个简单的步骤：

1. 首先，将经纬度从度数转换到弧度，因为Haversine公式需要的是弧度值。
2. 计算两个点之间的纬度差（y轴）和经度差（x轴）。
3. 使用反正切函数（atan2函数）计算出角度，该函数会返回一个介于-π和π之间的角度，代表两者之间的方向。

以下是一个基本的C代码示例：

#include <math.h>

double toRadians(double degrees) {
    return degrees * M_PI / 180.0;
}

double calculateAzimuth(double lat1, double lon1, double lat2, double lon2) {
    double dLat = toRadians(lat2 - lat1);
    double dLon = toRadians(lon2 - lon1);
    double a = sin(dLat / 2) * sin(dLat / 2) + cos(toRadians(lat1)) * cos(toRadians(lat2)) * sin(dLon / 2) * sin(dLon / 2);
    double c = 2 * atan2(sqrt(a), sqrt(1 - a));
    return (c * 180) / M_PI; // 返回角度范围在0-360之间
}

int main() {
    double lat1 = 39.9042; // 经纬度1
    double lon1 = 116.4074; // 经纬度2
    double lat2 = 51.5074; // 另一坐标点的经纬度
    double lon2 = -0.1278; // 另一坐标点的经纬度

    double azimuth = calculateAzimuth(lat1, lon1, lat2, lon2);

    printf("两地之间的角度是: %f°\n", azimuth);
    return 0;
}