遥感影像拼接分多种不同的情况,比如比较常用的是两幅影像拼接或者多个图像合并为一个图像。
GDAL中实现影像拼接的方式也有多种,比如常用可以通过numpy读取影像数组,并计算各个图像的numpy数组范围,合并多个numpy数组为一个最终的结果,并将其重新输出为一幅新的图像,即可实现图像的拼接。
另外比较方便的一种方法便是使用GDAL.warp()函数方法。这个函数的功能很多,可以实现多种不同的目的,比如用于坐标系转换、投影变换、图像合并与镶嵌、地理范围裁剪、更改分辨率、矢量裁剪等多种功能。这些功能可以通过设置Option参数来实现,如果不设置参数,一般按照原图像参数输出结果。
本文中图像拼接采用这个函数来实现。
案例一:两个文件拼接示例
通过以下函数实现:
outputfilePath:输出结果完整路径。
inputrasfile1,inputrasfile2:两个输入图像完整路径。
gdal.Warp(outputfilePath,[inputrasfile1,inputrasfile2],options=options)
import numpy as npfrom osgeo import gdal, gdalconstimport osinputfile_path1=r"E:\……1*.tif"inputfile_path2=r"E:\……2*.tif"outputfilePath =r"E:\……2*Mosaic.tif"print("output file name:",outputfilePath)print("图像拼接")inputr