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

实现这个功能，你需要用到以下库：

- GDAL（Geospatial Data Abstraction Library），用于读取和写入各种地理空间数据格式，包括TIF格式。
- Proj.4，用于地理空间数据的投影转换。

以下是C++代码示例：

#include <gdal_priv.h>
#include <gdalwarper.h>
#include <cpl_conv.h>
#include <proj_api.h>

int main() {
    // 输入文件路径
    const char* input_file = "input.tif";
    // 输出文件路径
    const char* output_file = "output.tif";
    // 输出文件的块大小
    const int block_size = 256;

    // 打开输入文件
    GDALDataset* input_ds = (GDALDataset*)GDALOpen(input_file, GA_ReadOnly);
    if (input_ds == NULL) {
        printf("Failed to open input file\n");
        return 1;
    }

    // 获取输入文件的地理参考信息
    OGRSpatialReference input_srs;
    input_srs.importFromWkt(input_ds->GetProjectionRef());

    // 定义输出文件的地理参考信息，这里假设输出文件的投影为 EPSG:4326
    OGRSpatialReference output_srs;
    output_srs.importFromEPSG(4326);

    // 创建投影转换器
    OGRCoordinateTransformation* transformer = OGRCreateCoordinateTransformation(&input_srs, &output_srs);
    if (transformer == NULL) {
        printf("Failed to create coordinate transformer\n");
        return 1;
    }

    // 获取输入文件的宽度和高度
    int input_width = input_ds->GetRasterXSize();
    int input_height = input_ds->GetRasterYSize();

    // 计算输出文件的宽度和高度
    double input_geo_transform[6];
    input_ds->GetGeoTransform(input_geo_transform);
    double x_min = input_geo_transform[0];
    double y_max = input_geo_transform[3];
    double x_max = x_min + input_width * input_geo_transform[1];
    double y_min = y_max + input_height * input_geo_transform[5];
    int output_width = (int)((x_max - x_min) / block_size + 0.5) * block_size;
    int output_height = (int)((y_max - y_min) / block_size + 0.5) * block_size;

    // 创建输出文件
    GDALDriver* driver = GetGDALDriverManager()->GetDriverByName("GTiff");
    char** options = NULL;
    GDALDataset* output_ds = driver->Create(output_file, output_width, output_height, input_ds->GetRasterCount(), GDT_Float32, options);
    if (output_ds == NULL) {
        printf("Failed to create output file\n");
        return 1;
    }

    // 设置输出文件的地理参考信息
    output_ds->SetProjection(output_srs.ExportToWkt());
    double output_geo_transform[6] = { x_min, input_geo_transform[1], 0, y_max, 0, input_geo_transform[5] };
    output_ds->SetGeoTransform(output_geo_transform);

    // 遍历输出文件的每个块，将坐标转换后写入输出文件
    for (int y = 0; y < output_height; y += block_size) {
        for (int x = 0; x < output_width; x += block_size) {
            // 计算当前块的大小
            int block_width = std::min(block_size, output_width - x);
            int block_height = std::min(block_size, output_height - y);

            // 创建输入块和输出块的内存缓存
            float* input_data = new float[block_width * block_height];
            float* output_data = new float[block_width * block_height];

            // 读取输入块的像素值
            CPLErr err = input_ds->RasterIO(GF_Read, x, y, block_width, block_height, input_data, block_width, block_height, GDT_Float32, input_ds->GetRasterCount(), NULL, 0, 0, 0);
            if (err != CE_None) {
                printf("Failed to read input block\n");
                return 1;
            }

            // 转换坐标
            for (int j = 0; j < block_height; j++) {
                for (int i = 0; i < block_width; i++) {
                    double in_x = x_min + (x + i) * input_geo_transform[1];
                    double in_y = y_max + (y + j) * input_geo_transform[5];
                    double out_x, out_y;
                    transformer->Transform(1, &in_x, &in_y, NULL);
                    out_x = (in_x - x_min) / output_geo_transform[1];
                    out_y = (y_max - in_y) / output_geo_transform[5];
                    if (out_x >= 0 && out_x < block_width && out_y >= 0 && out_y < block_height) {
                        int src_index = j * block_width + i;
                        int dst_index = (int)out_y * block_width + (int)out_x;
                        output_data[dst_index] = input_data[src_index];
                    }
                }
            }

            // 写入输出块的像素值
            err = output_ds->GetRasterBand(1)->WriteBlock(x / block_size, y / block_size, output_data);
            if (err != CE_None) {
                printf("Failed to write output block\n");
                return 1;
            }

            // 释放内存
            delete[] input_data;
            delete[] output_data;
        }
    }

    // 关闭文件和转换器
    GDALClose(input_ds);
    GDALClose(output_ds);
    OCTDestroyCoordinateTransformation(transformer);

    return 0;
}

这段代码的基本思路是先计算输出文件的大小，然后遍历输出文件的每个块，将坐标转换后写入输出文件。坐标转换的方法是先将输出坐标转换为输入坐标，然后根据输入坐标计算输出坐标。这个过程中需要用到Proj.4库的函数。最后要注意释放内存和关闭文件和转换器。