GDAL投影坐标系转大地坐标系python

在Python中，你可以使用GDAL库来实现投影坐标系到大地坐标系的转换。下面是一个示例代码：

from osgeo import gdal, osr

# 定义投影坐标系
source_projection = osr.SpatialReference()
source_projection.ImportFromEPSG(4326)  # EPSG代码，这里使用WGS84坐标系

# 定义大地坐标系
target_projection = osr.SpatialReference()
target_projection.ImportFromEPSG(3857)  # EPSG代码，这里使用Web墨卡托投影坐标系

# 创建坐标转换对象
transform = osr.CoordinateTransformation(source_projection, target_projection)

# 定义投影坐标点
x = 116.3975
y = 39.9085

# 进行坐标转换
transformed_point = transform.TransformPoint(x, y)

# 输出大地坐标点
print("转换后的坐标：", transformed_point)

这段代码中，首先定义了源投影坐标系和目标大地坐标系。然后，创建了一个坐标转换对象，并使用`TransformPoint`方法将投影坐标点转换为大地坐标点。最后，将转换后的大地坐标点打印出来。

请注意，你需要先安装GDAL库，可以使用pip进行安装：`pip install gdal`。另外，需要确保你已经正确安装了GDAL的依赖库和数据文件。