GDAL读取MODIS影像

int main(int argc, char* argv[])
{
    const char* pszFileName = argv[1];  // 从命令行传入MODIS数据文件路径
    GDALAllRegister();
    GDALDataset* poDataset = static_cast<GDALDataset*>(GDALOpen(pszFileName, GA_ReadOnly));
    if (!poDataset)
    {
        cout << "影像读取失败：" << pszFileName << "！" << '\n';
        return EXIT_FAILURE;
    }

    char** metadata = poDataset->GetMetadata("SUBDATASETS");
    vector<string> vtDatasets;
    // 我们可以将HDF数据看做一个容器，容器里面包含了很多具体的子Dataset，所以我们首先要取出子数据集的地址
    // 元数据metadata中保存的信息中，奇数索引的信息是子数据的地址，偶数索引的信息是该子数据的其它信息
    if (0 < CSLCount(const_cast<const char* const*>(metadata)))
    {
        for (int i = 0; metadata[i]; i += 2)
        {
            string strDSName(metadata[i]);
            cout << strDSName << '\n';
            strDSName = strDSName.substr(strDSName.find_first_of('=') + 1);
            vtDatasets.push_back(strDSName);
        }
    }
    GDALClose(poDataset);

    for (auto name : vtDatasets)
    {
        const char* pszDSName = name.c_str();
        cout << "数据路径：" << pszDSName << '\n';
        readGeoRaster(pszDSName);
    }

    return EXIT_SUCCESS;
}