大地坐标和空间直角坐标转换代码

以下是大地坐标和空间直角坐标间相互转换的Python代码：

import math

def geodetic_to_cartesian(latitude, longitude, height):
    a = 6378137.0
    b = 6356752.314245
    f = (a - b) / a
    e_sq = f * (2-f)
    N = a / math.sqrt(1 - e_sq * math.sin(latitude) ** 2)
    x = (N + height) * math.cos(latitude) * math.cos(longitude)
    y = (N + height) * math.cos(latitude) * math.sin(longitude)
    z = (N * (1 - e_sq) + height) * math.sin(latitude)
    return x, y, z

def cartesian_to_geodetic(x, y, z):
    a = 6378137.0
    b = 6356752.314245
    f = (a - b) / a
    e_sq = f * (2-f)
    p = math.sqrt(x ** 2 + y ** 2)
    theta = math.atan2(z * a, p * b)
    latitude = math.atan2(z + e_sq * b * math.sin(theta) ** 3, p - e_sq * a * math.cos(theta) ** 3)
    longitude = math.atan2(y, x)
    N = a / math.sqrt(1 - e_sq * math.sin(latitude) ** 2)
    height = p / math.cos(latitude) - N
    return latitude, longitude, height

其中，```geodetic_to_cartesian```函数接受一个纬度、经度和高度，返回一个空间直角坐标系下的x、y、z三个坐标值。```cartesian_to_geodetic```函数接受一个空间直角坐标系下的x、y、z三个坐标值，返回一个大地坐标系下的纬度、经度和高度。注意输入的角度应该是弧度制。