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;

 
};

} // namespace itk

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

#endif /* __itkDirectFourierReconstructionImageToImageFilter_h */