如何使用python拼接2张tif文件

在Python中，你可以使用GDAL库（Geospatial Data Abstraction Library）来合并两个TIFF（Tagged Image File Format）文件。下面是一个简单的步骤指南：

1. **安装GDAL**：首先确保已经安装了GDAL库，如果没有，可以使用pip安装：

   pip install gdal

2. **导入必要的库**：

   from osgeo import gdal

3. **打开两个TIFF文件**：

   dataset1 = gdal.Open('file1.tif')
   dataset2 = gdal.Open('file2.tif')

4. **获取每个文件的几何信息（例如宽、高、投影等）**：

   width1, height1 = dataset1.RasterXSize, dataset1.RasterYSize
   width2, height2 = dataset2.RasterXSize, dataset2.RasterYSize

5. **创建一个新的Raster bands数与原图相同的空白数据集**：

   combined_dataset = gdal.GetDriverByName('GTiff').Create('merged.tif', width1 + width2, max(height1, height2), dataset1.RasterCount)
   combined_dataset.SetGeoTransform(dataset1.GetGeoTransform())  # 使用第一个文件的坐标信息
   combined_dataset.SetProjection(dataset1.GetProjection())       # 设置投影

6. **复制数据块到新数据集中**：

   band_count = dataset1.RasterCount
   for i in range(band_count):
       src_band = dataset1.GetRasterBand(i+1)  # +1 因为索引从1开始
       dest_band = combined_dataset.GetRasterBand(i+1)
       dest_band.WriteArray(src_band.ReadAsArray(0, 0, width1, height1))  # 将第一个图象的数据复制过去

       if height1 < height2:  # 如果第二个图象更高，则填充剩余部分
           dest_band.WriteArray(src_band.ReadAsArray(width1, 0, width2, height2 - height1), width1, height1)

       else:  # 否则，如果第一个图象更高，则复制其余部分
           dest_band.WriteArray(src_band.ReadAsArray(0, height1, width2, height2 - height1))

7. **关闭所有打开的数据集**：

   dataset1 = None
   dataset2 = None
   combined_dataset = None

现在，你已经成功地合并了两个TIFF文件到一个新的文件`merged.tif`。