1 #include "ExtractSliceImageFilter.h"
2 #include <cpPlugins/Interface/Image.h>
4 #include <itkExtractImageFilter.h>
6 // -------------------------------------------------------------------------
7 cpPlugins::BasicFilters::ExtractSliceImageFilter::
8 ExtractSliceImageFilter( )
11 this->m_ClassName = "cpPlugins::BasicFilters::ExtractSliceImageFilter";
12 this->m_ClassCategory = "ImageToImageFilter";
13 this->SetNumberOfInputs( 1 );
14 this->SetNumberOfOutputs( 1 );
15 this->_MakeOutput< cpPlugins::Interface::Image >( 0 );
17 using namespace cpPlugins::Interface;
18 this->m_DefaultParameters.Configure( Parameters::Int, "Axis" );
19 this->m_DefaultParameters.Configure( Parameters::Int, "Slice" );
20 this->m_DefaultParameters.SetValueAsInt( "Axis", 0 );
21 this->m_DefaultParameters.SetValueAsInt( "Slice", 0 );
22 this->m_Parameters = this->m_DefaultParameters;
25 // -------------------------------------------------------------------------
26 cpPlugins::BasicFilters::ExtractSliceImageFilter::
27 ~ExtractSliceImageFilter( )
31 // -------------------------------------------------------------------------
32 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
35 cpPlugins::Interface::Image* image =
36 this->GetInput< cpPlugins::Interface::Image >( 0 );
38 return( "ExtractSliceImageFilter: No input image." );
40 itk::DataObject* itk_image = NULL;
43 cpPlugins_Image_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
44 else cpPlugins_VectorImage_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
45 else*/ r = "ExtractSliceImageFilter: Input image type not supported.";
49 // -------------------------------------------------------------------------
51 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
52 _GD0( itk::DataObject* image )
54 return( this->_RealGD< I, I >( image ) );
57 // -------------------------------------------------------------------------
58 template< class I, class O >
59 inline std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
60 _RealGD( itk::DataObject* image )
62 typedef itk::ExtractImageFilter< I, O > _F;
63 typedef typename O::PixelType _OP;
66 int axis = this->m_Parameters.GetValueAsInt( "Axis" );
67 int slice = this->m_Parameters.GetValueAsInt( "Slice" );
69 std::cout << "HOLA: " << slice << std::endl;
73 I* img = dynamic_cast< I* >( image );
74 typename I::RegionType region = img->GetRequestedRegion( );
75 typename I::SizeType size = region.GetSize( );
76 typename I::IndexType index = region.GetIndex( );
78 index[ axis ] = slice;
79 region.SetSize( size );
80 region.SetIndex( index );
82 std::cout << "HOLA-: " << region << std::endl;
85 _F* filter = this->_CreateITK< _F >( );
86 filter->SetInput( img );
87 filter->SetExtractionRegion( region );
88 filter->SetDirectionCollapseToIdentity( );
92 cpPlugins::Interface::Image* out =
93 this->GetOutput< cpPlugins::Interface::Image >( 0 );
96 out->SetITK< O >( filter->GetOutput( ) );
100 return( "ExtractSliceImageFilter: output not correctly created." );