itkPointSetToListSampleAdaptor.h

Go to the documentation of this file.

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

Program:   Insight Segmentation & Registration Toolkit
Module:    $RCSfile: itkPointSetToListSampleAdaptor.h,v $
Language:  C++
Date:      $Date: 2009-05-22 12:55:07 $
Version:   $Revision: 1.3 $

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 __itkPointSetToListSampleAdaptor_h
#define __itkPointSetToListSampleAdaptor_h

#include <typeinfo>

#include "itkPointSet.h"
#include "itkListSample.h"
#include "itkSmartPointer.h"

namespace itk { 
namespace Statistics {

template < class TPointSet >
class ITK_EXPORT PointSetToListSampleAdaptor :
    public ListSample< typename TPointSet::PointType >
{
public:
  typedef PointSetToListSampleAdaptor                           Self;
  typedef ListSample< typename TPointSet::PointType >           Superclass;
  typedef SmartPointer< Self >                                  Pointer;
  typedef SmartPointer<const Self>                              ConstPointer;

  itkTypeMacro(PointSetToListSampleAdaptor, ListSample);

  itkNewMacro(Self);

  itkStaticConstMacro(MeasurementVectorSize, unsigned int,
                      TPointSet::PointDimension);

  typedef TPointSet                                        PointSetType;
  typedef typename TPointSet::Pointer                      PointSetPointer;
  typedef typename TPointSet::ConstPointer                 PointSetConstPointer;
  typedef typename TPointSet::PointsContainer              PointsContainer;
  typedef typename TPointSet::PointsContainerPointer       PointsContainerPointer;
  typedef typename TPointSet::PointsContainerConstPointer  PointsContainerConstPointer;
  typedef typename TPointSet::PointsContainerIterator      PointsContainerIteratorType;
  typedef typename TPointSet::PointsContainerConstIterator PointsContainerConstIteratorType;
  typedef typename TPointSet::PointType                    PointType;


  itkSuperclassTraitMacro( MeasurementType )
  itkSuperclassTraitMacro( MeasurementVectorType )
  itkSuperclassTraitMacro( AbsoluteFrequencyType )
  itkSuperclassTraitMacro( TotalAbsoluteFrequencyType )
  itkSuperclassTraitMacro( MeasurementVectorSizeType )
  itkSuperclassTraitMacro( InstanceIdentifier )

  typedef MeasurementVectorType                           ValueType;

  void SetPointSet(const TPointSet* pointSet);

  const TPointSet* GetPointSet();

  InstanceIdentifier Size() const;

  const MeasurementVectorType & GetMeasurementVector( InstanceIdentifier id ) const;

  AbsoluteFrequencyType GetFrequency( InstanceIdentifier id ) const;

  TotalAbsoluteFrequencyType GetTotalFrequency() const;

  class ConstIterator
    {
    friend class PointSetToListSampleAdaptor;
    public:

    ConstIterator( const PointSetToListSampleAdaptor * adaptor )
      {
      *this = adaptor->Begin();
      }

    ConstIterator(const ConstIterator &iter)
      {
      m_Iter = iter.m_Iter;
      m_InstanceIdentifier = iter.m_InstanceIdentifier;
      }

    ConstIterator& operator=( const ConstIterator & iter )
      {
      m_Iter = iter.m_Iter;
      m_InstanceIdentifier = iter.m_InstanceIdentifier;
      return *this;
      }

    AbsoluteFrequencyType GetFrequency() const
      {
      return 1;
      }

    const MeasurementVectorType & GetMeasurementVector() const
      {
      return (const MeasurementVectorType&) m_Iter.Value();
      }

    InstanceIdentifier GetInstanceIdentifier() const
      {
      return m_InstanceIdentifier;
      }

    ConstIterator& operator++()
      {
      ++m_Iter;
      ++m_InstanceIdentifier;
      return *this;
      }

    bool operator!=(const ConstIterator &it)
      {
      return (m_Iter != it.m_Iter);
      }

    bool operator==(const ConstIterator &it)
      {
      return (m_Iter == it.m_Iter);
      }

#if !(defined(_MSC_VER) && (_MSC_VER <= 1200))
  protected:
#endif
    // This method should only be available to the ListSample class
    ConstIterator(
      PointsContainerConstIteratorType iter,
      InstanceIdentifier iid)
      {
      m_Iter = iter;
      m_InstanceIdentifier = iid;
      }

    // This method is purposely not implemented
    ConstIterator();

  private:
    PointsContainerConstIteratorType      m_Iter;
    InstanceIdentifier                    m_InstanceIdentifier;
    };

  class Iterator : public ConstIterator
  {

  friend class PointSetToListSampleAdaptor;

  public:

    Iterator(Self * adaptor):ConstIterator(adaptor)
      {
      }

    Iterator(const Iterator &iter):ConstIterator( iter )
      {
      }

    Iterator& operator =(const Iterator & iter)
      {
      this->ConstIterator::operator=( iter );
      return *this;
      }

#if !(defined(_MSC_VER) && (_MSC_VER <= 1200))
  protected:
#endif
    // To ensure const-correctness these method must not be in the public API.
    // The are purposly not implemented, since they should never be called.
    Iterator();
    Iterator(const Self * adaptor);
    Iterator(  PointsContainerConstIteratorType iter, InstanceIdentifier iid);
    Iterator(const ConstIterator & it);
    ConstIterator& operator=(const ConstIterator& it);
    Iterator(
      PointsContainerIteratorType iter,
      InstanceIdentifier iid):ConstIterator( iter, iid )
      {
      }

  private:
  };


  Iterator Begin()
    {
    PointsContainerPointer nonConstPointsDataContainer = 
                      const_cast< PointsContainer * > (m_PointsContainer.GetPointer());
    Iterator iter(nonConstPointsDataContainer->Begin(), 0);
    return iter;
    }

  Iterator End()
    {
    PointsContainerPointer nonConstPointsDataContainer = 
                      const_cast< PointsContainer * > (m_PointsContainer.GetPointer());

    Iterator iter(nonConstPointsDataContainer->End(), m_PointsContainer->Size());
    return iter;
    }

  ConstIterator Begin() const
    {
    ConstIterator iter(m_PointsContainer->Begin(), 0);
    return iter;
    }

  ConstIterator End() const
    {
    ConstIterator iter(m_PointsContainer->End(), m_PointsContainer->Size());
    return iter;
    }


protected:
  PointSetToListSampleAdaptor(); 
  
  virtual ~PointSetToListSampleAdaptor() {}
  void PrintSelf(std::ostream& os, Indent indent) const;  

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

  PointSetConstPointer                     m_PointSet;

  PointsContainerConstPointer              m_PointsContainer;

  mutable PointType                        m_TempPoint;

}; // end of class PointSetToListSampleAdaptor

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


#ifndef ITK_MANUAL_INSTANTIATION
#include "itkPointSetToListSampleAdaptor.txx"
#endif

#endif