lib/cpPlugins/Plugins/BasicFilters/ExtractSliceImageFilter.cxx

   1 #include "ExtractSliceImageFilter.h"
   2 #include <cpPlugins/Interface/Image.h>
   3
   4 #include <itkExtractImageFilter.h>
   5
   6 // -------------------------------------------------------------------------
   7 cpPlugins::BasicFilters::ExtractSliceImageFilter::
   8 ExtractSliceImageFilter( )
   9   : Superclass( )
  10 {
  11   this->_AddInput( "Input" );
  12   this->_AddOutput< cpPlugins::Interface::Image >( "Output" );
  13
  14   this->m_Parameters->ConfigureAsUint( "Axis" );
  15   this->m_Parameters->ConfigureAsInt( "Slice" );
  16
  17   this->m_Parameters->SetUint( "Axis", 0 );
  18   this->m_Parameters->SetInt( "Slice", 0 );
  19 }
  20
  21 // -------------------------------------------------------------------------
  22 cpPlugins::BasicFilters::ExtractSliceImageFilter::
  23 ~ExtractSliceImageFilter( )
  24 {
  25 }
  26
  27 // -------------------------------------------------------------------------
  28 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
  29 _GenerateData( )
  30 {
  31   auto image = this->GetInputData< cpPlugins::Interface::Image >( "Input" );
  32   if( image == NULL )
  33     return( "ExtractSliceImageFilter: No input image." );
  34
  35   itk::DataObject* itk_image = NULL;
  36   std::string r = "";
  37 /*
  38   cpPlugins_Image_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
  39   else cpPlugins_VectorImage_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
  40   else*/ r = "ExtractSliceImageFilter: Input image type not supported.";
  41   return( r );
  42 }
  43
  44 // -------------------------------------------------------------------------
  45 template< class I >
  46 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
  47 _GD0( itk::DataObject* image )
  48 {
  49   return( this->_RealGD< I, I >( image ) );
  50 }
  51
  52 // -------------------------------------------------------------------------
  53 template< class I, class O >
  54 inline std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
  55 _RealGD( itk::DataObject* image )
  56 {
  57   typedef itk::ExtractImageFilter< I, O > _F;
  58   typedef typename O::PixelType _OP;
  59
  60   // Get parameters
  61   int axis = this->m_Parameters->GetUint( "Axis" );
  62   int slice = this->m_Parameters->GetInt( "Slice" );
  63
  64   // Compute region
  65   I* img = dynamic_cast< I* >( image );
  66   typename I::RegionType region = img->GetRequestedRegion( );
  67   typename I::SizeType size = region.GetSize( );
  68   typename I::IndexType index = region.GetIndex( );
  69   size[ axis ] = 1;
  70   index[ axis ] = slice;
  71   region.SetSize( size );
  72   region.SetIndex( index );
  73
  74   // Configure filter
  75   _F* filter = this->_CreateITK< _F >( );
  76   filter->SetInput( img );
  77   filter->SetExtractionRegion( region );
  78   filter->SetDirectionCollapseToIdentity( );
  79   filter->Update( );
  80
  81   // Connect output
  82   auto out = this->GetOutputData< cpPlugins::Interface::Image >( "Output" );
  83   if( out != NULL )
  84   {
  85     out->SetITK< O >( filter->GetOutput( ) );
  86     return( "" );
  87   }
  88   else
  89     return( "ExtractSliceImageFilter: output not correctly created." );
  90 }
  91
  92 // eof - $RCSfile$