1 #include <ImageGradientFilters/MultiScaleGaussianImageFilter.h>
2 #include <cpPlugins/DataObjects/Image.h>
3 #include <cpPlugins/DataObjects/Image_Demanglers.h>
5 #include <cpExtensions/Algorithms/MultiScaleGaussianImageFilter.h>
8 #include <cpExtensions/Algorithms/MultiScaleGaussianImageFilter.hxx>
9 #include <itkGradientRecursiveGaussianImageFilter.hxx>
10 #include <itkRecursiveGaussianImageFilter.hxx>
11 #include <itkImageAdaptor.hxx>
12 #include <itkUnaryFunctorImageFilter.hxx>
13 #include <itkImageConstIteratorWithIndex.hxx>
14 #include <itkBinaryFunctorImageFilter.hxx>
15 #include <itkImageRegionIteratorWithIndex.hxx>
16 #include <itkInPlaceImageFilter.hxx>
17 #include <itkImageRegionConstIteratorWithIndex.hxx>
18 #include <itkSimpleDataObjectDecorator.hxx>
21 // -------------------------------------------------------------------------
22 cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
23 MultiScaleGaussianImageFilter( )
26 this->_ConfigureInput< cpPlugins::DataObjects::Image >( "Input", true, false );
27 this->_ConfigureOutput< cpPlugins::DataObjects::Image >( "Output" );
29 this->m_Parameters.ConfigureAsRealList( "Sigmas" );
31 std::vector< std::string > choices;
32 choices.push_back( "float" );
33 choices.push_back( "double" );
34 this->m_Parameters.ConfigureAsChoices( "ScalarType", choices );
35 this->m_Parameters.SetSelectedChoice( "ScalarType", "float" );
38 // -------------------------------------------------------------------------
39 cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
40 ~MultiScaleGaussianImageFilter( )
44 // -------------------------------------------------------------------------
45 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
48 auto o = this->GetInputData( "Input" );
49 cpPlugins_Demangle_Image_ScalarPixels_AllDims_1( o, _GD0 )
50 this->_Error( "Invalid input image." );
53 // -------------------------------------------------------------------------
54 template< class _TImage >
55 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
56 _GD0( _TImage* image )
58 auto choice = this->m_Parameters.GetSelectedChoice( "ScalarType" );
59 if( choice == "float" ) this->_GD1< _TImage, float >( image );
60 else if( choice == "double" ) this->_GD1< _TImage, double >( image );
61 else this->_Error( "Invalid scalar type." );
64 // -------------------------------------------------------------------------
65 template< class _TInputImage, class _TScalar >
66 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
67 _GD1( _TInputImage* image )
69 typedef itk::CovariantVector< _TScalar, _TInputImage::ImageDimension > _TGradient;
70 typedef itk::Image< _TGradient, _TInputImage::ImageDimension > _TOutputImage;
72 cpExtensions::Algorithms::MultiScaleGaussianImageFilter< _TInputImage, _TOutputImage >
75 auto sigmas = this->m_Parameters.GetRealList( "Sigmas" );
76 if( sigmas.size( ) == 0 )
77 this->_Error( "No given sigmas." );
80 _TFilter* filter = this->_CreateITK< _TFilter >( );
81 filter->SetInput( image );
82 for( auto sIt = sigmas.begin( ); sIt != sigmas.end( ); ++sIt )
83 filter->AddScale( *sIt );
87 this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );