1 #include "ExtractSliceImageFilter.h"
2 #include <cpPlugins/Interface/Image.h>
4 #include <itkExtractImageFilter.h>
6 // -------------------------------------------------------------------------
7 cpPlugins::BasicFilters::ExtractSliceImageFilter::
8 ExtractSliceImageFilter( )
11 this->SetNumberOfInputs( 1 );
12 this->SetNumberOfOutputs( 1 );
13 this->_MakeOutput< cpPlugins::Interface::Image >( 0 );
15 this->m_Parameters->ConfigureAsUint( "Axis", 0 );
16 this->m_Parameters->ConfigureAsInt( "Slice", 0 );
19 // -------------------------------------------------------------------------
20 cpPlugins::BasicFilters::ExtractSliceImageFilter::
21 ~ExtractSliceImageFilter( )
25 // -------------------------------------------------------------------------
26 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
29 cpPlugins::Interface::Image* image =
30 this->GetInput< cpPlugins::Interface::Image >( 0 );
32 return( "ExtractSliceImageFilter: No input image." );
34 itk::DataObject* itk_image = NULL;
37 cpPlugins_Image_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
38 else cpPlugins_VectorImage_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
39 else*/ r = "ExtractSliceImageFilter: Input image type not supported.";
43 // -------------------------------------------------------------------------
45 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
46 _GD0( itk::DataObject* image )
48 return( this->_RealGD< I, I >( image ) );
51 // -------------------------------------------------------------------------
52 template< class I, class O >
53 inline std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
54 _RealGD( itk::DataObject* image )
56 typedef itk::ExtractImageFilter< I, O > _F;
57 typedef typename O::PixelType _OP;
60 int axis = this->m_Parameters->GetUint( "Axis" );
61 int slice = this->m_Parameters->GetInt( "Slice" );
64 I* img = dynamic_cast< I* >( image );
65 typename I::RegionType region = img->GetRequestedRegion( );
66 typename I::SizeType size = region.GetSize( );
67 typename I::IndexType index = region.GetIndex( );
69 index[ axis ] = slice;
70 region.SetSize( size );
71 region.SetIndex( index );
74 _F* filter = this->_CreateITK< _F >( );
75 filter->SetInput( img );
76 filter->SetExtractionRegion( region );
77 filter->SetDirectionCollapseToIdentity( );
81 cpPlugins::Interface::Image* out =
82 this->GetOutput< cpPlugins::Interface::Image >( 0 );
85 out->SetITK< O >( filter->GetOutput( ) );
89 return( "ExtractSliceImageFilter: output not correctly created." );