大地坐标转地心直角坐标 python

可以使用Python中的geopy和geographiclib库来进行大地坐标和地心直角坐标的转换。以下是一个简单的示例代码：

from geopy import Point, distance
from geographiclib.geodesic import Geodesic

# 大地坐标
lat, lng, altitude = 39.9075, 116.39723, 50

# 计算地球半径
geod = Geodesic()
r = geod.a

# 转换为地心直角坐标
x = (r + altitude) * distance.distance(Point(lat, lng), Point(0, lng)).km
y = (r + altitude) * distance.distance(Point(lat, lng), Point(lat, 0)).km
z = (r + altitude) * (1 - geod.f) ** 2 * ((1 / geod.f ** 2) * ((1 - geod.f ** 2) * r ** 2 / (r + altitude) ** 2) + 1) ** 0.5

print("地心直角坐标：{:.2f}, {:.2f}, {:.2f}".format(x, y, z))

其中，`lat`和`lng`分别表示纬度和经度，`altitude`表示高度。`distance.distance`函数用于计算两个坐标点之间的距离，`geod.f`表示地球的扁率。最终输出的结果为地心直角坐标系下的x、y、z值。