plugins/cpPluginsImageFilters/MultiScaleGaussianImageFilter.cxx

   1 #include <cpPluginsImageFilters/MultiScaleGaussianImageFilter.h>
   2 #include <cpPlugins/Image.h>
   3 #include <cpPlugins_ITKInstances/ImageFilters.h>
   4
   5 #include <cpExtensions/Algorithms/MultiScaleGaussianImageFilter.h>
   6 #include <cpExtensions/Algorithms/MultiScaleGaussianImageFilter.hxx>
   7 #include <itkGradientRecursiveGaussianImageFilter.hxx>
   8 #include <itkImageAdaptor.hxx>
   9 #include <itkImageToImageFilter.hxx>
  10 #include <itkInPlaceImageFilter.hxx>
  11 #include <itkUnaryFunctorImageFilter.hxx>
  12 #include <itkRecursiveGaussianImageFilter.hxx>
  13 #include <itkImageConstIteratorWithIndex.hxx>
  14 #include <itkRecursiveSeparableImageFilter.hxx>
  15 #include <itkBinaryFunctorImageFilter.hxx>
  16 #include <itkImageRegionConstIterator.hxx>
  17 #include <itkImageScanlineIterator.hxx>
  18 #include <itkImageSource.hxx>
  19 #include <itkImageRegionIteratorWithIndex.hxx>
  20 #include <itkImageScanlineConstIterator.hxx>
  21 #include <itkImageRegionIterator.hxx>
  22 #include <itkSimpleDataObjectDecorator.hxx>
  23 #include <itkImageRegionConstIteratorWithIndex.hxx>
  24
  25 // -------------------------------------------------------------------------
  26 cpPluginsImageFilters::MultiScaleGaussianImageFilter::
  27 MultiScaleGaussianImageFilter( )
  28   : Superclass( )
  29 {
  30   this->_AddInput( "Input" );
  31   this->_AddOutput< cpPlugins::Image >( "Output" );
  32
  33   this->m_Parameters.ConfigureAsRealList( "Sigmas" );
  34
  35   std::vector< std::string > choices;
  36   choices.push_back( "float" );
  37   choices.push_back( "double" );
  38   this->m_Parameters.ConfigureAsChoices( "ScalarType", choices );
  39   this->m_Parameters.SetSelectedChoice( "ScalarType", "float" );
  40 }
  41
  42 // -------------------------------------------------------------------------
  43 cpPluginsImageFilters::MultiScaleGaussianImageFilter::
  44 ~MultiScaleGaussianImageFilter( )
  45 {
  46 }
  47
  48 // -------------------------------------------------------------------------
  49 std::string cpPluginsImageFilters::MultiScaleGaussianImageFilter::
  50 _GenerateData( )
  51 {
  52   auto image = this->GetInputData( "Input" )->GetITK< itk::DataObject >( );
  53   std::string   cpPlugin_Image_Demangle_Pixel_AllScalars( r, _GD0, image, 2 );
  54   if( r != "" ) cpPlugin_Image_Demangle_Pixel_AllScalars( r, _GD0, image, 3 );
  55   return( r );
  56 }
  57
  58 // -------------------------------------------------------------------------
  59 template< class _TImage >
  60 std::string cpPluginsImageFilters::MultiScaleGaussianImageFilter::
  61 _GD0( _TImage* image )
  62 {
  63   if( image != NULL )
  64   {
  65     auto choice = this->m_Parameters.GetSelectedChoice( "ScalarType" );
  66     if( choice == "float" )
  67       return( this->_GD1< _TImage, float >( image ) );
  68     else if( choice == "double" )
  69       return( this->_GD1< _TImage, double >( image ) );
  70     else return( "MultiScaleGaussianImageFilter: no valid scalar type." );
  71   }
  72   else
  73     return(
  74       "MultiScaleGaussianImageFilter: No valid input image."
  75       );
  76 }
  77
  78 // -------------------------------------------------------------------------
  79 template< class _TImage, class _TScalar >
  80 std::string cpPluginsImageFilters::MultiScaleGaussianImageFilter::
  81 _GD1( _TImage* image )
  82 {
  83   typedef itk::CovariantVector< _TScalar, _TImage::ImageDimension > _TVector;
  84   typedef itk::Image< _TVector, _TImage::ImageDimension > _TGradient;
  85   typedef
  86     cpExtensions::Algorithms::MultiScaleGaussianImageFilter< _TImage, _TGradient >
  87     _TFilter;
  88
  89   auto sigmas = this->m_Parameters.GetRealList( "Sigmas" );
  90   if( sigmas.size( ) == 0 )
  91     return(
  92       "MultiScaleGaussianImageFilter: No given sigmas."
  93       );
  94
  95   // Configure filter
  96   _TFilter* filter = this->_CreateITK< _TFilter >( );
  97   filter->SetInput( image );
  98   for( auto sIt = sigmas.begin( ); sIt != sigmas.end( ); ++sIt )
  99     filter->AddScale( *sIt );
 100   filter->Update( );
 101
 102   // Connect output
 103   this->GetOutputData( "Output" )->SetITK( filter->GetOutput( ) );
 104   return( "" );
 105 }
 106
 107 // eof - $RCSfile$