python代码实现CGCS2000大地坐标系XY值转化为对应经纬度

在Python中，转换CGCS2000大地坐标系下的XYZ值到经纬度，通常需要使用地理信息系统库如`pyproj`来进行坐标转换。这里是一个简单的例子：

from pyproj import Transformer

# 创建一个坐标投影器，源坐标系统是CGCS2000，目标是WGS84（经纬度常用的标准）
transformer = Transformer.from_crs('EPSG:15722', 'EPSG:4326', always_xy=True)

# 假设你有XYZ坐标的列表，例如 [(x1, y1, z1), (x2, y2, z2), ...]
coordinates_xyz = [(x1, y1, z1), (x2, y2, z2)]

# 对每个点进行转换
coordinates_latlon = []
for xyz in coordinates_xyz:
    xy = (xyz[0], xyz[1])  # 只取平面坐标，因为Z通常是海拔信息
    lat, lon = transformer.transform(xy[0], xy[1])
    coordinates_latlon.append((lat, lon))

# 结果存储为经纬度对列表
print(coordinates_latlon)

在这个例子中，你需要确保已经安装了`pyproj`库，如果没有，可以使用pip安装：`pip install pyproj`