#include <cpPlugins/Plugins/MedianImageFilter.h>
#include <cpPlugins/Interface/Image.h>

#include <itkMedianImageFilter.h>

// -------------------------------------------------------------------------
cpPlugins::Plugins::MedianImageFilter::
MedianImageFilter( )
  : Superclass( )
{
  this->m_ClassName = "cpPlugins::MedianImageFilter";
  this->m_ClassCategory = "ImageToImageFilter";
  this->SetNumberOfInputs( 1 );
  this->SetNumberOfOutputs( 1 );
  this->_MakeOutput< cpPlugins::Interface::Image >( 0 );

  using namespace cpPlugins::Interface;
  this->m_DefaultParameters.Configure(Parameters::Real, "Radius");
 
  this->m_DefaultParameters.SetValueAsReal("Radius", 3);
 
  this->m_Parameters = this->m_DefaultParameters;
}

// -------------------------------------------------------------------------
cpPlugins::Plugins::MedianImageFilter::
~MedianImageFilter( )
{
}

// -------------------------------------------------------------------------
std::string cpPlugins::Plugins::MedianImageFilter::
_GenerateData( )
{
  cpPlugins::Interface::Image* image =
    this->GetInput< cpPlugins::Interface::Image >( 0 );
  if( image == NULL )
    return( "MedianImageFilter: No input image." );

  itk::DataObject* itk_image = NULL;
  std::string r = "";
  cpPlugins_Image_Input_Demangle_Dimension_AllScalarTypes(
    2, image, itk_image, r, _DemangleOutput
    );
  else cpPlugins_Image_Input_Demangle_Dimension_AllScalarTypes(
    3, image, itk_image, r, _DemangleOutput
    );
  else cpPlugins_Image_Input_Demangle_Dimension_AllScalarTypes(
    4, image, itk_image, r, _DemangleOutput
    );
  else r = "MedianImageFilter: Input image type not supported.";
  return( r );
}

// -------------------------------------------------------------------------
template< class I >
std::string cpPlugins::Plugins::MedianImageFilter::
_DemangleOutput( itk::DataObject* image )
{
  return(
    this->_RealGD< I, itk::Image< unsigned char, I::ImageDimension > >(
      image
      )
    );
}

// -------------------------------------------------------------------------
template< class I, class O >
inline std::string cpPlugins::Plugins::MedianImageFilter::
_RealGD( itk::DataObject* image )
{
  typedef itk::MedianImageFilter< I, O > _F;
  typedef typename O::PixelType _OP;

  // Get parameters

  _OP rad_val = _OP(this->m_Parameters.GetValueAsReal("Radius"));
  

  // Configure filter
  _F* filter = dynamic_cast< _F* >( this->m_RealProcessObject.GetPointer( ) );
  if( filter == NULL )
  {
    this->m_RealProcessObject = _F::New( );
    filter = dynamic_cast< _F* >( this->m_RealProcessObject.GetPointer( ) );

  } // fi
  filter->SetInput( dynamic_cast< I* >( image ) );
  
  filter->SetRadius(rad_val);
 
  filter->Update( );

  // Connect output
  cpPlugins::Interface::Image* out =
    this->GetOutput< cpPlugins::Interface::Image >( 0 );
  if( out != NULL )
  {
    out->SetITKImage< O >( filter->GetOutput( ) );
    return( "" );
  }
  else
    return( "MedianImageFilter: output not correctly created." );
}

// eof - $RCSfile$