cgcs2000转换到ECEF，来个例程

CGCS2000是中国大地坐标系统2000版本的简称，而ECEF（地球中心惯性坐标系）是一种基于地球参考椭球体的地理坐标系。将CGCS2000转换到ECEF通常需要考虑地球椭球体模型和转换矩阵。

以下是一个简单的Python示例，假设我们已经有了CGCS2000的经度、纬度和高度（海拔），使用WGS84作为ECEF的基础椭球体模型，因为CGCS2000和WGS84在近似上相似：

from math import sin, cos, radians
from earthpy.core import ECEF_from_geodetic

# CGCS2000坐标（以弧度表示）
lat_radians = radians(39.9)  # 北纬
lon_radians = radians(116.4)  # 东经
height_m = 0  # 海拔（如果在地面则设为0）

# WGS84的半长轴和扁率（已知常量）
a_WGS84 = 6378137.0
b_WGS84 = a_WGS84 * (1 - 0.00335281066474)

# 转换为ECEF坐标
ecef_xyz = ECEF_from_geodetic(lat_radians, lon_radians, height_m, a_WGS84, f=WGS84_first_derivative(b_WGS84))

x, y, z = ecef_xyz  # 分别得到ECEF坐标的三个分量
print(f"转换后的ECEF坐标：({x:.2f}, {y:.2f}, {z:.2f})")

请注意，实际应用中可能需要精确的转换矩阵或库来进行更准确的转换，并且上述示例仅适用于大致计算，可能存在精度误差。