ITK  4.2.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 private:
96  VectorCurvatureNDAnisotropicDiffusionFunction(const Self &); //purposely not
97  // implemented
98  void operator=(const Self &); //purposely not
99 
100  // implemented
101 
104 
106  std::slice x_slice[ImageDimension];
107  std::slice xa_slice[ImageDimension][ImageDimension];
108  std::slice xd_slice[ImageDimension][ImageDimension];
109 
112 
114  double m_K;
115 
116  static double m_MIN_NORM;
118  SizeValueType m_Stride[ImageDimension];
119 };
120 } // end namespace itk
121 
122 #ifndef ITK_MANUAL_INSTANTIATION
123 #include "itkVectorCurvatureNDAnisotropicDiffusionFunction.hxx"
124 #endif
125 
126 #endif
127