python 空间直角坐标系转化为地心坐标系

可以通过以下代码实现：

import math

# 空间直角坐标系坐标值
x = 1
y = 2
z = 3

# 定义WGS84椭球体参数
a = 6378137
b = 6356752.314245

# 计算椭球面的第一偏心率平方值
e2 = math.pow(a, 2) - math.pow(b, 2)
e2 = e2 / math.pow(a, 2)

# 计算空间直角坐标系到地心坐标系的转换矩阵
p = math.sqrt(math.pow(x, 2) + math.pow(y, 2))
cos_theta = x / p
sin_theta = y / p
cos_phi = p / math.sqrt(math.pow(p, 2) + math.pow(z, 2))
sin_phi = z / math.sqrt(math.pow(p, 2) + math.pow(z, 2))
x1 = (a / math.sqrt(1 - e2 * math.pow(sin_phi, 2))) * cos_phi * cos_theta
y1 = (a / math.sqrt(1 - e2 * math.pow(sin_phi, 2))) * cos_phi * sin_theta
z1 = (a * (1 - e2) / math.sqrt(1 - e2 * math.pow(sin_phi, 2))) * sin_phi

# 输出结果
print("地心坐标系坐标值为：")
print("x: ", x1)
print("y: ", y1)
print("z: ", z1)