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
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkSymmetricEigenVectorAnalysisImageFilter.h,v $
Language: C++
Date: $Date: 2007/03/30 17:58:14 $
Version: $Revision: 1.2 $
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 __itkSymmetricEigenVectorAnalysisImageFilter_h
#define __itkSymmetricEigenVectorAnalysisImageFilter_h
#include "itkUnaryFunctorImageFilter.h"
#include "itkSymmetricEigenAnalysis.h"
namespace itk
{
// This functor class invokes the computation of Eigen Analysis for
// every pixel. The input pixel type must provide the API for the [][]
// operator, while the output pixel type must provide the API for the
// [] operator. Input pixel matrices should be symmetric.
//
// The default operation is to order eigen values in ascending order.
// You may also use OrderEigenValuesBy( ) to order eigen values by
// magnitude as is common with use of tensors in vessel extraction.
namespace Functor {
template< typename TInput, typename TOutput, typename TMatrix >
class SymmetricEigenVectorAnalysisFunction
{
public:
SymmetricEigenVectorAnalysisFunction() {}
~SymmetricEigenVectorAnalysisFunction() {}
typedef SymmetricEigenAnalysis< TInput, TOutput, TMatrix > CalculatorType;
bool operator!=( const SymmetricEigenVectorAnalysisFunction & ) const
{
return false;
}
bool operator==( const SymmetricEigenVectorAnalysisFunction & other ) const
{
return !(*this != other);
}
inline TMatrix operator()( const TInput & x )
{
TOutput eigenValues;
TMatrix eigenVectorMatrix;
m_Calculator.ComputeEigenValuesAndVectors( x, eigenValues, eigenVectorMatrix );
return eigenVectorMatrix;
}
/** Method to explicitly set the dimension of the matrix */
void SetDimension( unsigned int n )
{
m_Calculator.SetDimension(n);
}
/** Typdedefs to order eigen values.
* OrderByValue: lambda_1 < lambda_2 < ....
* OrderByMagnitude: |lambda_1| < |lambda_2| < .....
* DoNotOrder: Default order of eigen values obtained after QL method
*/
typedef enum {
OrderByValue=1,
OrderByMagnitude,
DoNotOrder
}EigenValueOrderType;
/** Order eigen values. Default is to OrderByValue: lambda_1 < lambda_2 < ....*/
void OrderEigenValuesBy( EigenValueOrderType order )
{
if( order == OrderByMagnitude )
{
m_Calculator.SetOrderEigenMagnitudes( true );
}
else if( order == DoNotOrder )
{
m_Calculator.SetOrderEigenValues( false );
}
}
private:
CalculatorType m_Calculator;
};
} // end namespace functor
/** \class SymmetricEigenVectorAnalysisImageFilter
*
* \ingroup IntensityImageFilters Multithreaded TensorObjects
*
*/
template <typename TInputImage, typename TOutputImage, typename TOutputMatrix>
class ITK_EXPORT SymmetricEigenVectorAnalysisImageFilter :
public
UnaryFunctorImageFilter<TInputImage,TOutputMatrix,
Functor::SymmetricEigenVectorAnalysisFunction<
typename TInputImage::PixelType,
typename TOutputImage::PixelType,
typename TOutputMatrix::PixelType > >
{
public:
/** Standard class typedefs. */
typedef SymmetricEigenVectorAnalysisImageFilter Self;
typedef UnaryFunctorImageFilter<TInputImage,TOutputMatrix,
Functor::SymmetricEigenVectorAnalysisFunction<
typename TInputImage::PixelType,
typename TOutputImage::PixelType,
typename TOutputMatrix::PixelType > >
Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
typedef typename Superclass::OutputImageType OutputImageType;
typedef typename TOutputImage::PixelType OutputPixelType;
typedef typename TInputImage::PixelType InputPixelType;
typedef typename Superclass::FunctorType FunctorType;
/** Typdedefs to order eigen values.
* OrderByValue: lambda_1 < lambda_2 < ....
* OrderByMagnitude: |lambda_1| < |lambda_2| < .....
* DoNotOrder: Default order of eigen values obtained after QL method
*/
typedef typename FunctorType::EigenValueOrderType EigenValueOrderType;
/** Order eigen values. Default is to OrderByValue: lambda_1 < lambda_2 < ....*/
void OrderEigenValuesBy( EigenValueOrderType order )
{
this->GetFunctor().OrderEigenValuesBy( order );
}
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Print internal ivars */
void PrintSelf(std::ostream& os, Indent indent) const
{ this->Superclass::PrintSelf( os, indent ); }
/** Set the dimension of the tensor. (For example the SymmetricSecondRankTensor
* is a pxp matrix) */
void SetDimension( unsigned int p )
{
this->GetFunctor().SetDimension(p);
}
protected:
SymmetricEigenVectorAnalysisImageFilter() {};
virtual ~SymmetricEigenVectorAnalysisImageFilter() {};
private:
SymmetricEigenVectorAnalysisImageFilter(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
};
} // end namespace itk
#endif
|
