利用VS软件做一个大地坐标与空间直角坐标相互转换算例（包括克拉索夫斯基椭球、IUGG11975椭球、CGCS2000椭球）的程序

以下是一个简单的示例程序，可以将大地坐标转换为空间直角坐标，支持克拉索夫斯基椭球、IUGG11975椭球和CGCS2000椭球。

#include <iostream>
#include <cmath>
using namespace std;

// 定义椭球参数
const double a_Kras = 6378245.0; // 克拉索夫斯基椭球长半轴
const double f_Kras = 1 / 298.3; // 克拉索夫斯基椭球扁率
const double e2_Kras = 2 * f_Kras - pow(f_Kras, 2); // 克拉索夫斯基椭球第一偏心率的平方
const double a_IUGG = 6378140.0; // IUGG1975椭球长半轴
const double f_IUGG = 1 / 298.257; // IUGG1975椭球扁率
const double e2_IUGG = 2 * f_IUGG - pow(f_IUGG, 2); // IUGG1975椭球第一偏心率的平方
const double a_CGCS = 6378137.0; // CGCS2000椭球长半轴
const double f_CGCS = 1 / 298.257222101; // CGCS2000椭球扁率
const double e2_CGCS = 2 * f_CGCS - pow(f_CGCS, 2); // CGCS2000椭球第一偏心率的平方

// 定义弧度转度数的函数
double rad2deg(double rad) {
    return rad * 180 / M_PI;
}

// 定义度数转弧度的函数
double deg2rad(double deg) {
    return deg * M_PI / 180;
}

// 定义大地坐标转空间直角坐标的函数
void geodetic2cartesian(double B, double L, double H, double a, double e2, double &X, double &Y, double &Z) {
    double N = a / sqrt(1 - e2 * pow(sin(deg2rad(B)), 2)); // 卯酉圈曲率半径
    X = (N + H) * cos(deg2rad(B)) * cos(deg2rad(L));
    Y = (N + H) * cos(deg2rad(B)) * sin(deg2rad(L));
    Z = (N * (1 - e2) + H) * sin(deg2rad(B));
}

int main() {
    // 输入大地坐标和椭球参数
    double B, L, H, a, e2;
    cout << "请输入纬度B（度）：" << endl;
    cin >> B;
    cout << "请输入经度L（度）：" << endl;
    cin >> L;
    cout << "请输入高程H（米）：" << endl;
    cin >> H;
    cout << "请选择椭球：" << endl;
    cout << "1. 克拉索夫斯基椭球" << endl;
    cout << "2. IUGG1975椭球" << endl;
    cout << "3. CGCS2000椭球" << endl;
    int choice;
    cin >> choice;
    switch (choice) {
        case 1:
            a = a_Kras;
            e2 = e2_Kras;
            break;
        case 2:
            a = a_IUGG;
            e2 = e2_IUGG;
            break;
        case 3:
            a = a_CGCS;
            e2 = e2_CGCS;
            break;
        default:
            cout << "输入错误！" << endl;
            return 0;
    }

    // 转换为空间直角坐标
    double X, Y, Z;
    geodetic2cartesian(B, L, H, a, e2, X, Y, Z);

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

    return 0;
}

这个程序可以根据用户输入的大地坐标和椭球参数，计算出相应的空间直角坐标。用户可以选择使用克拉索夫斯基椭球、IUGG1975椭球或CGCS2000椭球进行计算。