ITK  4.2.0
Insight Segmentation and Registration Toolkit
itkRecursiveGaussianImageFilter.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 __itkRecursiveGaussianImageFilter_h
19 #define __itkRecursiveGaussianImageFilter_h
20 
22 
23 namespace itk
24 {
61 template< typename TInputImage, typename TOutputImage = TInputImage >
63  public RecursiveSeparableImageFilter< TInputImage, TOutputImage >
64 {
65 public:
71 
72  typedef typename Superclass::RealType RealType;
73  typedef typename Superclass::ScalarRealType ScalarRealType;
74 
76  itkNewMacro(Self);
77 
80 
83  itkGetConstMacro(Sigma, ScalarRealType);
84  itkSetMacro(Sigma, ScalarRealType);
86 
90  typedef enum { ZeroOrder, FirstOrder, SecondOrder } OrderEnumType;
91 
93  typedef TOutputImage OutputImageType;
94 
134  itkSetMacro(NormalizeAcrossScale, bool);
135  itkGetConstMacro(NormalizeAcrossScale, bool);
136 
143  itkSetMacro(Order, OrderEnumType);
144  itkGetConstMacro(Order, OrderEnumType);
146 
148  void SetZeroOrder();
149 
151  void SetFirstOrder();
152 
154  void SetSecondOrder();
155 
156 protected:
159  void PrintSelf(std::ostream & os, Indent indent) const;
160 
165  virtual void SetUp(ScalarRealType spacing);
166 
167 private:
168  RecursiveGaussianImageFilter(const Self &); //purposely not implemented
169  void operator=(const Self &); //purposely not implemented
170 
172  void ComputeNCoefficients(ScalarRealType sigmad,
173  ScalarRealType A1, ScalarRealType B1, ScalarRealType W1, ScalarRealType L1,
174  ScalarRealType A2, ScalarRealType B2, ScalarRealType W2, ScalarRealType L2,
175  ScalarRealType & N0, ScalarRealType & N1,
176  ScalarRealType & N2, ScalarRealType & N3,
177  ScalarRealType & SN, ScalarRealType & DN, ScalarRealType & EN);
178 
180  void ComputeDCoefficients(ScalarRealType sigmad,
181  ScalarRealType W1, ScalarRealType L1, ScalarRealType W2, ScalarRealType L2,
182  ScalarRealType & SD, ScalarRealType & DD, ScalarRealType & ED);
183 
186  void ComputeRemainingCoefficients(bool symmetric);
187 
190 
193 
195 };
196 } // end namespace itk
197 
198 #ifndef ITK_MANUAL_INSTANTIATION
199 #include "itkRecursiveGaussianImageFilter.hxx"
200 #endif
201 
202 #endif
203