python 使用gdal转换影像投影

使用 GDAL 库来转换影像投影，需要按照以下步骤：

1. 安装 GDAL 库：可以使用 pip install gdal 命令来安装。

2. 打开源影像：使用 gdal.Open() 函数来打开源影像。

3. 指定目标投影：使用 gdal.WarpOptions() 函数来指定目标投影，可以设置目标投影的 EPSG 编码、影像分辨率等参数。

4. 进行投影转换：使用 gdal.Warp() 函数来进行投影转换，将源影像转换成目标投影下的影像。

下面是一个示例代码，演示如何将 WGS84 坐标系下的影像转换成 UTM 坐标系下的影像：

import gdal

# 定义源影像路径和目标影像路径
src_path = 'source.tif'
dst_path = 'target.tif'

# 打开源影像
src_ds = gdal.Open(src_path)

# 定义目标投影
dst_options = gdal.WarpOptions(format='GTiff', dstSRS='EPSG:32651')

# 进行投影转换
gdal.Warp(dst_path, src_ds, options=dst_options)

# 关闭源影像
src_ds = None

需要注意的是，在实际使用中，需要根据具体的业务需求和数据特点，来选择合适的投影转换方式和参数。

相关问题

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

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

1. 导入GDAL库

import gdal

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

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

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

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

4. 创建输出数据集

from osgeo import gdalconst

driver = gdal.GetDriverByName("GTiff")
output_dataset = driver.Create("output.tif", dataset.RasterXSize, dataset.RasterYSize, 1, gdalconst.GDT_Float32)
output_dataset.SetProjection(projection)
output_dataset.SetGeoTransform(geotransform)

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. 读取输入数据集的像素值，并转换为投影坐标

for y in range(dataset.RasterYSize):
    for x in range(dataset.RasterXSize):
        # 获取像素值
        pixel_value = dataset.GetRasterBand(1).ReadAsArray(x, y, 1, 1)[0, 0]
        # 将像素坐标转换为地理坐标
        geo_x = geotransform[0] + x * geotransform[1] + y * geotransform[2]
        geo_y = geotransform[3] + x * geotransform[4] + y * geotransform[5]
        # 将地理坐标转换为投影坐标
        point = coord_transform.TransformPoint(geo_x, geo_y)
        # 将投影坐标写入输出数据集
        output_dataset.GetRasterBand(1).WriteArray([[point[0]]], x, y)

7. 保存输出数据集

output_dataset.FlushCache()

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

python gdal 影像坐标转换

你可以使用 GDAL 库中的 Python 绑定来进行影像坐标转换。下面是一个简单的例子：

from osgeo import gdal, osr

# 打开影像文件
ds = gdal.Open('input.tif')

# 获取源影像的投影信息和地理信息
src_proj = osr.SpatialReference()
src_proj.ImportFromWkt(ds.GetProjectionRef())
src_geo = ds.GetGeoTransform()

# 定义目标投影
dst_proj = osr.SpatialReference()
dst_proj.SetWellKnownGeogCS('WGS84')  # 例如，设置目标坐标系为 WGS84

# 创建坐标转换对象
coord_trans = osr.CoordinateTransformation(src_proj, dst_proj)

# 获取影像的宽度和高度
width = ds.RasterXSize
height = ds.RasterYSize

# 计算转换后的坐标
ulx, uly, _ = coord_trans.TransformPoint(src_geo[0], src_geo[3])
lrx, lry, _ = coord_trans.TransformPoint(src_geo[0] + src_geo[1] * width, src_geo[3] + src_geo[5] * height)

# 输出转换后的坐标范围
print('Upper Left: ({}, {})'.format(ulx, uly))
print('Lower Right: ({}, {})'.format(lrx, lry))

# 关闭影像文件
ds = None

在上面的例子中，我们打开了一个名为 `input.tif` 的影像文件，并获取了其投影信息和地理信息。然后，我们定义了一个名为 `dst_proj` 的目标投影，创建了一个名为 `coord_trans` 的坐标转换对象，并使用 `coord_trans.TransformPoint` 方法计算了转换后的坐标范围。最后，我们输出了转换后的坐标范围，并关闭了影像文件。