C++中GDAL：用像素均为0的栅格填补时序遥感数据中缺少的时相

#include <iostream>
#include <iomanip>
#include <sstream>
#include "gdal_priv.h"
#include "cpl_conv.h"

using namespace std;

void create_missing_raster(string path);

int main() {
    string file_path = R"(E:\02_Project\TIFF\TEST)";
    create_missing_raster(file_path);
    return 0;
}

void create_missing_raster(string path)
{
    vector<string> all_file_path;
    for (int year = 2018; year <= 2022; year++)
    {
        for (int day = 1; day <= 361; day += 10)
        {
            ostringstream oss;
            oss << path << "/Albers_MuSyQ.NDVI.16m." << to_string(year) << setfill('0') << setw(3) << to_string(day) << "000000.49SDB.001.h5NDVI-NDVI.tif";
            all_file_path.push_back(oss.str());
        }
    }

    GDALAllRegister();
    int x_size, y_size;
    string one_actual_path;
    for (string& one_file_path : all_file_path)
    {
        GDALDataset* poDataset_actual;
        poDataset_actual = (GDALDataset*)GDALOpen(one_file_path.c_str(), GA_ReadOnly);
        if (poDataset_actual != NULL)
        {
            one_actual_path = one_file_path;
            x_size = poDataset_actual->GetRasterXSize();
            y_size = poDataset_actual->GetRasterYSize();
            GDALClose((GDALDatasetH)poDataset_actual);
            break;
        }
    }

    for (string& one_file_path : all_file_path)
    {
        if (CPLCheckForFile((char*)one_file_path.c_str(), NULL) == FALSE)
        {
            GDALDataset* poDataset;
            GDALDriver* poDriver;
            poDriver = GetGDALDriverManager()->GetDriverByName("GTiff");
            GDALDataset* poDataset_actual;
            poDataset_actual = (GDALDataset*)GDALOpen(one_actual_path.c_str(), GA_ReadOnly);
            poDataset = poDriver->CreateCopy(one_file_path.c_str(), poDataset_actual, FALSE, NULL, NULL, NULL);
            double* pafScanline = new double[x_size * y_size];
            for (int i = 0; i < x_size * y_size; i++)
            {
                pafScanline[i] = 0.0;
            }
            GDALRasterBand* poBand;
            poBand = poDataset->GetRasterBand(1);
            poBand->RasterIO(GF_Write, 0, 0, x_size, y_size, pafScanline, x_size, y_size, GDT_Float64, 0, 0);
            delete[] pafScanline;
            GDALClose((GDALDatasetH)poDataset);
            GDALClose((GDALDatasetH)poDataset_actual);
            cout << "New file is :" << one_file_path << endl;
        }
    }
    GDALDestroyDriverManager();
}