plugins/ImageArithmeticFilters/DivideOrZeroOutImageFilter.cxx

   1 #include <plugins/ImageArithmeticFilters/DivideOrZeroOutImageFilter.h>
   2 #include <cpPlugins/DataObjects/Image.h>
   3
   4 #include <itkDivideOrZeroOutImageFilter.h>
   5 #include <itkBinaryFunctorImageFilter.hxx>
   6
   7 // -------------------------------------------------------------------------
   8 cpPluginsImageArithmeticFilters::DivideOrZeroOutImageFilter::
   9 DivideOrZeroOutImageFilter( )
  10   : Superclass( )
  11 {
  12   this->_ConfigureInput< cpPlugins::DataObjects::Image >( "Input0", true, false );
  13   this->_ConfigureInput< cpPlugins::DataObjects::Image >( "Input1", true, false );
  14   this->_ConfigureOutput< cpPlugins::DataObjects::Image >( "Output" );
  15
  16   this->m_Parameters.ConfigureAsReal( "Threshold" );
  17   this->m_Parameters.SetReal( "Threshold", 1e-10 );
  18 }
  19
  20 // -------------------------------------------------------------------------
  21 cpPluginsImageArithmeticFilters::DivideOrZeroOutImageFilter::
  22 ~DivideOrZeroOutImageFilter( )
  23 {
  24 }
  25
  26 // -------------------------------------------------------------------------
  27 void cpPluginsImageArithmeticFilters::DivideOrZeroOutImageFilter::
  28 _GenerateData( )
  29 {
  30   auto o = this->GetInputData( "Input0" );
  31   cpPlugins_Demangle_ImageScalars_Dims( o, _GD0 );
  32   else this->_Error( "Invalid input image." );
  33 }
  34
  35 // -------------------------------------------------------------------------
  36 template< class _TImage >
  37 void cpPluginsImageArithmeticFilters::DivideOrZeroOutImageFilter::
  38 _GD0( _TImage* image0 )
  39 {
  40   typedef itk::DivideOrZeroOutImageFilter< _TImage, _TImage, _TImage > _TFilter;
  41
  42   auto image1 = this->GetInputData< _TImage >( "Input1" );
  43   if( image1 == NULL )
  44     this->_Error( "Incompatible second input image." );
  45
  46   // Configure filter
  47   auto filter = this->_CreateITK< _TFilter >( );
  48   filter->SetInput( 0, image0 );
  49   filter->SetInput( 1, image1 );
  50   filter->SetThreshold( this->m_Parameters.GetReal( "Threshold" ) );
  51   filter->Update( );
  52
  53   // Connect output
  54   this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );
  55 }
  56
  57 // eof - $RCSfile$