0

我正在尝试使用 HDF5 在磁盘上写出 32 位unsigned int到 64 位的向量；unsigned long然后将其读回unsigned int内存中的 32 位。为此，我的读写函数如下所示（我为我的 32 位容器实现了所有定义良好的函数push_back,等）：resizeunsigned intmy_container

bool write(const my_container<unsigned int>& p, const std::string& datapath, const H5::DateType& datatype):
    try {
       const hsize_t h5size[1] = { p.size() };
       const H5::DataSpace h5space(1, h5size);
       const H5::DataSet h5set = fileptr_->createDataSet(datapath, datatype, h5space);
       //-Tried using void* here for data, no difference
       const unsigned int* data = &(*p.begin()); 
       h5set.write(data, datatype);
    } 
    catch( H5::Exception& ) {
       //-Handle exception here
       return false;
    }
    return true;
}

read(my_container<unsigned int>& p, const H5::DataType& datatype, const std::string& datapath) {
    H5::DataType h5set_datatype = h5set.getDataType();
    const std::size_t size = (std::size_t)h5space.getSimpleExtentNpoints();
    try {
       if(h5set_datatype == H5::PredType::NATIVE_UINT64 && datatype == H5::PredType::NATIVE_UINT32 ) {
           typedef unsigned long long u64; 
           typedef std::vector<u64> u64vec;
           u64vec ivector;
           ivector.resize(size);
           void* data = (void*)(&(*ivector.begin()));
           h5set.read(data, h5set_datatype);
           p.resize(0);
           BOOST_FOREACH(const u64 &v, ivector) {
              //-I've handled the cast below using numeric cast separately
              p.push_back(v); 
           }
       } //-End compare datatypes on disk and memory
    } //-End try
    catch(const H5::Exception &e) {
        //-Handle exception
        return false;
    }
    return true;
}

我write用参数调用：const-reference to my_container,H5::Pred::NATIVE_UINT64和read参数：reference to my_containerand H5::Pred::NATIVE_UINT32。这可能是问题的根源之一。让我知道是否需要进一步澄清。基本上，当我重新读回它时，我得到了垃圾。感谢任何 HDF5 专家的建议。谢谢你的时间。

4

1 回答 1

0

解决方案在于更改写入函数以接受文件和内存数据类型：

bool write(const my_container<unsigned int>& p, const H5::DataType& file_datatype, const H5::DataType& mem_datatype, const std::string& datapath) const {

    try {
        const hsize_t h5size[1] = { p.size() };
        const H5::DataSpace h5space(1, h5size);

        const H5::DataSet h5set = fileptr_->createDataSet(datapath, file_datatype, h5space);
        const void* data = &(*p.begin());
        h5set.write(data, mem_datatype);
    } 
    catch( H5::Exception& ) {
        // Handle exception
        return false;
    }

    return true;
}

然后其余的按预期工作-读取功能基本上没有变化；即使对于成对的无符号整数等容器也能做到这一点。HTH。

于 2013-06-25T19:27:58.763 回答