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

以下是一个简单的示例代码，可以将空间直角坐标系（X，Y，Z）转换为大地坐标系（纬度，经度，高度）：

#include <iostream>
#include <cmath>

using namespace std;

const double PI = 3.14159265358979323846;
const double a = 6378137.0; // 地球长半轴
const double b = 6356752.314245; // 地球短半轴
const double f = (a - b) / a; // 扁率

double deg2rad(double deg) {
    return deg * PI / 180.0;
}

double rad2deg(double rad) {
    return rad * 180.0 / PI;
}

void cartesian2geodetic(double x, double y, double z, double &lat, double &lon, double &h) {
    double p = sqrt(x * x + y * y);
    double theta = atan2(z * a, p * b);

    lat = atan2(z + f * f * b * pow(sin(theta), 3), p - f * a * pow(cos(theta), 3));
    lon = atan2(y, x);
    double N = a / sqrt(1 - f * f * pow(sin(lat), 2));
    h = p / cos(lat) - N;

    lat = rad2deg(lat);
    lon = rad2deg(lon);
}

int main() {
    double x = 500000;
    double y = 500000;
    double z = 0;
    double lat, lon, h;

    cartesian2geodetic(x, y, z, lat, lon, h);

    cout << "X: " << x << endl;
    cout << "Y: " << y << endl;
    cout << "Z: " << z << endl;
    cout << "Latitude: " << lat << endl;
    cout << "Longitude: " << lon << endl;
    cout << "Height: " << h << endl;

    return 0;
}

该代码使用了椭球体参数WGS84，可以将地球看作一个旋转椭球体。该代码的输入为空间直角坐标系（X，Y，Z），输出为大地坐标系（纬度，经度，高度）。