Newer
Older
/* histogram.cc
Mark Woolrich, FMRIB Image Analysis Group
Copyright (C) 1999-2000 University of Oxford */
/* CCOPYRIGHT */
#include "miscmaths.h"
#include "histogram.h"
using namespace std;
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
#ifndef NO_NAMESPACE
namespace MISCMATHS {
#endif
void Histogram::generate()
{
Tracer ts("Histogram::generate");
int size = sourceData.Nrows();
if(calcRange)
{
// calculate range automatically
histMin=histMax=sourceData(1);
for(int i=1; i<=size; i++)
{
if (sourceData(i)>histMax)
histMax=sourceData(i);
if (sourceData(i)<histMin)
histMin=sourceData(i);
}
}
// zero histogram
histogram.ReSize(bins);
histogram=0;
// create histogram; the MIN is so that the maximum value falls in the
// last valid bin, not the (last+1) bin
for(int i=1; i<=size; i++)
{
histogram(getBin(sourceData(i)))++;
}
}
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
void Histogram::smooth()
{
Tracer ts("Histogram::smooth");
ColumnVector newhist=histogram;
// smooth in i direction
newhist=0;
for(int i=1; i<=bins; i++)
{
float val=0.5*histogram(i);
float norm=0.5;
if(i>1)
{
val+=0.2283*(histogram(i-1));
norm+=0.2283;
}
if(i>2)
{
val+=0.0219*(histogram(i-2));
norm+=0.0219;
}
if(i<bins)
{
val+=0.2283*(histogram(i+1));
norm+=0.2283;
}
if(i<bins-1)
{
val+=0.0219*(histogram(i+2));
norm+=0.0219;
}
val/=norm;
newhist(i)=val;
}
histogram=newhist;
}
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
int Histogram::integrate(float value1, float value2) const
{
int upperLimit = getBin(value2);
int sum = 0;
for(int i = getBin(value1)+1; i< upperLimit; i++)
{
sum += (int)histogram(i);
}
return sum;
}
float Histogram::mode() const
{
int maxbin = 0;
int maxnum = 0;
for(int i = 1; i< bins; i++)
{
if((int)histogram(i) > maxnum) {
maxnum = (int)histogram(i);
maxbin = i;
}
}
return getValue(maxbin);
}
#ifndef NO_NAMESPACE
}
#endif