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
|
/*=========================================================================
*
* Copyright Insight Software Consortium
*
* 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.
*
*=========================================================================*/
#ifndef itkPermuteAxesImageFilter_h
#define itkPermuteAxesImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkFixedArray.h"
namespace itk
{
/** \class PermuteAxesImageFilter
* \brief Permutes the image axes according to a user specified order.
*
* PermuateAxesImageFilter permutes the image axes according to a
* user specified order. The permutation order is set via method
* SetOrder( order ) where the input is an array of ImageDimension
* number of unsigned int. The elements of the array must be a rearrangment
* of the numbers from 0 to ImageDimension - 1.
*
* The i-th axis of the output image corresponds with the order[i]-th
* axis of the input image.
*
* The output meta image information (LargestPossibleRegion, spacing, origin)
* is computed by permuting the corresponding input meta information.
*
* \ingroup GeometricTransform
* \ingroup MultiThreaded
* \ingroup Streamed
* \ingroup ITKImageGrid
*
* \wiki
* \wikiexample{ImageProcessing/PermuteAxesImageFilter,Switch the axes of an image}
* \endwiki
*/
template< typename TImage >
class ITK_TEMPLATE_EXPORT PermuteAxesImageFilter:
public ImageToImageFilter< TImage, TImage >
{
public:
/** Standard class typedefs. */
typedef PermuteAxesImageFilter Self;
typedef ImageToImageFilter< TImage, TImage > Superclass;
typedef SmartPointer< Self > Pointer;
typedef SmartPointer< const Self > ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(PermuteAxesImageFilter, ImageToImageFilter);
/** ImageDimension enumeration */
itkStaticConstMacro(ImageDimension, unsigned int,
TImage::ImageDimension);
/** Inherited types */
typedef typename Superclass::InputImagePointer InputImagePointer;
typedef typename Superclass::OutputImagePointer OutputImagePointer;
typedef typename Superclass::OutputImageRegionType OutputImageRegionType;
/** PermuteOrderArray type. */
typedef FixedArray< unsigned int, itkGetStaticConstMacro(ImageDimension) > PermuteOrderArrayType;
/** Set the permutation order. The elements of order must be
* a rearrangement of the numbers from 0 to ImageDimension - 1. */
void SetOrder(const PermuteOrderArrayType & order);
/** Get the permutation order. */
itkGetConstReferenceMacro(Order, PermuteOrderArrayType);
/** Get the inverse permutation order. */
itkGetConstReferenceMacro(InverseOrder, PermuteOrderArrayType);
protected:
/** PermuteAxesImageFilter produces an image which is a different
* resolution and with a different pixel spacing than its input
* image. As such, PermuteAxesImageFilter 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() ITK_OVERRIDE;
/** PermuteAxesImageFilter needs different input requested region than the output
* requested region. As such, PermuteAxesImageFilter needs to provide an
* implementation for GenerateInputRequestedRegion() in order to inform the
* pipeline execution model.
* \sa ProcessObject::GenerateInputRequestedRegion() */
virtual void GenerateInputRequestedRegion() ITK_OVERRIDE;
PermuteAxesImageFilter();
~PermuteAxesImageFilter() ITK_OVERRIDE {}
void PrintSelf(std::ostream & os, Indent indent) const ITK_OVERRIDE;
/** PermuteAxesImageFilter can be implemented as a multithreaded filter.
* Therefore, this implementation provides a ThreadedGenerateData() routine
* which is called for each processing thread. The output image data is
* allocated automatically by the superclass prior to calling
* ThreadedGenerateData(). ThreadedGenerateData can only write to the
* portion of the output image specified by the parameter
* "outputRegionForThread"
*
* \sa ImageToImageFilter::ThreadedGenerateData(),
* ImageToImageFilter::GenerateData() */
void ThreadedGenerateData(const OutputImageRegionType & outputRegionForThread,
ThreadIdType threadId) ITK_OVERRIDE;
private:
ITK_DISALLOW_COPY_AND_ASSIGN(PermuteAxesImageFilter);
PermuteOrderArrayType m_Order;
PermuteOrderArrayType m_InverseOrder;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkPermuteAxesImageFilter.hxx"
#endif
#endif
|
