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
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkSparseFrequencyContainer.cxx,v $
Language: C++
Date: $Date: 2009-08-08 20:19:11 $
Version: $Revision: 1.9 $
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.
=========================================================================*/
#include "itkSparseFrequencyContainer.h"
namespace itk{
namespace Statistics{
SparseFrequencyContainer
::SparseFrequencyContainer()
{
m_TotalFrequency = NumericTraits< TotalFrequencyType >::Zero;
}
void
SparseFrequencyContainer
::Initialize(unsigned long)
{
this->SetToZero();
}
void
SparseFrequencyContainer
::SetToZero()
{
typedef FrequencyContainerType::iterator IteratorType;
IteratorType iter = m_FrequencyContainer.begin();
IteratorType end = m_FrequencyContainer.end();
while ( iter != end )
{
iter->second = NumericTraits< FrequencyType >::Zero;
++iter;
}
m_TotalFrequency = NumericTraits< TotalFrequencyType >::Zero;
}
bool
SparseFrequencyContainer
::SetFrequency(const InstanceIdentifier id, const FrequencyType value)
{
// No need to test for bounds because in a map container the
// element is allocated if the key doesn't exist yet
FrequencyType frequency = this->GetFrequency(id);
m_FrequencyContainer[id] = value;
m_TotalFrequency += (value - frequency);
return true;
}
SparseFrequencyContainer::FrequencyType
SparseFrequencyContainer
::GetFrequency(const InstanceIdentifier id) const
{
FrequencyContainerType::const_iterator iter = m_FrequencyContainer.find(id);
if ( iter != m_FrequencyContainer.end() )
{
return iter->second;
}
else
{
return 0;
}
}
bool
SparseFrequencyContainer
::IncreaseFrequency(const InstanceIdentifier id, const FrequencyType value)
{
// No need to test for bounds because in a map container the
// element is allocated if the key doesn't exist yet
FrequencyType frequency = this->GetFrequency(id);
const FrequencyType largestIntegerThatFitsInFloat = 16777216;
if( largestIntegerThatFitsInFloat - frequency < value )
{
itkExceptionMacro("Frequency container saturated for Instance ");
}
else
{
m_FrequencyContainer[id] = frequency + value;
}
m_TotalFrequency += value;
return true;
}
void
SparseFrequencyContainer
::PrintSelf(std::ostream& os, Indent indent) const
{
Superclass::PrintSelf(os,indent);
}
} // end of namespace Statistics
} // end of namespace itk
|
