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
|
/*=========================================================================
Program: ORFEO Toolbox
Language: C++
Date: $Date$
Version: $Revision$
Copyright (c) Centre National d'Etudes Spatiales. All rights reserved.
See OTBCopyright.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 otbMeanShiftSegmentationFilter_h
#define otbMeanShiftSegmentationFilter_h
#include "itkMacro.h"
#include "otbMeanShiftSmoothingImageFilter.h"
#include "otbLabelImageRegionMergingFilter.h"
#include "otbLabelImageRegionPruningFilter.h"
#include "itkRelabelComponentImageFilter.h"
#include "itkConnectedComponentFunctorImageFilter.h"
namespace otb {
namespace Functor
{
template<class TInput>
class ITK_EXPORT ConnectedLabelFunctor
{
public:
typedef ConnectedLabelFunctor Self;
std::string GetNameOfClass()
{
return "ConnectedLabelFunctor";
}
inline bool operator()( TInput &p1, TInput &p2)
{
//return static_cast<bool> (0);
return static_cast<bool> (p1==p2);
}
ConnectedLabelFunctor()
{
}
~ConnectedLabelFunctor()
{
}
private:
ConnectedLabelFunctor(const Self &); //purposely not implemented
void operator =(const Self &); //purposely not implemented
};
} // end of Functor namespace
/** \class MeanShiftSegmentationFilter
*
* Performs segmentation of an image by chaining a mean shift filter and region
* merging filter.
* This filter relies on MeanShift Labeled image created from clustered filtered output of meanshift filter.
* Mode search algorithm is not performed on tiles boundaries. Thus output depends on Thread numbers.
*
* \ingroup OTBMeanShift
*/
template <class TInputImage, class TOutputLabelImage, class TOutputClusteredImage = TInputImage, class TKernel = Meanshift::KernelUniform>
class ITK_EXPORT MeanShiftSegmentationFilter : public itk::ImageToImageFilter<TInputImage, TOutputLabelImage>
{
public:
/** Standard Self typedef */
typedef MeanShiftSegmentationFilter Self;
typedef itk::ImageToImageFilter<TInputImage, TOutputLabelImage> Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Some convenient typedefs. */
typedef TInputImage InputSpectralImageType;
typedef TOutputLabelImage OutputLabelImageType;
typedef TOutputClusteredImage OutputClusteredImageType;
typedef TKernel KernelType;
typedef double RealType;
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Runtime information support. */
itkTypeMacro(MeanShiftSegmentationFilter, ImageToImageFilter);
/** ImageDimension constants */
itkStaticConstMacro(ImageDimension, unsigned int, TInputImage::ImageDimension);
// Mean shift filter
typedef OutputClusteredImageType MeanShiftFilteredImageType;
typedef MeanShiftSmoothingImageFilter<InputSpectralImageType, MeanShiftFilteredImageType, KernelType> MeanShiftFilterType;
typedef typename MeanShiftFilterType::Pointer MeanShiftFilterPointerType;
// Region merging filter
typedef typename MeanShiftFilterType::OutputLabelImageType InputLabelImageType;
typedef typename MeanShiftFilterType::LabelType InputLabelPixelType;
typedef LabelImageRegionMergingFilter<InputLabelImageType, MeanShiftFilteredImageType,
OutputLabelImageType, OutputClusteredImageType> RegionMergingFilterType;
typedef typename RegionMergingFilterType::Pointer RegionMergingFilterPointerType;
typedef LabelImageRegionPruningFilter<OutputLabelImageType,OutputClusteredImageType,
OutputLabelImageType, OutputClusteredImageType> RegionPruningFilterType;
typedef typename RegionPruningFilterType::Pointer RegionPruningFilterPointerType;
typedef Functor::ConnectedLabelFunctor<InputLabelPixelType> LabelFunctorType;
typedef itk::ConnectedComponentFunctorImageFilter<InputLabelImageType,InputLabelImageType,LabelFunctorType> RelabelComponentFilterType;
typedef typename RelabelComponentFilterType::Pointer RelabelComponentFilterPointerType;
/** Sets the spatial bandwidth (or radius in the case of a uniform kernel)
* of the neighborhood for each pixel
*/
otbSetObjectMemberMacro(MeanShiftFilter,SpatialBandwidth,RealType);
otbGetObjectMemberConstReferenceMacro(MeanShiftFilter,SpatialBandwidth,RealType);
/** Sets the spectral bandwidth (or radius for a uniform kernel) for pixels
* to be included in the same mode
*/
otbSetObjectMemberMacro(MeanShiftFilter,RangeBandwidth,RealType);
otbGetObjectMemberConstReferenceMacro(MeanShiftFilter,RangeBandwidth,RealType);
/** Set the maximum number of iterations */
otbSetObjectMemberMacro(MeanShiftFilter,MaxIterationNumber,unsigned int);
otbGetObjectMemberConstReferenceMacro(MeanShiftFilter,MaxIterationNumber,unsigned int);
/** Set the convergence threshold */
otbSetObjectMemberMacro(MeanShiftFilter,Threshold,RealType);
otbGetObjectMemberConstReferenceMacro(MeanShiftFilter,Threshold,RealType);
/** Sets the minimum region size. (after merging step clustered regions, which size is under this
* threshold will be fused with adjacent region with smallest spectral distance).
*/
otbSetObjectMemberMacro(RegionPruningFilter,MinRegionSize,RealType);
otbGetObjectMemberMacro(RegionPruningFilter,MinRegionSize,RealType);
/** Returns the const image of region labels */
const OutputLabelImageType * GetLabelOutput() const;
/** Returns the image of region labels */
OutputLabelImageType * GetLabelOutput();
/** Returns the const clustered output image, with one spectral value per region */
const OutputClusteredImageType * GetClusteredOutput() const;
/** Returns the clustered output image, with one spectral value per region */
OutputClusteredImageType * GetClusteredOutput();
protected:
MeanShiftSegmentationFilter();
~MeanShiftSegmentationFilter() ITK_OVERRIDE;
// virtual void GenerateOutputInformation(void);
void GenerateData() ITK_OVERRIDE;
private:
MeanShiftFilterPointerType m_MeanShiftFilter;
RegionMergingFilterPointerType m_RegionMergingFilter;
RegionPruningFilterPointerType m_RegionPruningFilter;
RelabelComponentFilterPointerType m_RelabelFilter;
};
} // end namespace otb
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbMeanShiftSegmentationFilter.txx"
#endif
#endif
|
