#ifndef __fpa__Image__MoriRegionGrow__hxx__
#define __fpa__Image__MoriRegionGrow__hxx__

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
MoriRegionGrow( )
  : Superclass( ),
    m_Step( TPixel( 1 ) ),
    m_Sensitivity( double( 1 ) )
{
  this->m_Upper = std::numeric_limits< TPixel >::max( );
  if( std::numeric_limits< TPixel >::is_integer )
    this->m_Lower = std::numeric_limits< TPixel >::min( );
  else
    this->m_Lower = -this->m_Upper;
  typename TBinThresholdFunction::Pointer functor =
    TBinThresholdFunction::New( );
  this->SetGrowFunction( functor );
}

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

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
bool fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_ContinueGenerateData( )
{
  TBinThresholdFunction* functor =
    dynamic_cast< TBinThresholdFunction* >( this->GetGrowFunction( ) );
  TPixel u = functor->GetUpper( );

  // Analyze flooding
  bool stop = false;
  if( this->m_N > double( 10 ) )
  {
    double Q = double( this->m_PrevCount );
    Q += ( this->m_STD * this->m_Sensitivity );
    if( double( this->m_ActualCount ) > Q )
      stop = true;

  } // fi
  double v;
  if( this->m_PrevCount > 0 )
    v = double( this->m_ActualCount ) - double( this->m_PrevCount );
  else
    v = double( 0 );
  this->m_PrevCount = this->m_ActualCount;
  this->m_N += double( 1 );
  this->m_S1 += v;
  this->m_S2 += v * v;
  if( this->m_N > double( 1 ) )
    this->m_STD =
      ( this->m_S2 - ( ( this->m_S1 * this->m_S1 ) / this->m_N ) ) /
      ( this->m_N - double( 1 ) );
  else
    this->m_STD = double( 0 );
  this->m_STD = std::sqrt( this->m_STD );

  std::cout << u << " " << this->m_ActualCount << " " << this->m_STD << std::endl;

  if( u < this->m_Upper )
  {
    u += this->m_Step;
    if( u > this->m_Upper )
      u = this->m_Upper;
    functor->SetUpper( u );
    this->m_Queue = this->m_NextQueue;
    while( this->m_NextQueue.size( ) > 0 )
      this->m_NextQueue.pop( );
    return( !stop );
  }
  else
    return( false );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_BeforeGenerateData( )
{
  this->Superclass::_BeforeGenerateData( );

  while( this->m_NextQueue.size( ) > 0 )
    this->m_NextQueue.pop( );
  this->m_ActualCount = 0;
  this->m_PrevCount = 0;
  this->m_N   = double( 0 );
  this->m_S1  = double( 0 );
  this->m_S2  = double( 0 );
  this->m_STD = double( 0 );
  TBinThresholdFunction* functor =
    dynamic_cast< TBinThresholdFunction* >( this->GetGrowFunction( ) );
  functor->SetLower( this->m_Lower );
  functor->SetUpper( this->m_Lower + this->m_Step );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_AfterGenerateData( )
{
  this->Superclass::_AfterGenerateData( );
  while( this->m_NextQueue.size( ) > 0 )
    this->m_NextQueue.pop( );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_BeforeLoop( )
{
  this->Superclass::_BeforeLoop( );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_AfterLoop( )
{
  this->Superclass::_AfterLoop( );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
bool fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_UpdateValue( _TQueueNode& v, const _TQueueNode& p )
{
  bool ret = this->Superclass::_UpdateValue( v, p );
  if( !ret )
  {
    v.Result = this->m_InsideValue;
    this->m_NextQueue.push( v );

  } // fi
  return( ret );
}

// -------------------------------------------------------------------------
template< class _TInputImage, class _TOutputImage >
void fpa::Image::MoriRegionGrow< _TInputImage, _TOutputImage >::
_UpdateResult( const _TQueueNode& n )
{
  this->Superclass::_UpdateResult( n );
  this->m_ActualCount += 1;
}

#endif // __fpa__Image__MoriRegionGrow__hxx__

// eof - $RCSfile$