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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkComplexBSplineInterpolateImageFunction.txx
Language: C++
Date: $Date$
Version: $Revision$
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 __itkComplexBSplineInterpolateImageFunction_txx
#define __itkComplexBSplineInterpolateImageFunction_txx
#include "itkComplexBSplineInterpolateImageFunction.h"
namespace itk
{
/**
* Constructor
*/
template <class TImageType, class TCoordRep, class TCoefficientType>
ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>
::ComplexBSplineInterpolateImageFunction()
{
m_RealInterpolator = InterpolatorType::New();
m_ImaginaryInterpolator = InterpolatorType::New();
m_RealFilter = RealFilterType::New();
m_ImaginaryFilter = ImaginaryFilterType::New();
this->SetSplineOrder( 3 );
}
/**
* Standard "PrintSelf" method
*/
template <class TImageType, class TCoordRep, class TCoefficientType>
void ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>
::PrintSelf( std::ostream& os, Indent indent) const
{
Superclass::PrintSelf( os, indent );
os << indent << "Spline Order: " << m_SplineOrder << std::endl;
os << indent << "Real Interpolator: " << m_RealInterpolator << std::endl;
os << indent << "Imaginary Interpolator: " << m_ImaginaryInterpolator << std::endl;
os << indent << "Complex to Real Filter: " << m_RealFilter << std::endl;
os << indent << "Complex to Imaginary Filter: " << m_ImaginaryFilter << std::endl;
}
template <class TImageType, class TCoordRep, class TCoefficientType>
void ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>
::SetInputImage(const TImageType * inputData)
{
if ( inputData )
{
m_RealFilter->SetInput( inputData );
m_ImaginaryFilter->SetInput( inputData );
m_RealInterpolator->SetInputImage( m_RealFilter->GetOutput() );
m_ImaginaryInterpolator->SetInputImage( m_ImaginaryFilter->GetOutput() );
Superclass::SetInputImage( inputData );
}
}
template <class TImageType, class TCoordRep, class TCoefficientType>
void ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>
::SetSplineOrder(unsigned int SplineOrder)
{
m_SplineOrder = SplineOrder;
m_RealInterpolator->SetSplineOrder( SplineOrder );
m_ImaginaryInterpolator->SetSplineOrder( SplineOrder );
}
template <class TImageType, class TCoordRep, class TCoefficientType>
typename ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>::OutputType ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>
::EvaluateAtContinuousIndex( const ContinuousIndexType & x ) const
{
typename InterpolatorType::OutputType realPart =
m_RealInterpolator->EvaluateAtContinuousIndex( x );
typename InterpolatorType::OutputType imagPart =
m_ImaginaryInterpolator->EvaluateAtContinuousIndex( x );
typedef typename ComplexBSplineInterpolateImageFunction<TImageType,TCoordRep,TCoefficientType>::OutputType
ReturnType;
return ReturnType( realPart, imagPart );
}
} // namespace itk
#endif
|
