orthanc-stone: Resources/Orthanc/Core/Images/ImageProcessing.cpp annotate

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Fri, 23 Mar 2018 11:04:03 +0100
parents
children

rev	line source
200 03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Orthanc - A Lightweight, RESTful DICOM Store
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	4 * Department, University Hospital of Liege, Belgium
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	5 * Copyright (C) 2017-2018 Osimis S.A., Belgium
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * modify it under the terms of the GNU General Public License as
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9 * published by the Free Software Foundation, either version 3 of the
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * License, or (at your option) any later version.
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * In addition, as a special exception, the copyright holders of this
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * program give permission to link the code of its release with the
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 * OpenSSL project's "OpenSSL" library (or with modified versions of it
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	15 * that use the same license as the "OpenSSL" library), and distribute
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	16 * the linked executables. You must obey the GNU General Public License
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17 * in all respects for all of the code used other than "OpenSSL". If you
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * modify file(s) with this exception, you may extend this exception to
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	19 * your version of the file(s), but you are not obligated to do so. If
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20 * you do not wish to do so, delete this exception statement from your
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21 * version. If you delete this exception statement from all source files
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 * in the program, then also delete it here.
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 *
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24 * This program is distributed in the hope that it will be useful, but
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 * WITHOUT ANY WARRANTY; without even the implied warranty of
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 * General Public License for more details.
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 *
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29 * You should have received a copy of the GNU General Public License
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30 * along with this program. If not, see <http://www.gnu.org/licenses/>.
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 **/
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34 #include "../PrecompiledHeaders.h"
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35 #include "ImageProcessing.h"
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 #include "PixelTraits.h"
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 #include <boost/math/special_functions/round.hpp>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 #include <cassert>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 #include <string.h>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 #include <limits>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 #include <stdint.h>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	45
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46 namespace Orthanc
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	48 template <typename TargetType, typename SourceType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 static void ConvertInternal(ImageAccessor& target,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52 const TargetType minValue = std::numeric_limits<TargetType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53 const TargetType maxValue = std::numeric_limits<TargetType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	54
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	55 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	56 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62 if (static_cast<int32_t>(*s) < static_cast<int32_t>(minValue))
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	63 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	64 *t = minValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	65 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	66 else if (static_cast<int32_t>(*s) > static_cast<int32_t>(maxValue))
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 *t = maxValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72 t = static_cast<TargetType>(s);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	73 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	74 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	75 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	77
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79 template <typename SourceType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 static void ConvertGrayscaleToFloat(ImageAccessor& target,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	81 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	83 assert(sizeof(float) == 4);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	85 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	87 float* t = reinterpret_cast<float*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	88 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	89
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	90 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	91 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	92 t = static_cast<float>(s);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	93 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	94 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	95 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	96
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	97
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98 template <typename TargetType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99 static void ConvertColorToGrayscale(ImageAccessor& target,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	100 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	101 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 assert(source.GetFormat() == PixelFormat_RGB24);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 const TargetType minValue = std::numeric_limits<TargetType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105 const TargetType maxValue = std::numeric_limits<TargetType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110 const uint8_t* s = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	111
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	112 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s += 3)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	114 // Y = 0.2126 R + 0.7152 G + 0.0722 B
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	115 int32_t v = (2126 * static_cast<int32_t>(s[0]) +
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	116 7152 * static_cast<int32_t>(s[1]) +
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 0722 * static_cast<int32_t>(s[2])) / 10000;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 if (static_cast<int32_t>(v) < static_cast<int32_t>(minValue))
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	120 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	121 *t = minValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	123 else if (static_cast<int32_t>(v) > static_cast<int32_t>(maxValue))
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	124 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	125 *t = maxValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	126 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 *t = static_cast<TargetType>(v);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	134
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	135
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 template <typename PixelType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 static void SetInternal(ImageAccessor& image,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138 int64_t constant)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140 for (unsigned int y = 0; y < image.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	141 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	142 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	143
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	144 for (unsigned int x = 0; x < image.GetWidth(); x++, p++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 *p = static_cast<PixelType>(constant);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	150
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	151
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	152 template <typename PixelType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	153 static void GetMinMaxValueInternal(PixelType& minValue,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	154 PixelType& maxValue,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	155 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	156 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	157 // Deal with the special case of empty image
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	158 if (source.GetWidth() == 0 \|\|
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	159 source.GetHeight() == 0)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	160 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 minValue = 0;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162 maxValue = 0;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166 minValue = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	167 maxValue = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	168
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	169 const unsigned int width = source.GetWidth();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173 const PixelType* p = reinterpret_cast<const PixelType*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 for (unsigned int x = 0; x < width; x++, p++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	177 if (*p < minValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	178 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	179 minValue = *p;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	180 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182 if (*p > maxValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	183 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	184 maxValue = *p;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	188 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	189
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 template <typename PixelType>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 static void AddConstantInternal(ImageAccessor& image,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 int64_t constant)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196 if (constant == 0)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	198 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201 const int64_t minValue = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	202 const int64_t maxValue = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	203
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	204 for (unsigned int y = 0; y < image.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	208 for (unsigned int x = 0; x < image.GetWidth(); x++, p++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	209 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	210 int64_t v = static_cast<int64_t>(*p) + constant;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	211
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	212 if (v > maxValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	213 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	214 *p = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	216 else if (v < minValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218 *p = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222 *p = static_cast<PixelType>(v);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	223 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	224 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	230 template <typename PixelType,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	231 bool UseRound>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 static void MultiplyConstantInternal(ImageAccessor& image,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 float factor)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 if (std::abs(factor - 1.0f) <= std::numeric_limits<float>::epsilon())
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 const int64_t minValue = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241 const int64_t maxValue = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 const unsigned int width = image.GetWidth();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 for (unsigned int y = 0; y < image.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248 for (unsigned int x = 0; x < width; x++, p++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 int64_t v;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 if (UseRound)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 // The "round" operation is very costly
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 v = boost::math::llround(static_cast<float>(p) factor);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258 v = static_cast<int64_t>(static_cast<float>(p) factor);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 if (v > maxValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 *p = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	264 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265 else if (v < minValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267 *p = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 *p = static_cast<PixelType>(v);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	272 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	273 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	274 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 template <typename PixelType,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279 bool UseRound>
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	280 static void ShiftScaleInternal(ImageAccessor& image,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281 float offset,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282 float scaling)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	283 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 const float minFloatValue = static_cast<float>(std::numeric_limits<PixelType>::min());
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 const float maxFloatValue = static_cast<float>(std::numeric_limits<PixelType>::max());
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 const PixelType minPixelValue = std::numeric_limits<PixelType>::min();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 const PixelType maxPixelValue = std::numeric_limits<PixelType>::max();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289 const unsigned int height = image.GetHeight();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	290 const unsigned int width = image.GetWidth();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292 for (unsigned int y = 0; y < height; y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 for (unsigned int x = 0; x < width; x++, p++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298 float v = (static_cast<float>(p) + offset) scaling;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 if (v > maxFloatValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 *p = maxPixelValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304 else if (v < minFloatValue)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306 *p = minPixelValue;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	308 else if (UseRound)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	309 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 // The "round" operation is very costly
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	311 *p = static_cast<PixelType>(boost::math::iround(v));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 else
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315 *p = static_cast<PixelType>(v);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 void ImageProcessing::Copy(ImageAccessor& target,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 if (target.GetWidth() != source.GetWidth() \|\|
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 target.GetHeight() != source.GetHeight())
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 throw OrthancException(ErrorCode_IncompatibleImageSize);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 if (target.GetFormat() != source.GetFormat())
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 throw OrthancException(ErrorCode_IncompatibleImageFormat);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 unsigned int lineSize = GetBytesPerPixel(source.GetFormat()) * source.GetWidth();
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 assert(source.GetPitch() >= lineSize && target.GetPitch() >= lineSize);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 memcpy(target.GetRow(y), source.GetConstRow(y), lineSize);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 void ImageProcessing::Convert(ImageAccessor& target,
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 const ImageAccessor& source)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 if (target.GetWidth() != source.GetWidth() \|\|
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 target.GetHeight() != source.GetHeight())
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 throw OrthancException(ErrorCode_IncompatibleImageSize);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356 if (source.GetFormat() == target.GetFormat())
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 Copy(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 if (target.GetFormat() == PixelFormat_Grayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 source.GetFormat() == PixelFormat_Grayscale8)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	364 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 ConvertInternal<uint16_t, uint8_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370 source.GetFormat() == PixelFormat_Grayscale8)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	371 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	372 ConvertInternal<int16_t, uint8_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	373 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	374 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	375
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	376 if (target.GetFormat() == PixelFormat_Grayscale8 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 source.GetFormat() == PixelFormat_Grayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	379 ConvertInternal<uint8_t, uint16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	380 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	381 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	382
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	383 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384 source.GetFormat() == PixelFormat_Grayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 ConvertInternal<int16_t, uint16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 if (target.GetFormat() == PixelFormat_Grayscale8 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391 source.GetFormat() == PixelFormat_SignedGrayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 ConvertInternal<uint8_t, int16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 if (target.GetFormat() == PixelFormat_Grayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398 source.GetFormat() == PixelFormat_SignedGrayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 ConvertInternal<uint16_t, int16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 if (target.GetFormat() == PixelFormat_Grayscale8 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 source.GetFormat() == PixelFormat_RGB24)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 ConvertColorToGrayscale<uint8_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411 if (target.GetFormat() == PixelFormat_Grayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412 source.GetFormat() == PixelFormat_RGB24)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 ConvertColorToGrayscale<uint16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	418 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	419 source.GetFormat() == PixelFormat_RGB24)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	421 ConvertColorToGrayscale<int16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	422 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 if (target.GetFormat() == PixelFormat_Float32 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 source.GetFormat() == PixelFormat_Grayscale8)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 ConvertGrayscaleToFloat<uint8_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432 if (target.GetFormat() == PixelFormat_Float32 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433 source.GetFormat() == PixelFormat_Grayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435 ConvertGrayscaleToFloat<uint16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 if (target.GetFormat() == PixelFormat_Float32 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 source.GetFormat() == PixelFormat_Grayscale32)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 ConvertGrayscaleToFloat<uint32_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446 if (target.GetFormat() == PixelFormat_Float32 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447 source.GetFormat() == PixelFormat_SignedGrayscale16)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	448 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	449 ConvertGrayscaleToFloat<int16_t>(target, source);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	451 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 if (target.GetFormat() == PixelFormat_Grayscale8 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 source.GetFormat() == PixelFormat_RGBA32)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 for (unsigned int x = 0; x < source.GetWidth(); x++, q++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 q = static_cast<uint8_t>((2126 static_cast<uint32_t>(p[0]) +
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 7152 * static_cast<uint32_t>(p[1]) +
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464 0722 * static_cast<uint32_t>(p[2])) / 10000);
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 p += 4;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	466 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	469 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	470 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	471
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	472 if (target.GetFormat() == PixelFormat_RGB24 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473 source.GetFormat() == PixelFormat_RGBA32)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	476 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 for (unsigned int x = 0; x < source.GetWidth(); x++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 q[0] = p[0];
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 q[1] = p[1];
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 q[2] = p[2];
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 p += 4;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 q += 3;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 return;
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 }
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 if (target.GetFormat() == PixelFormat_RGB24 &&
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 source.GetFormat() == PixelFormat_BGRA32)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 for (unsigned int y = 0; y < source.GetHeight(); y++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496 {
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 for (unsigned int x = 0; x < source.GetWidth(); x++)
03afbee0cc7b integration of Orthanc core into Stone Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500 {

200

03afbee0cc7b integration of Orthanc core into Stone