orthanc-stone: Framework/Volumes/VolumeReslicer.cpp annotate

annotate Framework/Volumes/VolumeReslicer.cpp @ 186:fac46066b933 wasm

missing headers

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Thu, 15 Mar 2018 16:11:45 +0100
parents	d61224de4883
children	a12ca0f4aaaa

rev	line source
186 fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	1 /**
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	2 * Stone of Orthanc
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	4 * Department, University Hospital of Liege, Belgium
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	5 * Copyright (C) 2017-2018 Osimis S.A., Belgium
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	6 *
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	7 * This program is free software: you can redistribute it and/or
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	8 * modify it under the terms of the GNU Affero General Public License
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	9 * as published by the Free Software Foundation, either version 3 of
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	10 * the License, or (at your option) any later version.
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	11 *
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	12 * This program is distributed in the hope that it will be useful, but
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	13 * WITHOUT ANY WARRANTY; without even the implied warranty of
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	15 * Affero General Public License for more details.
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	16 *
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	17 * You should have received a copy of the GNU Affero General Public License
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	18 * along with this program. If not, see <http://www.gnu.org/licenses/>.
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	19 **/
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	20
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	21
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 #include "VolumeReslicer.h"
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23
159 0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	24 #include "../Toolbox/GeometryToolbox.h"
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	25 #include "../Toolbox/SubvoxelReader.h"
159 0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	26
179 db21c1810c89 moving PixelTraits to the Orthanc core Sebastien Jodogne <s.jodogne@gmail.com> parents: 177 diff changeset	27 #include <Core/Images/ImageTraits.h>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 #include <Core/Logging.h>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29 #include <Core/OrthancException.h>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31 #include <boost/math/special_functions/round.hpp>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	33 #include <stdio.h>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35 namespace OrthancStone
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 // Anonymous namespace to avoid clashes between compilation modules
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 namespace
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40 enum TransferFunction
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 TransferFunction_Copy,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 TransferFunction_Float,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 TransferFunction_Linear
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	45 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	48 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	49 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	50 ImageInterpolation Interpolation,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52 class PixelShader;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	54
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	55 template <Orthanc::PixelFormat Format>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	56 class PixelShader<Format,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	57 Format,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	58 ImageInterpolation_Nearest,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	59 TransferFunction_Copy>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	62 typedef SubvoxelReader<Format, ImageInterpolation_Nearest> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	63 typedef Orthanc::PixelTraits<Format> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	64
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	65 VoxelReader reader_;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	66
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	67 public:
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	68 PixelShader(const ImageBuffer3D& image,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	69 float /* scaling */,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	70 float /* offset */) :
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	71 reader_(image)
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	72 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	73 }
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	74
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	75 ORTHANC_FORCE_INLINE
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	76 void Apply(typename PixelWriter::PixelType* pixel,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	77 float volumeX,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	78 float volumeY,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	79 float volumeZ)
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	80 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	81 typename VoxelReader::PixelType value;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	82
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	83 if (!reader_.GetValue(value, volumeX, volumeY, volumeZ))
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	84 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	85 VoxelReader::Traits::SetMinValue(value);
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	86 }
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	87
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	88 *pixel = value;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	89 }
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	90 };
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	91
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	92
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	93 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	94 Orthanc::PixelFormat OutputFormat>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	95 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	96 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	97 ImageInterpolation_Nearest,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	98 TransferFunction_Copy>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	99 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	100 private:
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	101 typedef SubvoxelReader<InputFormat, ImageInterpolation_Nearest> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	102 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	103
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 float /* offset */) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110 reader_(image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	111 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	112 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	114 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	115 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	116 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	117 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	118 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	120 typename VoxelReader::PixelType value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	121
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	122 if (!reader_.GetValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	123 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	124 VoxelReader::Traits::SetMinValue(value);
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	125 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	126
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	127 PixelWriter::FloatToPixel(*pixel, VoxelReader::Traits::PixelToFloat(value));
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	129 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	130
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	132 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	133 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	134 ImageInterpolation Interpolation>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	135 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	136 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	137 Interpolation,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	138 TransferFunction_Float>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	141 typedef SubvoxelReader<InputFormat, Interpolation> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	142 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	143
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	144 VoxelReader reader_;
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	145 float outOfVolume_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	150 float /* offset */) :
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	151 reader_(image),
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	152 outOfVolume_(static_cast<float>(std::numeric_limits<typename VoxelReader::PixelType>::min()))
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	153 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	154 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	155
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	156 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	157 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	158 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	159 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	160 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	162 float value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	163
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	164 if (!reader_.GetFloatValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	165 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	166 value = outOfVolume_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	167 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	168
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	169 PixelWriter::FloatToPixel(*pixel, value);
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	170 }
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	174 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	175 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	176 ImageInterpolation Interpolation>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	177 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	178 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	179 Interpolation,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	180 TransferFunction_Linear>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	183 typedef SubvoxelReader<InputFormat, Interpolation> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	184 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	185
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 float scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	188 float offset_;
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	189 float outOfVolume_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 float offset) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 reader_(image),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196 scaling_(scaling),
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	197 offset_(offset),
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	198 outOfVolume_(static_cast<float>(std::numeric_limits<typename VoxelReader::PixelType>::min()))
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	202 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	203 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	204 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	205 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	206 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	208 float value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	209
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	210 if (reader_.GetFloatValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	211 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	212 value = scaling_ * value + offset_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	213 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	214 else
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	215 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	216 value = outOfVolume_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	217 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	218
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	219 PixelWriter::FloatToPixel(*pixel, value);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	223
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	224
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 class FastRowIterator : public boost::noncopyable
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228 float position_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229 float offset_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	230
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	231 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 FastRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236 unsigned int y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 const double width = static_cast<double>(slice.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239 const double height = static_cast<double>(slice.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240 assert(y < height);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 Vector q1 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	243 (extent.GetX1() + extent.GetWidth() * static_cast<double>(0) / static_cast<double>(width + 1),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	244 extent.GetY1() + extent.GetHeight() * static_cast<double>(y) / static_cast<double>(height + 1));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 Vector q2 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247 (extent.GetX1() + extent.GetWidth() * static_cast<double>(width - 1) / static_cast<double>(width + 1),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248 extent.GetY1() + extent.GetHeight() * static_cast<double>(y) / static_cast<double>(height + 1));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 Vector r1, r2;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 box.ToInternalCoordinates(r1, q1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252 box.ToInternalCoordinates(r2, q2);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 position_[0] = static_cast<float>(r1[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255 position_[1] = static_cast<float>(r1[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256 position_[2] = static_cast<float>(r1[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258 Vector tmp = (r2 - r1) / static_cast<double>(width - 1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259 offset_[0] = static_cast<float>(tmp[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260 offset_[1] = static_cast<float>(tmp[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 offset_[2] = static_cast<float>(tmp[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	264 ORTHANC_FORCE_INLINE
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265 void Next()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267 position_[0] += offset_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268 position_[1] += offset_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 position_[2] += offset_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	272 ORTHANC_FORCE_INLINE
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	273 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	274 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	275 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277 x = position_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 y = position_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279 z = position_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	280 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	281 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	282
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	283
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 class SlowRowIterator : public boost::noncopyable
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 const Orthanc::ImageAccessor& slice_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288 const Extent2D& extent_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289 const CoordinateSystem3D& plane_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	290 const OrientedBoundingBox& box_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291 unsigned int x_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292 unsigned int y_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 SlowRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	297 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299 unsigned int y) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 slice_(slice),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 extent_(extent),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 plane_(plane),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 box_(box),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304 x_(0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 y_(y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 assert(y_ < slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	308 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	309
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 void Next()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	311 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 x_++;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	314
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	315 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	316 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	317 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319 assert(x_ < slice_.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 const double width = static_cast<double>(slice_.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322 const double height = static_cast<double>(slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 Vector q = plane_.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 (extent_.GetX1() + extent_.GetWidth() * static_cast<double>(x_) / (width + 1.0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 extent_.GetY1() + extent_.GetHeight() * static_cast<double>(y_) / (height + 1.0));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 Vector r;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 box_.ToInternalCoordinates(r, q);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 x = static_cast<float>(r[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 y = static_cast<float>(r[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 z = static_cast<float>(r[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	336
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	337
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 Orthanc::PixelFormat OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 ImageInterpolation Interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350 {
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	351 typedef PixelShader<InputFormat, OutputFormat, Interpolation, Function> Shader;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 const unsigned int outputWidth = slice.GetWidth();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354 const unsigned int outputHeight = slice.GetHeight();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	355
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	356 const float sourceWidth = static_cast<float>(source.GetWidth());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	357 const float sourceHeight = static_cast<float>(source.GetHeight());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	358 const float sourceDepth = static_cast<float>(source.GetDepth());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	359
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360 Shader shader(source, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 for (unsigned int y = 0; y < outputHeight; y++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 {
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	364 typedef typename Orthanc::ImageTraits<OutputFormat>::PixelType PixelType;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	365 PixelType* p = reinterpret_cast<PixelType*>(slice.GetRow(y));
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367 RowIterator it(slice, extent, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 for (unsigned int x = 0; x < outputWidth; x++, p++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	371 float volumeX, volumeY, volumeZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	372 it.GetVolumeCoordinates(volumeX, volumeY, volumeZ);
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	373
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	374 shader.Apply(p,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	375 volumeX * sourceWidth,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	376 volumeY * sourceHeight,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	377 volumeZ * sourceDepth);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 it.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	379 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	380 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	381 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	382
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	383
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 Orthanc::PixelFormat OutputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 if (hasLinearFunction)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399 switch (interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 ImageInterpolation_Nearest, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409 ImageInterpolation_Bilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 ImageInterpolation_Trilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	418
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	419 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	421 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	422 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 switch (interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 ImageInterpolation_Nearest, TransferFunction_Copy>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435 ImageInterpolation_Bilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 ImageInterpolation_Trilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	447 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	448 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	449 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	450
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	451
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 template <typename RowIterator>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464 slice.GetFormat() == Orthanc::PixelFormat_Grayscale8)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	466 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468 Orthanc::PixelFormat_Grayscale8>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	469 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	470 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	471 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	472 slice.GetFormat() == Orthanc::PixelFormat_Grayscale16)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	476 Orthanc::PixelFormat_Grayscale16>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 else if (source.GetFormat() == Orthanc::PixelFormat_SignedGrayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 Orthanc::PixelFormat_SignedGrayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	496 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	497 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	501
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	502
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	503
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	504 void VolumeReslicer::CheckIterators(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 const OrientedBoundingBox& box) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	507 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	508 for (unsigned int y = 0; y < slice_->GetHeight(); y++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	509 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	510 FastRowIterator fast(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	511 SlowRowIterator slow(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	512
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	513 for (unsigned int x = 0; x < slice_->GetWidth(); x++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	515 float px, py, pz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	516 fast.GetVolumeCoordinates(px, py, pz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	518 float qx, qy, qz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	519 slow.GetVolumeCoordinates(qx, qy, qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	520
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	521 Vector d;
158 a053ca7fa5c6 LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 154 diff changeset	522 LinearAlgebra::AssignVector(d, px - qx, py - qy, pz - qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	523 double norm = boost::numeric::ublas::norm_2(d);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 if (norm > 0.0001)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	525 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	526 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	527 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	528
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	529 fast.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530 slow.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	532 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	536 void VolumeReslicer::Reset()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	537 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	538 success_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 extent_.Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 slice_.reset(NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	541 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	542
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	543
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	544 float VolumeReslicer::GetMinOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	545 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	546 switch (outputFormat_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	547 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	549 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 return 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	553
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	554 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	555 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	556 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	557 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	558
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 float VolumeReslicer::GetMaxOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 switch (outputFormat_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	564 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	566 return static_cast<float>(std::numeric_limits<uint8_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	568
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	569 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	570 return static_cast<float>(std::numeric_limits<uint16_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	571 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	572
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	573 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	574 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	575 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	576 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	577
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	578
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	579 VolumeReslicer::VolumeReslicer() :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 outputFormat_(Orthanc::PixelFormat_Grayscale8),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 interpolation_(ImageInterpolation_Nearest),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582 fastMode_(true),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 success_(false)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	584 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585 ResetLinearFunction();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	586 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	587
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	588
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	589 void VolumeReslicer::GetLinearFunction(float& scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590 float& offset) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 if (hasLinearFunction_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	593 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 scaling = scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 offset = offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	596 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	597 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	598 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	599 scaling = 1.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	600 offset = 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	601 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	602 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	603
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	604
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	605 void VolumeReslicer::ResetLinearFunction()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	606 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	607 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	608 hasLinearFunction_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	609 scaling_ = 1.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	610 offset_ = 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	611 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	612
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	613
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	614 void VolumeReslicer::SetLinearFunction(float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	615 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	616 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	617 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	618 hasLinearFunction_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 scaling_ = scaling;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 offset_ = offset;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	621 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	622
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	623
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	624 void VolumeReslicer::SetWindow(float low,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	625 float high)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	626 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	627 //printf("Range in pixel values: %f->%f\n", low, high);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 float scaling = (GetMaxOutputValue() - GetMinOutputValue()) / (high - low);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 float offset = GetMinOutputValue() - scaling * low;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 SetLinearFunction(scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	632
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	633 /float x = scaling_ low + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	634 float y = scaling_ * high + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	635 printf("%f %f (should be %f->%f)\n", x, y, GetMinOutputValue(), GetMaxOutputValue());*/
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	636 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	637
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 void VolumeReslicer::FitRange(const ImageBuffer3D& image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	640 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	641 float minInputValue, maxInputValue;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	642
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	643 if (!image.GetRange(minInputValue, maxInputValue) \|\|
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	644 maxInputValue < 1)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	646 ResetLinearFunction();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	647 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	648 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	649 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650 SetWindow(minInputValue, maxInputValue);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	652 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	653
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	654
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	655 void VolumeReslicer::SetWindowing(ImageWindowing windowing,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	656 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657 float rescaleSlope,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 float rescaleIntercept)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	659 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	660 if (windowing == ImageWindowing_Custom \|\|
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 windowing == ImageWindowing_Default)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	662 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	663 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	664 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	665
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	666 float center, width;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	667 ComputeWindowing(center, width, windowing, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	668
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	669 float a = (center - width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	670 float b = (center + width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	671 SetWindow(a, b);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	672 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	673
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	674
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	675 void VolumeReslicer::SetOutputFormat(Orthanc::PixelFormat format)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	676 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	677 if (format != Orthanc::PixelFormat_Grayscale8 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	678 format != Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	679 format != Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	682 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	683
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	684 if (hasLinearFunction_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	685 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	686 LOG(WARNING) << "Calls to VolumeReslicer::SetOutputFormat() should be done before VolumeReslicer::FitRange()";
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689 outputFormat_ = format;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	691 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	692
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	693
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 void VolumeReslicer::SetInterpolation(ImageInterpolation interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	696 if (interpolation != ImageInterpolation_Nearest &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	697 interpolation != ImageInterpolation_Bilinear &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	698 interpolation != ImageInterpolation_Trilinear)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	701 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	702
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	703 interpolation_ = interpolation;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708 const Extent2D& VolumeReslicer::GetOutputExtent() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	709 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	711 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	712 return extent_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	713 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	715 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	716 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	717 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	718 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	719
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	720
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	721 const Orthanc::ImageAccessor& VolumeReslicer::GetOutputSlice() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	722 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	723 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	725 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	726 return *slice_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	729 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	730 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	731 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	732 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	733
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	734
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735 Orthanc::ImageAccessor* VolumeReslicer::ReleaseOutputSlice()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	737 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	738 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	739 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740 success_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	741 return slice_.release();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	744 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	745 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	746 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	747 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	748
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	749
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	750 void VolumeReslicer::Apply(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	751 const CoordinateSystem3D& plane)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753 // Choose the default voxel size as the finest voxel dimension
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754 // of the source volumetric image
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755 const OrthancStone::Vector dim = source.GetVoxelDimensions(OrthancStone::VolumeProjection_Axial);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756 double voxelSize = dim[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	757
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 if (dim[1] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	759 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	760 voxelSize = dim[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	761 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	762
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	763 if (dim[2] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	764 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765 voxelSize = dim[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	766 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	767
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 if (voxelSize <= 0)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	769 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	770 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	771 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	772
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	773 Apply(source, plane, voxelSize);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	774 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	775
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	776
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	777 void VolumeReslicer::Apply(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779 double voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	780 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	781 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783 // Firstly, compute the intersection of the source volumetric
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 // image with the reslicing plane. This leads to a polygon with 3
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 // to 6 vertices. We compute the extent of the intersection
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	786 // polygon, with respect to the coordinate system of the reslicing
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	787 // plane.
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	788 OrientedBoundingBox box(source);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	789
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	790 if (!box.ComputeExtent(extent_, plane))
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	792 // The plane does not intersect with the bounding box of the volume
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793 slice_.reset(new Orthanc::Image(outputFormat_, 0, 0, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795 return;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	797
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798 // Secondly, the extent together with the voxel size gives the
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 // size of the output image
172 316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	800 unsigned int width = boost::math::iround(extent_.GetWidth() / voxelSize);
316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	801 unsigned int height = boost::math::iround(extent_.GetHeight() / voxelSize);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 slice_.reset(new Orthanc::Image(outputFormat_, width, height, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	804
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805 //CheckIterators(source, plane, box);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	806
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	807 if (fastMode_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	809 ProcessImage<FastRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	812 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	813 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	814 ProcessImage<SlowRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	816 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	817
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	818 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820 }

Mercurial > hg > orthanc-stone

annotate Framework/Volumes/VolumeReslicer.cpp @ 186:fac46066b933 wasm