用gdal正射校正遥感影像_gdal 正射校正

使用 gdal.Warp函数可以非常方便对遥感影像进行正射校正，这个过程需要我们确定目标影像的几何信息，包括坐标系、分辨率以及需要配准到的区域或基准影像

代码示例

以下是一个使用gdal.Warp配准影像的基本Python示例：

from osgeo import gdal

# 输入和输出文件路径
input_file_path = '/path/to/your_input_image.tif'  # 要配准的影像
output_file_path = '/path/to/your_output_image.tif'  # 配准后的输出影像
base_file_path = '/path/to/your_base_image.tif'  # 基准影像（如果需要配准到特定的基准影像）

# 打开基准影像文件获取其投影和地理变换信息
base_dataset = gdal.Open(base_file_path)
base_projection = base_dataset.GetProjection()
base_geo_transform = base_dataset.GetGeoTransform()
base_x_size = base_dataset.RasterXSize
base_y_size = base_dataset.RasterYSize
base_dataset = None  # 关闭基准影像

# 使用gdal.Warp函数进行配准
warp_options = gdal.WarpOptions(
    format='GTiff',
    outputBounds=base_geo_transform,  # 输出图像的边界（可选，如果您希望输出影像与基准影像边界一致）
    xRes=base_geo_transform[1],  # 基准影像的X方向分辨率
    yRes=-base_geo_transform[5],  # 基准影像的Y方向分辨率（注意负号，因为GeoTransform的第六个参数通常是负的）
    dstSRS=base_projection,  # 目标坐标系统
    resampleAlg=gdal.GRA_Bilinear  # 重采样算法，这里使用双线性插值
)

# 执行配准
gdal.Warp(destNameOrDestDS=output_file_path, srcDSOrSrcDSTab=input_file_path, options=warp_options)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

在这个脚本中，我们首先从基准影像中获取所需的空间参考信息，然后使用这些信息设置gdal.WarpOptions。gdal.Warp函数使用这些选项将输入影像配准到基准影像的空间参考下。

请确保将 /path/to/your_input_image.tif、/path/to/your_output_image.tif 和 /path/to/your_base_image.tif 替换为实际文件路径。

此脚本假设有一个基准影像来提供配准的参考。如果没有基准影像而是有目标的分辨率和坐标系，可以直接在gdal.WarpOptions中指定这些值。

有相应的RPC文件

如果您有相应的RPC文件，那么可以使用GDAL的gdalwarp命令或者Python API来进行基于RPC的正射校正。以下是如何使用GDAL Python进行操作的示例：

from osgeo import gdal

# 输入和输出文件路径
input_tif_path = '/path/to/your_input_image.tif'  # 影像文件
input_rpc_path = '/path/to/your_input_image_rpc.txt'  # RPC文件路径
output_tif_path = '/path/to/your_output_ortho_image.tif'  # 正射校正后的输出影像

# 加载原始影像和RPC信息
src_ds = gdal.Open(input_tif_path)

# 如果RPC信息不在TIFF的元数据中，可以从外部RPC文件加载
if not src_ds.GetMetadata('RPC'):
    src_ds.SetMetadata(gdal.Open(input_rpc_path).GetMetadata('RPC'), 'RPC')

# 设置正射校正选项
warp_options = gdal.WarpOptions(
    format='GTiff',
    rpc=True,
    resampleAlg=gdal.GRA_Bilinear,  # 可以选择不同的重采样算法
    dstSRS='EPSG:4326'  # 输出的SRS，根据需要更改，这里使用WGS84
)

# 执行正射校正
gdal.Warp(destNameOrDestDS=output_tif_path, srcDSOrSrcDSTab=src_ds, options=warp_options)

# 关闭数据集
src_ds = None
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

在这个脚本中，我们首先打开TIFF文件，并检查其元数据中是否包含RPC信息。如果不包含，我们尝试从外部RPC文件中加载它。然后我们设置gdal.WarpOptions，启用RPC并执行正射校正。

用gdal命令行校正

当然，也可以使用GDAL命令行工具gdalwarp来进行基于RPC的正射校正。如果您的TIFF影像和 RPC 文件是分开的，需要确保RPC文件与TIFF文件同名，并且在同一文件夹中，GDAL才能自动识别并应用RPC信息。

以下是使用gdalwarp进行RPC正射校正的命令行示例：

gdalwarp -rpc -t_srs EPSG:4326 input.tif output.tif
1

在这个命令中：

• rpc 标志告诉 GDAL 使用 TIFF 影像文件中的 RPC 信息。
• t_srs EPSG:4326 指定目标空间参考系统，这里使用的是 WGS84。如果想要投影到不同的坐标系，请更改为相应的 EPSG 代码。
• input.tif是原始影像文件。
• output.tif是正射校正后将生成的 TIFF 文件。
  确保替换 input.tif 和 output.tif 为实际的文件名。

如果RPC信息不在 TIFF 文件中，而是在单独的.rpc文件中，确保 .rpc 文件与 TIFF 文件同名并位于相同的目录下。GDAL 会自动查找并应用这些 RPC 信息。

在运行这个命令之前，请确保已经安装了 GDAL。在 Windows 中，可以使用或者通过Conda安装 GDAL ：

conda install -c conda-forge gdal
1

在运行命令之前，最好是先检查TIFF文件和RPC文件确保它们是匹配的，以及RPC数据是有效的。