python gdal 遥感影像配准

要在Python中使用GDAL库进行遥感影像配准，您可以使用以下代码示例：

from osgeo import gdal

def image_registration(input_image_path, reference_image_path, output_image_path):
    # 打开需要配准的影像
    src_ds = gdal.Open(input_image_path)

    # 打开参考影像
    ref_ds = gdal.Open(reference_image_path)

    # 获取需要配准影像的地理转换信息
    src_geo_transform = src_ds.GetGeoTransform()
    # 获取参考影像的地理转换信息
    ref_geo_transform = ref_ds.GetGeoTransform()

    # 创建一个空的输出影像，用于存储配准结果
    out_ds = gdal.GetDriverByName('GTiff').Create(output_image_path, src_ds.RasterXSize, src_ds.RasterYSize, src_ds.RasterCount, src_ds.GetRasterBand(1).DataType)

    # 设置输出影像的地理转换信息为参考影像的地理转换信息
    out_ds.SetGeoTransform(ref_geo_transform)

    # 设置输出影像的投影信息为参考影像的投影信息
    out_ds.SetProjection(ref_ds.GetProjection())

    # 进行影像配准
    gdal.ReprojectImage(src_ds, out_ds, src_ds.GetProjection(), ref_ds.GetProjection(), gdal.GRA_NearestNeighbour)

    # 关闭数据集
    src_ds = None
    ref_ds = None
    out_ds = None

# 调用函数进行影像配准
input_image_path = 'input_image.tif'
reference_image_path = 'reference_image.tif'
output_image_path = 'output_image.tif'
image_registration(input_image_path, reference_image_path, output_image_path)

请确保将`input_image.tif`替换为需要配准的影像路径，`reference_image.tif`替换为参考影像路径，`output_image.tif`替换为输出影像路径。

这个示例代码中的`image_registration`函数接受三个参数，分别是需要配准的影像路径、参考影像路径和输出影像路径。函数会打开需要配准的影像和参考影像，并根据参考影像的地理转换信息和投影信息创建一个空的输出影像。然后使用最近邻插值方法进行影像配准，最后关闭数据集。

希望这能对您有所帮助！如果您还有其他问题，请随时提问。