如何通过GDAL Python API实现GeoTIFF地理变换偏移（替代系统调用）

直接设置地理变换（等价于-a_gt参数） #

你不需要通过系统调用gdal_translate，GDAL Python API提供了直接修改数据集地理变换的方法，完全对应-a_gt的功能，代码更简洁高效：

from osgeo import gdal

# 打开源GeoTIFF文件
src_ds = gdal.Open(file_in, gdal.GA_ReadOnly)
if not src_ds:
    raise RuntimeError(f"无法读取文件: {file_in}")

# 获取原地理变换并计算偏移后的值
original_gt = src_ds.GetGeoTransform()
shifted_gt = (
    original_gt[0] + shift_x,  # 左上角X偏移
    original_gt[1],            # 像素宽度
    original_gt[2],            # 旋转参数
    original_gt[3] + shift_y,  # 左上角Y偏移
    original_gt[4],            # 旋转参数
    original_gt[5]             # 像素高度
)

# 创建输出文件（复制源文件的结构和数据）
driver = gdal.GetDriverByName("GTiff")
dst_ds = driver.CreateCopy(file_out, src_ds, strict=0)
if not dst_ds:
    raise RuntimeError(f"无法创建输出文件: {file_out}")

# 设置新的地理变换，这就是`-a_gt`参数的Python API实现
dst_ds.SetGeoTransform(shifted_gt)

# 可选：复制源文件的投影信息（gdal_translate默认会复制，除非指定-a_srs）
dst_ds.SetProjection(src_ds.GetProjection())

# 保存修改并释放资源
dst_ds.FlushCache()
del dst_ds, src_ds

关于gdal.Transformer的说明 #

gdal.Transformer的核心作用是坐标转换（比如将像素坐标转地理坐标，或不同投影之间的坐标转换），而不是直接修改数据集的地理变换。如果你的需求只是给输出文件设置新的地理变换，直接使用SetGeoTransform()是最直接的等价实现。

如果一定要结合Transformer使用（比如基于已有变换做复杂坐标映射），可以通过Transformer计算新的地理变换参数，但这属于冗余操作，不如直接设置高效。

内容的提问来源于stack exchange，提问作者David Boutelier