1 #include "ExtractSliceImageFilter.h"
2 #include <cpPlugins/Interface/Image.h>
4 #include <itkExtractImageFilter.h>
6 // -------------------------------------------------------------------------
7 cpPlugins::BasicFilters::ExtractSliceImageFilter::
8 ExtractSliceImageFilter( )
11 this->_AddInput( "Input" );
12 this->_AddOutput< cpPlugins::Interface::Image >( "Output" );
14 this->m_Parameters->ConfigureAsUint( "Axis" );
15 this->m_Parameters->ConfigureAsInt( "Slice" );
17 this->m_Parameters->SetUint( "Axis", 0 );
18 this->m_Parameters->SetInt( "Slice", 0 );
21 // -------------------------------------------------------------------------
22 cpPlugins::BasicFilters::ExtractSliceImageFilter::
23 ~ExtractSliceImageFilter( )
27 // -------------------------------------------------------------------------
28 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
31 auto image = this->GetInputData< cpPlugins::Interface::Image >( "Input" );
32 itk::DataObject* itk_image = NULL;
35 cpPlugins_Image_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
36 else cpPlugins_VectorImage_Demangle_AllTypes( 3, image, itk_image, r, _GD0 );
37 else*/ r = "ExtractSliceImageFilter: Input image type not supported.";
41 // -------------------------------------------------------------------------
43 std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
44 _GD0( itk::DataObject* image )
46 return( this->_RealGD< I, I >( image ) );
49 // -------------------------------------------------------------------------
50 template< class I, class O >
51 inline std::string cpPlugins::BasicFilters::ExtractSliceImageFilter::
52 _RealGD( itk::DataObject* image )
54 typedef itk::ExtractImageFilter< I, O > _F;
55 typedef typename O::PixelType _OP;
58 int axis = this->m_Parameters->GetUint( "Axis" );
59 int slice = this->m_Parameters->GetInt( "Slice" );
62 I* img = dynamic_cast< I* >( image );
63 typename I::RegionType region = img->GetRequestedRegion( );
64 typename I::SizeType size = region.GetSize( );
65 typename I::IndexType index = region.GetIndex( );
67 index[ axis ] = slice;
68 region.SetSize( size );
69 region.SetIndex( index );
72 _F* filter = this->_CreateITK< _F >( );
73 filter->SetInput( img );
74 filter->SetExtractionRegion( region );
75 filter->SetDirectionCollapseToIdentity( );
79 auto out = this->GetOutputData< cpPlugins::Interface::Image >( "Output" );
80 out->SetITK( filter->GetOutput( ) );