#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;
typedef itk::ImageFileWriter< TFilter::TAuxImage >   TAuxWriter;

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

  } // fi
  std::string in_fname = argv[ 1 ];
  std::string out_fname = argv[ 2 ];
  std::string out_aux_fname = 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->SetSeed( seed );

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

  TAuxWriter::Pointer aux_writer = TAuxWriter::New( );
  aux_writer->SetInput( filter->GetAuxiliaryImage( ) );
  aux_writer->SetFileName( out_aux_fname );
  aux_writer->Update( );
  

  return( 0 );
}

// eof - $RCSfile$