1 #include "SignedMaurerDistanceMapImageFilter.h"
2 #include <cpPlugins/Interface/Image.h>
4 #include <itkSignedMaurerDistanceMapImageFilter.h>
6 // -------------------------------------------------------------------------
7 cpPlugins::BasicFilters::SignedMaurerDistanceMapImageFilter::
8 SignedMaurerDistanceMapImageFilter( )
11 this->_AddInput( "Input" );
12 this->_MakeOutput< cpPlugins::Interface::Image >( "Output" );
14 this->m_Parameters->ConfigureAsReal( "BackgroundValue" );
15 this->m_Parameters->ConfigureAsBool( "InsideIsPositive" );
16 this->m_Parameters->ConfigureAsBool( "SquaredDistance" );
17 this->m_Parameters->ConfigureAsBool( "UseImageSpacing" );
19 std::vector< std::string > choices;
20 choices.push_back( "float" );
21 choices.push_back( "double" );
22 this->m_Parameters->ConfigureAsChoices( "OutputResolution", choices );
24 this->m_Parameters->SetReal( "BackgroundValue", 0 );
25 this->m_Parameters->SetBool( "InsideIsPositive", true );
26 this->m_Parameters->SetBool( "SquaredDistance", false );
27 this->m_Parameters->SetBool( "UseImageSpacing", true );
28 this->m_Parameters->SetSelectedChoice( "OutputResolution", "float" );
31 // -------------------------------------------------------------------------
32 cpPlugins::BasicFilters::SignedMaurerDistanceMapImageFilter::
33 ~SignedMaurerDistanceMapImageFilter( )
37 // -------------------------------------------------------------------------
38 std::string cpPlugins::BasicFilters::SignedMaurerDistanceMapImageFilter::
41 cpPlugins::Interface::Image* image =
42 this->GetInput< cpPlugins::Interface::Image >( "Input" );
44 return( "SignedMaurerDistanceMapImageFilter: No input image." );
46 itk::DataObject* itk_image = NULL;
48 cpPlugins_Image_Demangle_AllScalarTypes( 2, image, itk_image, r, _GD0 );
49 else cpPlugins_Image_Demangle_AllScalarTypes( 3, image, itk_image, r, _GD0 );
50 else r = "SignedMaurerDistanceMapImageFilter: Input image type not supported.";
54 // -------------------------------------------------------------------------
56 std::string cpPlugins::BasicFilters::SignedMaurerDistanceMapImageFilter::
57 _GD0( itk::DataObject* image )
60 this->m_Parameters->GetSelectedChoice( "OutputResolution" );
61 if( out_res == "float" )
63 this->_RealGD< I, itk::Image< float, I::ImageDimension > >(
67 else if( out_res == "double" )
69 this->_RealGD< I, itk::Image< double, I::ImageDimension > >(
74 return( "SignedMaurerDistanceMapImageFilter: Output resolution not supported." );
77 // -------------------------------------------------------------------------
78 template< class I, class O >
79 inline std::string cpPlugins::BasicFilters::SignedMaurerDistanceMapImageFilter::
80 _RealGD( itk::DataObject* image )
82 typedef itk::SignedMaurerDistanceMapImageFilter< I, O > _F;
85 double back_value = this->m_Parameters->GetReal( "BackgroundValue" );
86 bool pos_inside = this->m_Parameters->GetBool( "InsideIsPositive" );
87 bool sqr_dist = this->m_Parameters->GetBool( "SquaredDistance" );
88 bool use_spac = this->m_Parameters->GetBool( "UseImageSpacing" );
91 _F* filter = this->_CreateITK< _F >( );
92 filter->SetInput( dynamic_cast< I* >( image ) );
93 filter->SetBackgroundValue( ( typename I::PixelType )( back_value ) );
94 filter->SetInsideIsPositive( pos_inside );
95 filter->SetSquaredDistance( sqr_dist );
96 filter->SetUseImageSpacing( use_spac );
101 cpPlugins::Interface::Image* out =
102 this->GetOutput< cpPlugins::Interface::Image >( "Output" );
105 out->SetITK< O >( filter->GetOutput( ) );
109 return( "SignedMaurerDistanceMapImageFilter: output not correctly created." );