C++将h5转tif：支持高分数据等szip压缩的图像

2025-07-27 13:36:09

C++将h5转tif：支持高分数据等szip压缩的图像本文介绍基于C++ 语言的hdf5库与gdal库，将.h5格式的多波段HDF5图像批量转换为.tif格式的方法；其中，本方法支持对szip压缩的HDF5图像（例如高分一号卫星遥感影像）加以转换。将HDF5图像批量转换为.tif格式，在部分场景下操作并不难——在我们之前的文章ArcPy将HDF格式栅格文件批量转为TIFF格式中，就介绍过

C++将h5转tif：支持高分数据等szip压缩的图像

本文介绍基于C++ 语言的hdf5库与gdal库，将.h5格式的多波段HDF5图像批量转换为.tif格式的方法；其中，本方法支持对szip压缩的HDF5图像（例如高分一号卫星遥感影像）加以转换。

将HDF5图像批量转换为.tif格式，在部分场景下操作并不难——在我们之前的文章ArcPy将HDF格式栅格文件批量转为TIFF格式中，就介绍过基于Python中的arcpy模块实现这一需求的方法。但是，正如我们在文章HDF4与HDF5文件的打开方式：HDFView软件中提到的那样，由于szip这个压缩模块不再受到hdf5库的支持，导致用szip程序压缩的HDF5图像，比如高分系列遥感影像的.h5文件，就没办法在Windows中通过Python的h5py、gdal等库直接打开了。

那么在这里，我们就介绍一下基于C++ 语言的hdf5库，打开.h5格式图像（包括那些用到szip压缩程序的HDF5图像）的方法。不过需要注意，我这里是在Linux的Ubuntu系统中操作的，至少可以保证这个代码在Linux下可以正常运行；但能否在Windows中的C++ 环境下也正常运行，我暂时还没试过——按道理应该也是可行的，大家如果有需要的话可以尝试一下。

本文所用代码如下。

代码语言：javascript代码运行次数：0运行复制

#include <iostream>
#include <sstream>
#include <vector>
#include <filesystem>
#include <gdal.h>
#include <gdal_priv.h>
#include "hdf5.h"
#include "ogr_spatialref.h" 

int main(int argc, char *ar[]) {
    ct std::string h5_path = "/home/ctj/data/H5/";
    ct std::string tif_path = "/home/ctj/data/TIFF_48SUB/";
    // ct std::string h5_path = ar[1];
    // ct std::string tif_path = ar[2];
    ct char *dataset_0 = "/Cloud_Mask/cloudmask";
    ct char *dataset_1 = "/GeometricCorrection/DataSet_16_1";
    ct char *dataset_2 = "/GeometricCorrection/DataSet_16_2";
    ct char *dataset_ = "/GeometricCorrection/DataSet_16_";
    ct char *dataset_4 = "/GeometricCorrection/DataSet_16_4";
    ct char *projection_para = "ProjectionPara";
    ct char *projection_str = "ProjectionStr";

    hid_t file_id;
    hid_t dataset_id;
    hid_t attr_id;
    hid_t attr_dtype;
    herr_t status;
    hid_t mem_type_id = H5T_ATIVE_UIT16;
    int size = 686;
    int band_num = 5;

    // namespace fs = filesystem;

    status = H5open();
    GDALAllRegister();

    for (ct auto& entry : std::filesystem::directory_iterator(h5_path)) {
        if (entry.path().extension() == ".h5") {
            std::string filename = entry.path().filename().string();
            std::cout << filename << std::endl;
            std::string baseame = filename.substr(0, filename.find_last_of('.'));
            ct std::string output_filename = tif_path + baseame + ".tif";

            file_id = H5Fopen((h5_path + filename).c_str(), H5F_ACC_RDOLY, H5P_DEFAULT);

            attr_id = H5Aopen(file_id, projection_para, H5P_DEFAULT);
            attr_dtype = H5Aget_type(attr_id);
            size_t string_length = H5Tget_size(attr_dtype);
            char *attr_data = new char[1000];
            status = H5Aread(attr_id, attr_dtype, attr_data);

            std::istringstream iss(attr_data);
            std::vector<double> transform(6);
            int i = 0;
            std::string str;
            while (getline(iss, str, ',')) {
                    transform[i] = stod(str);
                    ++i;
                }

            attr_id = H5Aopen(file_id, projection_str, H5P_DEFAULT);
            attr_dtype = H5Aget_type(attr_id);
            char *attr_data_str = new char[1000];
            status = H5Aread(attr_id, attr_dtype, attr_data_str);

            dataset_id = H5Dopen1(file_id, dataset_0);
            std::vector<u_int16_t> data_0(size * size);
            status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_0.data());

            dataset_id = H5Dopen1(file_id, dataset_1);
            std::vector<u_int16_t> data_1(size * size);
            status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_1.data());

            dataset_id = H5Dopen1(file_id, dataset_2);
            std::vector<u_int16_t> data_2(size * size);
            status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_2.data());

            dataset_id = H5Dopen1(file_id, dataset_);
            std::vector<u_int16_t> data_(size * size);
            status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_.data());

            dataset_id = H5Dopen1(file_id, dataset_4);
            std::vector<u_int16_t> data_4(size * size);
            status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_4.data());

            status = H5Fclose(file_id);

            GDALDriver *poDriver = GetGDALDriverManager()->GetDriverByame("GTiff");

            GDALDataset *poDstDS = poDriver->Create(output__str(), size, size, band_num, GDT_UInt16, nullptr);
            u_int16_t *band_data_0 = &data_0[0];
            poDstDS->GetRasterBand(1)->RasterIO(GF_Write, 0, 0, size, size, band_data_0, size, size, GDT_UInt16, 0, 0);
            u_int16_t *band_data_1 = &data_1[0];
            poDstDS->GetRasterBand(2)->RasterIO(GF_Write, 0, 0, size, size, band_data_1, size, size, GDT_UInt16, 0, 0);
            u_int16_t *band_data_2 = &data_2[0];
            poDstDS->GetRasterBand()->RasterIO(GF_Write, 0, 0, size, size, band_data_2, size, size, GDT_UInt16, 0, 0);
            u_int16_t *band_data_ = &data_[0];
            poDstDS->GetRasterBand(4)->RasterIO(GF_Write, 0, 0, size, size, band_data_, size, size, GDT_UInt16, 0, 0);
            u_int16_t *band_data_4 = &data_4[0];
            poDstDS->GetRasterBand(5)->RasterIO(GF_Write, 0, 0, size, size, band_data_4, size, size, GDT_UInt16, 0, 0);

            for (int i = 1; i <= band_num; ++i) {
                GDALRasterBand *poBand = poDstDS->GetRasterBand(i);
                if (poBand != nullptr) {
                    poBand->SetoDataValue(0);
                }
            }

            poDstDS->SetGeoTransform(transform.data());
            poDstDS->SetProjection(attr_data_str);

            GDALClose(poDstDS);
        }
    }

    status = H5close();
    return 0;
}

上述是本文完整代码。接下来，就分段介绍一下每段代码的具体含义。

首先，需要包含必要的头文件。在这里，包括标准输入输出、字符串流、向量、文件系统等功能，以及hdf5库与gdal库。同时，定义了两个常量字符串h5_path与tif_path，分别指向转换前的HDF5图像和转换后的TIFF图像的目录。

代码语言：javascript代码运行次数：0运行复制

#include <iostream>
#include <sstream>
#include <vector>
#include <filesystem>
#include <gdal.h>
#include <gdal_priv.h>
#include "hdf5.h"
#include "ogr_spatialref.h" 

int main(int argc, char *ar[]) {
    ct std::string h5_path = "/home/ctj/data/H5/";
    ct std::string tif_path = "/home/ctj/data/TIFF_48SUB/";

随后，设定要读取的HDF5图像的数据集（波段）的路径，以及空间参考信息的属性名称；这些参数大家就按照自己HDF5图像的实际情况来修改即可。

接下来，初始化hdf5库的状态变量，这些变量是hdf5库操作需要的。同时，用size表示图像的宽度和高度，因为我这里HDF5图像是正方形，所以只需指定1个值。此外，band_num表示待转换遥感影像的波段数。

代码语言：javascript代码运行次数：0运行复制

ct char *dataset_0 = "/Cloud_Mask/cloudmask";
ct char *dataset_1 = "/GeometricCorrection/DataSet_16_1";
// ... 省略部分代码 ...
ct char *projection_para = "ProjectionPara";
ct char *projection_str = "ProjectionStr";

hid_t file_id;
hid_t dataset_id;
hid_t attr_id;
hid_t attr_dtype;
herr_t status;
hid_t mem_type_id = H5T_ATIVE_UIT16;
int size = 686;
int band_num = 5;

紧接着，初始化hdf5库，注册所有可用的GDAL驱动程序。

代码语言：javascript代码运行次数：0运行复制

status = H5open();
GDALAllRegister();

随后，使用std::filesystem::directory_iterator遍历指定目录中的所有文件，并只处理扩展名为.h5的文件；对于这些文件，构建输出文件名——基于原始文件名，去掉扩展名并添加.tif扩展名。

代码语言：javascript代码运行次数：0运行复制

for (ct auto& entry : std::filesystem::directory_iterator(h5_path)) {
    if (entry.path().extension() == ".h5") {
        std::string filename = entry.path().filename().string();
        std::cout << filename << std::endl;
        std::string baseame = filename.substr(0, filename.find_last_of('.'));
        ct std::string output_filename = tif_path + baseame + ".tif";

随后，使用H5Fopen打开HDF5图像，在这里选择以只读模式访问。

代码语言：javascript代码运行次数：0运行复制

file_id = H5Fopen((h5_path + filename).c_str(), H5F_ACC_RDOLY, H5P_DEFAULT);

随后，需要读取原本HDF5图像的空间参考信息。在这里，首先打开名为projection_para的属性，读取其内容到attr_data中；随后，解析attr_data为一个包含6个元素的double向量transform——这些元素用于地理变换。

代码语言：javascript代码运行次数：0运行复制

attr_id = H5Aopen(file_id, projection_para, H5P_DEFAULT);
attr_dtype = H5Aget_type(attr_id);
size_t string_length = H5Tget_size(attr_dtype);
char *attr_data = new char[1000];
status = H5Aread(attr_id, attr_dtype, attr_data);

std::istringstream iss(attr_data);
std::vector<double> transform(6);
int i = 0;
std::string str;
while (getline(iss, str, ',')) {
    transform[i] = stod(str);
    ++i;
}

类似地，读取名为projection_str的属性，该属性包含投影信息的WKT字符串。

代码语言：javascript代码运行次数：0运行复制

attr_id = H5Aopen(file_id, projection_str, H5P_DEFAULT);
attr_dtype = H5Aget_type(attr_id);
char *attr_data_str = new char[1000];
status = H5Aread(attr_id, attr_dtype, attr_data_str);

到这里，我们就可以对每个数据集调用H5Dopen1将其打开，并使用H5Dread将数据读入向量中

代码语言：javascript代码运行次数：0运行复制

dataset_id = H5Dopen1(file_id, dataset_0);
std::vector<u_int16_t> data_0(size * size);
status = H5Dread(dataset_id, mem_type_id, H5S_ALL, H5S_ALL, H5P_DEFAULT, data_0.data());

// ... 重复上述步骤读取其他数据集 ...

随后，记得关闭HDF5图像以释放资源。

代码语言：javascript代码运行次数：0运行复制

status = H5Fclose(file_id);

接下来，就该gdal库登场了。使用gdal库创建一个新的TIFF文件，并使用RasterIO方法将每个波段的数据写入到TIFF文件中。

代码语言：javascript代码运行次数：0运行复制

GDALDriver *poDriver = GetGDALDriverManager()->GetDriverByame("GTiff");

GDALDataset *poDstDS = poDriver->Create(output__str(), size, size, band_num, GDT_UInt16, nullptr);
u_int16_t *band_data_0 = &data_0[0];
poDstDS->GetRasterBand(1)->RasterIO(GF_Write, 0, 0, size, size, band_data_0, size, size, GDT_UInt16, 0, 0);
// ... 写入其他波段 ...

同时，设置每个波段的oData值为0，同时按照前述从HDF5图像中读取到的信息，设置TIFF图像的地理变换参数和投影信息。

代码语言：javascript代码运行次数：0运行复制

for (int i = 1; i <= band_num; ++i) {
    GDALRasterBand *poBand = poDstDS->GetRasterBand(i);
    if (poBand != nullptr) {
        poBand->SetoDataValue(0);
    }
}

poDstDS->SetGeoTransform(transform.data());
poDstDS->SetProjection(attr_data_str);

GDALClose(poDstDS);

最后，不要忘记关闭hdf5库以释放资源。

代码语言：javascript代码运行次数：0运行复制

status = H5close();

至此，大功告成。

本文参与腾讯云自媒体同步曝光计划，分享自。原始发表：2024-12-15，如有侵权请联系 cloudcommunity@tencent 删除c++datasize数据压缩

#感谢您对电脑配置推荐网 - 最新i3 i5 i7组装电脑配置单推荐报价格的认可，转载请说明来源于"电脑配置推荐网 - 最新i3 i5 i7组装电脑配置单推荐报价格

本文地址：http://www.dnpztj.cn/biancheng/1241926.html

本站网友红河网络营销	14分钟前发表
H5P_DEFAULT); attr_id = H5Aopen(file_id
本站网友常熟万达	30分钟前发表
0); // ... 写入其他波段 ... 同时
本站网友和女朋友聊天话题	29分钟前发表
本方法支持对szip压缩的HDF5图像（例如高分一号卫星遥感影像）加以转换
本站网友三一重工起诉奥巴马	15分钟前发表
band_data_0
本站网友土地出让	23分钟前发表
size
本站网友南阳不孕不育	5分钟前发表
attr_dtype
本站网友巴黎证券交易所	0秒前发表
初始化hdf5库的状态变量
本站网友青岛11中	5分钟前发表
size
本站网友梦喃	11分钟前发表

本站网友牡丹江酒店	18分钟前发表
directory_iterator(h5_path)) { if (entry.path().extension() == ".h5") { std
本站网友台州商业银行	8分钟前发表
band_data_0
本站网友义乌租房信息网	22分钟前发表
attr_data); std
本站网友易聚网	4分钟前发表
'
本站网友乐山二手房网	16分钟前发表
至少可以保证这个代码在Linux下可以正常运行；但能否在Windows中的C++ 环境下也正常运行

C++将h5转tif：支持高分数据等szip压缩的图像

C++将h5转tif：支持高分数据等szip压缩的图像

《C++解锁机器学习特征工程：构建智能数据基石》

《C++巧铸随机森林：开启智能决策新境界》

《智领未来：C++ 与遗传算法在 AI 模型参数优化中的深度融合》

《数据流驱动：C++构建 AI 模型持续学习新范式》