1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkLevelSetVelocityNeighborhoodExtractor.h,v $
Language: C++
Date: $Date: 2003-09-10 14:28:33 $
Version: $Revision: 1.21 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#ifndef _itkLevelSetVelocityNeighborhoodExtractor_h
#define _itkLevelSetVelocityNeighborhoodExtractor_h
#include "itkLevelSetNeighborhoodExtractor.h"
#include "itkLevelSet.h"
#include "itkIndex.h"
namespace itk
{
/** \class LevelSetVelocityNeighborhoodExtractor
* \brief Locate pixels of a particular level set.
*
* LevelSetVelocityNeighborhoodExtractor extends the functionality of
* LevelSetNeighborhoodExtractor by also extracting the values
* of velocity variables at the specified level set. Specifically,
* it populates two containers: one containing the value of velocity
* variables immediately inside the contour defined by the level set and the
* other containing values for velocity variables immediately outside.
*
* The containers AuxInsideValues() and AuxOutsideValues() can
* be used in conjunction with Superclass::InsidePoints() and
* Superclass::OutsidePoints() in FastMarchingExtensionImageFilter
* to produce images which extends the velocity variables smoothly
* from the specified level set.
*
* This class is templated over the image type representing
* the level set, the type of the auxiliary/velocity variables and the
* number of auxiliary/velocity variables.
*
* Implemenation of this class is based on Chapter 11 of
* "Level Set Methods and Fast Marching Methods", J.A. Sethian,
* Cambridge Press, Second edition, 1999.
*
* \ingroup LevelSetSegmentation
*
*/
template <
class TLevelSet,
class TAuxValue,
unsigned int VAuxDimension = 1
>
class ITK_EXPORT LevelSetVelocityNeighborhoodExtractor :
public LevelSetNeighborhoodExtractor<TLevelSet>
{
public:
/** Standard class typdedefs. */
typedef LevelSetVelocityNeighborhoodExtractor Self;
typedef LevelSetNeighborhoodExtractor<TLevelSet> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(LevelSetVelocityNeighborhoodExtractor,
LevelSetNeighborhoodExtractor);
/** The type of the level set. */
typedef LevelSetTypeDefault<TLevelSet> LevelSetType;
/** The dimension of the level set. */
itkStaticConstMacro(SetDimension, unsigned int,
LevelSetType::SetDimension);
/** Index typedef support */
typedef Index<itkGetStaticConstMacro(SetDimension)> Index;
/** AuxVarType typedef support. */
typedef AuxVarTypeDefault<TAuxValue,VAuxDimension,itkGetStaticConstMacro(SetDimension)> AuxVarType;
typedef typename AuxVarType::AuxValueType AuxValueType;
typedef typename AuxVarType::AuxValueVectorType AuxValueVectorType;
typedef typename AuxVarType::AuxValueContainer AuxValueContainer;
typedef typename AuxVarType::AuxImageType AuxImageType;
typedef typename AuxVarType::AuxImagePointer AuxImagePointer;
typedef typename AuxVarType::AuxImageConstPointer AuxImageConstPointer;
/** Set the auxiliary images. */
void SetAuxImage( const AuxImageType * ptr, unsigned int idx = 0 )
{
if ( idx < VAuxDimension && m_AuxImage[idx] != ptr )
{
m_AuxImage[idx] = ptr;
}
this->Modified();
}
/** Get the auxiliary images. */
AuxImageConstPointer GetAuxImage( unsigned int idx = 0 )
{
if ( idx >= VAuxDimension )
{
return NULL;
}
else
{
return m_AuxImage[idx];
}
}
/** Get the container of auxiliary values associated with the inside
* points. */
itkGetObjectMacro( AuxInsideValues, AuxValueContainer );
/** Get the container of auxiliary values associate with the outside
* points. */
itkGetObjectMacro( AuxOutsideValues, AuxValueContainer );
protected:
LevelSetVelocityNeighborhoodExtractor();
~LevelSetVelocityNeighborhoodExtractor(){};
void PrintSelf( std::ostream& os, Indent indent ) const;
virtual void Initialize();
virtual double CalculateDistance( Index & index );
private:
LevelSetVelocityNeighborhoodExtractor(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
typename AuxValueContainer::Pointer m_AuxInsideValues;
typename AuxValueContainer::Pointer m_AuxOutsideValues;
AuxImageConstPointer m_AuxImage[VAuxDimension];
};
} // namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkLevelSetVelocityNeighborhoodExtractor.txx"
#endif
#endif
|
