GNSS xyz转经纬高 Qt代码

以下是一个简单的Qt代码示例，用于将GNSS坐标的XYZ值转换为经度、纬度和高度：

#include <QtMath>

// 定义WGS84椭球参数
const double a = 6378137.0; // 长半轴
const double b = 6356752.3142; // 短半轴
const double f = (a - b) / a; // 扁率
const double e_sq = f * (2 - f); // 第一偏心率平方

// XYZ转经纬高
void xyz2llh(double x, double y, double z, double& lat, double& lon, double& alt)
{
    double p = qSqrt(qPow(x, 2) + qPow(y, 2)); // 平面投影距离
    double theta = qAtan2(z * a, p * b); // 地心纬度
    double sin_theta = qSin(theta);
    double cos_theta = qCos(theta);
    lat = qAtan2(z + e_sq * b * qPow(sin_theta, 3), p - e_sq * a * qPow(cos_theta, 3)); // 大地纬度
    double N = a / qSqrt(1 - e_sq * qPow(qSin(lat), 2)); // 卯酉圈半径
    lon = qAtan2(y, x); // 经度
    alt = p / qCos(lat) - N; // 高度
}

该代码使用WGS84椭球参数计算大地坐标系经度、纬度和高度。输入参数为XYZ坐标，输出参数为经度、纬度和高度。

注意：该代码仅为示例代码，未考虑异常情况和精度问题，具体使用时请根据实际情况进行调整。