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
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
|
/*=========================================================================
*
* Copyright UMC Utrecht and contributors
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0.txt
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
*
*=========================================================================*/
/** Parts of the code were taken from an ITK file.
* Original ITK copyright message, just for reference: */
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile$
Date: $Date: 2008-04-15 19:54:41 +0200 (Tue, 15 Apr 2008) $
Version: $Revision: 1573 $
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 itkMultiResolutionGaussianSmoothingPyramidImageFilter_h
#define itkMultiResolutionGaussianSmoothingPyramidImageFilter_h
#include "itkMultiResolutionPyramidImageFilter.h"
namespace itk
{
/** \class MultiResolutionGaussianSmoothingPyramidImageFilter
* \brief Framework for creating images in a multi-resolution
* pyramid.
*
* MultiResolutionGaussianSmoothingPyramidImageFilter creates
* an image pryamid according to a user defined multi-resolution schedule.
*
* This class inherits from the MultiResolutionPyramidImageFilter. It
* applies the same smoothing but does NOT do the downsampling.
*
* The multi-resolution schedule is still specified in terms for
* 'shrink factors' at each multi-resolution level for each dimension
* (although, actual shrinking is not performed).
*
* A user can either use the default schedules or specify
* each factor in the schedules directly.
*
* The schedule is stored as an unsigned int matrix.
* An element of the table can be access via the double bracket
* notation: table[resLevel][dimension]
*
* For example:
* 8 4 4
* 4 4 2
*
* is a schedule for two computation level. In the first (coarest)
* level the image is reduce by a factor of 8 in the column dimension,
* factor of 4 in the row dimension and factor of 4 in the slice dimension.
* In the second level, the image is reduce by a factor of 4 in the column
* dimension, 4 is the row dimension and 2 in the slice dimension.
*
* The method SetNumberOfLevels() set the number of
* computation levels in the pyramid. This method will
* allocate memory for the multi-resolution schedule table.
* This method generates defaults tables with the starting
* shrink factor for all dimension set to 2^(NumberOfLevel - 1).
* All factors are halved for all subsequent levels.
* For example if the number of levels was set to 4, the default table is:
*
* 8 8 8
* 4 4 4
* 2 2 2
* 1 1 1
*
* The user can get a copy of the schedule via GetSchedule()
* They may make alteration and reset it using SetSchedule().
*
* A user can create a default table by specifying the starting
* shrink factors via methods SetStartingShrinkFactors()
* The factors for subsequent level is generated by
* halving the factor or setting to one, depending on which is larger.
*
* For example, for 4 levels and starting factors of 8,8,4
* the default table would be:
*
* 8 8 4
* 4 4 2
* 2 2 1
* 1 1 1
*
* When this filter is updated, NumberOfLevels outputs are produced.
* The N'th output correspond to the N'th level of the pyramid.
*
* To generate each output image, Gaussian smoothing is first performed
* using a series of RecursiveGaussianImageFilter with standard deviation
* (shrink factor / 2)*imagespacing.
* The smoothed images are NOT downsampled, in contrast to the superclass's
* behaviour.
*
* This class is templated over the input image type and the output image
* type.
*
* This filter uses multithreaded filters to perform the smoothing.
*
* This filter supports streaming.
*
* \ingroup PyramidImageFilter Multithreaded Streamed
*/
template <class TInputImage, class TOutputImage>
class ITK_TEMPLATE_EXPORT MultiResolutionGaussianSmoothingPyramidImageFilter
: public MultiResolutionPyramidImageFilter<TInputImage, TOutputImage>
{
public:
ITK_DISALLOW_COPY_AND_MOVE(MultiResolutionGaussianSmoothingPyramidImageFilter);
/** Standard class typedefs. */
using Self = MultiResolutionGaussianSmoothingPyramidImageFilter;
using Superclass = MultiResolutionPyramidImageFilter<TInputImage, TOutputImage>;
using Pointer = SmartPointer<Self>;
using ConstPointer = SmartPointer<const Self>;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(MultiResolutionGaussianSmoothingPyramidImageFilter, MultiResolutionPyramidImageFilter);
/** ImageDimension enumeration. */
itkStaticConstMacro(ImageDimension, unsigned int, TInputImage::ImageDimension);
itkStaticConstMacro(OutputImageDimension, unsigned int, TOutputImage::ImageDimension);
/** Inherit types from Superclass. */
using typename Superclass::ScheduleType;
using typename Superclass::InputImageType;
using typename Superclass::OutputImageType;
using typename Superclass::InputImagePointer;
using typename Superclass::OutputImagePointer;
using typename Superclass::InputImageConstPointer;
/** Set a multi-resolution schedule. The input schedule must have only
* ImageDimension number of columns and NumberOfLevels number of rows. In
* contrast to the superclass, any schedule is allowed:
* - For each dimension, the shrink factor may be non-increasing with respect to
* subsequent levels.
* - shrink factors of 0 are allowed. This results in almost no smoothing.
* Because of lazy programming, the image is then smoothed with a gauss with sigma
* of 0.01*spacing...
*
* Note that the images are not actually shrunk by this class. They are
* only smoothed with the same standard deviation gaussian as used by
* the superclass.
*/
void
SetSchedule(const ScheduleType & schedule) override;
/** Set spacing etc. */
void
GenerateOutputInformation() override;
/** Given one output whose requested region has been set, this method sets
* the requested region for the remaining output images. The original
* documentation of this method is below. \sa
* ProcessObject::GenerateOutputRequestedRegion(); */
void
GenerateOutputRequestedRegion(DataObject * output) override;
/** MultiResolutionGaussianSmoothingPyramidImageFilter requires a larger input requested
* region than the output requested regions to accomdate the
* smoothing operations. As such, MultiResolutionGaussianSmoothingPyramidImageFilter needs
* to provide an implementation for GenerateInputRequestedRegion(). The
* original documentation of this method is below. \sa
* ProcessObject::GenerateInputRequestedRegion() */
void
GenerateInputRequestedRegion() override;
protected:
MultiResolutionGaussianSmoothingPyramidImageFilter() = default;
~MultiResolutionGaussianSmoothingPyramidImageFilter() override = default;
using Superclass::PrintSelf;
/** Generate the output data. */
void
GenerateData() override;
/** This filter by default generates the largest possible region,
* because it uses internally a filter that does this. */
void
EnlargeOutputRequestedRegion(DataObject * output) override;
};
} // namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
# include "itkMultiResolutionGaussianSmoothingPyramidImageFilter.hxx"
#endif
#endif
|
