ITKでの* .mhd/*。rawフォーマットの3D画像の読み込み

の読み込み方法.mhd/ .raw形式の3D画像はITKで表示されますか？次のコードを使用しようとしましたが、ロードされたイメージのディメンションが0,0,0として表示されるため、ロードされません。ITKでの* .mhd/*。rawフォーマットの3D画像の読み込み

私が作っている間違いを指摘できますか？

typedef float InputPixelType; 
const unsigned int DimensionOfRaw = 3; 
typedef itk::Image< InputPixelType, DimensionOfRaw > InputImageType; 

//typedef itk::RawImageIO<InputPixelType, DimensionOfRaw> ImageIOType; 
typedef itk::ImageFileReader<InputImageType > ReaderType; 

/* 
* --------------------Loader and saver of Raws, as well the function that takes a resulting (from inference matrix/vector) and creates a Raw out of it.----------------------- 
*/ 
InputImageType::Pointer loadRawImageItk(std::string RawFullFilepathname, ReaderType::Pointer & RawImageIO) { 
    //http://www.itk.org/Doxygen/html/classitk_1_1Image.html 
    //http://www.itk.org/Doxygen/html/classitk_1_1ImageFileReader.html 

    typedef itk::ImageFileReader<InputImageType> ReaderType; 

    ReaderType::Pointer reader = ReaderType::New(); 
    reader->SetFileName(RawFullFilepathname); 


    //ImageIOType::Pointer RawImageIO = ImageIOType::New(); 
    reader->SetImageIO(RawImageIO); 

    try { 
     reader->Update(); 
    } catch (itk::ExceptionObject& e) { 
     std::cerr << e.GetDescription() << std::endl; 
     exit(1); // You can choose to do something else, of course. 
    } 

    //InputImageType::Pointer inputImage = reader->GetOutput(); 
    InputImageType::Pointer inputImage = reader->GetOutput(); 

    return inputImage; 

} 


int saveRawImageItk(std::string RawFullFilepathname, InputImageType::Pointer & outputImageItkType , ImageIOType::Pointer & RawImageIO) { 
    std::cout << "Saving image to: " << RawFullFilepathname << "\n"; 

    typedef itk::ImageFileWriter<InputImageType> Writer1Type; 
    Writer1Type::Pointer writer1 = Writer1Type::New(); 

    writer1->SetInput(outputImageItkType); 
    writer1->SetFileName(RawFullFilepathname); 
    writer1->SetImageIO(RawImageIO); //seems like this is useless. 

    // Execution of the writer is triggered by invoking the \code{Update()} method. 
    try 
    { 
    writer1->Update(); 
    } 
    catch (itk::ExceptionObject & e) 
    { 
    std::cerr << "exception in file writer " << std::endl; 
    std::cerr << e.GetDescription() << std::endl; 
    std::cerr << e.GetLocation() << std::endl; 
    return 1; 
    } 

    return 0; 
}

出典

2016-06-14 Avijit Dasgupta

多分それはあなたのデータに問題があります。少なくとも、いくつかのサンプルデータを提供する必要があります。 – Trilarion

私は成功し、次のSimpleITKコードを使用してPythonでMHDとRAWファイルを読みました：

import SimpleITK as sitk 
import numpy as np 
def load_itk_image(filename): 
    itkimage = sitk.ReadImage(filename) 
    numpyImage = sitk.GetArrayFromImage(itkimage) 
    return numpyImage

たぶん、あなたはそれを参照として使用することができます。

ImageFileReaderではなくReadImage関数を使用する必要がありますか？試してみることができます。

出典

2016-07-22 02:51:10 Kongsea