Mercurial > hg > orthanc-stone

comparison UnitTestsSources/ImageToolboxTests.cpp @ 1286:ddb6676bbcbf bugs/2020-02-invisible-slice

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Added Histogram functions (new files: ImageToolbox...)
author Benjamin Golinvaux <bgo@osimis.io>
date Fri, 14 Feb 2020 14:59:32 +0100
parents
children 86400fa16091
comparison
equal deleted inserted replaced
1284:1bb54005e2bd 1286:ddb6676bbcbf
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-2020 Osimis S.A., Belgium
6 *
7 * This program is free software: you can redistribute it and/or
8 * modify it under the terms of the GNU Affero General Public License
9 * as published by the Free Software Foundation, either version 3 of
10 * the License, or (at your option) any later version.
11 *
12 * This program is distributed in the hope that it will be useful, but
13 * WITHOUT ANY WARRANTY; without even the implied warranty of
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 * Affero General Public License for more details.
16 *
17 * You should have received a copy of the GNU Affero General Public License
18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
19 **/
20
21 #include <Framework/Toolbox/ImageToolbox.h>
22
23 // #include <boost/chrono.hpp>
24 // #include <boost/lexical_cast.hpp>
25
26 #include <Core/Images/Image.h>
27 #include <Core/Images/PixelTraits.h>
28
29 #include "stdint.h"
30
31 #include "gtest/gtest.h"
32
33 #include <cmath>
34
35
36 TEST(ImageToolbox, SimpleHisto_Grayscale8_BinSize1)
37 {
38 using OrthancStone::HistogramData;
39 using OrthancStone::DumpHistogramResult;
40 using OrthancStone::ComputeHistogram;
41
42 const unsigned int W = 16;
43 const unsigned int H = 16;
44
45 // 256/17 = 15,...
46 // 256 % 17 = 1
47 // 0 will be 16 times
48 // 1 will be 15 times
49 // 2 will be 15 times
50 // ...
51 // 16 will be 15 times
52
53 size_t pixCounter = 0;
54
55 std::auto_ptr<Orthanc::Image> image(new Orthanc::Image(
56 Orthanc::PixelFormat_Grayscale8, W, H, false));
57
58 for (unsigned int y = 0; y < H; ++y)
59 {
60 uint8_t* buffer = reinterpret_cast<uint8_t*>(image->GetRow(y));
61 for (unsigned int x = 0; x < W; ++x, ++buffer, ++pixCounter)
62 {
63 *buffer = static_cast<uint8_t>(pixCounter % 17);
64 }
65 }
66
67 HistogramData hd;
68 ComputeHistogram(*image, hd, 1);
69 ASSERT_EQ(-0.5, hd.minValue);
70 ASSERT_EQ(17u, hd.bins.size());
71 ASSERT_EQ(16u, hd.bins[0]);
72 for (size_t i = 1; i < hd.bins.size(); ++i)
73 ASSERT_EQ(15u, hd.bins[i]);
74 }
75
76 TEST(ImageToolbox, SimpleHisto_Grayscale8_BinSize1_FormatString)
77 {
78 using OrthancStone::HistogramData;
79 using OrthancStone::DumpHistogramResult;
80 using OrthancStone::ComputeHistogram;
81
82 const unsigned int W = 16;
83 const unsigned int H = 16;
84
85 // 256/17 = 15,...
86 // 256 % 17 = 1
87 // 0 will be 16 times
88 // 1 will be 15 times
89 // 2 will be 15 times
90 // ...
91 // 16 will be 15 times
92
93 size_t pixCounter = 0;
94
95 std::auto_ptr<Orthanc::Image> image(new Orthanc::Image(
96 Orthanc::PixelFormat_Grayscale8, W, H, false));
97
98 for (unsigned int y = 0; y < H; ++y)
99 {
100 uint8_t* buffer = reinterpret_cast<uint8_t*>(image->GetRow(y));
101 for (unsigned int x = 0; x < W; ++x, ++buffer, ++pixCounter)
102 {
103 *buffer = static_cast<uint8_t>(pixCounter % 17);
104 }
105 }
106
107 HistogramData hd;
108 ComputeHistogram(*image, hd, 1);
109
110 // void DumpHistogramResult(std::string& s, const HistogramData& hd)
111 std::string s;
112 DumpHistogramResult(s, hd);
113 std::cout << s;
114 }
115
116 template<Orthanc::PixelFormat Format>
117 void SimpleHisto_T_BinSize1_2()
118 {
119 using OrthancStone::HistogramData;
120 using OrthancStone::DumpHistogramResult;
121 using OrthancStone::ComputeHistogram;
122
123 const unsigned int W = 16;
124 const unsigned int H = 16;
125
126 // 256/17 = 15,...
127 // 256 % 17 = 1
128 // 0 will be 16 times
129 // 1 will be 15 times
130 // 2 will be 15 times
131 // ...
132 // 16 will be 15 times
133
134 size_t pixCounter = 0;
135
136 std::auto_ptr<Orthanc::Image> image(new Orthanc::Image(
137 Format, W, H, false));
138
139 typedef Orthanc::PixelTraits<Format>::PixelType PixelType;
140
141 PixelType pixValue = 0;
142
143 for (unsigned int y = 0; y < H; ++y)
144 {
145 PixelType* buffer = reinterpret_cast<PixelType*>(image->GetRow(y));
146 for (unsigned int x = 0; x < W; ++x, ++buffer, ++pixCounter)
147 {
148 // 0..99 0..99 0..55
149 *buffer = pixValue;
150 pixValue++;
151 if (pixValue >= 100)
152 pixValue = 0;
153 }
154 }
155
156 HistogramData hd;
157 ComputeHistogram(*image, hd, 1);
158 ASSERT_EQ(-0.5, hd.minValue);
159 ASSERT_EQ(100u, hd.bins.size());
160 for (size_t i = 0; i <= 55; ++i)
161 ASSERT_EQ(3u, hd.bins[i]);
162 for (size_t i = 56; i <= 99; ++i)
163 ASSERT_EQ(2u, hd.bins[i]);
164 }
165
166 TEST(ImageToolbox, SimpleHisto_Grayscale8_BinSize1_2)
167 {
168 SimpleHisto_T_BinSize1_2<Orthanc::PixelFormat_Grayscale8>();
169 }
170
171 TEST(ImageToolbox, SimpleHisto_Grayscale16_BinSize1_2)
172 {
173 SimpleHisto_T_BinSize1_2<Orthanc::PixelFormat_Grayscale16>();
174 }
175
176 TEST(ImageToolbox, SimpleHisto_SignedGrayscale16_BinSize1_2)
177 {
178 SimpleHisto_T_BinSize1_2<Orthanc::PixelFormat_SignedGrayscale16>();
179 }
180
181 TEST(ImageToolbox, SimpleHisto_Grayscale32_BinSize1_2)
182 {
183 SimpleHisto_T_BinSize1_2<Orthanc::PixelFormat_Grayscale32>();
184 }
185
186 template<Orthanc::PixelFormat Format>
187 void SimpleHisto_T_BinSize10_2()
188 {
189 using OrthancStone::HistogramData;
190 using OrthancStone::DumpHistogramResult;
191 using OrthancStone::ComputeHistogram;
192
193 const unsigned int W = 16;
194 const unsigned int H = 16;
195
196 // 256/17 = 15,...
197 // 256 % 17 = 1
198 // 0 will be 16 times
199 // 1 will be 15 times
200 // 2 will be 15 times
201 // ...
202 // 16 will be 15 times
203
204 size_t pixCounter = 0;
205
206 std::auto_ptr<Orthanc::Image> image(new Orthanc::Image(
207 Format, W, H, false));
208
209 typedef Orthanc::PixelTraits<Format>::PixelType PixelType;
210
211 PixelType pixValue = 0;
212
213 for (unsigned int y = 0; y < H; ++y)
214 {
215 PixelType* buffer = reinterpret_cast<PixelType*>(image->GetRow(y));
216 for (unsigned int x = 0; x < W; ++x, ++buffer, ++pixCounter)
217 {
218 // 0..99 0..99 0..55
219 *buffer = pixValue;
220 pixValue++;
221 if (pixValue >= 100)
222 pixValue = 0;
223 }
224 }
225
226 HistogramData hd;
227 ComputeHistogram(*image, hd, 10);
228 ASSERT_EQ(-0.5, hd.minValue);
229 ASSERT_EQ(10u, hd.bins.size());
230
231 for (size_t i = 0; i <= 4; ++i)
232 ASSERT_EQ(30u, hd.bins[i]);
233
234 ASSERT_EQ(26u, hd.bins[5]);
235
236 for (size_t i = 6; i <= 9; ++i)
237 ASSERT_EQ(20u, hd.bins[i]);
238 }
239
240 TEST(ImageToolbox, SimpleHisto_Grayscale8_BinSize10_2)
241 {
242 SimpleHisto_T_BinSize10_2<Orthanc::PixelFormat_Grayscale8>();
243 }
244
245 TEST(ImageToolbox, SimpleHisto_Grayscale16_BinSize10_2)
246 {
247 SimpleHisto_T_BinSize10_2<Orthanc::PixelFormat_Grayscale16>();
248 }
249
250 TEST(ImageToolbox, SimpleHisto_SignedGrayscale16_BinSize10_2)
251 {
252 SimpleHisto_T_BinSize10_2<Orthanc::PixelFormat_SignedGrayscale16>();
253 }
254
255 TEST(ImageToolbox, SimpleHisto_Grayscale32_BinSize10_2)
256 {
257 SimpleHisto_T_BinSize10_2<Orthanc::PixelFormat_Grayscale32>();
258 }
259