itkListSampleBase.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkListSampleBase.h,v $
  Language:  C++
  Date:      $Date: 2010-02-02 13:45:16 $
  Version:   $Revision: 1.12 $

  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.

=========================================================================*/
#ifndef __itkListSampleBase_h
#define __itkListSampleBase_h

#include "itkSample.h"

#include <vector>

namespace itk { 
namespace Statistics {

template< class TMeasurementVector >
class ITK_EXPORT ListSampleBase : public Sample< TMeasurementVector >
{
public:
  typedef ListSampleBase               Self;
  typedef Sample< TMeasurementVector > Superclass;
  typedef SmartPointer< Self >         Pointer;
  typedef SmartPointer<const Self>     ConstPointer;

  itkTypeMacro(ListSampleBase, Sample);

  typedef typename Superclass::MeasurementVectorType     MeasurementVectorType;
  typedef typename Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;
  typedef typename Superclass::MeasurementType           MeasurementType;
  typedef typename Superclass::FrequencyType             FrequencyType;
  typedef typename Superclass::InstanceIdentifier        InstanceIdentifier;

  typedef std::vector< InstanceIdentifier > SearchResultVectorType;

  inline void Search(MeasurementVectorType center, double radius, 
                     SearchResultVectorType& result) const
    {
    if (radius == 0.0)
      {
      itkGenericExceptionMacro("Search radius should be greater than zero.");
      }

    unsigned int j;
    double squaredRadius;
    double distance;
    double coordinateDistance;
    
    MeasurementVectorType tempVector;
    
    squaredRadius = radius * radius;
    
    result.clear();
    for ( InstanceIdentifier identifier = 0; identifier < this->Size(); ++identifier )
      {
      distance = 0.0;
      tempVector = this->GetMeasurementVector( identifier );
      for (j = 0; j < this->GetMeasurementVectorSize() && distance < squaredRadius; j++)
        {
        coordinateDistance = (double)tempVector[j] - center[j];
        if (vnl_math_abs(coordinateDistance) > radius )
          {
          distance = squaredRadius;
          }
        }
      
      for (j = 0; j < this->GetMeasurementVectorSize() && distance < squaredRadius; j++)
        {
        coordinateDistance = (double)tempVector[j] - center[j];
        distance += coordinateDistance * coordinateDistance;
        }
      
      if (distance < squaredRadius)
        {
        result.push_back( identifier );
        }
      }
    }
  
protected:
  ListSampleBase() {}
  virtual ~ListSampleBase() {}
  void PrintSelf(std::ostream& os, Indent indent) const
    {
    Superclass::PrintSelf(os,indent);
    }

  
private:
  ListSampleBase(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

};

} // end of namespace Statistics 
} // end of namespace itk 

#endif