ITK  4.2.0
Insight Segmentation and Registration Toolkit
itkCurvatureNDAnisotropicDiffusionFunction.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 __itkCurvatureNDAnisotropicDiffusionFunction_h
19 #define __itkCurvatureNDAnisotropicDiffusionFunction_h
20 
24 #include "itkDerivativeOperator.h"
25 
26 namespace itk
27 {
74 template< class TImage >
76  public ScalarAnisotropicDiffusionFunction< TImage >
77 {
78 public:
79 
85 
87  itkNewMacro(Self);
88 
92 
94  typedef typename Superclass::ImageType ImageType;
95  typedef typename Superclass::PixelType PixelType;
100 
102 
104  itkStaticConstMacro(ImageDimension, unsigned int, Superclass::ImageDimension);
105 
107  virtual PixelType ComputeUpdate(const NeighborhoodType & neighborhood,
108  void *globalData,
109  const FloatOffsetType & offset = FloatOffsetType(0.0)
110  );
111 
113  virtual void InitializeIteration()
114  {
115  m_K = static_cast< PixelType >( this->GetAverageGradientMagnitudeSquared()
116  * this->GetConductanceParameter()
117  * this->GetConductanceParameter() * -2.0f );
118  }
119 
120 protected:
123 private:
124  CurvatureNDAnisotropicDiffusionFunction(const Self &); //purposely not
125  // implemented
126  void operator=(const Self &); //purposely not
127 
128  // implemented
129 
132 
134  std::slice x_slice[ImageDimension];
135  std::slice xa_slice[ImageDimension][ImageDimension];
136  std::slice xd_slice[ImageDimension][ImageDimension];
137 
140 
143 
145  static double m_MIN_NORM;
146 
148  NeighborhoodSizeValueType m_Stride[ImageDimension];
149 };
150 } // end namespace itk
151 
152 #ifndef ITK_MANUAL_INSTANTIATION
153 #include "itkCurvatureNDAnisotropicDiffusionFunction.hxx"
154 #endif
155 
156 #endif
157