#include <itkImage.h>
#include <itkImageFileReader.h>
#include <itkImageFileWriter.h>
#include <fpa/Image/MoriRegionGrow.h>

// -------------------------------------------------------------------------
typedef itk::Image< unsigned char, 2 >               TImage;
typedef itk::ImageFileReader< TImage >               TReader;
typedef itk::ImageFileWriter< TImage >               TWriter;
typedef fpa::Image::MoriRegionGrow< TImage, TImage > TFilter;

// -------------------------------------------------------------------------
int main( int argc, char* argv[] )
{
  if( argc < 4 )
  {
    std::cerr
      << "Usage: " << argv[ 0 ]
      << " input_filename output_filename sensitivity" << std::endl;
    return( 1 );

  } // fi
  std::string in_fname = argv[ 1 ];
  std::string out_fname = argv[ 2 ];
  double sensitivity = std::atof( argv[ 3 ] );
  int seed_x = 111;
  int seed_y = 91;

  TReader::Pointer reader = TReader::New( );
  reader->SetFileName( in_fname );

  TImage::IndexType seed;
  seed[ 0 ] = seed_x;
  seed[ 1 ] = seed_y;

  TFilter::Pointer filter = TFilter::New( );
  filter->SetInput( reader->GetOutput( ) );
  filter->SetLower( 0 );
  filter->SetUpper( 255 );
  filter->SetStep( 1 );
  filter->SetInsideValue( 255 );
  filter->SetOutsideValue( 0 );
  filter->SetSensitivity( sensitivity );
  filter->AddSeed( seed, filter->GetInsideValue( ) );

  TWriter::Pointer writer = TWriter::New( );
  writer->SetInput( filter->GetOutput( ) );
  writer->SetFileName( out_fname );
  writer->Update( );

  return( 0 );
}

// eof - $RCSfile$