用Visual Studio写一个大地坐标系转换为空间直角坐标系的代码

以下是C++语言的代码示例：

#include <iostream>
#include <cmath>

using namespace std;

// 定义常量
const double PI = 3.14159265358979323846;
const double a = 6378137.0;  // WGS84椭球长半轴
const double f = 1 / 298.257223563;  // WGS84椭球扁率

// 大地坐标系转空间直角坐标系函数
void geodetic2ecef(double B, double L, double H, double& X, double& Y, double& Z)
{
    // 弧度制转换为角度制
    double radB = B * PI / 180.0;
    double radL = L * PI / 180.0;

    // 计算第一偏心率的平方
    double e2 = 2 * f - f * f;

    // 计算卯酉圈曲率半径
    double N = a / sqrt(1 - e2 * sin(radB) * sin(radB));

    // 计算空间直角坐标系的坐标
    X = (N + H) * cos(radB) * cos(radL);
    Y = (N + H) * cos(radB) * sin(radL);
    Z = (N * (1 - e2) + H) * sin(radB);
}

int main()
{
    double B, L, H;  // 大地坐标系的纬度、经度、高程
    double X, Y, Z;  // 空间直角坐标系的坐标

    // 输入大地坐标系的坐标
    cout << "请输入大地坐标系的纬度（度分秒）：";
    cin >> B;
    cout << "请输入大地坐标系的经度（度分秒）：";
    cin >> L;
    cout << "请输入大地坐标系的高程（米）：";
    cin >> H;

    // 转换为十进制度数
    B = B / 3600.0;
    L = L / 3600.0;

    // 调用函数计算空间直角坐标系的坐标
    geodetic2ecef(B, L, H, X, Y, Z);

    // 输出结果
    cout << "空间直角坐标系的坐标为：" << endl;
    cout << "X = " << X << " m" << endl;
    cout << "Y = " << Y << " m" << endl;
    cout << "Z = " << Z << " m" << endl;

    return 0;
}

在上面的代码中，我们定义了一个 `geodetic2ecef` 函数，用来将大地坐标系的坐标转换为空间直角坐标系的坐标。该函数接受四个参数：大地纬度 B、大地经度 L、大地高程 H，以及空间直角坐标系的坐标 X、Y、Z。在函数中，我们首先将大地纬度和大地经度转换为弧度制，然后根据 WGS84 椭球的参数计算卯酉圈曲率半径 N，最后根据公式计算空间直角坐标系的坐标。

在 `main` 函数中，我们首先输入大地坐标系的坐标，然后调用 `geodetic2ecef` 函数计算空间直角坐标系的坐标，并输出结果。注意，在输入大地坐标系的纬度、经度和高程时，我们将其转换为十进制度数。