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

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

1. 使用GDAL库打开这4幅tif遥感影像。
2. 获取这4幅影像的地理坐标系和投影坐标系信息，并确保它们具有相同的参考系。
3. 获取这4幅影像的地理范围，也就是它们的左上角和右下角坐标。
4. 将这4幅影像的地理范围转换为像素坐标系，也就是它们在影像中的行列号。
5. 计算出这4幅影像在像素坐标系中的重叠区域，也就是它们在行列号上的交集。

以下是一个C++代码示例，用于获取4幅单波段tif遥感影像的重叠区域：

#include <gdal_priv.h>

int main()
{
    // 打开4幅tif遥感影像
    GDALDataset* ds1 = (GDALDataset*) GDALOpen("image1.tif", GA_ReadOnly);
    GDALDataset* ds2 = (GDALDataset*) GDALOpen("image2.tif", GA_ReadOnly);
    GDALDataset* ds3 = (GDALDataset*) GDALOpen("image3.tif", GA_ReadOnly);
    GDALDataset* ds4 = (GDALDataset*) GDALOpen("image4.tif", GA_ReadOnly);

    // 获取地理坐标系和投影坐标系信息
    const char* proj1 = ds1->GetProjectionRef();
    const char* proj2 = ds2->GetProjectionRef();
    const char* proj3 = ds3->GetProjectionRef();
    const char* proj4 = ds4->GetProjectionRef();
    double geoTransform1[6], geoTransform2[6], geoTransform3[6], geoTransform4[6];
    ds1->GetGeoTransform(geoTransform1);
    ds2->GetGeoTransform(geoTransform2);
    ds3->GetGeoTransform(geoTransform3);
    ds4->GetGeoTransform(geoTransform4);

    // 确保4幅影像具有相同的参考系
    if (strcmp(proj1, proj2) != 0 || strcmp(proj2, proj3) != 0 || strcmp(proj3, proj4) != 0)
    {
        printf("Error: the input images have different projections.\n");
        return 1;
    }

    // 获取地理范围
    double minX1 = geoTransform1[0];
    double maxX1 = geoTransform1[0] + geoTransform1[1] * ds1->GetRasterXSize();
    double minY1 = geoTransform1[3] + geoTransform1[5] * ds1->GetRasterYSize();
    double maxY1 = geoTransform1[3];
    double minX2 = geoTransform2[0];
    double maxX2 = geoTransform2[0] + geoTransform2[1] * ds2->GetRasterXSize();
    double minY2 = geoTransform2[3] + geoTransform2[5] * ds2->GetRasterYSize();
    double maxY2 = geoTransform2[3];
    double minX3 = geoTransform3[0];
    double maxX3 = geoTransform3[0] + geoTransform3[1] * ds3->GetRasterXSize();
    double minY3 = geoTransform3[3] + geoTransform3[5] * ds3->GetRasterYSize();
    double maxY3 = geoTransform3[3];
    double minX4 = geoTransform4[0];
    double maxX4 = geoTransform4[0] + geoTransform4[1] * ds4->GetRasterXSize();
    double minY4 = geoTransform4[3] + geoTransform4[5] * ds4->GetRasterYSize();
    double maxY4 = geoTransform4[3];

    // 转换为像素坐标系
    int minXPixel1, maxXPixel1, minYPixel1, maxYPixel1;
    ds1->GetGeoTransform(geoTransform1);
    ds1->GetRasterBand(1)->GetBlockSize(&minXPixel1, &minYPixel1);
    ds1->GetRasterBand(1)->GetBlockSize(&maxXPixel1, &maxYPixel1);
    minXPixel1 = (int) ((minX1 - geoTransform1[0]) / geoTransform1[1]);
    maxXPixel1 = (int) ((maxX1 - geoTransform1[0]) / geoTransform1[1]);
    minYPixel1 = (int) ((maxY1 - geoTransform1[3]) / geoTransform1[5]);
    maxYPixel1 = (int) ((minY1 - geoTransform1[3]) / geoTransform1[5]);
    int minXPixel2, maxXPixel2, minYPixel2, maxYPixel2;
    ds2->GetGeoTransform(geoTransform2);
    ds2->GetRasterBand(1)->GetBlockSize(&minXPixel2, &minYPixel2);
    ds2->GetRasterBand(1)->GetBlockSize(&maxXPixel2, &maxYPixel2);
    minXPixel2 = (int) ((minX2 - geoTransform2[0]) / geoTransform2[1]);
    maxXPixel2 = (int) ((maxX2 - geoTransform2[0]) / geoTransform2[1]);
    minYPixel2 = (int) ((maxY2 - geoTransform2[3]) / geoTransform2[5]);
    maxYPixel2 = (int) ((minY2 - geoTransform2[3]) / geoTransform2[5]);
    int minXPixel3, maxXPixel3, minYPixel3, maxYPixel3;
    ds3->GetGeoTransform(geoTransform3);
    ds3->GetRasterBand(1)->GetBlockSize(&minXPixel3, &minYPixel3);
    ds3->GetRasterBand(1)->GetBlockSize(&maxXPixel3, &maxYPixel3);
    minXPixel3 = (int) ((minX3 - geoTransform3[0]) / geoTransform3[1]);
    maxXPixel3 = (int) ((maxX3 - geoTransform3[0]) / geoTransform3[1]);
    minYPixel3 = (int) ((maxY3 - geoTransform3[3]) / geoTransform3[5]);
    maxYPixel3 = (int) ((minY3 - geoTransform3[3]) / geoTransform3[5]);
    int minXPixel4, maxXPixel4, minYPixel4, maxYPixel4;
    ds4->GetGeoTransform(geoTransform4);
    ds4->GetRasterBand(1)->GetBlockSize(&minXPixel4, &minYPixel4);
    ds4->GetRasterBand(1)->GetBlockSize(&maxXPixel4, &maxYPixel4);
    minXPixel4 = (int) ((minX4 - geoTransform4[0]) / geoTransform4[1]);
    maxXPixel4 = (int) ((maxX4 - geoTransform4[0]) / geoTransform4[1]);
    minYPixel4 = (int) ((maxY4 - geoTransform4[3]) / geoTransform4[5]);
    maxYPixel4 = (int) ((minY4 - geoTransform4[3]) / geoTransform4[5]);

    // 计算重叠区域
    int minXPixel = std::max(std::max(minXPixel1, minXPixel2), std::max(minXPixel3, minXPixel4));
    int maxXPixel = std::min(std::min(maxXPixel1, maxXPixel2), std::min(maxXPixel3, maxXPixel4));
    int minYPixel = std::max(std::max(minYPixel1, minYPixel2), std::max(minYPixel3, minYPixel4));
    int maxYPixel = std::min(std::min(maxYPixel1, maxYPixel2), std::min(maxYPixel3, maxYPixel4));

    // 输出结果
    printf("The overlap region in pixel coordinates is (%d,%d)-(%d,%d).\n", minXPixel, minYPixel, maxXPixel, maxYPixel);

    // 关闭数据集
    GDALClose(ds1);
    GDALClose(ds2);
    GDALClose(ds3);
    GDALClose(ds4);

    return 0;
}

注意，上述代码示例中使用了GDAL库来读取和处理tif遥感影像。如果您还没有安装GDAL库，请先安装。