WGS84坐标系转地心直角坐标系 代码

以下是Python代码示例，使用了pyproj库进行坐标转换：

import pyproj

# 定义WGS84坐标系和地心直角坐标系
wgs84 = pyproj.CRS('EPSG:4326')
ecef = pyproj.CRS('EPSG:4978')

# 定义经纬高坐标
lon, lat, height = 116.3975, 39.9086, 50

# 创建转换器并进行转换
converter = pyproj.Transformer.from_crs(wgs84, ecef, always_xy=True)
x, y, z = converter.transform(lon, lat, height)

print(f"ECEF坐标：({x:.2f}, {y:.2f}, {z:.2f})")

其中，EPSG:4326代表WGS84坐标系，EPSG:4978代表地心直角坐标系。经纬高坐标(lon, lat, height)可以根据需要进行修改。