Mercurial > hg > orthanc-webviewer
annotate Orthanc/Core/Images/ImageProcessing.cpp @ 159:5dc54316d68b
shared copyright with osimis
author | Sebastien Jodogne <s.jodogne@gmail.com> |
---|---|
date | Wed, 04 Jan 2017 16:41:43 +0100 |
parents | d850500b8ca6 |
children | f17ecfffdc75 |
rev | line source |
---|---|
78 | 1 /** |
2 * Orthanc - A Lightweight, RESTful DICOM Store | |
122 | 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics |
78 | 4 * Department, University Hospital of Liege, Belgium |
159
5dc54316d68b
shared copyright with osimis
Sebastien Jodogne <s.jodogne@gmail.com>
parents:
145
diff
changeset
|
5 * Copyright (C) 2017 Osimis, Belgium |
78 | 6 * |
7 * This program is free software: you can redistribute it and/or | |
8 * modify it under the terms of the GNU General Public License as | |
9 * published by the Free Software Foundation, either version 3 of the | |
10 * License, or (at your option) any later version. | |
11 * | |
12 * In addition, as a special exception, the copyright holders of this | |
13 * program give permission to link the code of its release with the | |
14 * OpenSSL project's "OpenSSL" library (or with modified versions of it | |
15 * that use the same license as the "OpenSSL" library), and distribute | |
16 * the linked executables. You must obey the GNU General Public License | |
17 * in all respects for all of the code used other than "OpenSSL". If you | |
18 * modify file(s) with this exception, you may extend this exception to | |
19 * your version of the file(s), but you are not obligated to do so. If | |
20 * you do not wish to do so, delete this exception statement from your | |
21 * version. If you delete this exception statement from all source files | |
22 * in the program, then also delete it here. | |
23 * | |
24 * This program is distributed in the hope that it will be useful, but | |
25 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
26 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
27 * General Public License for more details. | |
28 * | |
29 * You should have received a copy of the GNU General Public License | |
30 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
31 **/ | |
32 | |
33 | |
34 #include "../PrecompiledHeaders.h" | |
35 #include "ImageProcessing.h" | |
36 | |
37 #include "../OrthancException.h" | |
38 | |
39 #include <boost/math/special_functions/round.hpp> | |
40 | |
41 #include <cassert> | |
42 #include <string.h> | |
43 #include <limits> | |
44 #include <stdint.h> | |
45 | |
46 namespace Orthanc | |
47 { | |
48 template <typename TargetType, typename SourceType> | |
49 static void ConvertInternal(ImageAccessor& target, | |
50 const ImageAccessor& source) | |
51 { | |
52 const TargetType minValue = std::numeric_limits<TargetType>::min(); | |
53 const TargetType maxValue = std::numeric_limits<TargetType>::max(); | |
54 | |
55 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
56 { | |
57 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y)); | |
58 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y)); | |
59 | |
60 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++) | |
61 { | |
62 if (static_cast<int32_t>(*s) < static_cast<int32_t>(minValue)) | |
63 { | |
64 *t = minValue; | |
65 } | |
66 else if (static_cast<int32_t>(*s) > static_cast<int32_t>(maxValue)) | |
67 { | |
68 *t = maxValue; | |
69 } | |
70 else | |
71 { | |
72 *t = static_cast<TargetType>(*s); | |
73 } | |
74 } | |
75 } | |
76 } | |
77 | |
78 | |
129 | 79 template <typename SourceType> |
80 static void ConvertGrayscaleToFloat(ImageAccessor& target, | |
81 const ImageAccessor& source) | |
82 { | |
131 | 83 assert(sizeof(float) == 4); |
84 | |
129 | 85 for (unsigned int y = 0; y < source.GetHeight(); y++) |
86 { | |
87 float* t = reinterpret_cast<float*>(target.GetRow(y)); | |
88 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y)); | |
89 | |
90 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++) | |
91 { | |
92 *t = static_cast<float>(*s); | |
93 } | |
94 } | |
95 } | |
96 | |
97 | |
78 | 98 template <typename TargetType> |
99 static void ConvertColorToGrayscale(ImageAccessor& target, | |
129 | 100 const ImageAccessor& source) |
78 | 101 { |
102 assert(source.GetFormat() == PixelFormat_RGB24); | |
103 | |
104 const TargetType minValue = std::numeric_limits<TargetType>::min(); | |
105 const TargetType maxValue = std::numeric_limits<TargetType>::max(); | |
106 | |
107 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
108 { | |
109 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y)); | |
110 const uint8_t* s = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
111 | |
112 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s += 3) | |
113 { | |
114 // Y = 0.2126 R + 0.7152 G + 0.0722 B | |
115 int32_t v = (2126 * static_cast<int32_t>(s[0]) + | |
116 7152 * static_cast<int32_t>(s[1]) + | |
117 0722 * static_cast<int32_t>(s[2])) / 1000; | |
118 | |
119 if (static_cast<int32_t>(v) < static_cast<int32_t>(minValue)) | |
120 { | |
121 *t = minValue; | |
122 } | |
123 else if (static_cast<int32_t>(v) > static_cast<int32_t>(maxValue)) | |
124 { | |
125 *t = maxValue; | |
126 } | |
127 else | |
128 { | |
129 *t = static_cast<TargetType>(v); | |
130 } | |
131 } | |
132 } | |
133 } | |
134 | |
135 | |
136 template <typename PixelType> | |
137 static void SetInternal(ImageAccessor& image, | |
138 int64_t constant) | |
139 { | |
140 for (unsigned int y = 0; y < image.GetHeight(); y++) | |
141 { | |
142 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y)); | |
143 | |
144 for (unsigned int x = 0; x < image.GetWidth(); x++, p++) | |
145 { | |
146 *p = static_cast<PixelType>(constant); | |
147 } | |
148 } | |
149 } | |
150 | |
151 | |
152 template <typename PixelType> | |
153 static void GetMinMaxValueInternal(PixelType& minValue, | |
154 PixelType& maxValue, | |
155 const ImageAccessor& source) | |
156 { | |
157 // Deal with the special case of empty image | |
158 if (source.GetWidth() == 0 || | |
159 source.GetHeight() == 0) | |
160 { | |
161 minValue = 0; | |
162 maxValue = 0; | |
163 return; | |
164 } | |
165 | |
166 minValue = std::numeric_limits<PixelType>::max(); | |
167 maxValue = std::numeric_limits<PixelType>::min(); | |
168 | |
169 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
170 { | |
171 const PixelType* p = reinterpret_cast<const PixelType*>(source.GetConstRow(y)); | |
172 | |
173 for (unsigned int x = 0; x < source.GetWidth(); x++, p++) | |
174 { | |
175 if (*p < minValue) | |
176 { | |
177 minValue = *p; | |
178 } | |
179 | |
180 if (*p > maxValue) | |
181 { | |
182 maxValue = *p; | |
183 } | |
184 } | |
185 } | |
186 } | |
187 | |
188 | |
189 | |
190 template <typename PixelType> | |
191 static void AddConstantInternal(ImageAccessor& image, | |
192 int64_t constant) | |
193 { | |
194 if (constant == 0) | |
195 { | |
196 return; | |
197 } | |
198 | |
199 const int64_t minValue = std::numeric_limits<PixelType>::min(); | |
200 const int64_t maxValue = std::numeric_limits<PixelType>::max(); | |
201 | |
202 for (unsigned int y = 0; y < image.GetHeight(); y++) | |
203 { | |
204 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y)); | |
205 | |
206 for (unsigned int x = 0; x < image.GetWidth(); x++, p++) | |
207 { | |
208 int64_t v = static_cast<int64_t>(*p) + constant; | |
209 | |
210 if (v > maxValue) | |
211 { | |
212 *p = std::numeric_limits<PixelType>::max(); | |
213 } | |
214 else if (v < minValue) | |
215 { | |
216 *p = std::numeric_limits<PixelType>::min(); | |
217 } | |
218 else | |
219 { | |
220 *p = static_cast<PixelType>(v); | |
221 } | |
222 } | |
223 } | |
224 } | |
225 | |
226 | |
227 | |
228 template <typename PixelType> | |
229 void MultiplyConstantInternal(ImageAccessor& image, | |
230 float factor) | |
231 { | |
232 if (std::abs(factor - 1.0f) <= std::numeric_limits<float>::epsilon()) | |
233 { | |
234 return; | |
235 } | |
236 | |
237 const int64_t minValue = std::numeric_limits<PixelType>::min(); | |
238 const int64_t maxValue = std::numeric_limits<PixelType>::max(); | |
239 | |
240 for (unsigned int y = 0; y < image.GetHeight(); y++) | |
241 { | |
242 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y)); | |
243 | |
244 for (unsigned int x = 0; x < image.GetWidth(); x++, p++) | |
245 { | |
246 int64_t v = boost::math::llround(static_cast<float>(*p) * factor); | |
247 | |
248 if (v > maxValue) | |
249 { | |
250 *p = std::numeric_limits<PixelType>::max(); | |
251 } | |
252 else if (v < minValue) | |
253 { | |
254 *p = std::numeric_limits<PixelType>::min(); | |
255 } | |
256 else | |
257 { | |
258 *p = static_cast<PixelType>(v); | |
259 } | |
260 } | |
261 } | |
262 } | |
263 | |
264 | |
265 template <typename PixelType> | |
266 void ShiftScaleInternal(ImageAccessor& image, | |
267 float offset, | |
268 float scaling) | |
269 { | |
270 const float minValue = static_cast<float>(std::numeric_limits<PixelType>::min()); | |
271 const float maxValue = static_cast<float>(std::numeric_limits<PixelType>::max()); | |
272 | |
273 for (unsigned int y = 0; y < image.GetHeight(); y++) | |
274 { | |
275 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y)); | |
276 | |
277 for (unsigned int x = 0; x < image.GetWidth(); x++, p++) | |
278 { | |
279 float v = (static_cast<float>(*p) + offset) * scaling; | |
280 | |
281 if (v > maxValue) | |
282 { | |
283 *p = std::numeric_limits<PixelType>::max(); | |
284 } | |
285 else if (v < minValue) | |
286 { | |
287 *p = std::numeric_limits<PixelType>::min(); | |
288 } | |
289 else | |
290 { | |
291 *p = static_cast<PixelType>(boost::math::iround(v)); | |
292 } | |
293 } | |
294 } | |
295 } | |
296 | |
297 | |
298 void ImageProcessing::Copy(ImageAccessor& target, | |
299 const ImageAccessor& source) | |
300 { | |
301 if (target.GetWidth() != source.GetWidth() || | |
302 target.GetHeight() != source.GetHeight()) | |
303 { | |
304 throw OrthancException(ErrorCode_IncompatibleImageSize); | |
305 } | |
306 | |
307 if (target.GetFormat() != source.GetFormat()) | |
308 { | |
309 throw OrthancException(ErrorCode_IncompatibleImageFormat); | |
310 } | |
311 | |
312 unsigned int lineSize = GetBytesPerPixel(source.GetFormat()) * source.GetWidth(); | |
313 | |
314 assert(source.GetPitch() >= lineSize && target.GetPitch() >= lineSize); | |
315 | |
316 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
317 { | |
318 memcpy(target.GetRow(y), source.GetConstRow(y), lineSize); | |
319 } | |
320 } | |
321 | |
322 | |
323 void ImageProcessing::Convert(ImageAccessor& target, | |
324 const ImageAccessor& source) | |
325 { | |
326 if (target.GetWidth() != source.GetWidth() || | |
327 target.GetHeight() != source.GetHeight()) | |
328 { | |
329 throw OrthancException(ErrorCode_IncompatibleImageSize); | |
330 } | |
331 | |
332 if (source.GetFormat() == target.GetFormat()) | |
333 { | |
334 Copy(target, source); | |
335 return; | |
336 } | |
337 | |
338 if (target.GetFormat() == PixelFormat_Grayscale16 && | |
339 source.GetFormat() == PixelFormat_Grayscale8) | |
340 { | |
341 ConvertInternal<uint16_t, uint8_t>(target, source); | |
342 return; | |
343 } | |
344 | |
345 if (target.GetFormat() == PixelFormat_SignedGrayscale16 && | |
346 source.GetFormat() == PixelFormat_Grayscale8) | |
347 { | |
348 ConvertInternal<int16_t, uint8_t>(target, source); | |
349 return; | |
350 } | |
351 | |
352 if (target.GetFormat() == PixelFormat_Grayscale8 && | |
353 source.GetFormat() == PixelFormat_Grayscale16) | |
354 { | |
355 ConvertInternal<uint8_t, uint16_t>(target, source); | |
356 return; | |
357 } | |
358 | |
359 if (target.GetFormat() == PixelFormat_SignedGrayscale16 && | |
360 source.GetFormat() == PixelFormat_Grayscale16) | |
361 { | |
362 ConvertInternal<int16_t, uint16_t>(target, source); | |
363 return; | |
364 } | |
365 | |
366 if (target.GetFormat() == PixelFormat_Grayscale8 && | |
367 source.GetFormat() == PixelFormat_SignedGrayscale16) | |
368 { | |
369 ConvertInternal<uint8_t, int16_t>(target, source); | |
370 return; | |
371 } | |
372 | |
373 if (target.GetFormat() == PixelFormat_Grayscale16 && | |
374 source.GetFormat() == PixelFormat_SignedGrayscale16) | |
375 { | |
376 ConvertInternal<uint16_t, int16_t>(target, source); | |
377 return; | |
378 } | |
379 | |
380 if (target.GetFormat() == PixelFormat_Grayscale8 && | |
381 source.GetFormat() == PixelFormat_RGB24) | |
382 { | |
383 ConvertColorToGrayscale<uint8_t>(target, source); | |
384 return; | |
385 } | |
386 | |
387 if (target.GetFormat() == PixelFormat_Grayscale16 && | |
388 source.GetFormat() == PixelFormat_RGB24) | |
389 { | |
390 ConvertColorToGrayscale<uint16_t>(target, source); | |
391 return; | |
392 } | |
393 | |
394 if (target.GetFormat() == PixelFormat_SignedGrayscale16 && | |
395 source.GetFormat() == PixelFormat_RGB24) | |
396 { | |
397 ConvertColorToGrayscale<int16_t>(target, source); | |
398 return; | |
399 } | |
400 | |
129 | 401 if (target.GetFormat() == PixelFormat_Float32 && |
402 source.GetFormat() == PixelFormat_Grayscale8) | |
403 { | |
404 ConvertGrayscaleToFloat<uint8_t>(target, source); | |
405 return; | |
406 } | |
407 | |
408 if (target.GetFormat() == PixelFormat_Float32 && | |
409 source.GetFormat() == PixelFormat_Grayscale16) | |
410 { | |
411 ConvertGrayscaleToFloat<uint16_t>(target, source); | |
412 return; | |
413 } | |
414 | |
415 if (target.GetFormat() == PixelFormat_Float32 && | |
416 source.GetFormat() == PixelFormat_SignedGrayscale16) | |
417 { | |
418 ConvertGrayscaleToFloat<int16_t>(target, source); | |
419 return; | |
420 } | |
421 | |
78 | 422 if (target.GetFormat() == PixelFormat_Grayscale8 && |
423 source.GetFormat() == PixelFormat_RGBA32) | |
424 { | |
425 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
426 { | |
427 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
428 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
429 for (unsigned int x = 0; x < source.GetWidth(); x++, q++) | |
430 { | |
431 *q = static_cast<uint8_t>((2126 * static_cast<uint32_t>(p[0]) + | |
432 7152 * static_cast<uint32_t>(p[1]) + | |
433 0722 * static_cast<uint32_t>(p[2])) / 10000); | |
434 p += 4; | |
435 } | |
436 } | |
437 | |
438 return; | |
439 } | |
440 | |
441 if (target.GetFormat() == PixelFormat_RGB24 && | |
442 source.GetFormat() == PixelFormat_RGBA32) | |
443 { | |
444 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
445 { | |
446 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
447 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
448 for (unsigned int x = 0; x < source.GetWidth(); x++) | |
449 { | |
450 q[0] = p[0]; | |
451 q[1] = p[1]; | |
452 q[2] = p[2]; | |
453 p += 4; | |
454 q += 3; | |
455 } | |
456 } | |
457 | |
458 return; | |
459 } | |
460 | |
461 if (target.GetFormat() == PixelFormat_RGBA32 && | |
462 source.GetFormat() == PixelFormat_RGB24) | |
463 { | |
464 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
465 { | |
466 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
467 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
468 for (unsigned int x = 0; x < source.GetWidth(); x++) | |
469 { | |
470 q[0] = p[0]; | |
471 q[1] = p[1]; | |
472 q[2] = p[2]; | |
473 q[3] = 255; // Set the alpha channel to full opacity | |
474 p += 3; | |
475 q += 4; | |
476 } | |
477 } | |
478 | |
479 return; | |
480 } | |
481 | |
129 | 482 if (target.GetFormat() == PixelFormat_RGB24 && |
483 source.GetFormat() == PixelFormat_Grayscale8) | |
484 { | |
485 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
486 { | |
487 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
488 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
489 for (unsigned int x = 0; x < source.GetWidth(); x++) | |
490 { | |
491 q[0] = *p; | |
492 q[1] = *p; | |
493 q[2] = *p; | |
494 p += 1; | |
495 q += 3; | |
496 } | |
497 } | |
498 | |
499 return; | |
500 } | |
501 | |
502 if (target.GetFormat() == PixelFormat_RGBA32 && | |
503 source.GetFormat() == PixelFormat_Grayscale8) | |
504 { | |
505 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
506 { | |
507 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
508 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
509 for (unsigned int x = 0; x < source.GetWidth(); x++) | |
510 { | |
511 q[0] = *p; | |
512 q[1] = *p; | |
513 q[2] = *p; | |
514 q[3] = 255; | |
515 p += 1; | |
516 q += 4; | |
517 } | |
518 } | |
519 | |
520 return; | |
521 } | |
522 | |
145 | 523 if (target.GetFormat() == PixelFormat_BGRA32 && |
524 source.GetFormat() == PixelFormat_RGB24) | |
525 { | |
526 for (unsigned int y = 0; y < source.GetHeight(); y++) | |
527 { | |
528 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y)); | |
529 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y)); | |
530 for (unsigned int x = 0; x < source.GetWidth(); x++) | |
531 { | |
532 q[0] = p[2]; | |
533 q[1] = p[1]; | |
534 q[2] = p[0]; | |
535 q[3] = 255; | |
536 p += 3; | |
537 q += 4; | |
538 } | |
539 } | |
540 | |
541 return; | |
542 } | |
543 | |
78 | 544 throw OrthancException(ErrorCode_NotImplemented); |
545 } | |
546 | |
547 | |
548 | |
549 void ImageProcessing::Set(ImageAccessor& image, | |
550 int64_t value) | |
551 { | |
552 switch (image.GetFormat()) | |
553 { | |
554 case PixelFormat_Grayscale8: | |
555 SetInternal<uint8_t>(image, value); | |
556 return; | |
557 | |
558 case PixelFormat_Grayscale16: | |
559 SetInternal<uint16_t>(image, value); | |
560 return; | |
561 | |
562 case PixelFormat_SignedGrayscale16: | |
563 SetInternal<int16_t>(image, value); | |
564 return; | |
565 | |
131 | 566 case PixelFormat_Float32: |
567 assert(sizeof(float) == 4); | |
568 SetInternal<float>(image, value); | |
569 return; | |
570 | |
78 | 571 default: |
572 throw OrthancException(ErrorCode_NotImplemented); | |
573 } | |
574 } | |
575 | |
576 | |
145 | 577 void ImageProcessing::Set(ImageAccessor& image, |
578 uint8_t red, | |
579 uint8_t green, | |
580 uint8_t blue, | |
581 uint8_t alpha) | |
582 { | |
583 uint8_t p[4]; | |
584 unsigned int size; | |
585 | |
586 switch (image.GetFormat()) | |
587 { | |
588 case PixelFormat_RGBA32: | |
589 p[0] = red; | |
590 p[1] = green; | |
591 p[2] = blue; | |
592 p[3] = alpha; | |
593 size = 4; | |
594 break; | |
595 | |
596 case PixelFormat_BGRA32: | |
597 p[0] = blue; | |
598 p[1] = green; | |
599 p[2] = red; | |
600 p[3] = alpha; | |
601 size = 4; | |
602 break; | |
603 | |
604 case PixelFormat_RGB24: | |
605 p[0] = red; | |
606 p[1] = green; | |
607 p[2] = blue; | |
608 size = 3; | |
609 break; | |
610 | |
611 default: | |
612 throw OrthancException(ErrorCode_NotImplemented); | |
613 } | |
614 | |
615 for (unsigned int y = 0; y < image.GetHeight(); y++) | |
616 { | |
617 uint8_t* q = reinterpret_cast<uint8_t*>(image.GetRow(y)); | |
618 | |
619 for (unsigned int x = 0; x < image.GetWidth(); x++) | |
620 { | |
621 for (unsigned int i = 0; i < size; i++) | |
622 { | |
623 q[i] = p[i]; | |
624 } | |
625 | |
626 q += size; | |
627 } | |
628 } | |
629 } | |
630 | |
631 | |
78 | 632 void ImageProcessing::ShiftRight(ImageAccessor& image, |
633 unsigned int shift) | |
634 { | |
635 if (image.GetWidth() == 0 || | |
636 image.GetHeight() == 0 || | |
637 shift == 0) | |
638 { | |
639 // Nothing to do | |
640 return; | |
641 } | |
642 | |
643 throw OrthancException(ErrorCode_NotImplemented); | |
644 } | |
645 | |
646 | |
647 void ImageProcessing::GetMinMaxValue(int64_t& minValue, | |
648 int64_t& maxValue, | |
649 const ImageAccessor& image) | |
650 { | |
651 switch (image.GetFormat()) | |
652 { | |
653 case PixelFormat_Grayscale8: | |
654 { | |
655 uint8_t a, b; | |
656 GetMinMaxValueInternal<uint8_t>(a, b, image); | |
657 minValue = a; | |
658 maxValue = b; | |
659 break; | |
660 } | |
661 | |
662 case PixelFormat_Grayscale16: | |
663 { | |
664 uint16_t a, b; | |
665 GetMinMaxValueInternal<uint16_t>(a, b, image); | |
666 minValue = a; | |
667 maxValue = b; | |
668 break; | |
669 } | |
670 | |
671 case PixelFormat_SignedGrayscale16: | |
672 { | |
673 int16_t a, b; | |
674 GetMinMaxValueInternal<int16_t>(a, b, image); | |
675 minValue = a; | |
676 maxValue = b; | |
677 break; | |
678 } | |
679 | |
680 default: | |
681 throw OrthancException(ErrorCode_NotImplemented); | |
682 } | |
683 } | |
684 | |
685 | |
686 | |
687 void ImageProcessing::AddConstant(ImageAccessor& image, | |
688 int64_t value) | |
689 { | |
690 switch (image.GetFormat()) | |
691 { | |
692 case PixelFormat_Grayscale8: | |
693 AddConstantInternal<uint8_t>(image, value); | |
694 return; | |
695 | |
696 case PixelFormat_Grayscale16: | |
697 AddConstantInternal<uint16_t>(image, value); | |
698 return; | |
699 | |
700 case PixelFormat_SignedGrayscale16: | |
701 AddConstantInternal<int16_t>(image, value); | |
702 return; | |
703 | |
704 default: | |
705 throw OrthancException(ErrorCode_NotImplemented); | |
706 } | |
707 } | |
708 | |
709 | |
710 void ImageProcessing::MultiplyConstant(ImageAccessor& image, | |
711 float factor) | |
712 { | |
713 switch (image.GetFormat()) | |
714 { | |
715 case PixelFormat_Grayscale8: | |
716 MultiplyConstantInternal<uint8_t>(image, factor); | |
717 return; | |
718 | |
719 case PixelFormat_Grayscale16: | |
720 MultiplyConstantInternal<uint16_t>(image, factor); | |
721 return; | |
722 | |
723 case PixelFormat_SignedGrayscale16: | |
724 MultiplyConstantInternal<int16_t>(image, factor); | |
725 return; | |
726 | |
727 default: | |
728 throw OrthancException(ErrorCode_NotImplemented); | |
729 } | |
730 } | |
731 | |
732 | |
733 void ImageProcessing::ShiftScale(ImageAccessor& image, | |
734 float offset, | |
735 float scaling) | |
736 { | |
737 switch (image.GetFormat()) | |
738 { | |
739 case PixelFormat_Grayscale8: | |
740 ShiftScaleInternal<uint8_t>(image, offset, scaling); | |
741 return; | |
742 | |
743 case PixelFormat_Grayscale16: | |
744 ShiftScaleInternal<uint16_t>(image, offset, scaling); | |
745 return; | |
746 | |
747 case PixelFormat_SignedGrayscale16: | |
748 ShiftScaleInternal<int16_t>(image, offset, scaling); | |
749 return; | |
750 | |
751 default: | |
752 throw OrthancException(ErrorCode_NotImplemented); | |
753 } | |
754 } | |
755 } |