gps坐标转大地直角坐标代码c++

以下是C++实现的GPS坐标转大地直角坐标的代码：

#include <cmath>

const double a = 6378137.0; // 长半轴
const double b = 6356752.314245; // 短半轴
const double e = sqrt(1.0 - (b / a) * (b / a)); // 第一偏心率

// 将角度转换为弧度
double rad(double d)
{
    return d * 3.141592653589793 / 180.0;
}

// GPS坐标转大地直角坐标
void gps_to_xyz(double lat, double lng, double alt, double &x, double &y, double &z)
{
    double N = a / sqrt(1 - e * e * sin(rad(lat)) * sin(rad(lat)));
    x = (N + alt) * cos(rad(lat)) * cos(rad(lng));
    y = (N + alt) * cos(rad(lat)) * sin(rad(lng));
    z = (N * (1 - e * e) + alt) * sin(rad(lat));
}

其中，`lat`、`lng`和`alt`分别代表纬度、经度和海拔高度，`x`、`y`和`z`为大地直角坐标系下的坐标。

使用时，只需要调用`gps_to_xyz`函数，传入GPS坐标和空间坐标变量即可。例如：

double lat = 31.2304, lng = 121.4737, alt = 10.0;
double x, y, z;
gps_to_xyz(lat, lng, alt, x, y, z);