Python遥感图像处理应用篇(二十三)：Python+GDAL 批量拼接图像

Posted 2022-12-07 空中旋转篮球

遥感影像拼接分多种不同的情况，比如比较常用的是两幅影像拼接或者多个图像合并为一个图像。

GDAL中实现影像拼接的方式也有多种，比如常用可以通过numpy读取影像数组，并计算各个图像的numpy数组范围，合并多个numpy数组为一个最终的结果，并将其重新输出为一幅新的图像，即可实现图像的拼接。

另外比较方便的一种方法便是使用GDAL.warp()函数方法。这个函数的功能很多，可以实现多种不同的目的，比如用于坐标系转换、投影变换、图像合并与镶嵌、地理范围裁剪、更改分辨率、矢量裁剪等多种功能。这些功能可以通过设置Option参数来实现，如果不设置参数，一般按照原图像参数输出结果。

本文中图像拼接采用这个函数来实现。

案例一：两个文件拼接示例

通过以下函数实现：

outputfilePath：输出结果完整路径。

inputrasfile1,inputrasfile2：两个输入图像完整路径。

gdal.Warp(outputfilePath,[inputrasfile1,inputrasfile2],options=options)

import numpy as np
from osgeo import gdal, gdalconst
import os

inputfile_path1=r"E:\\……1*.tif"
inputfile_path2=r"E:\\……2*.tif"

outputfilePath =r"E:\\……2*Mosaic.tif"
print("output file name:",outputfilePath)

print("图像拼接")
inputr