#include "itkImageFileReader.h"
#include "itkImageFileWriter.h"

#include "itkBinaryWellComposed3DImageFilter.h"

int main( int argc, char * argv[])
{
  const unsigned int      ImageDimension = 3;
  typedef int             PixelType;
  typedef itk::Image<PixelType,  ImageDimension> ImageType;

  typedef itk::ImageFileReader<ImageType> ReaderType;
  ReaderType::Pointer reader = ReaderType::New();
  reader->SetFileName( "brain3d.hdr" );
  reader->Update();

  typedef itk::BinaryWellComposed3DImageFilter<ImageType> FilterType;
  FilterType::Pointer filter = FilterType::New();
  filter->SetInput( reader->GetOutput() );
  filter->SetBackgroundValue( static_cast<PixelType>( 0 ) );
  filter->SetForegroundValue( static_cast<PixelType>( 1 ) );
  filter->Update();

  typedef itk::ImageFileWriter<ImageType> WriterType;
  WriterType::Pointer writer = WriterType::New();
  writer->SetInput( filter->GetOutput() );
  writer->SetFileName( "brain3d_out.hdr" );
  writer->Update();
  
  return 0;
};