ITK  4.2.0
Insight Segmentation and Registration Toolkit
itkLevelSetEvolution.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 
19 
20 #ifndef __itkLevelSetEvolution_h
21 #define __itkLevelSetEvolution_h
22 
23 #include "itkLevelSetEvolutionBase.h"
24 #include "itkLevelSetDenseImage.h"
25 
26 #include "itkWhitakerSparseLevelSetImage.h"
27 #include "itkUpdateWhitakerSparseLevelSet.h"
28 
29 #include "itkShiSparseLevelSetImage.h"
30 #include "itkUpdateShiSparseLevelSet.h"
31 
32 #include "itkMalcolmSparseLevelSetImage.h"
33 #include "itkUpdateMalcolmSparseLevelSet.h"
34 
35 #include "itkLevelSetEvolutionComputeIterationThreader.h"
36 #include "itkLevelSetEvolutionUpdateLevelSetsThreader.h"
37 
38 namespace itk
39 {
49 template< class TEquationContainer, class TLevelSet >
50 class LevelSetEvolution
51 {};
52 
53 template< class TEquationContainer, class TImage >
54 class LevelSetEvolution< TEquationContainer,
55  LevelSetDenseImage< TImage > > :
56  public LevelSetEvolutionBase< TEquationContainer, LevelSetDenseImage< TImage > >
57 {
58 public:
59  typedef LevelSetDenseImage< TImage > LevelSetType;
60 
61  typedef LevelSetEvolution Self;
62  typedef SmartPointer< Self > Pointer;
63  typedef SmartPointer< const Self > ConstPointer;
64  typedef LevelSetEvolutionBase< TEquationContainer, LevelSetType > Superclass;
65 
67  itkNewMacro( Self );
68 
70  itkTypeMacro( LevelSetEvolution, LevelSetEvolutionBase );
71 
72  typedef typename Superclass::EquationContainerType EquationContainerType;
73  typedef typename Superclass::EquationContainerPointer EquationContainerPointer;
74  typedef typename Superclass::TermContainerType TermContainerType;
75  typedef typename Superclass::TermContainerPointer TermContainerPointer;
76 
77  typedef typename Superclass::TermType TermType;
78  typedef typename Superclass::TermPointer TermPointer;
79 
80  typedef typename Superclass::InputImageType InputImageType;
81  typedef typename Superclass::InputImagePixelType InputImagePixelType;
82  typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
83  typedef typename Superclass::InputImageRegionType InputImageRegionType;
84  typedef typename Superclass::InputPixelRealType InputPixelRealType;
85 
86  itkStaticConstMacro ( ImageDimension, unsigned int, Superclass::ImageDimension );
87 
88  typedef typename Superclass::LevelSetContainerType LevelSetContainerType;
89 
90  typedef typename Superclass::LevelSetIdentifierType LevelSetIdentifierType;
91 
92  typedef typename LevelSetType::ImageType LevelSetImageType;
93 
94  typedef typename Superclass::LevelSetOutputType LevelSetOutputType;
95  typedef typename Superclass::LevelSetOutputRealType LevelSetOutputRealType;
96  typedef typename Superclass::LevelSetDataType LevelSetDataType;
97 
98  typedef typename Superclass::IdListType IdListType;
99  typedef typename Superclass::IdListIterator IdListIterator;
100  typedef typename Superclass::IdListConstIterator IdListConstIterator;
101  typedef typename Superclass::IdListImageType IdListImageType;
102  typedef typename Superclass::CacheImageType CacheImageType;
103  typedef typename Superclass::DomainMapImageFilterType DomainMapImageFilterType;
104 
105  typedef typename Superclass::StoppingCriterionType StoppingCriterionType;
106  typedef typename Superclass::StoppingCriterionPointer StoppingCriterionPointer;
107 
108  typedef BinaryThresholdImageFilter< LevelSetImageType, LevelSetImageType > ThresholdFilterType;
109  typedef typename ThresholdFilterType::Pointer ThresholdFilterPointer;
110 
111  typedef SignedMaurerDistanceMapImageFilter< LevelSetImageType, LevelSetImageType > MaurerType;
112  typedef typename MaurerType::Pointer MaurerPointer;
113 
114  typedef ImageRegionIteratorWithIndex< LevelSetImageType > LevelSetImageIteratorType;
115 
116  typedef ImageRegionConstIteratorWithIndex< LevelSetImageType > LevelSetImageConstIteratorType;
117 
118  typedef ImageRegionConstIteratorWithIndex< InputImageType > InputImageConstIteratorType;
119 
120 protected:
121  LevelSetEvolution();
122  ~LevelSetEvolution();
123 
126  virtual void AllocateUpdateBuffer();
127 
129  virtual void ComputeIteration();
130 
132  virtual void ComputeTimeStepForNextIteration();
133 
135  virtual void UpdateLevelSets();
136 
138  virtual void UpdateEquations();
139 
141  void ReinitializeToSignedDistance();
142 
143  typename LevelSetContainerType::Pointer m_UpdateBuffer;
144 
145  friend class LevelSetEvolutionComputeIterationThreader< LevelSetType, ThreadedImageRegionPartitioner< TImage::ImageDimension >, Self >;
146  typedef LevelSetEvolutionComputeIterationThreader< LevelSetType, ThreadedImageRegionPartitioner< TImage::ImageDimension >, Self > SplitLevelSetComputeIterationThreaderType;
147  typename SplitLevelSetComputeIterationThreaderType::Pointer m_SplitLevelSetComputeIterationThreader;
148 
149  typedef typename DomainMapImageFilterType::DomainMapType::const_iterator DomainMapConstIteratorType;
150  typedef ThreadedIteratorRangePartitioner< DomainMapConstIteratorType > ThreadedDomainMapPartitionerType;
151  friend class LevelSetEvolutionComputeIterationThreader< LevelSetType, ThreadedDomainMapPartitionerType, Self >;
152  typedef LevelSetEvolutionComputeIterationThreader< LevelSetType, ThreadedDomainMapPartitionerType, Self > SplitDomainMapComputeIterationThreaderType;
153  typename SplitDomainMapComputeIterationThreaderType::Pointer m_SplitDomainMapComputeIterationThreader;
154 
155  friend class LevelSetEvolutionUpdateLevelSetsThreader< LevelSetType, ThreadedImageRegionPartitioner< TImage::ImageDimension >, Self >;
156  typedef LevelSetEvolutionUpdateLevelSetsThreader< LevelSetType, ThreadedImageRegionPartitioner< TImage::ImageDimension >, Self > SplitLevelSetUpdateLevelSetsThreaderType;
157  typename SplitLevelSetUpdateLevelSetsThreaderType::Pointer m_SplitLevelSetUpdateLevelSetsThreader;
158 
160  const IdListType * m_IdListToProcessWhenThreading;
161 };
162 
163 
164 template< class TEquationContainer, typename TOutput, unsigned int VDimension >
165 class LevelSetEvolution< TEquationContainer, WhitakerSparseLevelSetImage< TOutput, VDimension > > :
166  public LevelSetEvolutionBase< TEquationContainer, WhitakerSparseLevelSetImage< TOutput, VDimension > >
167 {
168 public:
169  typedef WhitakerSparseLevelSetImage< TOutput, VDimension > LevelSetType;
170 
171  typedef LevelSetEvolution Self;
172  typedef SmartPointer< Self > Pointer;
173  typedef SmartPointer< const Self > ConstPointer;
174  typedef LevelSetEvolutionBase< TEquationContainer, LevelSetType > Superclass;
175 
177  itkNewMacro( Self );
178 
180  itkTypeMacro( LevelSetEvolution, LevelSetEvolutionBase );
181 
182  typedef typename Superclass::EquationContainerType EquationContainerType;
183  typedef typename Superclass::EquationContainerPointer EquationContainerPointer;
184  typedef typename Superclass::TermContainerType TermContainerType;
185  typedef typename Superclass::TermContainerPointer TermContainerPointer;
186 
187  typedef typename Superclass::TermType TermType;
188  typedef typename Superclass::TermPointer TermPointer;
189 
190  typedef typename Superclass::InputImageType InputImageType;
191  typedef typename Superclass::InputImagePixelType InputImagePixelType;
192  typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
193  typedef typename Superclass::InputImageRegionType InputImageRegionType;
194  typedef typename Superclass::InputPixelRealType InputPixelRealType;
195 
196  itkStaticConstMacro ( ImageDimension, unsigned int, Superclass::ImageDimension );
197 
198  typedef typename Superclass::LevelSetContainerType LevelSetContainerType;
199  typedef typename Superclass::LevelSetIdentifierType LevelSetIdentifierType;
200 
201  typedef typename Superclass::LevelSetInputType LevelSetInputType;
202  typedef typename Superclass::LevelSetOutputType LevelSetOutputType;
203  typedef typename Superclass::LevelSetOutputRealType LevelSetOutputRealType;
204  typedef typename Superclass::LevelSetDataType LevelSetDataType;
205 
206  typedef typename LevelSetType::LayerType LevelSetLayerType;
207 
208  typedef typename LevelSetType::LabelMapType LevelSetLabelMapType;
209  typedef typename LevelSetType::LabelMapPointer LevelSetLabelMapPointer;
210 
211 
212  typedef typename Superclass::IdListType IdListType;
213  typedef typename Superclass::IdListIterator IdListIterator;
214  typedef typename Superclass::IdListImageType IdListImageType;
215  typedef typename Superclass::CacheImageType CacheImageType;
216  typedef typename Superclass::DomainMapImageFilterType DomainMapImageFilterType;
217 
218  typedef typename Superclass::StoppingCriterionType StoppingCriterionType;
219  typedef typename Superclass::StoppingCriterionPointer StoppingCriterionPointer;
220 
221  typedef ImageRegionConstIteratorWithIndex< InputImageType > InputImageConstIteratorType;
222 
223  typedef UpdateWhitakerSparseLevelSet< ImageDimension, LevelSetOutputType, EquationContainerType > UpdateLevelSetFilterType;
224  typedef typename UpdateLevelSetFilterType::Pointer UpdateLevelSetFilterPointer;
225 
226 protected:
227  LevelSetEvolution();
228  ~LevelSetEvolution();
229 
230  typedef std::pair< LevelSetInputType, LevelSetOutputType > NodePairType;
231 
232  // For sparse case, the update buffer needs to be the size of the active layer
233  std::map< IdentifierType, LevelSetLayerType* > m_UpdateBuffer;
234 
237  virtual void AllocateUpdateBuffer();
238 
240  virtual void ComputeIteration();
241 
243  virtual void ComputeTimeStepForNextIteration();
244 
246  virtual void UpdateLevelSets();
247 
249  virtual void UpdateEquations();
250 
251  typedef ThreadedIteratorRangePartitioner< typename LevelSetType::LayerConstIterator > SplitLevelSetPartitionerType;
252  friend class LevelSetEvolutionComputeIterationThreader< LevelSetType, SplitLevelSetPartitionerType, Self >;
253  typedef LevelSetEvolutionComputeIterationThreader< LevelSetType, SplitLevelSetPartitionerType, Self > SplitLevelSetComputeIterationThreaderType;
254  typename SplitLevelSetComputeIterationThreaderType::Pointer m_SplitLevelSetComputeIterationThreader;
255 
256 private:
257  LevelSetEvolution( const Self& );
258  void operator = ( const Self& );
259 };
260 
261 
262 // Shi
263 template< class TEquationContainer, unsigned int VDimension >
264 class LevelSetEvolution<
265  TEquationContainer,
266  ShiSparseLevelSetImage< VDimension > > :
267 public LevelSetEvolutionBase< TEquationContainer, ShiSparseLevelSetImage< VDimension > >
268 {
269 public:
270  typedef ShiSparseLevelSetImage< VDimension > LevelSetType;
271 
272  typedef LevelSetEvolution Self;
273  typedef SmartPointer< Self > Pointer;
274  typedef SmartPointer< const Self > ConstPointer;
275  typedef LevelSetEvolutionBase< TEquationContainer, LevelSetType > Superclass;
276 
278  itkNewMacro( Self );
279 
281  itkTypeMacro( LevelSetEvolution, LevelSetEvolutionBase );
282 
283  typedef typename Superclass::EquationContainerType EquationContainerType;
284  typedef typename Superclass::EquationContainerPointer EquationContainerPointer;
285  typedef typename Superclass::TermContainerType TermContainerType;
286  typedef typename Superclass::TermContainerPointer TermContainerPointer;
287 
288  typedef typename Superclass::TermType TermType;
289  typedef typename Superclass::TermPointer TermPointer;
290 
291  typedef typename Superclass::InputImageType InputImageType;
292  typedef typename Superclass::InputImagePixelType InputImagePixelType;
293  typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
294  typedef typename Superclass::InputImageRegionType InputImageRegionType;
295  typedef typename Superclass::InputPixelRealType InputPixelRealType;
296 
297  itkStaticConstMacro ( ImageDimension, unsigned int, Superclass::ImageDimension );
298 
299  typedef typename Superclass::LevelSetContainerType LevelSetContainerType;
300  typedef typename Superclass::LevelSetIdentifierType LevelSetIdentifierType;
301 
302  typedef typename Superclass::LevelSetInputType LevelSetInputType;
303  typedef typename Superclass::LevelSetOutputType LevelSetOutputType;
304  typedef typename Superclass::LevelSetOutputRealType LevelSetOutputRealType;
305  typedef typename Superclass::LevelSetDataType LevelSetDataType;
306 
307  typedef typename LevelSetType::LayerType LevelSetLayerType;
308 
309  typedef typename LevelSetType::LabelMapType LevelSetLabelMapType;
310  typedef typename LevelSetType::LabelMapPointer LevelSetLabelMapPointer;
311 
312 
313  typedef typename Superclass::IdListType IdListType;
314  typedef typename Superclass::IdListIterator IdListIterator;
315  typedef typename Superclass::IdListImageType IdListImageType;
316  typedef typename Superclass::CacheImageType CacheImageType;
317  typedef typename Superclass::DomainMapImageFilterType DomainMapImageFilterType;
318 
319  typedef typename Superclass::StoppingCriterionType StoppingCriterionType;
320  typedef typename Superclass::StoppingCriterionPointer StoppingCriterionPointer;
321 
322  typedef ImageRegionConstIteratorWithIndex< InputImageType > InputImageConstIteratorType;
323 
324  typedef UpdateShiSparseLevelSet< ImageDimension, EquationContainerType > UpdateLevelSetFilterType;
325  typedef typename UpdateLevelSetFilterType::Pointer UpdateLevelSetFilterPointer;
326 
327 protected:
328  LevelSetEvolution();
329  ~LevelSetEvolution();
330 
332  virtual void UpdateLevelSets();
333 
335  virtual void UpdateEquations();
336 
337 private:
338  LevelSetEvolution( const Self& );
339  void operator = ( const Self& );
340 };
341 
342 // Malcolm
343 template< class TEquationContainer, unsigned int VDimension >
344 class LevelSetEvolution< TEquationContainer,
345  MalcolmSparseLevelSetImage< VDimension > > :
346 public LevelSetEvolutionBase< TEquationContainer, MalcolmSparseLevelSetImage< VDimension > >
347 {
348 public:
349  typedef MalcolmSparseLevelSetImage< VDimension > LevelSetType;
350 
351  typedef LevelSetEvolution Self;
352  typedef SmartPointer< Self > Pointer;
353  typedef SmartPointer< const Self > ConstPointer;
354  typedef LevelSetEvolutionBase< TEquationContainer, LevelSetType > Superclass;
355 
357  itkNewMacro( Self );
358 
360  itkTypeMacro( LevelSetEvolution, LevelSetEvolutionBase );
361 
362  typedef typename Superclass::EquationContainerType EquationContainerType;
363  typedef typename Superclass::EquationContainerPointer EquationContainerPointer;
364  typedef typename Superclass::TermContainerType TermContainerType;
365  typedef typename Superclass::TermContainerPointer TermContainerPointer;
366 
367  typedef typename Superclass::TermType TermType;
368  typedef typename Superclass::TermPointer TermPointer;
369 
370  typedef typename Superclass::InputImageType InputImageType;
371  typedef typename Superclass::InputImagePixelType InputImagePixelType;
372  typedef typename Superclass::InputImageConstPointer InputImageConstPointer;
373  typedef typename Superclass::InputImageRegionType InputImageRegionType;
374  typedef typename Superclass::InputPixelRealType InputPixelRealType;
375 
376  itkStaticConstMacro ( ImageDimension, unsigned int, Superclass::ImageDimension );
377 
378  typedef typename Superclass::LevelSetContainerType LevelSetContainerType;
379  typedef typename Superclass::LevelSetIdentifierType LevelSetIdentifierType;
380 
381  typedef typename Superclass::LevelSetInputType LevelSetInputType;
382  typedef typename Superclass::LevelSetOutputType LevelSetOutputType;
383  typedef typename Superclass::LevelSetOutputRealType LevelSetOutputRealType;
384  typedef typename Superclass::LevelSetDataType LevelSetDataType;
385 
386  typedef typename LevelSetType::LayerType LevelSetLayerType;
387  typedef typename LevelSetType::LayerIterator LevelSetLayerIterator;
388 
389  typedef typename LevelSetType::LabelMapType LevelSetLabelMapType;
390  typedef typename LevelSetType::LabelMapPointer LevelSetLabelMapPointer;
391 
392 
393  typedef typename Superclass::IdListType IdListType;
394  typedef typename Superclass::IdListIterator IdListIterator;
395  typedef typename Superclass::IdListImageType IdListImageType;
396  typedef typename Superclass::CacheImageType CacheImageType;
397  typedef typename Superclass::DomainMapImageFilterType DomainMapImageFilterType;
398 
399  typedef typename Superclass::StoppingCriterionType StoppingCriterionType;
400  typedef typename Superclass::StoppingCriterionPointer StoppingCriterionPointer;
401 
402  typedef ImageRegionConstIteratorWithIndex< InputImageType > InputImageConstIteratorType;
403 
404  typedef UpdateMalcolmSparseLevelSet< ImageDimension, EquationContainerType > UpdateLevelSetFilterType;
405  typedef typename UpdateLevelSetFilterType::Pointer UpdateLevelSetFilterPointer;
406 
407 protected:
408  LevelSetEvolution();
409  virtual ~LevelSetEvolution();
410 
411  virtual void UpdateLevelSets();
412 
413  virtual void UpdateEquations();
414 
415 private:
416  LevelSetEvolution( const Self& ); // purposely not implemented
417  void operator = ( const Self& ); // purposely not implemented
418 };
419 }
420 
421 #ifndef ITK_MANUAL_INSTANTIATION
422 #include "itkLevelSetEvolution.hxx"
423 #endif
424 
425 #endif // __itkLevelSetEvolution_h
426 
Tarballs of the nightly generated Doxygen documentation are available for the html, xml, and tag file.
Generated on Tue Jul 10 2012 23:35:13 for ITK by doxygen 1.8.1