itkDirectFourierReconstructionImageToImageFilter.h

Go to the documentation of this file.

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

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkDirectFourierReconstructionImageToImageFilter.h,v $
  Language:  C++
  Date:      $Date: 2010-06-14 15:01:05 $
  Version:   $Revision: 1.4 $

  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 __itkDirectFourierReconstructionImageToImageFilter_h
#define __itkDirectFourierReconstructionImageToImageFilter_h

#include "itkImageToImageFilter.h"
#include "itkImage.h"

#include "itkVnlFFTRealToComplexConjugateImageFilter.h"
#include "itkVnlFFTComplexConjugateToRealImageFilter.h"

#include "itkImageRegionIteratorWithIndex.h"
#include "itkImageSliceConstIteratorWithIndex.h"

#include "itkComplexBSplineInterpolateImageFunction.h"

#include <math.h>

namespace itk
{
template< class TInputPixelType = double, class TOutputPixelType = double >
 class ITK_EXPORT DirectFourierReconstructionImageToImageFilter :
  public ImageToImageFilter< Image< TInputPixelType, 3 >, Image< TOutputPixelType, 3 > >
 {
public:
  typedef DirectFourierReconstructionImageToImageFilter Self;
  typedef Image< TInputPixelType, 3 > InputImageType;
  typedef Image< TOutputPixelType, 3 > OutputImageType;
  
  typedef ImageToImageFilter< InputImageType, OutputImageType > Superclass;
  
  typedef SmartPointer< Self > Pointer;
  typedef SmartPointer< const Self > ConstPointer;
  
  itkNewMacro( Self );
  itkTypeMacro( DirectFourierReconstructionImageToImageFilter, ImageToImageFilter );
  
  typedef typename InputImageType::RegionType RegionType;
  typedef typename InputImageType::IndexType IndexType;
  typedef typename InputImageType::SizeType SizeType;
  typedef typename InputImageType::PointType PointType;
  typedef typename InputImageType::SpacingType SpacingType;
  
  typedef typename InputImageType::ConstPointer ConstInputImagePointer;
  typedef typename InputImageType::Pointer InputImagePointer;
  typedef typename OutputImageType::Pointer OutputImagePointer;
  
  itkSetMacro( ZeroPadding, unsigned short int );
  itkGetConstMacro( ZeroPadding, unsigned short int );
  itkSetMacro( OverSampling, unsigned short int );
  itkGetConstMacro( OverSampling, unsigned short int );
  itkSetMacro( Cutoff, double );
  itkGetConstMacro( Cutoff, double );
  itkSetMacro( AlphaRange, double );
  itkGetConstMacro( AlphaRange, double );
  itkSetMacro( AlphaDirection, unsigned short int );
  itkGetConstMacro( AlphaDirection, unsigned short int );
  itkSetMacro( ZDirection, unsigned short int );
  itkGetConstMacro( ZDirection, unsigned short int );
  itkSetMacro( RDirection, unsigned short int );
  itkGetConstMacro( RDirection, unsigned short int );
  itkSetMacro( RadialSplineOrder, unsigned short int );
  itkGetConstMacro( RadialSplineOrder, unsigned short int ); 
  
protected:
  DirectFourierReconstructionImageToImageFilter();
  ~DirectFourierReconstructionImageToImageFilter() {};
  
  void PrintSelf( std::ostream &os, Indent indent ) const;
  
  void GenerateOutputInformation();
  void GenerateInputRequestedRegion();
  void GenerateData();
  
private:
  typedef ImageSliceConstIteratorWithIndex< InputImageType > InputSliceIteratorType;

  typedef VnlFFTRealToComplexConjugateImageFilter< double, 1 > FFTLineFilterType;
  typedef FFTLineFilterType::OutputImageType FFTLineType;
  typedef FFTLineFilterType::InputImageType ProjectionLineType;
  typedef ImageRegionIteratorWithIndex< FFTLineType > FFTLineIteratorType;
  typedef ComplexBSplineInterpolateImageFunction< FFTLineType, double, double > FFTLineInterpolatorType;
  
  typedef VnlFFTComplexConjugateToRealImageFilter< double, 2 > IFFTSliceFilterType;
  typedef IFFTSliceFilterType::InputImageType FFTSliceType;
  typedef IFFTSliceFilterType::OutputImageType OutputSliceType;
  typedef ImageRegionIteratorWithIndex< FFTSliceType > FFTSliceIteratorType;
  typedef ImageRegionIteratorWithIndex< OutputSliceType > OutputSliceIteratorType;
 
 
  unsigned short int m_ZeroPadding;       
  unsigned short int m_OverSampling;      
  double             m_Cutoff;            
  double             m_AlphaRange;        
  unsigned short int m_ZDirection;        
  unsigned short int m_AlphaDirection;    
  unsigned short int m_RDirection;        
  unsigned short int m_RadialSplineOrder; 
  double m_PI; 
  RegionType m_InputRequestedRegion; 
  DirectFourierReconstructionImageToImageFilter(const Self&); 
  void operator=(const Self&); 
};

} // namespace itk

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

#endif /* __itkDirectFourierReconstructionImageToImageFilter_h */