#include "DiscontinuityMapImageFilter.h"
#include <cpInstances/DataObjects/Image.h>

#include <cpExtensions/Algorithms/DiscontinuityMapImageFilter.h>

// -------------------------------------------------------------------------
cpPluginscpExtensions::DiscontinuityMapImageFilter::
DiscontinuityMapImageFilter( )
  : Superclass( )
{
  typedef cpInstances::DataObjects::Image _TImage;

  this->_ConfigureInput< _TImage >( "Input", true, false );
  this->_ConfigureOutput< _TImage >( "Output" );
  this->m_Parameters.ConfigureAsInt( "Radius", 1 );
  this->m_Parameters.ConfigureAsBool( "UseSquareRoot", false );
}

// -------------------------------------------------------------------------
cpPluginscpExtensions::DiscontinuityMapImageFilter::
~DiscontinuityMapImageFilter( )
{
}

// -------------------------------------------------------------------------
void cpPluginscpExtensions::DiscontinuityMapImageFilter::
_GenerateData( )
{
  auto o = this->GetInputData( "Input" );
  cpPlugins_Demangle_Image_ScalarPixels_AllDims_1( o, _GD0 )
    this->_Error( "Invalid input image pixel type." );
}

// -------------------------------------------------------------------------
template< class _TInputImage >
void cpPluginscpExtensions::DiscontinuityMapImageFilter::
_GD0( _TInputImage* input )
{
  typedef cpExtensions::Algorithms::DiscontinuityMapImageFilter< _TInputImage > _TFilter;

  auto filter = this->_CreateITK< _TFilter >( );
  filter->SetInput( input );
  filter->SetRadius( this->m_Parameters.GetInt( "Radius" ) );
  filter->SetUseSquareRoot( this->m_Parameters.GetBool( "UseSquareRoot" ) );
  filter->Update( );
  this->GetOutput( "Output" )->SetITK( filter->GetOutput( ) );
}

// eof - $RCSfile$