#include <cstdlib>
#include <iostream>
#include <string>

#include <itkImage.h>
#include <itkImageFileReader.h>
#include <itkRGBPixel.h>

#include <cpPlugins/Extensions/Algorithms/ParallelImageMean.h>

// -------------------------------------------------------------------------
const unsigned int Dim = 2;
typedef itk::RGBPixel< unsigned char > TPixel;
typedef itk::Image< TPixel, Dim >      TImage;

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

  } // fi
  std::string image_fn = argv[ 1 ];

  // Read image
  itk::ImageFileReader< TImage >::Pointer image_reader =
    itk::ImageFileReader< TImage >::New( );
  image_reader->SetFileName( image_fn );
  try
  {
    image_reader->Update( );
  }
  catch( itk::ExceptionObject& err )
  {
    std::cerr << "Error caught: " << err << std::endl;
    return( 1 );

  } // yrt
  TImage::Pointer image = image_reader->GetOutput( );

  // Compute mean
  cpPlugins::Extensions::Algorithms::ParallelImageMean< TImage >::Pointer mean =
    cpPlugins::Extensions::Algorithms::ParallelImageMean< TImage >::New( );
  mean->Execute( image, image->GetRequestedRegion( ) );
  std::cout << "Mean: " << mean->GetMean( ) << std::endl;

  return( 0 );
}

// eof - $RCSfile$