已知地心地固坐标系的vx和vy，如何转换成地平坐标系下的速度，给我一个C++代码

在地理信息系统（GIS）中，从地球的地心地固坐标（如地心坐标系的水平分量vx和vy）转换到地平坐标系（如WGS84坐标系下的速度）通常涉及到经纬度的转换以及地球半径的影响。这涉及到了一些复杂的数学计算，包括球面三角学和投影变换。

以下是一个简化的C++代码示例，它假设了你已经有了一个转换函数用于经纬度的转换，这里仅展示基本步骤，实际应用中需要考虑精度和浮点数处理：

#include <cmath>

// 定义地球平均半径 (米)
const double EARTH_RADIUS = 6371e3;

// 地心地固坐标系的速度
double vx, vy;

// 转换为弧度
double phi_rad = atan2(vy, vx);
double lambda_rad = asin(sin(phi_rad) / cos(EarthRadius / distance));

// WGS84转换（这里仅作演示，实际应用需精确的经纬度转换）
double lat_deg = degrees(phi_rad); // 经度（纬度不变）
double lon_deg = degrees(lambda_rad);

// 如果需要转换速度方向，例如通过球面余弦定律计算偏航角
double bearing_rad = acos(cos(lat_rad_1) * cos(lat_rad_2) + sin(lat_rad_1) * sin(lat_rad_2) * cos(lon_deg_1 - lon_deg_2));
double bearing_degrees = degrees(bearing_rad);

// 返回地平坐标系的速度矢量（假设水平速度不变）
std::pair<double, double> velocity_on_plane(double speed_mps) {
    return {speed_mps * cos(bearing_degrees), speed_mps * sin(bearing_degrees)};
}

// 使用方法
std::pair<double, double> new_velocity = velocity_on_plane(std::hypot(vx, vy)); // 取水平速度