1 #include <ImageArithmeticFilters/PowImageFilter.h>
2 #include <cpPlugins/DataObjects/Image.h>
3 #include <cpPlugins/DataObjects/Image_Demanglers.h>
5 #include <itkPowImageFilter.h>
7 // -------------------------------------------------------------------------
8 cpPluginsImageArithmeticFilters::PowImageFilter::
12 typedef cpPlugins::DataObjects::Image _TImage;
13 this->_ConfigureInput< _TImage >( "Input1", true, false );
14 this->_ConfigureInput< _TImage >( "Input2", true, false );
15 this->_ConfigureOutput< _TImage >( "Output" );
17 this->m_Parameters.ConfigureAsReal( "Constant", 1 );
20 // -------------------------------------------------------------------------
21 cpPluginsImageArithmeticFilters::PowImageFilter::
26 // -------------------------------------------------------------------------
27 void cpPluginsImageArithmeticFilters::PowImageFilter::
30 auto o = this->GetInputData( "Input1" );
31 cpPlugins_Demangle_Image_ScalarPixels_AllDims_1( o, _GD0 )
32 this->_Error( "Invalid input image (0)." );
35 // -------------------------------------------------------------------------
36 template< class _TInput1 >
37 void cpPluginsImageArithmeticFilters::PowImageFilter::
38 _GD0( _TInput1* input1 )
40 auto input2 = this->GetInputData< _TInput1 >( "Input2" );
42 this->_Error( "Incompatible second input image." );
43 this->_GD1( input2, input1 );
46 // -------------------------------------------------------------------------
47 template< class _TInput2, class _TInput1 >
48 void cpPluginsImageArithmeticFilters::PowImageFilter::
49 _GD1( _TInput2* input2, _TInput1* input1 )
51 typedef itk::PowImageFilter< _TInput1, _TInput2, _TInput1 > _TFilter;
54 auto filter = this->_CreateITK< _TFilter >( );
55 filter->SetInput1( input1 );
56 filter->SetInput2( input2 );
60 this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );
66 filter->SetInput1( image0 );
67 auto image1 = this->GetInputData< _TImage >( "Input1" );
69 filter->SetConstant( this->m_Parameters.GetReal( "Constant" ) );
71 filter->SetInput2( image1 );
75 this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );