-#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$