itkSimultaneousReconstructionRegistrationMetric.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 itkSimultaneousReconstructionRegistrationMetric_h
#define itkSimultaneousReconstructionRegistrationMetric_h

#include <itkImageToImageFilter.h>
#include <itkConceptChecking.h>
#include <itkSingleValuedCostFunction.h>

#include <itkForwardProjectionWithAffineTransformDifferenceFilter.h>


namespace itk
{
  
template <class IntensityType = float>
class ITK_EXPORT SimultaneousReconstructionRegistrationMetric : public SingleValuedCostFunction
{
public:

  typedef SimultaneousReconstructionRegistrationMetric    Self;
  typedef SingleValuedCostFunction                        Superclass;
  typedef SmartPointer<Self>                              Pointer;
  typedef SmartPointer<const Self>                        ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(SimultaneousReconstructionRegistrationMetric, SingleValuedCostFunction);

  typedef itk::ForwardProjectionWithAffineTransformDifferenceFilter<IntensityType>    ForwardProjectionWithAffineTransformDifferenceFilterType;
  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::Pointer  ForwardProjectionWithAffineTransformDifferenceFilterPointerType;

  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::InputVolumeType    InputVolumeType;
  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::InputVolumePointer InputVolumePointer;

  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::InputProjectionVolumeType    InputProjectionVolumeType;
  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::InputProjectionVolumePointer InputProjectionVolumePointer;

  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::ProjectionGeometryType     ProjectionGeometryType;
  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::ProjectionGeometryPointer  ProjectionGeometryPointer;

  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::OutputBackProjectedDifferencesType     OutputBackProjectedDifferencesType;
  typedef typename ForwardProjectionWithAffineTransformDifferenceFilterType::OutputBackProjectedDifferencesPointer  OutputBackProjectedDifferencesPointer;
  

  typedef typename SingleValuedCostFunction::ParametersType ParametersType;
  typedef typename SingleValuedCostFunction::MeasureType    MeasureType;
  typedef typename SingleValuedCostFunction::DerivativeType DerivativeType;

  void SetInputVolume( InputVolumeType *im3D );

  void SetInputProjectionVolumeOne( InputProjectionVolumeType *im2DNumberOne );
  void SetInputProjectionVolumeTwo( InputProjectionVolumeType *im2DNumberTwo );

  void SetProjectionGeometry( ProjectionGeometryType *projGeometry );

  unsigned int GetNumberOfParameters(void) const;

  const InputVolumeType *GetReconstructedVolume(void) {
    return m_FwdAndBackProjDiffFilter->GetPointerToInputVolume();
  }

  void SetParametersAddress( const ParametersType & parameters ) const;
  void SetDerivativesAddress( const DerivativeType & derivatives ) const;

  void SetParameters( const ParametersType & parameters ) const;

  void Initialise(void) {m_FwdAndBackProjDiffFilter->Initialise();}

  virtual void GetValueAndDerivative( const ParametersType & parameters,
                                      MeasureType & value,
                                      DerivativeType & derivative ) const;

  virtual MeasureType GetValue( const ParametersType & parameters ) const;

  virtual void GetDerivative( const ParametersType & parameters,
                              DerivativeType & derivative ) const;


protected:

  SimultaneousReconstructionRegistrationMetric();
  virtual ~SimultaneousReconstructionRegistrationMetric() {};

  void PrintSelf(std::ostream& os, Indent indent) const;
  
  ForwardProjectionWithAffineTransformDifferenceFilterPointerType m_FwdAndBackProjDiffFilter;


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

} // end namespace itk

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

#endif