itkHistogramImageToImageMetric.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkHistogramImageToImageMetric.h,v $
  Language:  C++
  Date:      $Date: 2003/09/10 14:28:31 $
  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.

=========================================================================*/
#ifndef __itkHistogramImageToImageMetric_h
#define __itkHistogramImageToImageMetric_h

#include "itkHistogram.h"
#include "itkImageToImageMetric.h"

namespace itk
{
template <class TFixedImage, class TMovingImage>
class ITK_EXPORT HistogramImageToImageMetric : 
public ImageToImageMetric<TFixedImage, TMovingImage>
{
public:
  typedef HistogramImageToImageMetric Self;
  typedef ImageToImageMetric<TFixedImage, TMovingImage> Superclass;
  typedef SmartPointer<Self>                            Pointer;
  typedef SmartPointer<const Self>                      ConstPointer;

  itkNewMacro(Self);
 
  itkTypeMacro(HistogramImageToImageMetric, ImageToImageMetric);
 
  typedef typename Superclass::RealType                   RealType;
  typedef typename Superclass::TransformType              TransformType;
  typedef typename Superclass::TransformPointer           TransformPointer;
  typedef typename Superclass::TransformParametersType
    TransformParametersType;
  typedef typename Superclass::TransformJacobianType
    TransformJacobianType;
  typedef typename Superclass::GradientPixelType          GradientPixelType;
  typedef typename Superclass::InputPointType             InputPointType;
  typedef typename Superclass::OutputPointType            OutputPointType;
  typedef typename Superclass::MeasureType                MeasureType;
  typedef typename Superclass::DerivativeType             DerivativeType;
  typedef typename Superclass::FixedImageType             FixedImageType;
  typedef typename Superclass::FixedImageType::PixelType  FixedImagePixelType;
  typedef typename Superclass::MovingImageType            MovingImageType;
  typedef typename Superclass::MovingImageType::PixelType MovingImagePixelType;
  typedef typename Superclass::FixedImageConstPointer
    FixedImageConstPointerType;
  typedef typename Superclass::MovingImageConstPointer
    MovingImageConstPointerType;

  typedef Statistics::Histogram<double, 2> HistogramType;
  typedef typename HistogramType::MeasurementVectorType  MeasurementVectorType;
  typedef typename HistogramType::SizeType               HistogramSizeType;

  void Initialize() throw (ExceptionObject);

  inline void SetHistogramSize(const HistogramSizeType & size)
    {
      m_HistogramSize = size;
    }

  inline const HistogramSizeType & GetHistogramSize() const
    {
      return m_HistogramSize;
    }

  inline void SetPaddingValue(FixedImagePixelType const padding)
    {
      m_PaddingValue = padding;
    }

  inline FixedImagePixelType GetPaddingValue() const
    {
      return m_PaddingValue;
    }

  inline void UsePaddingValue(bool use = true)
    {
      m_UsePaddingValue = use;
    }

  inline void SetDerivativeStepLength(double length)
    {
      m_DerivativeStepLength = length;
    }

  inline double GetDerivativeStepLength()
    {
      return m_DerivativeStepLength;
    }

  typedef Array<double> ScalesType;

  inline void SetDerivativeStepLengthScales(ScalesType const& scales)
    {
      m_DerivativeStepLengthScales = scales;
    }

  itkGetConstReferenceMacro(DerivativeStepLengthScales, ScalesType);

  MeasureType GetValue(const TransformParametersType& parameters) const;

  void GetDerivative(const TransformParametersType & parameters,
                     DerivativeType & derivative) const;

  void GetValueAndDerivative(const TransformParametersType & parameters,
                             MeasureType& Value,
                             DerivativeType& Derivative) const;

protected:
  HistogramImageToImageMetric();
  virtual ~HistogramImageToImageMetric() {};

  HistogramSizeType m_HistogramSize;
  mutable MeasurementVectorType m_LowerBound;
  mutable MeasurementVectorType m_UpperBound;

  void ComputeHistogram(const TransformParametersType & parameters,
                        HistogramType& histogram) const;
  void ComputeHistogram(const TransformParametersType & parameters,
                        unsigned int parameter,
                        double step,
                        HistogramType& histogram) const;
  void CopyHistogram(HistogramType& target, HistogramType& source) const;

  virtual MeasureType EvaluateMeasure(HistogramType& histogram) const = 0;
  
  void PrintSelf(std::ostream& os, Indent indent) const;

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

  FixedImagePixelType m_PaddingValue;

  bool m_UsePaddingValue;

  double m_DerivativeStepLength;

  ScalesType m_DerivativeStepLengthScales;
};

} // end namespace itk

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

#endif // __itkHistogramImageToImageMetric_h

---