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
|
/*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.txt for details.
Copyright (c) Institut Mines-Telecom. All rights reserved.
See IMTCopyright.txt 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 otbSubsampleImageFilter_h
#define otbSubsampleImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkImage.h"
//Just to get the forward/inverse definitions
#include "otbWaveletOperatorBase.h"
namespace otb {
/** \class SubsampleImageFilter
* \brief Performs a down sampling of an image
*
* This class performs a simple resampling image filtering.
*
* The third template is a Wavelet::FORWARD/Wavelet::INVERSE Wavelet::WaveletDirection
* enumeration type. With Wavelet::FORWARD, the output image is of the same size of the one of
* the input image. But the output image contains only co-located subsample, the rest being 0.
* When Wavelet::INVERSE is used, the input image is considered as the sub-sampled. Then, the output
* image is larger. Initial pixel values are preserved but the output image is interleaved with
* 0.
*
* \sa ResampleImageFilter
* \sa SubsampleImageRegionConstIterator
* \sa DecimateImageFilter
*
* \ingroup OTBWavelet
*/
template <class TInputImage, class TOutputImage,
Wavelet::WaveletDirection TDirectionOfTransformation>
class ITK_EXPORT SubsampleImageFilter :
public itk::ImageToImageFilter<TInputImage, TOutputImage>
{
public:
/** Standard class typedefs. */
typedef SubsampleImageFilter Self;
typedef itk::ImageToImageFilter<TInputImage, TOutputImage> Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Run-time type information (and related methods). */
itkTypeMacro(SubsampleImageFilter, ImageToImageFilter);
/** Extract dimension from input and output image. */
itkStaticConstMacro(InputImageDimension, unsigned int,
TInputImage::ImageDimension);
itkStaticConstMacro(OutputImageDimension, unsigned int,
TOutputImage::ImageDimension);
/** Direction definition */
typedef Wavelet::WaveletDirection DirectionOfTransformationEnumType;
itkStaticConstMacro(DirectionOfTransformation, DirectionOfTransformationEnumType, TDirectionOfTransformation);
/** Image typedef support. */
typedef TInputImage InputImageType;
typedef typename InputImageType::RegionType InputImageRegionType;
typedef typename InputImageType::IndexType InputImageIndexType;
typedef typename InputImageIndexType::IndexValueType InputImageIndexValueType;
typedef TOutputImage OutputImageType;
typedef typename OutputImageType::Pointer OutputImagePointerType;
typedef typename OutputImageType::RegionType OutputImageRegionType;
typedef typename OutputImageType::SizeType OutputImageSizeType;
typedef typename OutputImageType::IndexType OutputImageIndexType;
typedef typename OutputImageType::PixelType OutputPixelType;
/** Set/Get the SubsampleFactor */
itkGetMacro(SubsampleFactor, const InputImageIndexType &);
itkSetMacro(SubsampleFactor, InputImageIndexType &);
void SetSubSampleFactor(InputImageIndexValueType factor)
{
InputImageIndexType indexFactor;
indexFactor.Fill(factor);
SetSubSampleFactor(indexFactor);
}
protected:
SubsampleImageFilter ()
{
m_SubsampleFactor.Fill(1);
//this->SetNumberOfThreads(1);
}
~SubsampleImageFilter() ITK_OVERRIDE {}
/** Internal test function to check if there is any direction to subsample */
bool IsSubsampleFactorOne() const;
/** Since input and output image are very likely to be of different size.
* Region estimation functions has to be reimplemented
*/
void CallCopyOutputRegionToInputRegion
(InputImageRegionType& destRegion, const OutputImageRegionType& srcRegion) ITK_OVERRIDE;
void CallCopyInputRegionToOutputRegion
(OutputImageRegionType& destRegion, const InputImageRegionType& srcRegion) ITK_OVERRIDE;
/** Output image region size is not of the same dimension as the input.
* That is why GenerateOutputInformation has to be redefined.
*/
void GenerateOutputInformation() ITK_OVERRIDE;
/** Set output image to 0 before processing */
void BeforeThreadedGenerateData() ITK_OVERRIDE;
/** Allows multithreading */
void ThreadedGenerateData
(const OutputImageRegionType& outputRegionForThread, itk::ThreadIdType threadId) ITK_OVERRIDE;
void PrintSelf(std::ostream& os, itk::Indent indent) const ITK_OVERRIDE;
private:
SubsampleImageFilter (const Self &); // purposely not implemented
void operator =(const Self&); // purposely not implemented
InputImageIndexType m_SubsampleFactor;
}; // end of class
} // end of namespace otb
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbSubsampleImageFilter.txx"
#endif
#endif
|
