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
|
/*
//
// Copyright 1997-2009 Torsten Rohlfing
//
// Copyright 2004-2012 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 __cmtkLogHistogram_h_included_
#define __cmtkLogHistogram_h_included_
#include <cmtkconfig.h>
#include <Base/cmtkHistogram.h>
#include <System/cmtkSmartPtr.h>
#include <math.h>
namespace
cmtk
{
/** \addtogroup Base */
//@{
/** Histogram of log intensities.
*/
template<class T>
class LogHistogram :
/// Inherit from non-log histogram..
public Histogram<T>
{
public:
/// This class.
typedef LogHistogram<T> Self;
/// Smart pointer.
typedef SmartPointer<Self> SmartPtr;
/// Parent class.
typedef Histogram<T> Superclass;
/** Constructor.
*/
LogHistogram ( const size_t numBins = 0 ) : Superclass( numBins ), m_LogNumBins( log( static_cast<double>( numBins ) ) ) {}
/** Destructor.
* All bin arrays and the precomputed data bin index arrays are
* de-allocated.
*/
virtual ~LogHistogram() {}
/// Resize and allocate histogram bins.
virtual void Resize( const size_t numberOfBins, const bool reset = true )
{
this->Superclass::Resize( numberOfBins, reset );
this->m_LogNumBins = log( static_cast<double>( numberOfBins ) );
}
/// Make an identical copy of this object.
typename Self::SmartPtr Clone () const
{
return typename Self::SmartPtr( this->CloneVirtual() );
}
/** Return bin corresponding to a certain value of the distribution.
*\param value A value from the distribution.
*\return The index of the bin corresponding to the given value.
*/
virtual size_t ValueToBin ( const Types::DataItem value ) const
{
return static_cast<size_t>( this->ValueToBinFractional( value ) );
}
/** Return fractional bin corresponding to a value of the distribution.
*\param value A value from the distribution.
*\return The index of the fractional bin index corresponding to the given
* value. This value is an integer if and only if the given value is
* identical to the lower bound of a bin.
*/
virtual Types::DataItem ValueToBinFractional ( const Types::DataItem value ) const
{
const Types::DataItem binIndex = this->Superclass::ValueToBinFractional( value );
return (this->GetNumberOfBins()-1) * std::max<Types::DataItem>( 0.0, std::min<Types::DataItem>( 1.0, log( 1+binIndex ) / this->m_LogNumBins ) );
}
/** Get value range of a given bin.
*/
virtual const Types::DataItemRange GetRangeBin( const size_t bin ) const
{
return Types::DataItemRange( this->BinToValue( bin ), this->BinToValue( bin+1 ) );
}
/** Return center of values represented by a certain bin.
*\param bin Index of a bin from the distribution.
*\return Average of upper and lower margin values of the given bin.
*/
virtual Types::DataItem BinToValue ( const size_t bin ) const
{
return this->Superclass::BinToValue( static_cast<size_t>( exp( static_cast<Types::DataItem>( bin ) / (this->GetNumberOfBins()-1) * this->m_LogNumBins ) - 1 ) );
}
protected:
/// Make an identical copy of this object including derived class objects
virtual Self* CloneVirtual() const
{
return new Self( *this );
}
private:
/// Pre-computed log of number of bins.
double m_LogNumBins;
};
//@}
} // namespace cmtk
#endif // #ifndef __cmtkLogHistogram_h_included_
|
