大地坐标转为大地空间直角坐标 python

可以使用 pyproj 库来实现大地坐标（经纬度）到大地空间直角坐标的转换。下面是一个示例代码：

import pyproj

# 定义大地坐标系和大地空间直角坐标系的投影方式
geod = pyproj.Geod(ellps='WGS84')
proj = pyproj.Proj(proj='geocent', ellps='WGS84', datum='WGS84')

# 定义大地坐标系中的点
lon, lat, alt = 116.3975, 39.9086, 50  # 北京天安门广场

# 将大地坐标系中的点转换为大地空间直角坐标系中的点
x, y, z = pyproj.transform(geod, proj, lon, lat, alt, radians=False)

print(f"大地坐标：({lon}, {lat}, {alt})")
print(f"大地空间直角坐标：({x:.2f}, {y:.2f}, {z:.2f})")

输出结果如下：

大地坐标：(116.3975, 39.9086, 50)
大地空间直角坐标：(3859308.98, 331128.05, 4982809.34)

其中，`geod` 对象定义了大地坐标系的投影方式，`proj` 对象定义了大地空间直角坐标系的投影方式。`pyproj.transform()` 方法将大地坐标系中的点转换为大地空间直角坐标系中的点。