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
|
/*
//
// Copyright 2010-2011, 2013 SRI International
//
// This file is part of the Computational Morphometry Toolkit.
//
// http://www.nitrc.org/projects/cmtk/
//
// The Computational Morphometry Toolkit is free software: you can
// redistribute it and/or modify it under the terms of the GNU General Public
// License as published by the Free Software Foundation, either version 3 of
// the License, or (at your option) any later version.
//
// The Computational Morphometry Toolkit is distributed in the hope that it
// will be useful, but WITHOUT ANY WARRANTY; without even the implied
// warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License along
// with the Computational Morphometry Toolkit. If not, see
// <http://www.gnu.org/licenses/>.
//
// $Revision: 5436 $
//
// $LastChangedDate: 2018-12-10 19:01:20 -0800 (Mon, 10 Dec 2018) $
//
// $LastChangedBy: torstenrohlfing $
//
*/
#ifndef __cmtkSimpleLevelset_h_included_
#define __cmtkSimpleLevelset_h_included_
#include <cmtkconfig.h>
#include <Base/cmtkUniformVolume.h>
#include <Base/cmtkUnits.h>
#include <System/cmtkException.h>
namespace
cmtk
{
/** \addtogroup Segmentation */
//@{
/** Class for computing a simple two-phase levelset evolution.
*/
class SimpleLevelset
{
public:
/// This class.
typedef SimpleLevelset Self;
/// Exception that is thrown for degenerate levelsets, i.e., all foreground or all background.
class DegenerateLevelsetException : public Exception {};
/// Constructor.
SimpleLevelset( UniformVolume::SmartConstPtr& volume ) : m_Volume( volume ), m_ScaleInitialSphere( 1.0 ), m_TimeDelta( 1.0 ), m_LevelsetThreshold( 0.0 ) {}
/// Set initial sphere scale factor.
void SetScaleInitialSphere( const Types::Coordinate scale )
{
this->m_ScaleInitialSphere = scale;
}
/// Set filter sigma parameter.
void SetFilterSigma( const Units::GaussianSigma filterSigma )
{
this->m_FilterSigma = filterSigma;
}
/// Set evolution time delta.
void SetTimeDelta( const Types::Coordinate timeDelta )
{
this->m_TimeDelta = timeDelta;
}
/// Set levelset threshold.
void SetLevelsetThreshold( const Types::Coordinate levelsetThreshold )
{
this->m_LevelsetThreshold = levelsetThreshold;
}
/// Initialize levelset with a centered sphere.
void InitializeCenteredSphere();
/// Levelset evolution.
virtual void Evolve( const int numberOfIterations /*!< Number of iterations */,
const bool forceIterations = false /*!< If this is set, evolution continues until maximum iteration count is reached, even when convergence is detected */ );
/** Return levelset, optionally converting to a binarized byte pixel representation.
*\warning If the levelset is retrieved with the "binarize" flag set, then the
* levelset stored in this object will remain binarized after the call and
* should be re-initialized before calling "Evolve" again.
*/
UniformVolume::SmartPtr& GetLevelset( const bool binarize = false, /*!< If set, levelset is binarized and converted to byte data */ const float threshold = 0.5 /*!< Threshold for optional binarization */ );
protected:
/// The volume to compute a levelset segmentation for.
UniformVolume::SmartConstPtr m_Volume;
/// The evolving levelset.
UniformVolume::SmartPtr m_Levelset;
/// Initial sphere scale factor.
Types::Coordinate m_ScaleInitialSphere;
/// Sigma parameter of the Gaussian filter kernel.
Units::GaussianSigma m_FilterSigma;
/// Delta time constant.
Types::Coordinate m_TimeDelta;
/// Levelset threshold.
Types::Coordinate m_LevelsetThreshold;
};
} // namespace cmtk
#endif // #ifndef __cmtkSimpleLevelset_h_included_
|
