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
|
/*=========================================================================
Program: ITK-SNAP
Module: $RCSfile: SNAPAdvectionFieldImageFilter.h,v $
Language: C++
Date: $Date: 2009/01/23 20:09:38 $
Version: $Revision: 1.3 $
Copyright (c) 2007 Paul A. Yushkevich
This file is part of ITK-SNAP
ITK-SNAP is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program. If not, see <http://www.gnu.org/licenses/>.
-----
Copyright (c) 2003 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 __SNAPAdvectionFieldImageFilter_h_
#define __SNAPAdvectionFieldImageFilter_h_
#include "itkCovariantVector.h"
#include "itkImage.h"
#include "itkImageToImageFilter.h"
/**
* \class SNAPAdvectionFieldImageFilter
* \brief A filter used to compute the advection field in the SNAP level set
* equation.
*/
template <class TInputImage, class TOutputValueType=float>
class SNAPAdvectionFieldImageFilter:
public itk::ImageToImageFilter<
TInputImage,
itk::Image<itk::CovariantVector<TOutputValueType, TInputImage::ImageDimension>,
TInputImage::ImageDimension> >
{
public:
/** Input image types */
typedef TInputImage InputImageType;
typedef itk::SmartPointer<InputImageType> InputImagePointer;
/** Standard class typedefs. */
typedef SNAPAdvectionFieldImageFilter Self;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Image dimension. */
itkStaticConstMacro(ImageDimension, unsigned int,
TInputImage::ImageDimension);
/** Output image types */
typedef itk::CovariantVector<
TOutputValueType,
itkGetStaticConstMacro(ImageDimension)> VectorType;
typedef itk::Image<
VectorType,
itkGetStaticConstMacro(ImageDimension)> OutputImageType;
typedef itk::SmartPointer<OutputImageType> OutputImagePointer;
typedef itk::ImageToImageFilter<InputImageType,OutputImageType> Superclass;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Set the power of g() by which the gradient is scaled */
itkSetMacro(Exponent,unsigned int);
/** Get the power of g() by which the gradient is scaled */
itkGetMacro(Exponent,unsigned int);
protected:
SNAPAdvectionFieldImageFilter();
virtual ~SNAPAdvectionFieldImageFilter() {};
void PrintSelf(std::ostream& os, itk::Indent indent) const;
/** Generate Data */
void GenerateData( void );
private:
/** The g-scaling exponent */
unsigned int m_Exponent;
};
#ifndef ITK_MANUAL_INSTANTIATION
#include "SNAPAdvectionFieldImageFilter.txx"
#endif
#endif // __SNAPAdvectionFieldImageFilter_h_
|
