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
|
/*=========================================================================
*
* 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 itkHardConnectedComponentImageFilter_h
#define itkHardConnectedComponentImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkImage.h"
namespace itk
{
/**
* \class HardConnectedComponentImageFilter
* The purpose of this program is to produce the connected components
* for any input binary image of dimensionality n.
*
* The program does a forward pass line by line through the entire image.
* Each cell in the foreground is assigned the same label value as cells
* in its neighborhood. If there is no label among the cells in its neighborhood,
* a new label value is assigned to the cell. This means that this cell belongs
* to a different connected component. We set up an equivalence table for each
* label to indicate the equivalence of the labels stored in the table. After
* the forward pass goes through the entire image, we merge the different
* connected components corresponding to the equivalence labels in the table.
* We implement this strategy in the function GenerateData().
*
* There are two options in the program.
* 1. Take an nD binary image as input, and produce an nD gray image, where intensity indicates label assigned to a connected component.
* 2. Take an nD binary image and a set of seed points as input, and output an nD binary image containing the cells connected to the seeds.
* For option 2, users need to assign the member variable std::list<IndexType> m_Seeds before calling function GenerateData().
* \sa ImageToImageFilter
* \ingroup ITKConnectedComponents
*/
template< typename TInputImage, typename TOutputImage >
class ITK_TEMPLATE_EXPORT HardConnectedComponentImageFilter:
public ImageToImageFilter< TInputImage, TOutputImage >
{
public:
/**
* Standard class typedef's
*/
typedef HardConnectedComponentImageFilter Self;
typedef ImageToImageFilter< TInputImage, TOutputImage > 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;
itkStaticConstMacro(InputImageDimension, unsigned int,
TInputImage::ImageDimension);
itkStaticConstMacro(ImageDimension, unsigned int,
TOutputImage::ImageDimension);
/**
* Image typedef support
*/
typedef TInputImage InputImageType;
typedef TOutputImage OutputImageType;
typedef typename TInputImage::IndexType IndexType;
typedef typename TInputImage::SizeType SizeType;
typedef typename TOutputImage::RegionType RegionType;
typedef std::list< IndexType > ListType;
/**
* Run-time type information (and related methods)
*/
itkTypeMacro(HardConnectedComponentImageFilter, ImageToImageFilter);
/**
* Method for creation through the object factory.
*/
itkNewMacro(Self);
/** Setting the seed points for specified object. */
void SetObjectSeed(const IndexType & seed)
{ m_Seeds.push_front(seed); }
#ifdef ITK_USE_CONCEPT_CHECKING
// Begin concept checking
itkConceptMacro( SameDimensionCheck,
( Concept::SameDimension< InputImageDimension, ImageDimension > ) );
itkConceptMacro( IntConvertibleToOutputCheck,
( Concept::Convertible< int, OutputPixelType > ) );
itkConceptMacro( UnsignedShortConvertibleToOutputCheck,
( Concept::Convertible< unsigned short, OutputPixelType > ) );
itkConceptMacro( OutputEqualityComparableCheck,
( Concept::EqualityComparable< OutputPixelType > ) );
itkConceptMacro( UnsignedCharConvertibleToOutputCheck,
( Concept::Convertible< unsigned char, OutputPixelType > ) );
itkConceptMacro( OutputIncrementDecrementOperatorsCheck,
( Concept::IncrementDecrementOperators< OutputPixelType > ) );
// End concept checking
#endif
protected:
HardConnectedComponentImageFilter() {}
virtual ~HardConnectedComponentImageFilter() ITK_OVERRIDE {}
/**
* Standard pipeline method.
*/
void GenerateData() ITK_OVERRIDE;
void PrintSelf(std::ostream & os, Indent indent) const ITK_OVERRIDE
{ Superclass::PrintSelf(os, indent); }
private:
ITK_DISALLOW_COPY_AND_ASSIGN(HardConnectedComponentImageFilter);
ListType m_Seeds;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkHardConnectedComponentImageFilter.hxx"
#endif
#endif
|
