#include <cpPluginsImageFilters/SignedMaurerDistanceMapImageFilter.h>
#include <cpPlugins/Image.h>
#include <cpPlugins_Instances/DistanceMapFilters.h>

// -------------------------------------------------------------------------
cpPluginsImageFilters::SignedMaurerDistanceMapImageFilter::
SignedMaurerDistanceMapImageFilter( )
  : Superclass( )
{
  this->_AddInput( "Input" );
  this->_AddOutput< cpPlugins::Image >( "Output" );

  this->m_Parameters.ConfigureAsReal( "BackgroundValue" );
  this->m_Parameters.ConfigureAsBool( "InsideIsPositive" );
  this->m_Parameters.ConfigureAsBool( "SquaredDistance" );
  this->m_Parameters.ConfigureAsBool( "UseImageSpacing" );

  std::vector< std::string > choices;
  choices.push_back( "float" );
  choices.push_back( "double" );
  this->m_Parameters.ConfigureAsChoices( "OutputResolution", choices );

  this->m_Parameters.SetReal( "BackgroundValue", 0 );
  this->m_Parameters.SetBool( "InsideIsPositive", true );
  this->m_Parameters.SetBool( "SquaredDistance", false );
  this->m_Parameters.SetBool( "UseImageSpacing", true );
  this->m_Parameters.SetSelectedChoice( "OutputResolution", "float" );
}

// -------------------------------------------------------------------------
cpPluginsImageFilters::SignedMaurerDistanceMapImageFilter::
~SignedMaurerDistanceMapImageFilter( )
{
}

// -------------------------------------------------------------------------
std::string cpPluginsImageFilters::SignedMaurerDistanceMapImageFilter::
_GenerateData( )
{
  auto image = this->GetInputData( "Input" )->GetITK< itk::DataObject >( );
  std::string   cpPlugin_Image_Demangle_Pixel_AllScalars( r, _GD0, image, 2 );
  if( r != "" ) cpPlugin_Image_Demangle_Pixel_AllScalars( r, _GD0, image, 3 );
  return( r );
}

// -------------------------------------------------------------------------
template< class _TImage >
std::string cpPluginsImageFilters::SignedMaurerDistanceMapImageFilter::
_GD0( _TImage* image )
{
  if( image != NULL )
  {
    std::string out_res =
      this->m_Parameters.GetSelectedChoice( "OutputResolution" );
    if( out_res == "float" )
      return(
        this->_GD1< _TImage, itk::Image< float, _TImage::ImageDimension > >( image )
        );
  else if( out_res == "double" )
    return(
      this->_GD1< _TImage, itk::Image< double, _TImage::ImageDimension > >( image )
      );
  else
    return( "ImageFilters::SignedMaurerDistanceMapImageFilter: Output resolution not supported." );
  }
  else
    return(
      "ImageFilters::SignedMaurerDistanceMapImageFilter: No valid input image."
      );
}

// -------------------------------------------------------------------------
template< class _TImage, class _TDMap >
std::string cpPluginsImageFilters::SignedMaurerDistanceMapImageFilter::
_GD1( _TImage* image )
{
  typedef itk::SignedMaurerDistanceMapImageFilter< _TImage, _TDMap > _F;

  // Get parameters
  double back_value = this->m_Parameters.GetReal( "BackgroundValue" );
  bool pos_inside = this->m_Parameters.GetBool( "InsideIsPositive" );
  bool sqr_dist = this->m_Parameters.GetBool( "SquaredDistance" );
  bool use_spac = this->m_Parameters.GetBool( "UseImageSpacing" );

  // Configure filter
  _F* filter = this->_CreateITK< _F >( );
  filter->SetInput( image );
  filter->SetBackgroundValue( ( typename _TImage::PixelType )( back_value ) );
  filter->SetInsideIsPositive( pos_inside );
  filter->SetSquaredDistance( sqr_dist );
  filter->SetUseImageSpacing( use_spac );
  filter->Update( );

  // Connect output
  this->GetOutputData( "Output" )->SetITK( filter->GetOutput( ) );
  return( "" );
}

// eof - $RCSfile$