ITK  4.3.0
Insight Segmentation and Registration Toolkit
itkVectorCurvatureNDAnisotropicDiffusionFunction.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 __itkVectorCurvatureNDAnisotropicDiffusionFunction_h
19 #define __itkVectorCurvatureNDAnisotropicDiffusionFunction_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  typedef typename PixelType::ValueType ScalarValueType;
69 
71  itkStaticConstMacro(ImageDimension, unsigned int,
72  Superclass::ImageDimension);
73  itkStaticConstMacro(VectorDimension, unsigned int,
74  Superclass::VectorDimension);
76 
78  virtual PixelType ComputeUpdate(const NeighborhoodType & neighborhood,
79  void *globalData,
80  const FloatOffsetType & offset = FloatOffsetType(0.0)
81  );
82 
84  virtual void InitializeIteration()
85  {
86  m_K = this->GetAverageGradientMagnitudeSquared() * this->GetConductanceParameter()
87  * this->GetConductanceParameter() * -2.0f;
88  }
89 
90 protected:
93  void PrintSelf(std::ostream & os, Indent indent) const
94  { Superclass::PrintSelf(os, indent); }
95 
96 private:
97  VectorCurvatureNDAnisotropicDiffusionFunction(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 
115  double m_K;
116 
117  static double m_MIN_NORM;
119  SizeValueType m_Stride[ImageDimension];
120 };
121 } // end namespace itk
122 
123 #ifndef ITK_MANUAL_INSTANTIATION
124 #include "itkVectorCurvatureNDAnisotropicDiffusionFunction.hxx"
125 #endif
126 
127 #endif
128