1 #include "OtsuThresholdImageFilter.h"
2 #include <cpPlugins/Interface/Image.h>
4 #include <itkOtsuThresholdImageFilter.h>
6 // -------------------------------------------------------------------------
7 cpPlugins::BasicFilters::OtsuThresholdImageFilter::
8 OtsuThresholdImageFilter( )
11 this->_AddInput( "Input" );
12 this->_AddOutput< cpPlugins::Interface::Image >( "Output" );
14 this->m_Parameters->ConfigureAsUint( "NumberOfHistogramBins" );
15 this->m_Parameters->ConfigureAsUint( "InsideValue" );
16 this->m_Parameters->ConfigureAsUint( "OutsideValue" );
18 this->m_Parameters->SetUint( "NumberOfHistogramBins", 100 );
19 this->m_Parameters->SetUint( "InsideValue", 255 );
20 this->m_Parameters->SetUint( "OutsideValue", 0 );
23 // -------------------------------------------------------------------------
24 cpPlugins::BasicFilters::OtsuThresholdImageFilter::
25 ~OtsuThresholdImageFilter( )
29 // -------------------------------------------------------------------------
30 std::string cpPlugins::BasicFilters::OtsuThresholdImageFilter::
33 cpPlugins::Interface::Image* image =
34 this->GetInput< cpPlugins::Interface::Image >( "Input" );
36 return( "OtsuThresholdImageFilter: No input image." );
38 itk::DataObject* itk_image = NULL;
40 cpPlugins_Image_Demangle_AllScalarTypes( 2, image, itk_image, r, _GD0 );
41 else cpPlugins_Image_Demangle_AllScalarTypes( 3, image, itk_image, r, _GD0 );
42 else cpPlugins_Image_Demangle_AllScalarTypes( 4, image, itk_image, r, _GD0 );
43 else r = "OtsuThresholdImageFilter: Input image type not supported.";
47 // -------------------------------------------------------------------------
49 std::string cpPlugins::BasicFilters::OtsuThresholdImageFilter::
50 _GD0( itk::DataObject* image )
53 this->_RealGD< I, itk::Image< unsigned char, I::ImageDimension > >(
59 // -------------------------------------------------------------------------
60 template< class I, class O >
61 inline std::string cpPlugins::BasicFilters::OtsuThresholdImageFilter::
62 _RealGD( itk::DataObject* image )
64 typedef itk::OtsuThresholdImageFilter< I, O > _F;
65 typedef typename O::PixelType _OP;
68 unsigned int bins = this->m_Parameters->GetUint( "NumberOfHistogramBins" );
69 _OP in_val = _OP( this->m_Parameters->GetUint( "InsideValue" ) );
70 _OP out_val = _OP( this->m_Parameters->GetUint( "OutsideValue" ) );
73 _F* filter = this->_CreateITK< _F >( );
74 filter->SetInput( dynamic_cast< I* >( image ) );
75 filter->SetNumberOfHistogramBins( bins );
76 filter->SetInsideValue( out_val ); // WARNING: these are inverted
77 filter->SetOutsideValue( in_val );
81 cpPlugins::Interface::Image* out =
82 this->GetOutput< cpPlugins::Interface::Image >( "Output" );
85 out->SetITK< O >( filter->GetOutput( ) );
89 return( "OtsuThresholdImageFilter: output not correctly created." );