ITK  4.4.0
Insight Segmentation and Registration Toolkit
itkCurvatureAnisotropicDiffusionImageFilter.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 __itkCurvatureAnisotropicDiffusionImageFilter_h
19 #define __itkCurvatureAnisotropicDiffusionImageFilter_h
20 
23 #include "itkMacro.h"
24 namespace itk
25 {
57 template< class TInputImage, class TOutputImage >
59  public AnisotropicDiffusionImageFilter< TInputImage, TOutputImage >
60 {
61 public:
62 
69 
71  itkNewMacro(Self);
72 
76 
78  typedef typename Superclass::UpdateBufferType UpdateBufferType;
79 
81  itkStaticConstMacro(ImageDimension, unsigned int,
82  Superclass::ImageDimension);
83 
84 #ifdef ITK_USE_CONCEPT_CHECKING
85 
86  itkConceptMacro( OutputHasNumericTraitsCheck,
88 
90 #endif
91 
92 protected:
94  {
97  this->SetDifferenceFunction(q);
98  }
99 
101 
102  virtual void InitializeIteration()
103  {
104  Superclass::InitializeIteration();
105  if ( this->GetTimeStep() > 0.5 / vcl_pow( 2.0, static_cast< double >( ImageDimension ) ) )
106  {
107  itkWarningMacro(
108  << "Anisotropic diffusion is using a time step which may introduce instability into the solution.");
109  }
110  }
111 
112 private:
113  CurvatureAnisotropicDiffusionImageFilter(const Self &); //purposely not
114  // implemented
115  void operator=(const Self &); //purposely not
116 
117  // implemented
118 };
119 } // end namspace itk
120 
121 #endif
122