ITK／VTK对DICOM文件的读取

   typedef signed short InputPixelType;

    
   const unsigned int   InputDimension = 2;
    
   typedef itk::Image< InputPixelType, InputDimension > InputImageType;
    
   typedef itk::ImageFileReader< InputImageType > ReaderType;
    
   ReaderType::Pointer reader = ReaderType::New();
    
   reader->SetFileName( argv[1] );
    
   typedef itk::GDCMImageIO           ImageIOType;
    
   ImageIOType::Pointer gdcmImageIO = ImageIOType::New();
    
   reader->SetImageIO( gdcmImageIO );
    
   try
    
     {
    
     reader->Update();
    
     }
    
   catch (itk::ExceptionObject & e)
    
     {
    
     std::cerr << "exception in file reader " << std::endl;
    
     std::cerr << e << std::endl;
    
     return EXIT_FAILURE;
    
     }

  const unsigned int InputDimension = 3;

    
   const unsigned int OutputDimension = 2;
    
  
    
   typedef signed short PixelType;
    
  
    
   typedef itk::Image< PixelType, InputDimension >
    
     InputImageType;
    
   typedef itk::Image< PixelType, OutputDimension >
    
     OutputImageType;
    
   typedef itk::ImageSeriesReader< InputImageType >
    
     ReaderType;
    
   typedef itk::GDCMImageIO
    
     ImageIOType;
    
   typedef itk::GDCMSeriesFileNames
    
     InputNamesGeneratorType;
    
   typedef itk::NumericSeriesFileNames
    
     OutputNamesGeneratorType;
    
   typedef itk::IdentityTransform< double, InputDimension >
    
     TransformType;
    
   typedef itk::LinearInterpolateImageFunction< InputImageType, double >
    
     InterpolatorType;
    
   typedef itk::ResampleImageFilter< InputImageType, InputImageType >
    
     ResampleFilterType;
    
 #if ( ( ITK_VERSION_MAJOR == 4 ) && ( ITK_VERSION_MINOR < 6 ) )
    
   typedef itk::ShiftScaleImageFilter< InputImageType, InputImageType >
    
     ShiftScaleType;
    
 #endif
    
   typedef itk::ImageSeriesWriter< InputImageType, OutputImageType >
    
     SeriesWriterType;
    
  
    
  
    
 // Read the input series
    
  
    
   ImageIOType::Pointer gdcmIO = ImageIOType::New();
    
   InputNamesGeneratorType::Pointer inputNames = InputNamesGeneratorType::New();
    
   inputNames->SetInputDirectory( argv[1] );
    
  
    
   const ReaderType::FileNamesContainer & filenames =
    
                         inputNames->GetInputFileNames();
    
  
    
   ReaderType::Pointer reader = ReaderType::New();
    
  
    
   reader->SetImageIO( gdcmIO );
    
   reader->SetFileNames( filenames );
    
   try
    
     {
    
     reader->Update();
    
     }
    
   catch (itk::ExceptionObject &excp)
    
     {
    
     std::cerr << "Exception thrown while reading the series" << std::endl;
    
     std::cerr << excp << std::endl;
    
     return EXIT_FAILURE;
    
     }

 <pre name="code" class="cpp">   // Read all the DICOM files in the specified directory.

    
    vtkSmartPointer<vtkDICOMImageReader> reader =
    
       vtkSmartPointer<vtkDICOMImageReader>::New();
    
    reader->SetDirectoryName(folder.c_str());
    
    reader->Update();
    
  
    
    // Visualize
    
    vtkSmartPointer<vtkImageViewer2> imageViewer =
    
       vtkSmartPointer<vtkImageViewer2>::New();
    
    imageViewer->SetInputConnection(reader->GetOutputPort());
    
    imageViewer->SetSliceIndex(1);