orthanc: Core/Images/ImageProcessing.cpp annotate

annotate Core/Images/ImageProcessing.cpp @ 2495:cd7b854dbc05

convert grayscale16 to bgra32

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Tue, 20 Mar 2018 13:26:00 +0100
parents	a9459c7d4bc2
children	3d65adee289a

rev	line source
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 /**
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2 * Orthanc - A Lightweight, RESTful DICOM Store
1900 b1291df2f780 2016 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1612 diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
1288 6e7e5ed91c2d upgrade to year 2015 Sebastien Jodogne <s.jodogne@gmail.com> parents: 993 diff changeset	4 * Department, University Hospital of Liege, Belgium
2447 878b59270859 upgrade to year 2018 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2423 diff changeset	5 * Copyright (C) 2017-2018 Osimis S.A., Belgium
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 *
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7 * This program is free software: you can redistribute it and/or
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 * modify it under the terms of the GNU General Public License as
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9 * published by the Free Software Foundation, either version 3 of the
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 * License, or (at your option) any later version.
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 *
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 * In addition, as a special exception, the copyright holders of this
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 * program give permission to link the code of its release with the
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 * OpenSSL project's "OpenSSL" library (or with modified versions of it
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	15 * that use the same license as the "OpenSSL" library), and distribute
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	16 * the linked executables. You must obey the GNU General Public License
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17 * in all respects for all of the code used other than "OpenSSL". If you
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18 * modify file(s) with this exception, you may extend this exception to
839be3022203 DicomImageInformation 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
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20 * you do not wish to do so, delete this exception statement from your
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21 * version. If you delete this exception statement from all source files
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 * in the program, then also delete it here.
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 *
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24 * This program is distributed in the hope that it will be useful, but
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 * WITHOUT ANY WARRANTY; without even the implied warranty of
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 * General Public License for more details.
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 *
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29 * You should have received a copy of the GNU General Public License
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30 * along with this program. If not, see <http://www.gnu.org/licenses/>.
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 **/
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34 #include "../PrecompiledHeaders.h"
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35 #include "ImageProcessing.h"
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36
2489 e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	37 #include "PixelTraits.h"
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	39 #include <boost/math/special_functions/round.hpp>
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	40
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 #include <cassert>
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 #include <string.h>
859 610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	43 #include <limits>
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	44 #include <stdint.h>
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	45
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46 namespace Orthanc
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47 {
859 610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	48 template <typename TargetType, typename SourceType>
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	49 static void ConvertInternal(ImageAccessor& target,
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	50 const ImageAccessor& source)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	51 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	52 const TargetType minValue = std::numeric_limits<TargetType>::min();
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	53 const TargetType maxValue = std::numeric_limits<TargetType>::max();
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	54
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	55 for (unsigned int y = 0; y < source.GetHeight(); y++)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	56 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	57 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y));
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	58 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y));
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	59
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	60 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	61 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	62 if (static_cast<int32_t>(*s) < static_cast<int32_t>(minValue))
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	63 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	64 *t = minValue;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	65 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	66 else if (static_cast<int32_t>(*s) > static_cast<int32_t>(maxValue))
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	67 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	68 *t = maxValue;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	69 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	70 else
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	71 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	72 t = static_cast<TargetType>(s);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	73 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	74 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	75 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	76 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	77
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	78
1993 e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	79 template <typename SourceType>
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	80 static void ConvertGrayscaleToFloat(ImageAccessor& target,
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	81 const ImageAccessor& source)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	82 {
1994 4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	83 assert(sizeof(float) == 4);
4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	84
1993 e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	85 for (unsigned int y = 0; y < source.GetHeight(); y++)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	86 {
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	87 float* t = reinterpret_cast<float*>(target.GetRow(y));
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	88 const SourceType* s = reinterpret_cast<const SourceType*>(source.GetConstRow(y));
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	89
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	90 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s++)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	91 {
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	92 t = static_cast<float>(s);
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	93 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	94 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	95 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	96
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	97
993 501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	98 template <typename TargetType>
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	99 static void ConvertColorToGrayscale(ImageAccessor& target,
1993 e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	100 const ImageAccessor& source)
993 501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	101 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	102 assert(source.GetFormat() == PixelFormat_RGB24);
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	103
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	104 const TargetType minValue = std::numeric_limits<TargetType>::min();
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	105 const TargetType maxValue = std::numeric_limits<TargetType>::max();
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	106
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	107 for (unsigned int y = 0; y < source.GetHeight(); y++)
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	108 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	109 TargetType* t = reinterpret_cast<TargetType*>(target.GetRow(y));
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	110 const uint8_t* s = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	111
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	112 for (unsigned int x = 0; x < source.GetWidth(); x++, t++, s += 3)
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	113 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	114 // Y = 0.2126 R + 0.7152 G + 0.0722 B
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	115 int32_t v = (2126 * static_cast<int32_t>(s[0]) +
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	116 7152 * static_cast<int32_t>(s[1]) +
2491 a9459c7d4bc2 fix color-to-grayscale conversion Sebastien Jodogne <s.jodogne@gmail.com> parents: 2490 diff changeset	117 0722 * static_cast<int32_t>(s[2])) / 10000;
993 501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	118
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	119 if (static_cast<int32_t>(v) < static_cast<int32_t>(minValue))
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	120 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	121 *t = minValue;
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	122 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	123 else if (static_cast<int32_t>(v) > static_cast<int32_t>(maxValue))
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	124 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	125 *t = maxValue;
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	126 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	127 else
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	128 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	129 *t = static_cast<TargetType>(v);
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	130 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	131 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	132 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	133 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	134
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	135
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	136 template <typename PixelType>
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	137 static void SetInternal(ImageAccessor& image,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	138 int64_t constant)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	139 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	140 for (unsigned int y = 0; y < image.GetHeight(); y++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	141 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	142 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	143
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	144 for (unsigned int x = 0; x < image.GetWidth(); x++, p++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	145 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	146 *p = static_cast<PixelType>(constant);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	147 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	148 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	149 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	150
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	151
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	152 template <typename PixelType>
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	153 static void GetMinMaxValueInternal(PixelType& minValue,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	154 PixelType& maxValue,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	155 const ImageAccessor& source)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	156 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	157 // Deal with the special case of empty image
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	158 if (source.GetWidth() == 0 \|\|
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	159 source.GetHeight() == 0)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	160 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	161 minValue = 0;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	162 maxValue = 0;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	163 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	164 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	165
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	166 minValue = std::numeric_limits<PixelType>::max();
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	167 maxValue = std::numeric_limits<PixelType>::min();
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	168
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	169 const unsigned int width = source.GetWidth();
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	170
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	171 for (unsigned int y = 0; y < source.GetHeight(); y++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	172 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	173 const PixelType* p = reinterpret_cast<const PixelType*>(source.GetConstRow(y));
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	174
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	175 for (unsigned int x = 0; x < width; x++, p++)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	176 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	177 if (*p < minValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	178 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	179 minValue = *p;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	180 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	181
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	182 if (*p > maxValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	183 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	184 maxValue = *p;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	185 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	186 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	187 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	188 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	189
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	190
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	191
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	192 template <typename PixelType>
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	193 static void AddConstantInternal(ImageAccessor& image,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	194 int64_t constant)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	195 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	196 if (constant == 0)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	197 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	198 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	199 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	200
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	201 const int64_t minValue = std::numeric_limits<PixelType>::min();
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	202 const int64_t maxValue = std::numeric_limits<PixelType>::max();
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	203
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	204 for (unsigned int y = 0; y < image.GetHeight(); y++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	205 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	206 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	207
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	208 for (unsigned int x = 0; x < image.GetWidth(); x++, p++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	209 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	210 int64_t v = static_cast<int64_t>(*p) + constant;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	211
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	212 if (v > maxValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	213 {
876 e21d1a5f5934 fix msvc warnings jodogne parents: 863 diff changeset	214 *p = std::numeric_limits<PixelType>::max();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	215 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	216 else if (v < minValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	217 {
876 e21d1a5f5934 fix msvc warnings jodogne parents: 863 diff changeset	218 *p = std::numeric_limits<PixelType>::min();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	219 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	220 else
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	221 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	222 *p = static_cast<PixelType>(v);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	223 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	224 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	225 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	226 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	227
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	228
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	229
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	230 template <typename PixelType,
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	231 bool UseRound>
2489 e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	232 static void MultiplyConstantInternal(ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	233 float factor)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	234 {
1334 7f0aa3c0f659 code improvements Sebastien Jodogne <s.jodogne@gmail.com> parents: 1288 diff changeset	235 if (std::abs(factor - 1.0f) <= std::numeric_limits<float>::epsilon())
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	236 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	237 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	238 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	239
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	240 const int64_t minValue = std::numeric_limits<PixelType>::min();
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	241 const int64_t maxValue = std::numeric_limits<PixelType>::max();
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	242 const unsigned int width = image.GetWidth();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	243
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	244 for (unsigned int y = 0; y < image.GetHeight(); y++)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	245 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	246 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	247
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	248 for (unsigned int x = 0; x < width; x++, p++)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	249 {
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	250 int64_t v;
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	251 if (UseRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	252 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	253 // The "round" operation is very costly
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	254 v = boost::math::llround(static_cast<float>(p) factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	255 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	256 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	257 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	258 v = static_cast<int64_t>(static_cast<float>(p) factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	259 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	260
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	261 if (v > maxValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	262 {
876 e21d1a5f5934 fix msvc warnings jodogne parents: 863 diff changeset	263 *p = std::numeric_limits<PixelType>::max();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	264 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	265 else if (v < minValue)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	266 {
876 e21d1a5f5934 fix msvc warnings jodogne parents: 863 diff changeset	267 *p = std::numeric_limits<PixelType>::min();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	268 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	269 else
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	270 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	271 *p = static_cast<PixelType>(v);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	272 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	273 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	274 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	275 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	276
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	277
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	278 template <typename PixelType,
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	279 bool UseRound>
2489 e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	280 static void ShiftScaleInternal(ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	281 float offset,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	282 float scaling)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	283 {
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	284 const float minFloatValue = static_cast<float>(std::numeric_limits<PixelType>::min());
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	285 const float maxFloatValue = static_cast<float>(std::numeric_limits<PixelType>::max());
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	286 const PixelType minPixelValue = std::numeric_limits<PixelType>::min();
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	287 const PixelType maxPixelValue = std::numeric_limits<PixelType>::max();
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	288
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	289 const unsigned int height = image.GetHeight();
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	290 const unsigned int width = image.GetWidth();
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	291
509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	292 for (unsigned int y = 0; y < height; y++)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	293 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	294 PixelType* p = reinterpret_cast<PixelType*>(image.GetRow(y));
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	295
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	296 for (unsigned int x = 0; x < width; x++, p++)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	297 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	298 float v = (static_cast<float>(p) + offset) scaling;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	299
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	300 if (v > maxFloatValue)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	301 {
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	302 *p = maxPixelValue;
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	303 }
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	304 else if (v < minFloatValue)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	305 {
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	306 *p = minPixelValue;
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	307 }
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	308 else if (UseRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	309 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	310 // The "round" operation is very costly
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	311 *p = static_cast<PixelType>(boost::math::iround(v));
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	312 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	313 else
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	314 {
2482 509041cb57db speedup by truncating instead of rounding Sebastien Jodogne <s.jodogne@gmail.com> parents: 2447 diff changeset	315 *p = static_cast<PixelType>(v);
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	316 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	317 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	318 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	319 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	320
859 610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	321
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 void ImageProcessing::Copy(ImageAccessor& target,
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 const ImageAccessor& source)
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 if (target.GetWidth() != source.GetWidth() \|\|
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 target.GetHeight() != source.GetHeight())
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 throw OrthancException(ErrorCode_IncompatibleImageSize);
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 }
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 if (target.GetFormat() != source.GetFormat())
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 throw OrthancException(ErrorCode_IncompatibleImageFormat);
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 }
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336 unsigned int lineSize = GetBytesPerPixel(source.GetFormat()) * source.GetWidth();
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 assert(source.GetPitch() >= lineSize && target.GetPitch() >= lineSize);
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 for (unsigned int y = 0; y < source.GetHeight(); y++)
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 {
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	342 memcpy(target.GetRow(y), source.GetConstRow(y), lineSize);
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 }
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 }
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	346
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 void ImageProcessing::Convert(ImageAccessor& target,
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 const ImageAccessor& source)
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 if (target.GetWidth() != source.GetWidth() \|\|
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 target.GetHeight() != source.GetHeight())
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 throw OrthancException(ErrorCode_IncompatibleImageSize);
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354 }
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356 if (source.GetFormat() == target.GetFormat())
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	357 {
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	358 Copy(target, source);
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 return;
839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 }
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	361
859 610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	362 if (target.GetFormat() == PixelFormat_Grayscale16 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	363 source.GetFormat() == PixelFormat_Grayscale8)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	364 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	365 ConvertInternal<uint16_t, uint8_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	366 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	367 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	368
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	369 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	370 source.GetFormat() == PixelFormat_Grayscale8)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	371 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	372 ConvertInternal<int16_t, uint8_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	373 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	374 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	375
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	376 if (target.GetFormat() == PixelFormat_Grayscale8 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	377 source.GetFormat() == PixelFormat_Grayscale16)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	378 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	379 ConvertInternal<uint8_t, uint16_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	380 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	381 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	382
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	383 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	384 source.GetFormat() == PixelFormat_Grayscale16)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	385 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	386 ConvertInternal<int16_t, uint16_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	387 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	388 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	389
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	390 if (target.GetFormat() == PixelFormat_Grayscale8 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	391 source.GetFormat() == PixelFormat_SignedGrayscale16)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	392 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	393 ConvertInternal<uint8_t, int16_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	394 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	395 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	396
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	397 if (target.GetFormat() == PixelFormat_Grayscale16 &&
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	398 source.GetFormat() == PixelFormat_SignedGrayscale16)
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	399 {
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	400 ConvertInternal<uint16_t, int16_t>(target, source);
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	401 return;
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	402 }
610a9a1ed855 ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 854 diff changeset	403
993 501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	404 if (target.GetFormat() == PixelFormat_Grayscale8 &&
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	405 source.GetFormat() == PixelFormat_RGB24)
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	406 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	407 ConvertColorToGrayscale<uint8_t>(target, source);
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	408 return;
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	409 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	410
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	411 if (target.GetFormat() == PixelFormat_Grayscale16 &&
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	412 source.GetFormat() == PixelFormat_RGB24)
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	413 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	414 ConvertColorToGrayscale<uint16_t>(target, source);
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	415 return;
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	416 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	417
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	418 if (target.GetFormat() == PixelFormat_SignedGrayscale16 &&
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	419 source.GetFormat() == PixelFormat_RGB24)
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	420 {
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	421 ConvertColorToGrayscale<int16_t>(target, source);
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	422 return;
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	423 }
501880d76474 improvements to GDCM plugin Sebastien Jodogne <s.jodogne@gmail.com> parents: 876 diff changeset	424
1993 e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	425 if (target.GetFormat() == PixelFormat_Float32 &&
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	426 source.GetFormat() == PixelFormat_Grayscale8)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	427 {
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	428 ConvertGrayscaleToFloat<uint8_t>(target, source);
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	429 return;
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	430 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	431
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	432 if (target.GetFormat() == PixelFormat_Float32 &&
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	433 source.GetFormat() == PixelFormat_Grayscale16)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	434 {
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	435 ConvertGrayscaleToFloat<uint16_t>(target, source);
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	436 return;
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	437 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	438
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	439 if (target.GetFormat() == PixelFormat_Float32 &&
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	440 source.GetFormat() == PixelFormat_Grayscale32)
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	441 {
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	442 ConvertGrayscaleToFloat<uint32_t>(target, source);
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	443 return;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	444 }
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	445
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	446 if (target.GetFormat() == PixelFormat_Float32 &&
1993 e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	447 source.GetFormat() == PixelFormat_SignedGrayscale16)
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	448 {
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	449 ConvertGrayscaleToFloat<int16_t>(target, source);
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	450 return;
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	451 }
e2a3ff770b48 introducing float32 images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1992 diff changeset	452
1610 2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	453 if (target.GetFormat() == PixelFormat_Grayscale8 &&
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	454 source.GetFormat() == PixelFormat_RGBA32)
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	455 {
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	456 for (unsigned int y = 0; y < source.GetHeight(); y++)
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	457 {
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	458 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	459 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	460 for (unsigned int x = 0; x < source.GetWidth(); x++, q++)
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	461 {
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	462 q = static_cast<uint8_t>((2126 static_cast<uint32_t>(p[0]) +
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	463 7152 * static_cast<uint32_t>(p[1]) +
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	464 0722 * static_cast<uint32_t>(p[2])) / 10000);
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	465 p += 4;
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	466 }
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	467 }
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	468
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	469 return;
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	470 }
2dff2bdffdb8 font support within Orthanc Sebastien Jodogne <s.jodogne@gmail.com> parents: 1608 diff changeset	471
1608 adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	472 if (target.GetFormat() == PixelFormat_RGB24 &&
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	473 source.GetFormat() == PixelFormat_RGBA32)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	474 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	475 for (unsigned int y = 0; y < source.GetHeight(); y++)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	476 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	477 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	478 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	479 for (unsigned int x = 0; x < source.GetWidth(); x++)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	480 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	481 q[0] = p[0];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	482 q[1] = p[1];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	483 q[2] = p[2];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	484 p += 4;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	485 q += 3;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	486 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	487 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	488
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	489 return;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	490 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	491
2252 002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	492 if (target.GetFormat() == PixelFormat_RGB24 &&
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	493 source.GetFormat() == PixelFormat_BGRA32)
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	494 {
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	495 for (unsigned int y = 0; y < source.GetHeight(); y++)
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	496 {
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	497 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	498 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	499 for (unsigned int x = 0; x < source.GetWidth(); x++)
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	500 {
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	501 q[0] = p[2];
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	502 q[1] = p[1];
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	503 q[2] = p[0];
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	504 p += 4;
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	505 q += 3;
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	506 }
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	507 }
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	508
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	509 return;
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	510 }
002b94046c69 colorspace conversion from BGRA32 to RGB24 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2244 diff changeset	511
1608 adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	512 if (target.GetFormat() == PixelFormat_RGBA32 &&
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	513 source.GetFormat() == PixelFormat_RGB24)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	514 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	515 for (unsigned int y = 0; y < source.GetHeight(); y++)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	516 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	517 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	518 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	519 for (unsigned int x = 0; x < source.GetWidth(); x++)
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	520 {
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	521 q[0] = p[0];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	522 q[1] = p[1];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	523 q[2] = p[2];
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	524 q[3] = 255; // Set the alpha channel to full opacity
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	525 p += 3;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	526 q += 4;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	527 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	528 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	529
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	530 return;
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	531 }
adc6a5704cdb OrthancPluginConvertPixelFormat Sebastien Jodogne <s.jodogne@gmail.com> parents: 1334 diff changeset	532
1992 9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	533 if (target.GetFormat() == PixelFormat_RGB24 &&
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	534 source.GetFormat() == PixelFormat_Grayscale8)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	535 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	536 for (unsigned int y = 0; y < source.GetHeight(); y++)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	537 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	538 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	539 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	540 for (unsigned int x = 0; x < source.GetWidth(); x++)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	541 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	542 q[0] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	543 q[1] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	544 q[2] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	545 p += 1;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	546 q += 3;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	547 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	548 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	549
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	550 return;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	551 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	552
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	553 if (target.GetFormat() == PixelFormat_RGBA32 &&
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	554 source.GetFormat() == PixelFormat_Grayscale8)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	555 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	556 for (unsigned int y = 0; y < source.GetHeight(); y++)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	557 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	558 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	559 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	560 for (unsigned int x = 0; x < source.GetWidth(); x++)
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	561 {
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	562 q[0] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	563 q[1] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	564 q[2] = *p;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	565 q[3] = 255;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	566 p += 1;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	567 q += 4;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	568 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	569 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	570
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	571 return;
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	572 }
9161e3ef0d17 new conversions in ImageProcessing::Convert Sebastien Jodogne <s.jodogne@gmail.com> parents: 1900 diff changeset	573
2100 1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	574 if (target.GetFormat() == PixelFormat_BGRA32 &&
2495 cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	575 source.GetFormat() == PixelFormat_Grayscale16)
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	576 {
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	577 for (unsigned int y = 0; y < source.GetHeight(); y++)
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	578 {
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	579 const uint16_t* p = reinterpret_cast<const uint16_t*>(source.GetConstRow(y));
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	580 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	581 for (unsigned int x = 0; x < source.GetWidth(); x++)
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	582 {
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	583 uint8_t value = (p < 256 ? p : 255);
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	584 q[0] = value;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	585 q[1] = value;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	586 q[2] = value;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	587 q[3] = 255;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	588 p += 1;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	589 q += 4;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	590 }
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	591 }
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	592
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	593 return;
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	594 }
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	595
cd7b854dbc05 convert grayscale16 to bgra32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2491 diff changeset	596 if (target.GetFormat() == PixelFormat_BGRA32 &&
2100 1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	597 source.GetFormat() == PixelFormat_RGB24)
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	598 {
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	599 for (unsigned int y = 0; y < source.GetHeight(); y++)
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	600 {
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	601 const uint8_t* p = reinterpret_cast<const uint8_t*>(source.GetConstRow(y));
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	602 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	603 for (unsigned int x = 0; x < source.GetWidth(); x++)
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	604 {
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	605 q[0] = p[2];
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	606 q[1] = p[1];
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	607 q[2] = p[0];
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	608 q[3] = 255;
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	609 p += 3;
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	610 q += 4;
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	611 }
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	612 }
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	613
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	614 return;
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	615 }
1554fc153a93 conversion RGB24 to BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2089 diff changeset	616
2423 5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	617 if (target.GetFormat() == PixelFormat_RGB24 &&
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	618 source.GetFormat() == PixelFormat_RGB48)
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	619 {
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	620 for (unsigned int y = 0; y < source.GetHeight(); y++)
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	621 {
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	622 const uint16_t* p = reinterpret_cast<const uint16_t*>(source.GetConstRow(y));
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	623 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	624 for (unsigned int x = 0; x < source.GetWidth(); x++)
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	625 {
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	626 q[0] = p[0] >> 8;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	627 q[1] = p[1] >> 8;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	628 q[2] = p[2] >> 8;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	629 p += 3;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	630 q += 3;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	631 }
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	632 }
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	633
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	634 return;
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	635 }
5a7c5c541a1d Built-in decoding of palette images Sebastien Jodogne <s.jodogne@gmail.com> parents: 2415 diff changeset	636
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	637 throw OrthancException(ErrorCode_NotImplemented);
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638 }
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	639
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	640
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	641
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	642 void ImageProcessing::Set(ImageAccessor& image,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	643 int64_t value)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	644 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	645 switch (image.GetFormat())
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	646 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	647 case PixelFormat_Grayscale8:
2487 be1dbb1dcdd6 PixelTraits and ImageTraits Sebastien Jodogne <s.jodogne@gmail.com> parents: 2482 diff changeset	648 memset(image.GetBuffer(), static_cast<uint8_t>(value), image.GetPitch() * image.GetHeight());
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	649 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	650
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	651 case PixelFormat_Grayscale16:
2490 0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	652 if (value == 0)
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	653 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	654 memset(image.GetBuffer(), 0, image.GetPitch() * image.GetHeight());
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	655 }
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	656 else
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	657 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	658 SetInternal<uint16_t>(image, value);
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	659 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	660 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	661
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	662 case PixelFormat_Grayscale32:
2490 0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	663 if (value == 0)
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	664 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	665 memset(image.GetBuffer(), 0, image.GetPitch() * image.GetHeight());
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	666 }
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	667 else
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	668 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	669 SetInternal<uint32_t>(image, value);
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	670 }
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	671 return;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	672
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	673 case PixelFormat_SignedGrayscale16:
2490 0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	674 if (value == 0)
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	675 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	676 memset(image.GetBuffer(), 0, image.GetPitch() * image.GetHeight());
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	677 }
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	678 else
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	679 {
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	680 SetInternal<int16_t>(image, value);
0cacb38297a9 speed up ImageProcessing::Set Sebastien Jodogne <s.jodogne@gmail.com> parents: 2489 diff changeset	681 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	682 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	683
1994 4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	684 case PixelFormat_Float32:
4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	685 assert(sizeof(float) == 4);
4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	686 SetInternal<float>(image, value);
4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	687 return;
4d099fee5eca ImageProcessing::Set for float images Sebastien Jodogne <s.jodogne@gmail.com> parents: 1993 diff changeset	688
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	689 default:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	690 throw OrthancException(ErrorCode_NotImplemented);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	691 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	692 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	693
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	694
2089 7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	695 void ImageProcessing::Set(ImageAccessor& image,
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	696 uint8_t red,
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	697 uint8_t green,
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	698 uint8_t blue,
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	699 uint8_t alpha)
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	700 {
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	701 uint8_t p[4];
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	702 unsigned int size;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	703
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	704 switch (image.GetFormat())
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	705 {
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	706 case PixelFormat_RGBA32:
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	707 p[0] = red;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	708 p[1] = green;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	709 p[2] = blue;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	710 p[3] = alpha;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	711 size = 4;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	712 break;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	713
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	714 case PixelFormat_BGRA32:
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	715 p[0] = blue;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	716 p[1] = green;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	717 p[2] = red;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	718 p[3] = alpha;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	719 size = 4;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	720 break;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	721
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	722 case PixelFormat_RGB24:
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	723 p[0] = red;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	724 p[1] = green;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	725 p[2] = blue;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	726 size = 3;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	727 break;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	728
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	729 default:
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	730 throw OrthancException(ErrorCode_NotImplemented);
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	731 }
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	732
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	733 for (unsigned int y = 0; y < image.GetHeight(); y++)
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	734 {
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	735 uint8_t* q = reinterpret_cast<uint8_t*>(image.GetRow(y));
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	736
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	737 for (unsigned int x = 0; x < image.GetWidth(); x++)
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	738 {
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	739 for (unsigned int i = 0; i < size; i++)
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	740 {
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	741 q[i] = p[i];
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	742 }
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	743
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	744 q += size;
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	745 }
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	746 }
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	747 }
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	748
7a969f235adf PixelFormat_BGRA32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1994 diff changeset	749
854 ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	750 void ImageProcessing::ShiftRight(ImageAccessor& image,
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	751 unsigned int shift)
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	752 {
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	753 if (image.GetWidth() == 0 \|\|
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	754 image.GetHeight() == 0 \|\|
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	755 shift == 0)
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	756 {
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	757 // Nothing to do
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	758 return;
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	759 }
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	760
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	761 throw OrthancException(ErrorCode_NotImplemented);
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	762 }
ff530685e46a fast version of image copy Sebastien Jodogne <s.jodogne@gmail.com> parents: 853 diff changeset	763
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	764
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	765 void ImageProcessing::GetMinMaxIntegerValue(int64_t& minValue,
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	766 int64_t& maxValue,
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	767 const ImageAccessor& image)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	768 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	769 switch (image.GetFormat())
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	770 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	771 case PixelFormat_Grayscale8:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	772 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	773 uint8_t a, b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	774 GetMinMaxValueInternal<uint8_t>(a, b, image);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	775 minValue = a;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	776 maxValue = b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	777 break;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	778 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	779
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	780 case PixelFormat_Grayscale16:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	781 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	782 uint16_t a, b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	783 GetMinMaxValueInternal<uint16_t>(a, b, image);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	784 minValue = a;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	785 maxValue = b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	786 break;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	787 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	788
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	789 case PixelFormat_Grayscale32:
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	790 {
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	791 uint32_t a, b;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	792 GetMinMaxValueInternal<uint32_t>(a, b, image);
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	793 minValue = a;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	794 maxValue = b;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	795 break;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	796 }
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	797
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	798 case PixelFormat_SignedGrayscale16:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	799 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	800 int16_t a, b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	801 GetMinMaxValueInternal<int16_t>(a, b, image);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	802 minValue = a;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	803 maxValue = b;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	804 break;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	805 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	806
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	807 default:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	808 throw OrthancException(ErrorCode_NotImplemented);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	809 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	810 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	811
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	812
2415 7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	813 void ImageProcessing::GetMinMaxFloatValue(float& minValue,
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	814 float& maxValue,
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	815 const ImageAccessor& image)
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	816 {
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	817 switch (image.GetFormat())
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	818 {
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	819 case PixelFormat_Float32:
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	820 {
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	821 assert(sizeof(float) == 32);
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	822 float a, b;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	823 GetMinMaxValueInternal<float>(a, b, image);
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	824 minValue = a;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	825 maxValue = b;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	826 break;
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	827 }
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	828
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	829 default:
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	830 throw OrthancException(ErrorCode_NotImplemented);
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	831 }
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	832 }
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	833
7e217a1cc63f PixelFormat_Float32 Sebastien Jodogne <s.jodogne@gmail.com> parents: 2281 diff changeset	834
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	835
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	836 void ImageProcessing::AddConstant(ImageAccessor& image,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	837 int64_t value)
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	838 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	839 switch (image.GetFormat())
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	840 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	841 case PixelFormat_Grayscale8:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	842 AddConstantInternal<uint8_t>(image, value);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	843 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	844
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	845 case PixelFormat_Grayscale16:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	846 AddConstantInternal<uint16_t>(image, value);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	847 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	848
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	849 case PixelFormat_SignedGrayscale16:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	850 AddConstantInternal<int16_t>(image, value);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	851 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	852
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	853 default:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	854 throw OrthancException(ErrorCode_NotImplemented);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	855 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	856 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	857
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	858
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	859 void ImageProcessing::MultiplyConstant(ImageAccessor& image,
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	860 float factor,
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	861 bool useRound)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	862 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	863 switch (image.GetFormat())
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	864 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	865 case PixelFormat_Grayscale8:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	866 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	867 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	868 MultiplyConstantInternal<uint8_t, true>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	869 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	870 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	871 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	872 MultiplyConstantInternal<uint8_t, false>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	873 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	874 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	875
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	876 case PixelFormat_Grayscale16:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	877 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	878 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	879 MultiplyConstantInternal<uint16_t, true>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	880 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	881 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	882 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	883 MultiplyConstantInternal<uint16_t, false>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	884 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	885 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	886
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	887 case PixelFormat_SignedGrayscale16:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	888 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	889 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	890 MultiplyConstantInternal<int16_t, true>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	891 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	892 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	893 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	894 MultiplyConstantInternal<int16_t, false>(image, factor);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	895 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	896 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	897
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	898 default:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	899 throw OrthancException(ErrorCode_NotImplemented);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	900 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	901 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	902
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	903
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	904 void ImageProcessing::ShiftScale(ImageAccessor& image,
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	905 float offset,
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	906 float scaling,
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	907 bool useRound)
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	908 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	909 switch (image.GetFormat())
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	910 {
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	911 case PixelFormat_Grayscale8:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	912 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	913 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	914 ShiftScaleInternal<uint8_t, true>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	915 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	916 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	917 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	918 ShiftScaleInternal<uint8_t, false>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	919 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	920 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	921
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	922 case PixelFormat_Grayscale16:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	923 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	924 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	925 ShiftScaleInternal<uint16_t, true>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	926 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	927 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	928 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	929 ShiftScaleInternal<uint16_t, false>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	930 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	931 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	932
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	933 case PixelFormat_SignedGrayscale16:
2488 345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	934 if (useRound)
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	935 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	936 ShiftScaleInternal<int16_t, true>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	937 }
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	938 else
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	939 {
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	940 ShiftScaleInternal<int16_t, false>(image, offset, scaling);
345725b9350c back to rounding to fix integration tests Sebastien Jodogne <s.jodogne@gmail.com> parents: 2487 diff changeset	941 }
863 3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	942 return;
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	943
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	944 default:
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	945 throw OrthancException(ErrorCode_NotImplemented);
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	946 }
3c0d0836f704 refactoring Sebastien Jodogne <s.jodogne@gmail.com> parents: 862 diff changeset	947 }
2281 e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	948
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	949
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	950 void ImageProcessing::Invert(ImageAccessor& image)
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	951 {
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	952 switch (image.GetFormat())
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	953 {
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	954 case PixelFormat_Grayscale8:
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	955 {
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	956 for (unsigned int y = 0; y < image.GetHeight(); y++)
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	957 {
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	958 uint8_t* p = reinterpret_cast<uint8_t*>(image.GetRow(y));
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	959
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	960 for (unsigned int x = 0; x < image.GetWidth(); x++, p++)
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	961 {
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	962 p = 255 - (p);
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	963 }
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	964 }
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	965
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	966 return;
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	967 }
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	968
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	969 default:
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	970 throw OrthancException(ErrorCode_NotImplemented);
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	971 }
e002430baa41 Fix issue #44 (Bad interpretation of photometric interpretation MONOCHROME1) Sebastien Jodogne <s.jodogne@gmail.com> parents: 2252 diff changeset	972 }
2489 e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	973
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	974
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	975
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	976 namespace
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	977 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	978 template <Orthanc::PixelFormat Format>
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	979 class BresenhamPixelWriter
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	980 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	981 private:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	982 typedef typename PixelTraits<Format>::PixelType PixelType;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	983
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	984 Orthanc::ImageAccessor& image_;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	985 PixelType value_;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	986
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	987 void PlotLineLow(int x0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	988 int y0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	989 int x1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	990 int y1)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	991 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	992 int dx = x1 - x0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	993 int dy = y1 - y0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	994 int yi = 1;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	995
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	996 if (dy < 0)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	997 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	998 yi = -1;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	999 dy = -dy;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1000 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1001
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1002 int d = 2 * dy - dx;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1003 int y = y0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1004
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1005 for (int x = x0; x <= x1; x++)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1006 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1007 Write(x, y);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1008
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1009 if (d > 0)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1010 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1011 y = y + yi;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1012 d = d - 2 * dx;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1013 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1014
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1015 d = d + 2*dy;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1016 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1017 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1018
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1019 void PlotLineHigh(int x0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1020 int y0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1021 int x1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1022 int y1)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1023 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1024 int dx = x1 - x0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1025 int dy = y1 - y0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1026 int xi = 1;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1027
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1028 if (dx < 0)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1029 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1030 xi = -1;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1031 dx = -dx;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1032 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1033
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1034 int d = 2 * dx - dy;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1035 int x = x0;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1036
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1037 for (int y = y0; y <= y1; y++)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1038 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1039 Write(x, y);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1040
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1041 if (d > 0)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1042 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1043 x = x + xi;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1044 d = d - 2 * dy;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1045 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1046
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1047 d = d + 2 * dx;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1048 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1049 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1050
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1051 public:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1052 BresenhamPixelWriter(Orthanc::ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1053 int64_t value) :
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1054 image_(image),
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1055 value_(PixelTraits<Format>::IntegerToPixel(value))
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1056 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1057 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1058
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1059 BresenhamPixelWriter(Orthanc::ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1060 const PixelType& value) :
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1061 image_(image),
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1062 value_(value)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1063 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1064 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1065
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1066 void Write(int x,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1067 int y)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1068 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1069 if (x >= 0 &&
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1070 y >= 0 &&
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1071 static_cast<unsigned int>(x) < image_.GetWidth() &&
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1072 static_cast<unsigned int>(y) < image_.GetHeight())
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1073 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1074 PixelType* p = reinterpret_cast<PixelType*>(image_.GetRow(y));
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1075 p[x] = value_;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1076 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1077 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1078
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1079 void DrawSegment(int x0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1080 int y0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1081 int x1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1082 int y1)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1083 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1084 // This is an implementation of Bresenham's line algorithm
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1085 // https://en.wikipedia.org/wiki/Bresenham%27s_line_algorithm#All_cases
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1086
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1087 if (abs(y1 - y0) < abs(x1 - x0))
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1088 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1089 if (x0 > x1)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1090 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1091 PlotLineLow(x1, y1, x0, y0);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1092 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1093 else
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1094 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1095 PlotLineLow(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1096 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1097 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1098 else
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1099 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1100 if (y0 > y1)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1101 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1102 PlotLineHigh(x1, y1, x0, y0);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1103 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1104 else
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1105 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1106 PlotLineHigh(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1107 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1108 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1109 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1110 };
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1111 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1112
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1113
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1114 void ImageProcessing::DrawLineSegment(ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1115 int x0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1116 int y0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1117 int x1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1118 int y1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1119 int64_t value)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1120 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1121 switch (image.GetFormat())
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1122 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1123 case Orthanc::PixelFormat_Grayscale8:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1124 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1125 BresenhamPixelWriter<Orthanc::PixelFormat_Grayscale8> writer(image, value);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1126 writer.DrawSegment(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1127 break;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1128 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1129
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1130 case Orthanc::PixelFormat_Grayscale16:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1131 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1132 BresenhamPixelWriter<Orthanc::PixelFormat_Grayscale16> writer(image, value);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1133 writer.DrawSegment(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1134 break;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1135 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1136
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1137 case Orthanc::PixelFormat_SignedGrayscale16:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1138 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1139 BresenhamPixelWriter<Orthanc::PixelFormat_SignedGrayscale16> writer(image, value);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1140 writer.DrawSegment(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1141 break;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1142 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1143
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1144 default:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1145 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1146 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1147 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1148
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1149
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1150 void ImageProcessing::DrawLineSegment(ImageAccessor& image,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1151 int x0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1152 int y0,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1153 int x1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1154 int y1,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1155 uint8_t red,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1156 uint8_t green,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1157 uint8_t blue,
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1158 uint8_t alpha)
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1159 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1160 switch (image.GetFormat())
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1161 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1162 case Orthanc::PixelFormat_BGRA32:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1163 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1164 PixelTraits<Orthanc::PixelFormat_BGRA32>::PixelType pixel;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1165 pixel.red_ = red;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1166 pixel.green_ = green;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1167 pixel.blue_ = blue;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1168 pixel.alpha_ = alpha;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1169
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1170 BresenhamPixelWriter<Orthanc::PixelFormat_BGRA32> writer(image, pixel);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1171 writer.DrawSegment(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1172 break;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1173 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1174
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1175 case Orthanc::PixelFormat_RGB24:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1176 {
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1177 PixelTraits<Orthanc::PixelFormat_RGB24>::PixelType pixel;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1178 pixel.red_ = red;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1179 pixel.green_ = green;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1180 pixel.blue_ = blue;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1181
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1182 BresenhamPixelWriter<Orthanc::PixelFormat_RGB24> writer(image, pixel);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1183 writer.DrawSegment(x0, y0, x1, y1);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1184 break;
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1185 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1186
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1187 default:
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1188 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1189 }
e91bab2d8c75 Bresenham's line algorithm Sebastien Jodogne <s.jodogne@gmail.com> parents: 2488 diff changeset	1190 }
853 839be3022203 DicomImageInformation Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1191 }

Mercurial > hg > orthanc

annotate Core/Images/ImageProcessing.cpp @ 2495:cd7b854dbc05