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