给定两个经纬度，求相对位方位算法 c++

在C++中，计算两个经纬度之间的相对方位通常涉及到地理编码和角度计算。一个常见的做法是使用Haversine公式来确定两点之间的大圆距离，然后通过地图投影（如Mercator、UTM等）转换为平面坐标上角度差。这里是一个简单的示例：

#include <cmath>

// Haversine公式计算弧度距离
const double earthRadius = 6371; // 地球平均半径（公里）

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

double calculateDistance(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 earthRadius * c;
}

// 根据距离判断方位
std::string getDirection(double distance) {
    if (distance < 111.32 * earthRadius / 60) { // 约1°
        return "几乎正北";
    } else if (distance > 111.32 * earthRadius) { // 约180°
        return "几乎正南";
    } else {
        return "东北、东南、西北或西南方向"; // 余下的大致方位范围
    }
}

int main() {
    double lat1 = 39.9042; // 上海经度
    double lon1 = 121.4737;
    double lat2 = 40.7128; // 纽约经度
    double lon2 = -74.0060;

    double distance = calculateDistance(lat1, lon1, lat2, lon2);
    std::cout << "两地距离: " << distance << " km" << std::endl;
    std::cout << "方位: " << getDirection(distance) << std::endl;

    return 0;
}

请注意，这只是一个简化的例子，实际应用中可能会需要更精确的地图投影和考虑地球曲率的影响。此外，计算角度差和方位时可能需要根据需求调整精度。