itkLocalHistogramDerivativeForceFilter.h

Go to the documentation of this file.

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

  NifTK: A software platform for medical image computing.

  Copyright (c) University College London (UCL). All rights reserved.

  This software is distributed WITHOUT ANY WARRANTY; without even
  the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  PURPOSE.

  See LICENSE.txt in the top level directory for details.

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

#ifndef itkLocalHistogramDerivativeForceFilter_h
#define itkLocalHistogramDerivativeForceFilter_h

#include "itkRegistrationForceFilter.h"

namespace itk {
template< class TFixedImage, class TMovingImage, class TScalar >
class ITK_EXPORT LocalHistogramDerivativeForceFilter :
  public RegistrationForceFilter<TFixedImage, TMovingImage, TScalar>
{
public:

  typedef LocalHistogramDerivativeForceFilter                Self;
  typedef RegistrationForceFilter<TFixedImage, TMovingImage, TScalar> Superclass;
  typedef SmartPointer<Self>                                 Pointer;
  typedef SmartPointer<const Self>                           ConstPointer;
  
  itkTypeMacro(LocalHistogramDerivativeForceFilter, RegistrationForceFilter);

  itkStaticConstMacro(Dimension, unsigned int, TFixedImage::ImageDimension);

  typedef typename Superclass::OutputPixelType                OutputPixelType;
  typedef typename Superclass::OutputImageType                OutputImageType;
  typedef typename OutputImageType::SpacingType               OutputImageSpacingType;
  typedef typename Superclass::InputImageType                 InputImageType;
  typedef typename Superclass::InputImageRegionType           RegionType;
  typedef typename Superclass::MetricType                     MetricType;
  typedef typename Superclass::MetricPointer                  MetricPointer;
  typedef typename Superclass::HistogramType                  HistogramType;
  typedef typename Superclass::HistogramPointer               HistogramPointer;
  typedef typename Superclass::HistogramSizeType              HistogramSizeType;
  typedef typename Superclass::HistogramMeasurementVectorType HistogramMeasurementVectorType;
  typedef typename Superclass::HistogramFrequencyType         HistogramFrequencyType;
  typedef typename Superclass::HistogramIteratorType          HistogramIteratorType;
  typedef typename HistogramType::ConstPointer                HistogramConstPointer;
  typedef typename Superclass::MeasureType                    MeasureType;
  
protected:
  
  LocalHistogramDerivativeForceFilter();
  ~LocalHistogramDerivativeForceFilter() {};
  void PrintSelf(std::ostream& os, Indent indent) const;

  virtual void ThreadedGenerateData( const RegionType &outputRegionForThread, ThreadIdType threadId);

  virtual MeasureType ComputeForcePerVoxel(double totalFrequency, 
                                           double jointEntropy, 
                                           double fixedImageEntropy, 
                                           double movingImageEntropy, 
                                           double transformedMovingImageMinusHistogramIndexJointFrequency,
                                           double transformedMovingImagePlusHistogramIndexJointFrequency,
                                           double transformedMovingImageMinusHistogramIndexFrequency,
                                           double transformedMovingImagePlusHistogramIndexFrequency) const = 0;

private:

  LocalHistogramDerivativeForceFilter(const Self&); 
  void operator=(const Self&); 

};

}

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

#endif