写个大地坐标系转空间直角坐标系的Python代码

以下是将大地坐标系转换为空间直角坐标系的Python代码：

import math

def geodetic_to_ecef(lat, lon, alt):
    # WGS84椭球体参数
    a = 6378137.0
    b = 6356752.314245
    f = (a - b) / a
    e_sq = f * (2-f)

    # 转换为弧度
    lat_rad = math.radians(lat)
    lon_rad = math.radians(lon)

    # 计算参数
    N = a / math.sqrt(1 - e_sq * math.sin(lat_rad)**2)
    x = (N + alt) * math.cos(lat_rad) * math.cos(lon_rad)
    y = (N + alt) * math.cos(lat_rad) * math.sin(lon_rad)
    z = (N * (1 - e_sq) + alt) * math.sin(lat_rad)

    return x, y, z

其中，输入参数 `lat`, `lon`, `alt` 分别代表大地坐标系下的纬度、经度和高程，输出结果为空间直角坐标系下的 `x`, `y`, `z` 坐标。