Mercurial > hg > orthanc-stone

comparison OrthancStone/Sources/Toolbox/Internals/BucketMapper.cpp @ 1891:3716d72161d2

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
reorganization
author Sebastien Jodogne <s.jodogne@gmail.com>
date Wed, 19 Jan 2022 12:32:15 +0100
parents
children bd527bbc34df
comparison
equal deleted inserted replaced
1890:6ce81914f7e4 1891:3716d72161d2
1 /**
2 * Stone of Orthanc
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
4 * Department, University Hospital of Liege, Belgium
5 * Copyright (C) 2017-2022 Osimis S.A., Belgium
6 * Copyright (C) 2021-2022 Sebastien Jodogne, ICTEAM UCLouvain, Belgium
7 *
8 * This program is free software: you can redistribute it and/or
9 * modify it under the terms of the GNU Lesser General Public License
10 * as published by the Free Software Foundation, either version 3 of
11 * the License, or (at your option) any later version.
12 *
13 * This program is distributed in the hope that it will be useful, but
14 * WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
19 * License along with this program. If not, see
20 * <http://www.gnu.org/licenses/>.
21 **/
22
23
24 #include "BucketMapper.h"
25
26 #include <OrthancException.h>
27
28 #include <cassert>
29 #include <cmath>
30
31
32 namespace OrthancStone
33 {
34 namespace Internals
35 {
36 BucketMapper::BucketMapper(double minValue,
37 double maxValue,
38 size_t bucketsCount) :
39 minValue_(minValue),
40 maxValue_(maxValue),
41 bucketsCount_(bucketsCount)
42 {
43 if (minValue >= maxValue ||
44 bucketsCount <= 0)
45 {
46 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
47 }
48 }
49
50
51 void BucketMapper::CheckIndex(size_t i) const
52 {
53 if (i >= bucketsCount_)
54 {
55 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
56 }
57 }
58
59
60 double BucketMapper::GetBucketLow(size_t i) const
61 {
62 CheckIndex(i);
63 double alpha = static_cast<double>(i) / static_cast<double>(bucketsCount_);
64 return (1.0 - alpha) * minValue_ + alpha * maxValue_;
65 }
66
67
68 double BucketMapper::GetBucketHigh(size_t i) const
69 {
70 CheckIndex(i);
71 double alpha = static_cast<double>(i + 1) / static_cast<double>(bucketsCount_);
72 return (1.0 - alpha) * minValue_ + alpha * maxValue_;
73 }
74
75
76 double BucketMapper::GetBucketCenter(size_t i) const
77 {
78 CheckIndex(i);
79 double alpha = (static_cast<double>(i) + 0.5) / static_cast<double>(bucketsCount_);
80 return (1.0 - alpha) * minValue_ + alpha * maxValue_;
81 }
82
83
84 size_t BucketMapper::GetBucketIndex(double value) const
85 {
86 if (value < minValue_ ||
87 value > maxValue_)
88 {
89 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
90 }
91 else
92 {
93 double tmp = (value - minValue_) / (maxValue_ - minValue_) * static_cast<double>(bucketsCount_);
94 assert(tmp >= 0 && tmp <= static_cast<double>(bucketsCount_));
95
96 size_t bucket = static_cast<unsigned int>(std::floor(tmp));
97 if (bucket == bucketsCount_) // This is the case if "value == maxValue_"
98 {
99 return bucketsCount_ - 1;
100 }
101 else
102 {
103 return bucket;
104 }
105 }
106 }
107 }
108 }