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
|
/*=========================================================================
*
* 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 itkAreaOpeningImageFilter_h
#define itkAreaOpeningImageFilter_h
#include "itkAttributeMorphologyBaseImageFilter.h"
#include <functional>
namespace itk
{
/**
* \class AreaOpeningImageFilter
* \brief Morphological opening by attributes
*
* This is the base class for morphology attribute
* operations. Attribute openings remove blobs according to criteria
* such as area. When applied to grayscale images it has the effect of
* trimming peaks based on area while leaving the rest of the image
* unchanged. It is possible to use attributes besides area, but no
* others are implemented yet. This filter uses some dodgy coding
* practices - most notably copying the image data to a linear buffer
* to allow direct implementation of the published algorithm. It
* should therefore be quite a good candidate to carry out tests of
* itk iterator performance with randomish access patterns.
*
* This filter is implemented using the method of Wilkinson, "A
* comparison of algorithms for Connected set openings and Closings",
* A. Meijster and M. H. Wilkinson, PAMI, vol 24, no. 4, April 2002.
* Attempts at implementing the method from ISMM 2000 are also
* included, but operation appears incorrect. Check the ifdefs if you
* are interested.
*
* This code was contributed in the Insight Journal paper
*
* "Grayscale morphological attribute operations"
* by Beare R.
* https://hdl.handle.net/1926/1316
* http://www.insight-journal.org/browse/publication/203
*
*
* \author Richard Beare. Department of Medicine, Monash University, Melbourne, Australia.
*
* \ingroup ITKReview
*/
template< typename TInputImage, typename TOutputImage, typename TAttribute = typename TInputImage::SpacingType::ValueType >
class AreaOpeningImageFilter:
public AttributeMorphologyBaseImageFilter< TInputImage, TOutputImage, TAttribute,
std::greater< typename TInputImage::PixelType > >
{
public:
typedef AreaOpeningImageFilter Self;
typedef AttributeMorphologyBaseImageFilter< TInputImage, TOutputImage, TAttribute,
std::greater< typename TInputImage::PixelType > >
Superclass;
typedef SmartPointer< Self > Pointer;
typedef SmartPointer< const Self > ConstPointer;
/**
* Extract some information from the image types. Dimensionality
* of the two images is assumed to be the same.
*/
typedef typename TOutputImage::PixelType OutputPixelType;
typedef typename TOutputImage::InternalPixelType OutputInternalPixelType;
typedef typename TInputImage::PixelType InputPixelType;
typedef typename TInputImage::InternalPixelType InputInternalPixelType;
typedef typename TInputImage::IndexType IndexType;
typedef typename TInputImage::OffsetType OffsetType;
typedef typename TInputImage::SizeType SizeType;
typedef TAttribute AttributeType;
itkStaticConstMacro(ImageDimension, unsigned int,
TOutputImage::ImageDimension);
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Runtime information support. */
itkTypeMacro(AreaOpeningImageFilter,
AttributeMorphologyBaseImageFilter);
/**
* Set/Get whether the image spacing is used or not - defaults to true.
*/
itkSetMacro(UseImageSpacing, bool);
itkGetConstReferenceMacro(UseImageSpacing, bool);
itkBooleanMacro(UseImageSpacing);
protected:
AreaOpeningImageFilter()
{
m_UseImageSpacing = true;
}
virtual ~AreaOpeningImageFilter() {}
void GenerateData() ITK_OVERRIDE
{
this->m_AttributeValuePerPixel = 1;
if ( m_UseImageSpacing )
{
// compute pixel size
double psize = 1.0;
for ( unsigned i = 0; i < ImageDimension; i++ )
{
psize *= this->GetInput()->GetSpacing()[i];
}
this->m_AttributeValuePerPixel = static_cast< AttributeType >( psize );
// std::cout << "m_AttributeValuePerPixel: " <<
// this->m_AttributeValuePerPixel << std::endl;
// and call superclass implementation of GenerateData()
}
Superclass::GenerateData();
}
void PrintSelf(std::ostream & os, Indent indent) const ITK_OVERRIDE
{
Superclass::PrintSelf(os, indent);
os << indent << "UseImageSpacing: " << m_UseImageSpacing << std::endl;
}
private:
AreaOpeningImageFilter(const Self &) ITK_DELETE_FUNCTION;
void operator=(const Self &) ITK_DELETE_FUNCTION;
bool m_UseImageSpacing;
};
} // namespace itk
#endif
|
