c++给定一张TIF格式的影像，给定一套格网重投影坐标，输出一张新的TIF影像（要求是分块存储格式Tiled TIF）；

实现这个功能，需要用到一些库，包括`GDAL`和`OpenCV`。

1. 首先需要读取输入的TIF影像，可以使用`GDAL`库中的`GDALOpen()`函数打开TIF文件。具体的代码如下：

   GDALDataset* inputDataset = (GDALDataset*)GDALOpen(inputTifPath.c_str(), GA_ReadOnly);
   if (inputDataset == nullptr)
   {
       std::cout << "Failed to open input TIF file!" << std::endl;
       return -1;
   }

2. 接下来需要获取输入影像的一些信息，包括分辨率、坐标系、行列数等等。可以使用`GDAL`库中的相关函数来获取这些信息，具体代码如下：

   double adfGeoTransform[6];
   inputDataset->GetGeoTransform(adfGeoTransform);    // 获取影像的六参数
   int xSize = inputDataset->GetRasterXSize();        // 获取影像的列数
   int ySize = inputDataset->GetRasterYSize();        // 获取影像的行数
   const char* projectionRef = inputDataset->GetProjectionRef();    // 获取影像的投影信息

3. 接下来就是将影像进行重投影，并输出到新的TIF文件中。具体步骤如下：

1. 定义输出影像的一些参数，包括影像大小、坐标系、分块大小等等。这里以输出为EPSG:4326坐标系的影像为例：

      int outXSize = 1000;    // 输出影像的列数
      int outYSize = 1000;    // 输出影像的行数
      const char* outProjectionRef = "EPSG:4326";    // 输出影像的坐标系
      int tileWidth = 256;    // 分块的宽度
      int tileHeight = 256;   // 分块的高度

2. 创建输出影像，并设置其一些参数，包括大小、坐标系、分块大小等等。具体的代码如下：

      GDALDriver* driver = GetGDALDriverManager()->GetDriverByName("GTiff");
      if (driver == nullptr)
      {
          std::cout << "Failed to get GDAL driver!" << std::endl;
          return -1;
      }

      // 创建输出影像
      GDALDataset* outputDataset = driver->Create(outputTifPath.c_str(), outXSize, outYSize, inputDataset->GetRasterCount(), GDT_Byte, nullptr);
      if (outputDataset == nullptr)
      {
          std::cout << "Failed to create output TIF file!" << std::endl;
          return -1;
      }

      // 设置输出影像的一些参数
      double outGeoTransform[6];
      outGeoTransform[0] = adfGeoTransform[0];   // 左上角x坐标
      outGeoTransform[1] = adfGeoTransform[1];   // 东西方向像素分辨率
      outGeoTransform[2] = adfGeoTransform[2];   // 旋转角度，通常为0
      outGeoTransform[3] = adfGeoTransform[3];   // 左上角y坐标
      outGeoTransform[4] = adfGeoTransform[4];   // 旋转角度，通常为0
      outGeoTransform[5] = adfGeoTransform[5];   // 南北方向像素分辨率
      outputDataset->SetGeoTransform(outGeoTransform);    // 设置输出影像的六参数
      outputDataset->SetProjection(outProjectionRef);     // 设置输出影像的投影信息
      outputDataset->SetMetadataItem("TIFFTAG_TILEWIDTH", std::to_string(tileWidth).c_str());     // 设置输出影像的分块大小
      outputDataset->SetMetadataItem("TIFFTAG_TILELENGTH", std::to_string(tileHeight).c_str());

3. 对于每一个分块，将其重投影并写入输出影像中。具体的代码如下：

      for (int y = 0; y < outYSize; y += tileHeight)
      {
          for (int x = 0; x < outXSize; x += tileWidth)
          {
              // 计算当前分块在输入影像中的位置和大小
              int inX = (int)(adfGeoTransform[0] + x * adfGeoTransform[1] + y * adfGeoTransform[2]);
              int inY = (int)(adfGeoTransform[3] + x * adfGeoTransform[4] + y * adfGeoTransform[5]);
              int inWidth = std::min(tileWidth, outXSize - x);
              int inHeight = std::min(tileHeight, outYSize - y);

              // 创建输入影像的ROI
              cv::Rect roi(inX, inY, inWidth, inHeight);

              // 读取输入影像的ROI
              cv::Mat inputImage;
              cv::Mat inputImageROI = cv::imread(inputTifPath, cv::IMREAD_UNCHANGED);
              inputImageROI(roi).copyTo(inputImage);

              // 对输入影像进行重投影
              cv::Mat outputImage;
              cv::Size outputSize(inWidth, inHeight);
              cv::Mat mapX, mapY;
              cv::fisheye::initUndistortRectifyMap(intrinsic, distortion, cv::Mat(), intrinsic, outputSize, CV_32FC1, mapX, mapY);
              cv::remap(inputImage, outputImage, mapX, mapY, cv::INTER_LINEAR);

              // 将重投影后的分块写入输出影像中
              GDALRasterBand* outputBand = outputDataset->GetRasterBand(1);
              outputBand->RasterIO(GF_Write, x, y, inWidth, inHeight, outputImage.data, inWidth, inHeight, GDT_Byte, 0, 0);
          }
      }

4. 最后需要释放一些资源，包括输入影像和输出影像。具体的代码如下：

   GDALClose(inputDataset);
   GDALClose(outputDataset);

完整的代码示例可以参考下面的代码：