python写东北天坐标系转大地坐标系

可以借助第三方库 `pyproj` 来实现东北天坐标系（高斯-克吕格投影）到大地坐标系（经纬度）的转换。以下是一个简单的示例代码：

import pyproj

# 定义投影坐标系
crs_bj54 = pyproj.Proj(proj='tmerc', lon_0=120, lat_0=0, k=1, x_0=500000, y_0=0, ellps='intl', units='m')

# 定义地理坐标系
crs_wgs84 = pyproj.Proj(proj='latlong', datum='WGS84', ellps='WGS84')

# 定义东北天坐标系坐标
x, y = 579898.23, 4558078.55

# 投影转换
lon, lat = pyproj.transform(crs_bj54, crs_wgs84, x, y)

# 输出结果
print(lon, lat)

其中，`crs_bj54` 是北京54坐标系（即东北天坐标系）的投影坐标系参数，`crs_wgs84` 是WGS84坐标系（即地理坐标系）的参数。`x` 和 `y` 分别是东北天坐标系下的横纵坐标，`lon` 和 `lat` 分别是大地坐标系下的经纬度。