c++版本itk读取dicom序列

最新推荐文章于 2023-11-09 16:51:49 发布

吴天德少侠

最新推荐文章于 2023-11-09 16:51:49 发布

阅读量180

点赞数 1

文章标签： c++ 开发语言

本文链接：https://blog.csdn.net/sdhdsf132452/article/details/132560955

版权

医学图像处理1 专栏收录该内容

37 篇文章 23 订阅

订阅专栏

// 读取dicom序列.hablee
void readDicomSeries(const std::string& dicom_dir, itk::Image<float, 3>::Pointer& out_image)
{
    using ImageIOType = itk::GDCMImageIO;
    ImageIOType::Pointer DicomIO = ImageIOType::New();

    using ImageReaderType3D = itk::ImageSeriesReader<itk::Image<float, 3>>;
    ImageReaderType3D::Pointer ImageReader3D = ImageReaderType3D::New();

    ImageReader3D->SetImageIO(DicomIO);

    using NamesGeneratorType = itk::GDCMSeriesFileNames;
    NamesGeneratorType::Pointer NameGenerator = NamesGeneratorType::New();
    NameGenerator->SetDirectory(dicom_dir);

    using SeriesIdContainer = std::vector<std::string>;
    const SeriesIdContainer& SeriesUID = NameGenerator->GetSeriesUIDs();
    SeriesIdContainer::const_iterator SeriesItr_begin = SeriesUID.begin();
    SeriesIdContainer::const_iterator SeriesItr_end = SeriesUID.end();

    while (SeriesItr_begin != SeriesItr_end)
    {
        std::cout << SeriesItr_begin->c_str() << std::endl;
        SeriesItr_begin++;
    }

    std::string SeriesIdentifier;
    SeriesIdentifier = SeriesUID.begin()->c_str();
    using FileNamesContainer = std::vector<std::string>;
    FileNamesContainer filenames;
    filenames = NameGenerator->GetFileNames(SeriesIdentifier);

    ImageReader3D->SetFileNames(filenames);
    try
    {
        ImageReader3D->Update();
    }
    catch (itk::ExceptionObject& err)
    {
        std::cerr << "error, exception object caught!" << std::endl;
    }
    out_image = ImageReader3D->GetOutput();
}