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
172
173
174
175
176
177
178
179
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkDeformationFieldSource.h,v $
Language: C++
Date: $Date: 2005-07-26 15:06:23 $
Version: $Revision: 1.3 $
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 __itkDeformationFieldSource_h
#define __itkDeformationFieldSource_h
#include "itkImageSource.h"
#include "itkKernelTransform.h"
#include "itkImageRegionIteratorWithIndex.h"
namespace itk
{
/** \class DeformationFieldSource
* \brief Computes a deformation field from two sets of landmarks.
*
* DeformationFieldSource produces a deformation field from two set of input
* landmarks. One set of landmarks are associated to the input space while the
* second set of landmarks is associated with the output space.
*
* A KernelBase spline is used to interpolate the deformations and produce
* deformation values for all the nodes of the image grid that will be produced
* as output.
*
* The number of landmarks in the KernelBased spline will have a dramatic
* effect on both the precision of output deformation field and the
* computational time required for the filter to complete the estimation.
*
*
* This source object expects the image to be of pixel type Vector.
*
* \ingroup ImageSource
*/
template <class TOutputImage>
class ITK_EXPORT DeformationFieldSource:
public ImageSource<TOutputImage>
{
public:
/** Standard class typedefs. */
typedef DeformationFieldSource Self;
typedef ImageSource<TOutputImage> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
typedef TOutputImage OutputImageType;
typedef typename OutputImageType::Pointer OutputImagePointer;
typedef typename OutputImageType::RegionType OutputImageRegionType;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(DeformationFieldSource, ImageToImageFilter);
/** Number of dimensions. */
itkStaticConstMacro(ImageDimension, unsigned int,
TOutputImage::ImageDimension);
/** Transform typedef.
*
* The KernelBased spline transform types are defined here.
*/
typedef KernelTransform<double, itkGetStaticConstMacro(ImageDimension)> KernelTransformType;
typedef typename KernelTransformType::PointSetType LandmarkPointSetType;
typedef typename LandmarkPointSetType::PointType LandmarkPointType;
typedef typename KernelTransformType::Pointer KernelTransformPointerType;
typedef typename KernelTransformType::PointsContainer LandmarkContainer;
typedef typename LandmarkContainer::ConstPointer LandmarkContainerPointer;
/** Image size typedef. */
typedef typename OutputImageType::SizeType OutputSizeType;
/** Image index typedef. */
typedef typename OutputImageType::IndexType OutputIndexType;
/** Image pixel value typedef. */
typedef typename TOutputImage::PixelType OutputPixelType;
typedef typename OutputPixelType::ValueType OutputPixelComponentType;
/** Image spacing typedef */
typedef typename TOutputImage::SpacingType SpacingType;
typedef typename TOutputImage::PointType OriginPointType;
/** Set the coordinate transformation.
* Set the KernelBase spline used for resampling the deformation grid.
* */
itkSetObjectMacro( KernelTransform, KernelTransformType );
/** Get a pointer to the coordinate transform. */
itkGetObjectMacro( KernelTransform, KernelTransformType );
/** Set the size of the output image. */
itkSetMacro( OutputRegion, OutputImageRegionType );
/** Get the size of the output image. */
itkGetConstReferenceMacro( OutputRegion, OutputImageRegionType );
/** Set the output image spacing. */
itkSetMacro(OutputSpacing, SpacingType);
virtual void SetOutputSpacing( const double* values);
/** Get the output image spacing. */
itkGetConstReferenceMacro( OutputSpacing, SpacingType );
/** Set the output image origin. */
itkSetMacro(OutputOrigin, OriginPointType);
virtual void SetOutputOrigin( const double* values);
/** Get the output image origin. */
itkGetConstReferenceMacro( OutputOrigin, OriginPointType );
/** Set the list of source landmarks */
itkSetConstObjectMacro( SourceLandmarks, LandmarkContainer );
itkSetConstObjectMacro( TargetLandmarks, LandmarkContainer );
/** DeformationFieldSource produces an image which is a different size
* than its input. As such, it needs to provide an implementation
* for GenerateOutputInformation() in order to inform the pipeline
* execution model. The original documentation of this method is
* below. \sa ProcessObject::GenerateOutputInformaton() */
virtual void GenerateOutputInformation();
/** Method Compute the Modified Time based on changed to the components. */
unsigned long GetMTime( void ) const;
protected:
DeformationFieldSource();
~DeformationFieldSource() {};
void PrintSelf(std::ostream& os, Indent indent) const;
/**
* GenerateData() computes the internal KernelBase spline and resamples
* the deformation field.
*/
void GenerateData();
/** Subsample the input deformation field and generate the
* landmarks for the kernel base spline
*/
void PrepareKernelBaseSpline();
private:
DeformationFieldSource(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
KernelTransformPointerType m_KernelTransform; // Coordinate transform to use
OutputImageRegionType m_OutputRegion; // Region of the output image
SpacingType m_OutputSpacing; // output image spacing
OriginPointType m_OutputOrigin; // output image origin
LandmarkContainerPointer m_SourceLandmarks; // List of source landmarks
LandmarkContainerPointer m_TargetLandmarks; // List of target landmarks
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkDeformationFieldSource.txx"
#endif
#endif
|
