itkParzenWindowNMIDerivativeForceGenerator.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 itkParzenWindowNMIDerivativeForceGenerator_h
#define itkParzenWindowNMIDerivativeForceGenerator_h

#include "itkRegistrationForceFilter.h"
#include "itkLinearlyInterpolatedDerivativeFilter.h"

namespace itk {
template< class TFixedImage, class TMovingImage, class TScalarType, class TDeformationScalar>
class ITK_EXPORT ParzenWindowNMIDerivativeForceGenerator :
  public RegistrationForceFilter<TFixedImage, TMovingImage, TDeformationScalar>
{
public:

  typedef ParzenWindowNMIDerivativeForceGenerator                         Self;
  typedef RegistrationForceFilter<TFixedImage, TMovingImage, TDeformationScalar> Superclass;
  typedef SmartPointer<Self>                                              Pointer;
  typedef SmartPointer<const Self>                                        ConstPointer;

  itkNewMacro(Self);
  
  itkTypeMacro(ParzenWindowNMIDerivativeForceGenerator, RegistrationForceFilter);

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

  typedef typename Superclass::OutputDataType                 OutputDataType;
  typedef typename Superclass::OutputPixelType                OutputPixelType;
  typedef typename Superclass::OutputImageType                OutputImageType;
  typedef typename OutputImageType::SpacingType               OutputImageSpacingType;
  typedef typename Superclass::InputImageType                 InputImageType;
  typedef typename InputImageType::PixelType                  InputPixelType;
  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;
  typedef LinearlyInterpolatedDerivativeFilter<TFixedImage, TMovingImage, TScalarType, TDeformationScalar>
    ScalarImageGradientFilterType;
  typedef typename ScalarImageGradientFilterType::Pointer 
    ScalarImageGradientFilterPointer;

  itkSetObjectMacro( ScalarImageGradientFilter, ScalarImageGradientFilterType );

  itkGetConstObjectMacro( ScalarImageGradientFilter, ScalarImageGradientFilterType );

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

  ScalarImageGradientFilterPointer m_ScalarImageGradientFilter;

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

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

};

}

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

#endif