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
|
#ifndef __itkGradientVectorFlowImageFilter_h
#define __itkGradientVectorFlowImageFilter_h
#include "itkImageToImageFilter.h"
#include "vnl/vnl_matrix_fixed.h"
#include "vnl/vnl_math.h"
#include "itkImage.h"
#include "itkVector.h"
#include "itkLaplacianImageFilter.h"
#include "itkImageRegionIterator.h"
#include "itkImageRegionConstIteratorWithIndex.h"
//#include "itkImageRegionConstIterator.h"
namespace itk
{
/** \class GradientVectorFlowImageFilter
* \brief
* This class computes a diffusion of the gradient vectors for graylevel or binary
* edge map derive from the image. It enlarges the capture range of the gradient
* force and make external force derived from the gradient work effectively in the
* framework of deformable model.
*
* \ingroup ImageFilters
* \ingroup ImageSegmentation */
template <class TInputImage, class TOutputImage>
class ITK_EXPORT GradientVectorFlowImageFilter : public ImageToImageFilter<TInputImage, TOutputImage>
{
public:
/** Standard "Self" typedef. */
typedef GradientVectorFlowImageFilter Self;
/** Standard "Superclass" typedef. */
typedef ImageToImageFilter<TInputImage, TOutputImage> Superclass;
/** Smart pointer typedef support */
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Method of creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(GradientVectorFlowImageFilter, ImageToImageFilter);
/** Some typedefs. */
typedef TInputImage InputImageType;
typedef TOutputImage OutputImageType;
typedef typename TInputImage::IndexType IndexType;
typedef typename TInputImage::SizeType SizeType;
typedef typename TInputImage::PixelType PixelType;
typedef typename OutputImageType::Pointer OutputImagePointer;
typedef typename OutputImageType::RegionType RegionType;
/** Image and Image iterator definition. */
// typedef ImageRegionConstIterator<InputImageType> InputImageConstIterator;
typedef ImageRegionIterator<InputImageType> InputImageIterator;
typedef ImageRegionConstIteratorWithIndex<InputImageType> InputImageConstIterator;
typedef ImageRegionIterator<OutputImageType> OutputImageIterator;
/** Image dimension. */
itkStaticConstMacro(ImageDimension, unsigned int,
TInputImage::ImageDimension);
itkStaticConstMacro(OutputImageDimension, unsigned int,
TOutputImage::ImageDimension);
typedef itk::Image<double, itkGetStaticConstMacro(ImageDimension)> InternalImageType;
typedef typename InternalImageType::Pointer InternalImagePointer;
typedef ImageRegionIterator<InternalImageType> InternalImageIterator;
typedef LaplacianImageFilter<InternalImageType, InternalImageType> LaplacianFilterType;
typedef typename LaplacianFilterType::Pointer LaplacianFilterPointer;
/** Routines. */
itkSetMacro(LaplacianFilter, LaplacianFilterPointer);
itkSetMacro(TimeStep, double);
itkGetMacro(TimeStep, double);
itkSetMacro(NoiseLevel, double);
itkGetMacro(NoiseLevel, double);
itkSetMacro(IterationNum, int);
itkGetMacro(IterationNum, int);
#ifdef ITK_USE_CONCEPT_CHECKING
/** Begin concept checking */
itkConceptMacro(SameDimensionCheck,
(Concept::SameDimension<ImageDimension, OutputImageDimension>));
itkConceptMacro(InputHasNumericTraitsCheck,
(Concept::HasNumericTraits<typename PixelType::ValueType>));
itkConceptMacro(OutputHasNumericTraitsCheck,
(Concept::HasNumericTraits<typename TOutputImage::PixelType::ValueType>));
/** End concept checking */
#endif
protected:
GradientVectorFlowImageFilter();
~GradientVectorFlowImageFilter() {}
GradientVectorFlowImageFilter(const Self&) {}
void operator=(const Self&) {}
void PrintSelf(std::ostream& os, Indent indent) const;
virtual void GenerateData();
void InitInterImage();
void UpdateInterImage();
void UpdatePixels();
private:
// parameters;
double m_TimeStep; //the timestep of each iteration
double m_Steps[Superclass::InputImageDimension]; // set to be 1 in all directions in most cases
double m_NoiseLevel; // the noise level of the image
int m_IterationNum; // the iteration number
LaplacianFilterPointer m_LaplacianFilter;
typename Superclass::InputImagePointer m_IntermediateImage;
InternalImagePointer m_InternalImages[Superclass::InputImageDimension];
InternalImagePointer m_BImage; // store the "b" value for every pixel
typename Superclass::InputImagePointer m_CImage; // store the $c_i$ value for every pixel
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkGradientVectorFlowImageFilter.txx"
#endif
#endif
|
