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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkShapeLabelMapFilter.h
Language: C++
Date: $Date$
Version: $Revision$
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 __itkShapeLabelMapFilter_h
#define __itkShapeLabelMapFilter_h
#include "itkLabelPerimeterEstimationCalculator.h"
#include "itkInPlaceLabelMapFilter.h"
namespace itk {
/** \class ShapeLabelMapFilter
* \brief The valuator class for the ShapeLabelObject
*
* ShapeLabelMapFilter can be used to set the attributes values of the
* ShapeLabelObject in a LabelMap.
*
* ShapeLabelMapFilter takes an optional parameter, used only to
* optimize the computation time and the memory usage when the
* perimeter or the feret diameter is used: the exact copy of the
* input LabelMap is stored in an Image.
* It can be set with SetLabelImage(). It is cleared at the end of the
* computation, so must be reset before running Update() again. It is
* not part of the pipeline management design, to let the subclasses
* of ShapeLabelMapFilter use the pipeline design to specify truly
* required inputs.
*
* \author Gaetan Lehmann. Biologie du Developpement et de la Reproduction, INRA de Jouy-en-Josas, France.
*
* This implementation was taken from the Insight Journal paper:
* http://hdl.handle.net/1926/584 or
* http://www.insight-journal.org/browse/publication/176
*
* \ingroup ImageEnhancement MathematicalMorphologyImageFilters
*/
template<class TImage, class TLabelImage=Image< ITK_TYPENAME TImage::PixelType, ::itk::GetImageDimension<TImage>::ImageDimension> >
class ITK_EXPORT ShapeLabelMapFilter :
public
InPlaceLabelMapFilter<TImage>
{
public:
/** Standard class typedefs. */
typedef ShapeLabelMapFilter Self;
typedef InPlaceLabelMapFilter<TImage> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Some convenient typedefs. */
typedef TImage ImageType;
typedef typename ImageType::Pointer ImagePointer;
typedef typename ImageType::ConstPointer ImageConstPointer;
typedef typename ImageType::PixelType PixelType;
typedef typename ImageType::IndexType IndexType;
typedef typename ImageType::SizeType SizeType;
typedef typename ImageType::LabelObjectType LabelObjectType;
typedef typename LabelObjectType::MatrixType MatrixType;
typedef typename LabelObjectType::VectorType VectorType;
typedef TLabelImage LabelImageType;
typedef typename LabelImageType::Pointer LabelImagePointer;
typedef typename LabelImageType::ConstPointer LabelImageConstPointer;
typedef typename LabelImageType::PixelType LabelPixelType;
/** ImageDimension constants */
itkStaticConstMacro(ImageDimension, unsigned int, TImage::ImageDimension );
typedef LabelPerimeterEstimationCalculator< LabelImageType > PerimeterCalculatorType;
/** Standard New method. */
itkNewMacro(Self);
/** Runtime information support. */
itkTypeMacro(ShapeLabelMapFilter, InPlaceLabelMapFilter );
#ifdef ITK_USE_CONCEPT_CHECKING
/** Begin concept checking */
/* itkConceptMacro(InputEqualityComparableCheck,
(Concept::EqualityComparable<InputImagePixelType>));
itkConceptMacro(IntConvertibleToInputCheck,
(Concept::Convertible<int, InputImagePixelType>));
itkConceptMacro(InputOStreamWritableCheck,
(Concept::OStreamWritable<InputImagePixelType>));*/
/** End concept checking */
#endif
/**
* Set/Get whether the maximum Feret diameter should be computed or not.
* Default value is false because of the high computation time required.
*/
itkSetMacro( ComputeFeretDiameter, bool );
itkGetConstReferenceMacro( ComputeFeretDiameter, bool );
itkBooleanMacro( ComputeFeretDiameter );
/**
* Set/Get whether the perimeter should be computed or not.
* Default value is false because of the high computation time required.
*/
itkSetMacro( ComputePerimeter, bool );
itkGetConstReferenceMacro( ComputePerimeter, bool );
itkBooleanMacro( ComputePerimeter );
/** Set the label image */
void SetLabelImage( const TLabelImage *input )
{
m_LabelImage = input;
}
protected:
ShapeLabelMapFilter();
~ShapeLabelMapFilter() {};
virtual void ThreadedProcessLabelObject( LabelObjectType * labelObject );
virtual void BeforeThreadedGenerateData();
virtual void AfterThreadedGenerateData();
void PrintSelf(std::ostream& os, Indent indent) const;
private:
ShapeLabelMapFilter(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
bool m_ComputeFeretDiameter;
bool m_ComputePerimeter;
LabelImageConstPointer m_LabelImage;
typename PerimeterCalculatorType::Pointer m_PerimeterCalculator;
void ComputeFeretDiameter( LabelObjectType * labelObject );
/** utilities */
static long Factorial( long n );
static long DoubleFactorial( long n );
static double GammaN2p1( long n );
static double HyperSphereVolume( double radius );
static double HyperSpherePerimeter( double radius );
static double HyperSphereRadiusFromVolume( double volume );
}; // end of class
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkShapeLabelMapFilter.txx"
#endif
#endif
|
