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
|
/*
* Copyright (C) 2005-2017 Centre National d'Etudes Spatiales (CNES)
*
* This file is part of Orfeo Toolbox
*
* https://www.orfeo-toolbox.org/
*
* 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
*
* 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 otbPixelSuppressionByDirectionImageFilter_h
#define otbPixelSuppressionByDirectionImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkImage.h"
#include "itkNumericTraits.h"
namespace otb
{
/** \class PixelSuppressionByDirectionImageFilter
* \brief Application of a filter of suppression of isolated pixels,
* not belonging to a line, strating from the directions of pixels.
*
* This class implements an image filter which detects isolated pixels
* that have little chance of belonging to a raod and performs a pixel
* suppression. For each pixel kept with direction \f$ \theta_{i} \f$, we look
* for other pixels with a direction close to \f$ \theta_i \f$ in an angular
* beam around it. If none is found, the pixel is suppressed.
*
* This filter is the first step to generate an image of segments
* primitives. It starts from the two output images of the line detector
* image filters which are the image of intensity of detection and
* the image of direction.
*
*
* \ingroup OTBEdge
*/
template <class TInputImage, class TOutputImage>
class ITK_EXPORT PixelSuppressionByDirectionImageFilter : public itk::ImageToImageFilter<TInputImage, TOutputImage>
{
public:
/** Extract input and output image dimensions */
itkStaticConstMacro(InputImageDimension,
unsigned int,
TInputImage::ImageDimension);
itkStaticConstMacro(OutputImageDimension,
unsigned int,
TOutputImage::ImageDimension);
typedef TInputImage InputImageType;
typedef TOutputImage OutputImageType;
/** typedef for the classes standards. */
typedef PixelSuppressionByDirectionImageFilter Self;
typedef itk::ImageToImageFilter<InputImageType, OutputImageType> Superclass;
typedef itk::SmartPointer<Self> Pointer;
typedef itk::SmartPointer<const Self> ConstPointer;
/** Method for management of the "object factory". */
itkNewMacro(Self);
/** Return the nale of the class. */
itkTypeMacro(PixelSuppressionByDirectionImageFilter, ImageToImageFilter);
/** Definition of the input and output images */
typedef typename InputImageType::PixelType InputPixelType;
typedef typename OutputImageType::PixelType OutputPixelType;
typedef typename itk::NumericTraits<InputPixelType>::RealType InputRealType;
typedef typename InputImageType::RegionType InputImageRegionType;
typedef typename OutputImageType::RegionType OutputImageRegionType;
typedef typename InputImageType::SizeType SizeType;
/** Set the radius of one zone. */
itkSetMacro(Radius, SizeType);
/** Get the radius of one zone. */
itkGetConstReferenceMacro(Radius, SizeType);
/** Set the angular beam. */
itkSetMacro(AngularBeam, InputRealType);
/** Get the angular beam. */
itkGetConstReferenceMacro(AngularBeam, InputRealType);
/** Set/Get the image input of this process object. */
void SetInputImage(const InputImageType *image);
const InputImageType * GetInputImage(void);
void SetInputImageDirection(const InputImageType *image);
const InputImageType * GetInputImageDirection(void);
void GenerateInputRequestedRegion()
throw(itk::InvalidRequestedRegionError) override;
protected:
PixelSuppressionByDirectionImageFilter();
~PixelSuppressionByDirectionImageFilter() override {}
void PrintSelf(std::ostream& os, itk::Indent indent) const override;
void ThreadedGenerateData(const OutputImageRegionType& outputRegionForThread,
itk::ThreadIdType threadId) override;
private:
PixelSuppressionByDirectionImageFilter(const Self &); //purposely not implemented
void operator =(const Self&); //purposely not implemented
// Radius of the region
SizeType m_Radius;
// Angular Accuracy on the direction of the central pixel
InputRealType m_AngularBeam;
};
} // end namespace otb
#ifndef OTB_MANUAL_INSTANTIATION
#include "otbPixelSuppressionByDirectionImageFilter.txx"
#endif
#endif
|
