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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkHardConnectedComponentImageFilter.h,v $
Language: C++
Date: $Date: 2006-03-28 19:59:05 $
Version: $Revision: 1.6 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm 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 __itkHardConnectedComponentImageFilter_h
#define __itkHardConnectedComponentImageFilter_h
#include "itkImageToImageFilter.h"
#include "itkImage.h"
namespace itk
{/**
* 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 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
*
*
*
*/
template <class TInputImage, class TOutputImage>
class ITK_EXPORT HardConnectedComponentImageFilter :
public ImageToImageFilter< TInputImage, TOutputImage >
{
public:
/**
* Standard "Self" & Superclass typedef.
*/
typedef HardConnectedComponentImageFilter Self;
typedef ImageToImageFilter< TInputImage, TOutputImage > Superclass;
/**
* 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;
/**
* Smart pointer typedef support
*/
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/**
* 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() {}
HardConnectedComponentImageFilter(const Self&) {}
/**
* Standard pipeline method.
*/
void GenerateData();
void PrintSelf(std::ostream& os, Indent indent) const
{ Superclass::PrintSelf(os,indent); }
private:
ListType m_Seeds;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkHardConnectedComponentImageFilter.txx"
#endif
#endif
|
