C++ 实现地理坐标系转直角坐标系

C++ 实现地理坐标系转直角坐标系需要进行以下步骤：

1. 引入数学库，如 `<cmath>`。因为需要进行一些三角函数的计算。

2. 定义基准点的经纬度坐标系坐标和直角坐标系坐标。例如：

   double base_lon = 116.403963; // 基准点经度
   double base_lat = 39.915119;  // 基准点纬度
   double base_x = 0;            // 基准点直角坐标系X轴坐标
   double base_y = 0;            // 基准点直角坐标系Y轴坐标

3. 定义需要转换的点的经纬度坐标系坐标。例如：

   double point_lon = 116.404269; // 需要转换点的经度
   double point_lat = 39.914714;  // 需要转换点的纬度
   double point_x = 0;            // 需要转换点的直角坐标系X轴坐标
   double point_y = 0;            // 需要转换点的直角坐标系Y轴坐标

4. 进行经纬度坐标系到直角坐标系的转换。其中，需要先将经纬度坐标系的经度和纬度转换为弧度，然后进行一些三角函数的计算。例如：

   double PI = 3.14159265358979323846; // 圆周率
   double EARTH_RADIUS = 6378137.0;    // 地球半径，单位为米

   double rad_lon = base_lon * PI / 180.0; // 将经度转换为弧度
   double rad_lat = base_lat * PI / 180.0; // 将纬度转换为弧度
   double ec = EARTH_RADIUS * cos(rad_lat); // 小圆的半径
   double ed = ec * cos(rad_lon); // 大圆的半径

   base_x = ed * sin(rad_lon); // 计算X轴坐标
   base_y = ec * sin(rad_lat); // 计算Y轴坐标

   rad_lon = point_lon * PI / 180.0;
   rad_lat = point_lat * PI / 180.0;
   ec = EARTH_RADIUS * cos(rad_lat);
   ed = ec * cos(rad_lon);

   point_x = ed * sin(rad_lon) - base_x;
   point_y = ec * sin(rad_lat) - base_y;

在上面的代码中，`PI` 表示圆周率，`EARTH_RADIUS` 表示地球半径（单位为米）。计算过程中，需要先计算出小圆的半径 `ec` 和大圆的半径 `ed`，然后根据公式计算出直角坐标系的坐标。

5. 最后，将计算出的直角坐标系的坐标赋值给相应的变量。例如：

   std::cout << "经纬度坐标系坐标：(" << point_lon << ", " << point_lat << ")" << std::endl;
   std::cout << "直角坐标系坐标：(" << point_x << ", " << point_y << ")" << std::endl;

完整代码如下：