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