c++获取多幅单波段tif遥感影像的重叠区域

要获取多幅单波段tif遥感影像的重叠区域，可以按照以下步骤：

1. 读取多幅tif影像，使用GDAL开源库可以方便地读取tif影像数据。

2. 获取每幅tif影像的地理坐标信息，包括左上角和右下角的经纬度坐标，以及像元宽度和高度。

3. 将每幅tif影像的地理坐标信息转换为像元坐标信息，即行列号，可以使用GDAL库中的GDALRasterBand::GetGeoTransform函数。

4. 找到多幅tif影像的重叠区域，可以通过比较每幅tif影像的像元坐标信息，找到它们的重叠部分。

5. 计算重叠区域的像元坐标范围，即左上角和右下角的行列号，可以使用std::min和std::max函数。

6. 利用计算出的像元坐标范围，裁剪每幅tif影像，得到它们的重叠部分。

下面是一个简单的示例代码：

#include <iostream>
#include "gdal/gdal.h"
#include "gdal/gdal_priv.h"

int main()
{
    // 读取多幅tif影像
    GDALAllRegister();
    GDALDataset* dataset1 = (GDALDataset*)GDALOpen("image1.tif", GA_ReadOnly);
    GDALDataset* dataset2 = (GDALDataset*)GDALOpen("image2.tif", GA_ReadOnly);

    // 获取每幅tif影像的地理坐标信息
    double adfGeoTransform[6];
    dataset1->GetGeoTransform(adfGeoTransform);
    double x1 = adfGeoTransform[0];
    double y1 = adfGeoTransform[3];
    double pixelWidth1 = adfGeoTransform[1];
    double pixelHeight1 = adfGeoTransform[5];

    dataset2->GetGeoTransform(adfGeoTransform);
    double x2 = adfGeoTransform[0];
    double y2 = adfGeoTransform[3];
    double pixelWidth2 = adfGeoTransform[1];
    double pixelHeight2 = adfGeoTransform[5];

    // 将每幅tif影像的地理坐标信息转换为像元坐标信息
    int x1_start, x1_end, y1_start, y1_end;
    double x1_min = std::min(x1, x2);
    double y1_max = std::max(y1, y2);
    x1_start = (int)((x1_min - x1) / pixelWidth1 + 0.5);
    y1_start = (int)((y1_max - y1) / pixelHeight1 + 0.5);

    int x2_start, x2_end, y2_start, y2_end;
    double x2_min = std::max(x1, x2);
    double y2_min = std::min(y1, y2);
    x2_start = (int)((x2_min - x2) / pixelWidth2 + 0.5);
    y2_start = (int)((y2_min - y2) / pixelHeight2 + 0.5);

    int x1_size = dataset1->GetRasterXSize();
    int y1_size = dataset1->GetRasterYSize();
    int x2_size = dataset2->GetRasterXSize();
    int y2_size = dataset2->GetRasterYSize();

    x1_end = std::min(x1_start + x2_size - x2_start, x1_size);
    y1_end = std::min(y1_start + y2_size - y2_start, y1_size);
    x2_end = std::min(x2_start + x1_size - x1_start, x2_size);
    y2_end = std::min(y2_start + y1_size - y1_start, y2_size);

    // 输出重叠区域的像元坐标范围
    std::cout << "overlap region: (" << x1_start << ", " << y1_start << ") - (" << x1_end << ", " << y1_end << ")\n";

    // 裁剪每幅tif影像，得到它们的重叠部分
    GDALDataset* overlap1 = dataset1->GetRasterBand(1)->ReadAsArray(x1_start, y1_start, x1_end - x1_start, y1_end - y1_start);
    GDALDataset* overlap2 = dataset2->GetRasterBand(1)->ReadAsArray(x2_start, y2_start, x2_end - x2_start, y2_end - y2_start);

    // 释放资源
    GDALClose(overlap1);
    GDALClose(overlap2);
    GDALClose(dataset1);
    GDALClose(dataset2);

    return 0;
}

注意，上面的代码只是一个简单的示例，实际应用中可能需要处理更多的异常情况。