#include <cpPluginsImageFilters/MaskedImageToHistogramFilter.h>
#include <cpPlugins/Image.h>
#include <cpPlugins_ITKInstances/ImageFilters.h>

#include <itkMaskedImageToHistogramFilter.h>
#include <itkMaskedImageToHistogramFilter.hxx>

// -------------------------------------------------------------------------
cpPluginsImageFilters::MaskedImageToHistogramFilter::
MaskedImageToHistogramFilter( )
  : Superclass( )
{
  this->_AddInput( "Input" );
  this->_AddInput( "Mask" );
  this->_AddOutput< cpPlugins::DataObject >( "Output" );
}

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

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

// -------------------------------------------------------------------------
template< class _TImage >
std::string cpPluginsImageFilters::MaskedImageToHistogramFilter::
_GD0( _TImage* image )
{
  if( image != NULL )
    return(
      this->_GD1< _TImage, itk::Image< unsigned char, _TImage::ImageDimension > >( image )
      );
  else
    return(
      "ImageFilters::MaskedImageToHistogramFilter: No valid input image."
      );
}

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

  // Get parameters

  // Configure filter
  _F* filter = this->_CreateITK< _F >( );
  /*
    filter->SetInput( image );
    filter->SetLowerThreshold( lower_val );
    filter->SetUpperThreshold( upper_val );
    filter->SetInsideValue( in_val );
    filter->SetOutsideValue( out_val );
    filter->Update( );

    // Connect output
    this->GetOutputData( "Output" )->SetITK( filter->GetOutput( ) );
  */
  return( "MaskedImageToHistogramFilter: not ready yet." );
}

// eof - $RCSfile$