ITK  4.3.0
Insight Segmentation and Registration Toolkit
itkCurvesLevelSetFunction.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 __itkCurvesLevelSetFunction_h
19 #define __itkCurvesLevelSetFunction_h
20 
22 
23 namespace itk
24 {
61 template< class TImageType, class TFeatureImageType = TImageType >
62 class ITK_EXPORT CurvesLevelSetFunction:
63  public SegmentationLevelSetFunction< TImageType, TFeatureImageType >
64 {
65 public:
72  typedef TFeatureImageType FeatureImageType;
73 
75  itkNewMacro(Self);
76 
79 
81  typedef typename SuperSuperclass::PixelType PixelType;
82  typedef typename Superclass::ImageType ImageType;
83  typedef typename Superclass::NeighborhoodType NeighborhoodType;
84  typedef typename Superclass::ScalarValueType ScalarValueType;
85  typedef typename Superclass::FeatureScalarType FeatureScalarType;
86  typedef typename Superclass::RadiusType RadiusType;
89  typedef typename Superclass::VectorImageType VectorImageType;
90 
92  itkStaticConstMacro(ImageDimension, unsigned int,
93  Superclass::ImageDimension);
94 
96  virtual void CalculateSpeedImage();
97 
99  virtual void CalculateAdvectionImage();
100 
102  virtual ScalarValueType CurvatureSpeed(const NeighborhoodType & neighborhood,
103  const FloatOffsetType & offset, GlobalDataStruct *gd) const
104  { return this->PropagationSpeed(neighborhood, offset, gd); }
105 
108  void SetDerivativeSigma(const double v)
109  { m_DerivativeSigma = v; }
110  double GetDerivativeSigma()
111  { return m_DerivativeSigma; }
113 
114  virtual void Initialize(const RadiusType & r);
115 
116 protected:
118  {
119  //Curvature term is the minimal curvature.
120  this->UseMinimalCurvatureOn();
121  this->SetAdvectionWeight(NumericTraits< ScalarValueType >::One);
122  this->SetPropagationWeight(NumericTraits< ScalarValueType >::One);
123  this->SetCurvatureWeight(NumericTraits< ScalarValueType >::One);
124 
125  m_DerivativeSigma = 1.0;
126  }
127 
129 
130  CurvesLevelSetFunction(const Self &); //purposely not implemented
131  void operator=(const Self &); //purposely not implemented
132 
133  void PrintSelf(std::ostream & os, Indent indent) const
134  {
135  Superclass::PrintSelf(os, indent);
136  os << indent << "DerivativeSigma: " << m_DerivativeSigma << std::endl;
137  }
138 
139 private:
140 
142  std::slice x_slice[ImageDimension];
143 
146 
148  OffsetValueType m_xStride[ImageDimension];
149 
151 };
152 } // end namespace itk
153 
154 #ifndef ITK_MANUAL_INSTANTIATION
155 #include "itkCurvesLevelSetFunction.hxx"
156 #endif
157 
158 #endif
159