/*=========================================================================

Program:   Insight Segmentation & Registration Toolkit
Module:    $RCSfile: itkNeighborhoodSamplerTest.cxx,v $
Language:  C++
Date:      $Date: 2005/07/26 23:36:40 $
Version:   $Revision: 1.7 $

Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

This software is distributed WITHOUT ANY WARRANTY; without even 
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include "itkWin32Header.h"

#include <fstream>

#include "itkArray.h"

#include "itkListSample.h"
#include "itkSubsample.h"
#include "itkNeighborhoodSampler.h"

#define NO_OF_DIMENSIONS 2
typedef double ValueType ;

int itkNeighborhoodSamplerTest(int , char* [])
{
  namespace stat = itk::Statistics ;

  typedef itk::Vector< ValueType, NO_OF_DIMENSIONS > MeasurementVectorType ; 
  typedef stat::ListSample< MeasurementVectorType > SampleType ;

  SampleType::Pointer sample = SampleType::New() ;
  sample->SetMeasurementVectorSize( NO_OF_DIMENSIONS );

  double radius = 1.0 ;

  // add four vectors within the radius
  MeasurementVectorType tempVector ;
  tempVector[0] = 0.0 ;
  tempVector[1] = 0.5 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = 0.0 ;
  tempVector[1] = -0.5 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = 0.5 ;
  tempVector[1] = 0.0 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = -0.5 ;
  tempVector[1] = 0.0 ;
  sample->PushBack(tempVector) ;

  // add four vectors outside of the radius
  tempVector[0] = 0.0 ;
  tempVector[1] = 2.0 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = 0.0 ;
  tempVector[1] = -2.0 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = 2.0 ;
  tempVector[1] = 0.0 ;
  sample->PushBack(tempVector) ;

  tempVector[0] = -2.0 ;
  tempVector[1] = 0.0 ;
  sample->PushBack(tempVector) ;

  typedef stat::NeighborhoodSampler< SampleType > SamplerType ;
  typedef SamplerType::CenterType CenterType;
  CenterType center( NO_OF_DIMENSIONS );
  center.Fill( 0.0 );
  
  SamplerType::Pointer sampler = SamplerType::New() ;

  sampler->SetInputSample(sample) ;
  sampler->SetCenter(&center) ;
  sampler->SetRadius(&radius) ;
  sampler->Update() ;

  SamplerType::OutputType::Pointer output = sampler->GetOutput() ;
  
  if ( output->Size() != 4 )
    {
      std::cout 
        << "Test failed: The output size should be 4. The actual size =" 
        << output->Size() << std::endl;
      return EXIT_FAILURE;
    }

  for ( unsigned int i = 0 ; i < 4 ; i++ )
    {
      if ( output->GetInstanceIdentifier(i) != i )
        {
          std::cout << "Test failed." << std::endl;
          return EXIT_FAILURE;
        }
    }

  std::cout << "Test succeeded." << std::endl;
  return EXIT_SUCCESS ;
}
