python下gdal 地理坐标转投影坐标

在Python下使用GDAL库进行地理坐标转投影坐标，可以通过以下步骤实现：

1. 导入GDAL库

import gdal

2. 打开需要处理的栅格数据

dataset = gdal.Open("input.tif")

3. 获取数据集的投影和地理参考信息

projection = dataset.GetProjection()
geotransform = dataset.GetGeoTransform()

4. 定义需要进行转换的点的地理坐标

x = 116.0
y = 39.0

5. 创建坐标转换对象

from osgeo import osr

src_srs = osr.SpatialReference()
src_srs.ImportFromWkt(projection)

dst_srs = osr.SpatialReference()
dst_srs.SetWellKnownGeogCS("WGS84")

coord_transform = osr.CoordinateTransformation(src_srs, dst_srs)

6. 进行坐标转换

point = coord_transform.TransformPoint(x, y)

7. 输出投影坐标

print(point[0], point[1])

注意，以上代码仅适用于将地理坐标转换为WGS84投影坐标，如果需要将地理坐标转换为其他投影坐标系，需要修改第5步中的dst_srs对象。