// =========================================================================
// @author Leonardo Florez Valencia
// @email florez-l@javeriana.edu.co
// =========================================================================
#ifndef __fpa__Filters__Image__Mori__hxx__
#define __fpa__Filters__Image__Mori__hxx__

#include <itkBinaryThresholdImageFilter.h>

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
fpa::Filters::Image::Mori< _TInputImage, _TOutputImage >::
Mori( )
  : Superclass( )
{
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
fpa::Filters::Image::Mori< _TInputImage, _TOutputImage >::
~Mori( )
{
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Filters::Image::Mori< _TInputImage, _TOutputImage >::
_AfterGenerateData( )
{
  this->Superclass::_AfterGenerateData( );

  // Extract correct mark
  typedef typename Superclass::TMarksImage _TMarks;
  typedef itk::BinaryThresholdImageFilter< _TMarks, _TOutputImage > _TFilter;
  typename _TFilter::Pointer filter = _TFilter::New( );
  filter->SetInput( this->GetMarks( ) );
  filter->SetLowerThreshold( 1 );
  filter->SetUpperThreshold( this->m_PeakDetector.GetNumberOfSamples( ) - 2 );
  filter->SetInsideValue( this->GetInsideValue( ) );
  filter->SetOutsideValue( this->GetOutsideValue( ) );
  filter->Update( );
  this->GraftOutput( filter->GetOutput( ) );
}

#endif // __fpa__Filters__Image__Mori__hxx__
// eof - $RCSfile$