ITK  4.6.0
Insight Segmentation and Registration Toolkit
itkVectorCurvatureAnisotropicDiffusionImageFilter.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 __itkVectorCurvatureAnisotropicDiffusionImageFilter_h
19 #define __itkVectorCurvatureAnisotropicDiffusionImageFilter_h
20 
21 #include "itkMacro.h"
24 
25 namespace itk
26 {
62 template< typename TInputImage, typename TOutputImage >
64  public AnisotropicDiffusionImageFilter< TInputImage, TOutputImage >
65 {
66 public:
73 
75  itkNewMacro(Self);
76 
80 
83 
85  itkStaticConstMacro(ImageDimension, unsigned int,
87 
88 #ifdef ITK_USE_CONCEPT_CHECKING
89  // Begin concept checking
90  itkConceptMacro( InputHasNumericTraitsCheck,
92  itkConceptMacro( OutputHasNumericTraitsCheck,
94  // End concept checking
95 #endif
96 
97 protected:
99  {
102  this->SetDifferenceFunction(q);
103  }
104 
106 
107  virtual void InitializeIteration()
108  {
110  if ( this->GetTimeStep() > 0.5 / std::pow( 2.0, static_cast< double >( ImageDimension ) ) )
111  {
112  itkWarningMacro(
113  << "Anisotropic diffusion has attempted to use a time step which may introduce instability into the solution.");
114  }
115  }
116 
117 private:
118  VectorCurvatureAnisotropicDiffusionImageFilter(const Self &); //purposely not
119  // implemented
120  void operator=(const Self &); //purposely not
121 
122  // implemented
123 };
124 } // end namspace itk
125 
126 #endif
Base class for all process objects that output image data.
virtual void SetDifferenceFunction(FiniteDifferenceFunctionType *_arg)
AnisotropicDiffusionImageFilter< TInputImage, TOutputImage > Superclass
#define itkConceptMacro(name, concept)
virtual TimeStepType GetTimeStep() const