ITK  4.4.0
Insight Segmentation and Registration Toolkit
itkVectorGradientNDAnisotropicDiffusionFunction.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright Insight Software Consortium
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 #ifndef __itkVectorGradientNDAnisotropicDiffusionFunction_h
19 #define __itkVectorGradientNDAnisotropicDiffusionFunction_h
20 
24 #include "itkDerivativeOperator.h"
25 
26 namespace itk
27 {
43 template< class TImage >
45  public VectorAnisotropicDiffusionFunction< TImage >
46 {
47 public:
53 
55  itkNewMacro(Self);
56 
60 
62  typedef typename Superclass::ImageType ImageType;
63  typedef typename Superclass::PixelType PixelType;
68 
70  itkStaticConstMacro(ImageDimension, unsigned int,
71  Superclass::ImageDimension);
72  itkStaticConstMacro(VectorDimension, unsigned int,
73  Superclass::VectorDimension);
75 
77  typedef typename PixelType::ValueType ScalarValueType;
78 
80  virtual PixelType ComputeUpdate(const NeighborhoodType & neighborhood,
81  void *globalData,
82  const FloatOffsetType & offset = FloatOffsetType(0.0)
83  );
84 
86  virtual void InitializeIteration()
87  {
88  m_K = this->GetAverageGradientMagnitudeSquared() * this->GetConductanceParameter()
89  * this->GetConductanceParameter() * -2.0f;
90  }
91 
92 protected:
95 
96 private:
97  VectorGradientNDAnisotropicDiffusionFunction(const Self &); //purposely not
98  // implemented
99  void operator=(const Self &); //purposely not
100 
101  // implemented
102 
105 
107  std::slice x_slice[ImageDimension];
108  std::slice xa_slice[ImageDimension][ImageDimension];
109  std::slice xd_slice[ImageDimension][ImageDimension];
110 
113  itkGetStaticConstMacro(ImageDimension) > dx_op;
114 
117 
118  static double m_MIN_NORM;
119 
120  SizeValueType m_Stride[ImageDimension];
122 };
123 } // end namespace itk
124 
125 #ifndef ITK_MANUAL_INSTANTIATION
126 #include "itkVectorGradientNDAnisotropicDiffusionFunction.hxx"
127 #endif
128 
129 #endif
130