00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkAnisotropicDiffusionImageFilter_h_
00018 #define __itkAnisotropicDiffusionImageFilter_h_
00019
00020 #include "itkDenseFiniteDifferenceImageFilter.h"
00021 #include "itkAnisotropicDiffusionFunction.h"
00022 #include "itkNumericTraits.h"
00023
00024 namespace itk {
00025
00069 template <class TInputImage, class TOutputImage>
00070 class ITK_EXPORT AnisotropicDiffusionImageFilter
00071 : public DenseFiniteDifferenceImageFilter<TInputImage, TOutputImage>
00072 {
00073 public:
00075 typedef AnisotropicDiffusionImageFilter Self;
00076 typedef DenseFiniteDifferenceImageFilter<TInputImage, TOutputImage>
00077 Superclass;
00078 typedef SmartPointer<Self> Pointer;
00079 typedef SmartPointer<const Self> ConstPointer;
00080
00082 itkTypeMacro(AnisotropicDiffusionImageFilter,
00083 DenseFiniteDifferenceImageFilter);
00084
00086 typedef typename Superclass::InputImageType InputImageType;
00087 typedef typename Superclass::OutputImageType OutputImageType;
00088 typedef typename Superclass::UpdateBufferType UpdateBufferType;
00089
00092 itkStaticConstMacro(ImageDimension, unsigned int,Superclass::ImageDimension);
00093
00096 typedef typename Superclass::PixelType PixelType;
00097 typedef typename Superclass::TimeStepType TimeStepType;
00098
00100 itkSetMacro(TimeStep, TimeStepType);
00101 itkGetMacro(TimeStep, TimeStepType);
00102
00105 itkSetMacro(ConductanceParameter, double);
00106 itkGetMacro(ConductanceParameter, double);
00107
00110 itkSetMacro(ConductanceScalingUpdateInterval, unsigned int);
00111 itkGetMacro(ConductanceScalingUpdateInterval, unsigned int);
00112
00115 itkSetMacro(ConductanceScalingParameter, double);
00116 itkGetMacro(ConductanceScalingParameter, double);
00117
00125 void SetFixedAverageGradientMagnitude(double a)
00126 {
00127 m_FixedAverageGradientMagnitude= a;
00128 this->Modified();
00129 m_GradientMagnitudeIsFixed = true;
00130 }
00131 itkGetMacro(FixedAverageGradientMagnitude, double);
00132
00133 protected:
00134 AnisotropicDiffusionImageFilter();
00135 ~AnisotropicDiffusionImageFilter() {}
00136 void PrintSelf(std::ostream& os, Indent indent) const;
00137
00140
00141
00143 virtual void InitializeIteration();
00144
00145 bool m_GradientMagnitudeIsFixed;
00146
00147 private:
00148 AnisotropicDiffusionImageFilter(const Self&);
00149 void operator=(const Self&);
00150
00151 double m_ConductanceParameter;
00152 double m_ConductanceScalingParameter;
00153 unsigned int m_ConductanceScalingUpdateInterval;
00154 double m_FixedAverageGradientMagnitude;
00155
00156 TimeStepType m_TimeStep;
00157
00158 };
00159
00160 }
00161
00162 #ifndef ITK_MANUAL_INSTANTIATION
00163 #include "itkAnisotropicDiffusionImageFilter.txx"
00164 #endif
00165
00166 #endif