itkSHCoefficientsToGFAImageFilter.h

Go to the documentation of this file.

#ifndef __itkSHCoefficientsToGFAImageFilter_h
#define __itkSHCoefficientsToGFAImageFilter_h

#include "itkUnaryFunctorImageFilter.h"
#include <cmath>

namespace itk
{
namespace Functor
{
template <class TInput, class TOutput=double>
class SHCoefficientsToGFA
{
public:

  SHCoefficientsToGFA()
  {
  }

  ~SHCoefficientsToGFA()
  {
  }

  bool operator!=(const SHCoefficientsToGFA & other) const
  {
    return !( *this == other );
  }

  bool operator==(const SHCoefficientsToGFA & other) const
  {
    return false;
  }

  inline TOutput operator()(const TInput & A)
  {
  double norm2 = A.GetSquaredNorm();
  if (norm2>0)
    return std::sqrt(1 - A[0]*A[0]/norm2 );
  else
    return 0;
  }


};
}

template <class TInputImage, class TOutputImage=itk::Image<double,3> >
class ITK_EXPORT SHCoefficientsToGFAImageFilter :
  public
  UnaryFunctorImageFilter<TInputImage, TOutputImage,
                          Functor::SHCoefficientsToGFA<
                            typename TInputImage::PixelType,
                            typename TOutputImage::PixelType> >
{
public:
  typedef SHCoefficientsToGFAImageFilter Self;
  typedef UnaryFunctorImageFilter<
      TInputImage, TOutputImage,
      Functor::SHCoefficientsToGFA<
        typename TInputImage::PixelType,  typename TOutputImage::PixelType> 
      >             Superclass;

  typedef SmartPointer<Self>       Pointer;
  typedef SmartPointer<const Self> ConstPointer;

  itkNewMacro( Self );

  itkTypeMacro( SHCoefficientsToGFAImageFilter, UnaryFunctorImageFilter );

protected:
  SHCoefficientsToGFAImageFilter()
  {
  }

  virtual ~SHCoefficientsToGFAImageFilter()
  {
  }

private:
  SHCoefficientsToGFAImageFilter(const Self &); // purposely not implemented
  void operator=(const Self &);                      // purposely not implemented
};
} // end namespace itk


#endif