用c语言编写大地坐标系和直角坐标系的转化

下面是一个简单的C语言程序，用于将大地坐标系的经纬度转换为直角坐标系的X、Y、Z坐标。其中，WGS84椭球体参数被用作参考。

#include <stdio.h>
#include <math.h>

#define a 6378137.0 //WGS84椭球体长半轴
#define b 6356752.314245 //WGS84椭球体短半轴
#define f (1.0/298.257223563) //WGS84椭球体扁率
#define e2 (2*f-f*f) //WGS84椭球体第一偏心率的平方

double rad(double d) {
    return d * M_PI / 180.0;
}

void geodetic2rectangular(double lat, double lon, double height, double *x, double *y, double *z) {
    double N, X, Y, Z;

    //转换为弧度
    lat = rad(lat);
    lon = rad(lon);

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

    //计算直角坐标系坐标
    X = (N + height) * cos(lat) * cos(lon);
    Y = (N + height) * cos(lat) * sin(lon);
    Z = (N * (1 - e2) + height) * sin(lat);

    //返回坐标值
    *x = X;
    *y = Y;
    *z = Z;
}

int main() {
    double lat = 39.9042; //纬度
    double lon = 116.4074; //经度
    double height = 0; //高度
    double x, y, z;

    geodetic2rectangular(lat, lon, height, &x, &y, &z);

    printf("X = %.3f\nY = %.3f\nZ = %.3f\n", x, y, z);

    return 0;
}

在此示例中，我们使用WGS84椭球体长半轴a、短半轴b和扁率f来计算大地坐标系中的卯酉圈曲率半径N。然后，我们可以使用以下公式将大地坐标系转换为直角坐标系：

X = (N + height) * cos(lat) * cos(lon)
Y = (N + height) * cos(lat) * sin(lon)
Z = (N * (1 - e2) + height) * sin(lat)

最后，我们将计算出的X、Y和Z坐标打印到控制台上。