AnomalousFiltersExtension/IADImageFilter/itkIsotropicAnomalousDiffusionImageFilter.h at a613296dd05cfe1d1394864f12857d7d8fc4c5a3 · CSIM-Toolkits/AnomalousFiltersExtension · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
#ifndef __itkIsotropicAnomalousDiffusionImageFilter_h
#define __itkIsotropicAnomalousDiffusionImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkImage.h"
#include "itkNumericTraits.h"

namespace itk
{

template< typename TInputImage, typename TOutputImage >
class ITK_EXPORT IsotropicAnomalousDiffusionImageFilter:
  public ImageToImageFilter< TInputImage, TOutputImage >
{
public:
  /** Extract dimension from input and output image. */
  itkStaticConstMacro(InputImageDimension, unsigned int,
                      TInputImage::ImageDimension);
  itkStaticConstMacro(OutputImageDimension, unsigned int,
                      TOutputImage::ImageDimension);

  /** Convenient typedefs for simplifying declarations. */
  typedef TInputImage  InputImageType;
  typedef TOutputImage OutputImageType;

  /** Standard class typedefs. */
  typedef IsotropicAnomalousDiffusionImageFilter                             Self;
  typedef ImageToImageFilter< InputImageType, OutputImageType >                 Superclass;
  typedef SmartPointer< Self >                                                  Pointer;
  typedef SmartPointer< const Self >                                            ConstPointer;

  /** Method for creation through the object factory. */
  itkNewMacro(Self);

  /** Run-time type information (and related methods). */
  itkTypeMacro(IsotropicAnomalousDiffusionImageFilter, ImageToImageFilter);

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

  typedef typename Superclass::OutputImageRegionType OutputImageRegionType;

  itkSetMacro(GeneralizedDiffusion, double);
  itkSetMacro(Iterations, int);
  itkSetMacro(TimeStep, double);
  itkSetMacro(Q, double);
  itkGetMacro(GeneralizedDiffusion, double);
  itkGetMacro(Iterations, int);
  itkGetMacro(TimeStep, double);
  itkGetMacro(Q, double);

#ifdef ITK_USE_CONCEPT_CHECKING
  // Begin concept checking
  itkConceptMacro( InputHasNumericTraitsCheck,
                   ( Concept::HasNumericTraits< InputPixelType > ) );
  itkConceptMacro( SameDimensionCheck,
                   ( Concept::SameDimension< InputImageDimension, OutputImageDimension > ) );
#endif

protected:
  IsotropicAnomalousDiffusionImageFilter();
  virtual ~IsotropicAnomalousDiffusionImageFilter() {}
    double m_GeneralizedDiffusion;
    int m_Iterations;
    double m_TimeStep;
    double m_Q;
    virtual void ThreadedGenerateData(const OutputImageRegionType &, ThreadIdType);

private:
  IsotropicAnomalousDiffusionImageFilter(const Self &); //purposely not implemented
  void operator=(const Self &);  //purposely not implemented
  void TimeStepTestStability();
  double GeneralizedDiffCurve();
};
} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkIsotropicAnomalousDiffusionImageFilter.hxx"
#endif

#endif