itkForwardAndBackProjectionDifferenceFilter.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 itkForwardAndBackProjectionDifferenceFilter_h
#define itkForwardAndBackProjectionDifferenceFilter_h

#include <itkImageToImageFilter.h>
#include <itkConceptChecking.h>
#include "itkForwardImageProjector3Dto2D.h"
#include "itkBackwardImageProjector2Dto3D.h"
#include <itkPerspectiveProjectionTransform.h>
#include <itkEulerAffineTransform.h>
#include "itkProjectionGeometry.h"
#include "itkSubtract2DImageFromVolumeSliceFilter.h"

namespace itk
{
  
template <class IntensityType = float>
class ITK_EXPORT ForwardAndBackProjectionDifferenceFilter : 
  public ImageToImageFilter<Image<IntensityType, 3>,
          Image<IntensityType, 3> > 
{
public:

  typedef ForwardAndBackProjectionDifferenceFilter      Self;
  typedef ImageToImageFilter<Image< IntensityType, 3>, 
           Image< IntensityType, 3> > Superclass;
  typedef SmartPointer<Self>                            Pointer;
  typedef SmartPointer<const Self>                      ConstPointer;

  itkNewMacro(Self);

  itkTypeMacro(ForwardAndBackProjectionDifferenceFilter, ImageToImageFilter);

  typedef Image<IntensityType, 3>                   InputVolumeType;
  typedef typename    InputVolumeType::Pointer      InputVolumePointer;
  typedef typename    InputVolumeType::ConstPointer InputVolumeConstPointer;
  typedef typename    InputVolumeType::RegionType   InputVolumeRegionType;
  typedef typename    InputVolumeType::PixelType    InputVolumePixelType;
  typedef typename    InputVolumeType::SizeType     InputVolumeSizeType;
  typedef typename    InputVolumeType::SpacingType  InputVolumeSpacingType;
  typedef typename    InputVolumeType::PointType    InputVolumePointType;

  typedef Image<IntensityType, 3>                             InputProjectionVolumeType;
  typedef typename    InputProjectionVolumeType::Pointer      InputProjectionVolumePointer;
  typedef typename    InputProjectionVolumeType::RegionType   InputProjectionVolumeRegionType;
  typedef typename    InputProjectionVolumeType::PixelType    InputProjectionVolumePixelType;
  typedef typename    InputProjectionVolumeType::SizeType     InputProjectionVolumeSizeType;
  typedef typename    InputProjectionVolumeType::SpacingType  InputProjectionVolumeSpacingType;
  typedef typename    InputProjectionVolumeType::PointType    InputProjectionVolumePointType;

  typedef Image<IntensityType, 3>                                      OutputBackProjectedDifferencesType;
  typedef typename     OutputBackProjectedDifferencesType::Pointer     OutputBackProjectedDifferencesPointer;
  typedef typename     OutputBackProjectedDifferencesType::RegionType  OutputBackProjectedDifferencesRegionType;
  typedef typename     OutputBackProjectedDifferencesType::PixelType   OutputBackProjectedDifferencesPixelType;
  typedef typename     OutputBackProjectedDifferencesType::SizeType    OutputBackProjectedDifferencesSizeType;
  typedef typename     OutputBackProjectedDifferencesType::SpacingType OutputBackProjectedDifferencesSpacingType;
  typedef typename     OutputBackProjectedDifferencesType::PointType   OutputBackProjectedDifferencesPointType;

  typedef itk::EulerAffineTransform<double, 3, 3> EulerAffineTransformType;
  typedef typename EulerAffineTransformType::Pointer EulerAffineTransformPointer;

  typedef itk::PerspectiveProjectionTransform<double> PerspectiveProjectionTransformType;
  typedef typename PerspectiveProjectionTransformType::Pointer PerspectiveProjectionTransformPointer;

  typedef itk::ForwardImageProjector3Dto2D<IntensityType> ForwardImageProjector3Dto2DType;
  typedef typename ForwardImageProjector3Dto2DType::Pointer ForwardImageProjector3Dto2DPointer;

  typedef typename ForwardImageProjector3Dto2DType::OutputImageType ForwardProjectorOutputImageType;
  typedef typename ForwardProjectorOutputImageType::Pointer         ForwardProjectorOutputImagePointer;
  typedef typename ForwardProjectorOutputImageType::RegionType      ForwardProjectorOutputImageRegionType;
  typedef typename ForwardProjectorOutputImageType::SizeType        ForwardProjectorOutputImageSizeType;
  typedef typename ForwardProjectorOutputImageType::SpacingType     ForwardProjectorOutputImageSpacingType;
  typedef typename ForwardProjectorOutputImageType::PointType       ForwardProjectorOutputImagePointType;
  typedef typename ForwardProjectorOutputImageType::PixelType       ForwardProjectorOutputImagePixelType;
  typedef typename ForwardProjectorOutputImageType::IndexType       ForwardProjectorOutputImageIndexType;

  typedef itk::Subtract2DImageFromVolumeSliceFilter<IntensityType> Subtract2DImageFromVolumeSliceFilterType;
  typedef typename Subtract2DImageFromVolumeSliceFilterType::Pointer Subtract2DImageFromVolumeSliceFilterPointer;

  typedef itk::BackwardImageProjector2Dto3D<IntensityType> BackwardImageProjector2Dto3DType;
  typedef typename BackwardImageProjector2Dto3DType::Pointer BackwardImageProjector2Dto3DPointer;

  typedef itk::ProjectionGeometry<IntensityType> ProjectionGeometryType;
  typedef typename ProjectionGeometryType::Pointer ProjectionGeometryPointer;

  void SetInputVolume( InputVolumeType *im3D);

  void SetInputProjectionVolume( InputProjectionVolumeType *im2D);

  itkSetObjectMacro( ProjectionGeometry, ProjectionGeometryType );
  itkGetObjectMacro( ProjectionGeometry, ProjectionGeometryType );

  InputVolumePointer GetPointerToInputVolume(void);

  itkStaticConstMacro(InputVolumeDimension, unsigned int,
                      InputVolumeType::ImageDimension);
  itkStaticConstMacro(OutputBackProjectedDifferencesDimension, unsigned int,
                      OutputBackProjectedDifferencesType::ImageDimension);

  virtual void GenerateInputRequestedRegion();
  virtual void EnlargeOutputRequestedRegion(DataObject *output); 

  void Initialise(void);
  
  void ClearVolumePriorToNextBackProjection(void) {m_BackProjector->ClearVolumePriorToNextBackProjection();}

protected:

  ForwardAndBackProjectionDifferenceFilter();
  virtual ~ForwardAndBackProjectionDifferenceFilter() {};

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

  void GenerateData();

  bool m_FlagPipelineInitialised;

  unsigned int m_NumberOfProjections;

  ForwardImageProjector3Dto2DPointer m_ForwardProjector;

  Subtract2DImageFromVolumeSliceFilterPointer m_SubtractProjectionFromEstimate;

  BackwardImageProjector2Dto3DPointer m_BackProjector;

  ProjectionGeometryPointer m_ProjectionGeometry;

  virtual void VerifyInputInformation() {};

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

} // end namespace itk

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

#endif