itkRegularizedIVIMLocalVariationImageFilter.h

Go to the documentation of this file.

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

The Medical Imaging Interaction Toolkit (MITK)

Copyright (c) German Cancer Research Center,
Division of Medical and Biological Informatics.
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 or http://www.mitk.org for details.

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

#ifndef __itkRegularizedIVIMLocalVariationImageFilter_h
#define __itkRegularizedIVIMLocalVariationImageFilter_h

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

namespace itk
{

  template<class TPixelType>
  class IVIMSquaredEuclideanMetric
  {
  public:
    static double Calc(TPixelType p)
    {
      return p*p;
    }
  };

  template<>
  class IVIMSquaredEuclideanMetric<itk::Vector<double,3> >
  {
  public:
    static double Calc(itk::Vector<double,3> p)
    {
      return p[1]*p[1];
    }
  };

  template<>
  class IVIMSquaredEuclideanMetric<itk::VariableLengthVector<float> >
  {
  public:
    static double Calc(itk::VariableLengthVector<float> p)
    {
      return p.GetSquaredNorm();
    }
  };

  template<>
  class IVIMSquaredEuclideanMetric<itk::VariableLengthVector<double> >
  {
  public:
    static double Calc(itk::VariableLengthVector<double> p)
    {
      return p.GetSquaredNorm();
    }
  };

template <class TInputImage, class TOutputImage>
class RegularizedIVIMLocalVariationImageFilter :
    public ImageToImageFilter< TInputImage, TOutputImage >
{
public:
  itkStaticConstMacro(InputImageDimension, unsigned int,
                      TInputImage::ImageDimension);
  itkStaticConstMacro(OutputImageDimension, unsigned int,
                      TOutputImage::ImageDimension);

  typedef TInputImage InputImageType;
  typedef TOutputImage OutputImageType;

  typedef RegularizedIVIMLocalVariationImageFilter Self;
  typedef ImageToImageFilter< InputImageType, OutputImageType> Superclass;
  typedef SmartPointer<Self> Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  itkFactorylessNewMacro(Self)
  itkCloneMacro(Self)

  itkTypeMacro(RegularizedIVIMLocalVariationImageFilter, ImageToImageFilter);

  typedef typename InputImageType::PixelType InputPixelType;
  typedef typename OutputImageType::PixelType OutputPixelType;

  typedef typename InputImageType::RegionType InputImageRegionType;
  typedef typename OutputImageType::RegionType OutputImageRegionType;

  typedef typename InputImageType::SizeType InputSizeType;

  virtual void GenerateInputRequestedRegion()
    throw(InvalidRequestedRegionError);

protected:
  RegularizedIVIMLocalVariationImageFilter();
  virtual ~RegularizedIVIMLocalVariationImageFilter() {}
  void PrintSelf(std::ostream& os, Indent indent) const;

  void ThreadedGenerateData(const OutputImageRegionType& outputRegionForThread,
                            ThreadIdType threadId );

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

} // end namespace itk

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

#endif //RegularizedIVIMLocalVariationImageFilter