plugins/ImageGradientFilters/MultiScaleGaussianImageFilter.cxx

   1 #include <ImageGradientFilters/MultiScaleGaussianImageFilter.h>
   2 #include <cpPlugins/DataObjects/Image.h>
   3 #include <cpPlugins/DataObjects/Image_Demanglers.h>
   4
   5 #include <cpExtensions/Algorithms/MultiScaleGaussianImageFilter.h>
   6
   7 /* TODO
   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>
  19 */
  20
  21 // -------------------------------------------------------------------------
  22 cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
  23 MultiScaleGaussianImageFilter( )
  24   : Superclass( )
  25 {
  26   this->_ConfigureInput< cpPlugins::DataObjects::Image >( "Input", true, false );
  27   this->_ConfigureOutput< cpPlugins::DataObjects::Image >( "Output" );
  28
  29   this->m_Parameters.ConfigureAsRealList( "Sigmas" );
  30
  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" );
  36 }
  37
  38 // -------------------------------------------------------------------------
  39 cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
  40 ~MultiScaleGaussianImageFilter( )
  41 {
  42 }
  43
  44 // -------------------------------------------------------------------------
  45 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
  46 _GenerateData( )
  47 {
  48   auto o = this->GetInputData( "Input" );
  49   cpPlugins_Demangle_Image_ScalarPixels_AllDims_1( o, _GD0 )
  50     this->_Error( "Invalid input image." );
  51 }
  52
  53 // -------------------------------------------------------------------------
  54 template< class _TImage >
  55 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
  56 _GD0( _TImage* image )
  57 {
  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." );
  62 }
  63
  64 // -------------------------------------------------------------------------
  65 template< class _TInputImage, class _TScalar >
  66 void cpPluginsImageGradientFilters::MultiScaleGaussianImageFilter::
  67 _GD1( _TInputImage* image )
  68 {
  69   typedef itk::CovariantVector< _TScalar, _TInputImage::ImageDimension > _TGradient;
  70   typedef itk::Image< _TGradient, _TInputImage::ImageDimension > _TOutputImage;
  71   typedef
  72     cpExtensions::Algorithms::MultiScaleGaussianImageFilter< _TInputImage, _TOutputImage >
  73     _TFilter;
  74
  75   auto sigmas = this->m_Parameters.GetRealList( "Sigmas" );
  76   if( sigmas.size( ) == 0 )
  77     this->_Error( "No given sigmas." );
  78
  79   // Configure filter
  80   _TFilter* filter = this->_CreateITK< _TFilter >( );
  81   filter->SetInput( image );
  82   for( auto sIt = sigmas.begin( ); sIt != sigmas.end( ); ++sIt )
  83     filter->AddScale( *sIt );
  84   filter->Update( );
  85
  86   // Connect output
  87   this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );
  88 }
  89
  90 // eof - $RCSfile$