orthanc-stone: OrthancStone/Sources/Volumes/VolumeReslicer.cpp annotate

annotate OrthancStone/Sources/Volumes/VolumeReslicer.cpp @ 1772:6c246f862b00

unit test VolumeRendering.Basic

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Wed, 12 May 2021 17:24:24 +0200
parents	9ac2a65d4172
children	95ece40bb298

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
1739 9ac2a65d4172 upgrade to year 2021 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1640 diff changeset	5 * Copyright (C) 2017-2021 Osimis S.A., Belgium
186 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
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	8 * modify it under the terms of the GNU Lesser General Public License
186 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
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	15 * Lesser General Public License for more details.
1596 4fb8fdf03314 removed annoying whitespace Sebastien Jodogne <s.jodogne@gmail.com> parents: 1512 diff changeset	16 *
1598 8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	17 * You should have received a copy of the GNU Lesser General Public
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	18 * License along with this program. If not, see
8563ea5d8ae4 relicensing some files, cf. osimis bm26 and chu agreement on 2020-05-20 Sebastien Jodogne <s.jodogne@gmail.com> parents: 1596 diff changeset	19 * <http://www.gnu.org/licenses/>.
186 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
fac46066b933 missing headers Sebastien Jodogne <s.jodogne@gmail.com> parents: 185 diff changeset	22
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 #include "VolumeReslicer.h"
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24
159 0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	25 #include "../Toolbox/GeometryToolbox.h"
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	26 #include "../Toolbox/SubvoxelReader.h"
159 0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	27
1455 30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1270 diff changeset	28 #include <Images/ImageTraits.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1270 diff changeset	29 #include <Logging.h>
30deba7bc8e2 simplifying include_directories Sebastien Jodogne <s.jodogne@gmail.com> parents: 1270 diff changeset	30 #include <OrthancException.h>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32 #include <boost/math/special_functions/round.hpp>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	34 namespace OrthancStone
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	35 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36 // Anonymous namespace to avoid clashes between compilation modules
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 namespace
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	39 enum TransferFunction
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	40 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	41 TransferFunction_Copy,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	42 TransferFunction_Float,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 TransferFunction_Linear
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44 };
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
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	47 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	48 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	49 ImageInterpolation Interpolation,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 class PixelShader;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	53
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	54 template <Orthanc::PixelFormat Format>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	55 class PixelShader<Format,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	56 Format,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	57 ImageInterpolation_Nearest,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	58 TransferFunction_Copy>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	61 typedef SubvoxelReader<Format, ImageInterpolation_Nearest> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	62 typedef Orthanc::PixelTraits<Format> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	63
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	64 VoxelReader reader_;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	65
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	66 public:
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	67 PixelShader(const ImageBuffer3D& image,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	68 float /* scaling */,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	69 float /* offset */) :
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	70 reader_(image)
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	71 {
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 ORTHANC_FORCE_INLINE
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	75 void Apply(typename PixelWriter::PixelType* pixel,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	76 float volumeX,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	77 float volumeY,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	78 float volumeZ)
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	79 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	80 typename VoxelReader::PixelType value;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	81
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	82 if (!reader_.GetValue(value, volumeX, volumeY, volumeZ))
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	83 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	84 VoxelReader::Traits::SetMinValue(value);
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	85 }
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 *pixel = value;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	88 }
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 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	93 Orthanc::PixelFormat OutputFormat>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	94 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	95 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	96 ImageInterpolation_Nearest,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	97 TransferFunction_Copy>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	98 {
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	99 private:
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	100 typedef SubvoxelReader<InputFormat, ImageInterpolation_Nearest> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	101 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	102
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 float /* offset */) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 reader_(image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110 {
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
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	113 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	114 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	115 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	116 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	117 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	119 typename VoxelReader::PixelType value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	120
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	121 if (!reader_.GetValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	122 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	123 VoxelReader::Traits::SetMinValue(value);
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	124 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	125
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	126 PixelWriter::FloatToPixel(*pixel, VoxelReader::Traits::PixelToFloat(value));
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 }
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
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	131 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	132 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	133 ImageInterpolation Interpolation>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	134 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	135 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	136 Interpolation,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	137 TransferFunction_Float>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	140 typedef SubvoxelReader<InputFormat, Interpolation> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	141 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	142
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	143 VoxelReader reader_;
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	144 float outOfVolume_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149 float /* offset */) :
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	150 reader_(image),
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	151 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	152 {
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
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	155 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	156 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	157 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	158 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	159 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	160 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	161 float value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	162
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	163 if (!reader_.GetFloatValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	164 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	165 value = outOfVolume_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	166 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	167
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	168 PixelWriter::FloatToPixel(*pixel, value);
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	169 }
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170 };
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
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	173 template <Orthanc::PixelFormat InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	174 Orthanc::PixelFormat OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	175 ImageInterpolation Interpolation>
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	176 class PixelShader<InputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	177 OutputFormat,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	178 Interpolation,
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	179 TransferFunction_Linear>
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	180 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 private:
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	182 typedef SubvoxelReader<InputFormat, Interpolation> VoxelReader;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	183 typedef Orthanc::PixelTraits<OutputFormat> PixelWriter;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	184
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 float scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 float offset_;
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	188 float outOfVolume_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	189
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193 float offset) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 reader_(image),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 scaling_(scaling),
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	196 offset_(offset),
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	197 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	198 {
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
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	201 ORTHANC_FORCE_INLINE
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	202 void Apply(typename PixelWriter::PixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	203 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	204 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	205 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 {
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	207 float value;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	208
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	209 if (reader_.GetFloatValue(value, volumeX, volumeY, volumeZ))
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	210 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	211 value = scaling_ * value + offset_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	212 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	213 else
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	214 {
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	215 value = outOfVolume_;
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	216 }
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	217
184 9523ce4f44cc removing PixelWriter Sebastien Jodogne <s.jodogne@gmail.com> parents: 183 diff changeset	218 PixelWriter::FloatToPixel(*pixel, value);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 }
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 class FastRowIterator : public boost::noncopyable
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227 float position_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228 float offset_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	230 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	231 FastRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	234 const OrientedVolumeBoundingBox& box,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235 unsigned int y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	237 const double width = static_cast<double>(slice.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 const double height = static_cast<double>(slice.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	239 assert(y < height);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	240
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	241 Vector q1 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	242 (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	243 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	244
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	245 Vector q2 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	246 (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	247 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	248
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	249 Vector r1, r2;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	250 box.ToInternalCoordinates(r1, q1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	251 box.ToInternalCoordinates(r2, q2);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	252
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	253 position_[0] = static_cast<float>(r1[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	254 position_[1] = static_cast<float>(r1[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	255 position_[2] = static_cast<float>(r1[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	256
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	257 Vector tmp = (r2 - r1) / static_cast<double>(width - 1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	258 offset_[0] = static_cast<float>(tmp[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259 offset_[1] = static_cast<float>(tmp[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	260 offset_[2] = static_cast<float>(tmp[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	262
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	263 ORTHANC_FORCE_INLINE
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	264 void Next()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	265 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 position_[0] += offset_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	267 position_[1] += offset_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268 position_[2] += offset_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	270
180 4da803580da9 remove macro ORTHANC_STONE_FORCE_INLINE Sebastien Jodogne <s.jodogne@gmail.com> parents: 179 diff changeset	271 ORTHANC_FORCE_INLINE
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	272 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	273 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	274 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	275 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	276 x = position_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	277 y = position_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 z = position_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279 }
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 class SlowRowIterator : public boost::noncopyable
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	284 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	285 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	286 const Orthanc::ImageAccessor& slice_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	287 const Extent2D& extent_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	288 const CoordinateSystem3D& plane_;
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	289 const OrientedVolumeBoundingBox& box_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	290 unsigned int x_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	291 unsigned int y_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	292
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	294 SlowRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	295 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	296 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	297 const OrientedVolumeBoundingBox& box,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	298 unsigned int y) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	299 slice_(slice),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	300 extent_(extent),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301 plane_(plane),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	302 box_(box),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	303 x_(0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	304 y_(y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	305 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	306 assert(y_ < slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	307 }
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 void Next()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	310 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	311 x_++;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	312 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	313
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	314 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	315 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	316 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 assert(x_ < slice_.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320 const double width = static_cast<double>(slice_.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 const double height = static_cast<double>(slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323 Vector q = plane_.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 (extent_.GetX1() + extent_.GetWidth() * static_cast<double>(x_) / (width + 1.0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 extent_.GetY1() + extent_.GetHeight() * static_cast<double>(y_) / (height + 1.0));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	327 Vector r;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	328 box_.ToInternalCoordinates(r, q);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 x = static_cast<float>(r[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 y = static_cast<float>(r[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 z = static_cast<float>(r[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 }
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 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339 Orthanc::PixelFormat OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340 ImageInterpolation Interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	342 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	343 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	344 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	346 const OrientedVolumeBoundingBox& box,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	347 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 {
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	350 typedef PixelShader<InputFormat, OutputFormat, Interpolation, Function> Shader;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 const unsigned int outputWidth = slice.GetWidth();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	353 const unsigned int outputHeight = slice.GetHeight();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	354
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	355 const float sourceWidth = static_cast<float>(source.GetWidth());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	356 const float sourceHeight = static_cast<float>(source.GetHeight());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	357 const float sourceDepth = static_cast<float>(source.GetDepth());
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	358
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359 Shader shader(source, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	360
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	361 for (unsigned int y = 0; y < outputHeight; y++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	362 {
185 d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	363 typedef typename Orthanc::ImageTraits<OutputFormat>::PixelType PixelType;
d61224de4883 simplification Sebastien Jodogne <s.jodogne@gmail.com> parents: 184 diff changeset	364 PixelType* p = reinterpret_cast<PixelType*>(slice.GetRow(y));
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 RowIterator it(slice, extent, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 for (unsigned int x = 0; x < outputWidth; x++, p++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	370 float volumeX, volumeY, volumeZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	371 it.GetVolumeCoordinates(volumeX, volumeY, volumeZ);
183 98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	372
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	373 shader.Apply(p,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	374 volumeX * sourceWidth,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	375 volumeY * sourceHeight,
98da3a8d4820 SubvoxelReader Sebastien Jodogne <s.jodogne@gmail.com> parents: 180 diff changeset	376 volumeZ * sourceDepth);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 it.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 }
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 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 Orthanc::PixelFormat OutputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	390 const OrientedVolumeBoundingBox& box,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	391 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	394 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	395 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	396 if (hasLinearFunction)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	398 switch (interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	400 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	402 ImageInterpolation_Nearest, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408 ImageInterpolation_Bilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	409 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	411
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	412 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	413 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 ImageInterpolation_Trilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	418 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	419 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 }
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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	423 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	424 switch (interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	425 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	426 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	427 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428 ImageInterpolation_Nearest, TransferFunction_Copy>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	432 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	433 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	434 ImageInterpolation_Bilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	435 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	436 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	438 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	439 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440 ImageInterpolation_Trilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	441 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	442 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	443
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	445 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	446 }
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 template <typename RowIterator>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	452 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	456 const OrientedVolumeBoundingBox& box,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	459 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	460 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 {
1772 6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	462 if (source.GetFormat() == Orthanc::PixelFormat_Grayscale8 &&
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 slice.GetFormat() == Orthanc::PixelFormat_Grayscale8)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 ProcessImage<RowIterator,
1772 6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	466 Orthanc::PixelFormat_Grayscale8,
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	467 Orthanc::PixelFormat_Grayscale8>
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	468 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	469 }
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	470 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	471 slice.GetFormat() == Orthanc::PixelFormat_Grayscale8)
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	472 {
6c246f862b00 unit test VolumeRendering.Basic Sebastien Jodogne <s.jodogne@gmail.com> parents: 1739 diff changeset	473 ProcessImage<RowIterator,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475 Orthanc::PixelFormat_Grayscale8>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	476 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	477 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	478 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 slice.GetFormat() == Orthanc::PixelFormat_Grayscale16)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	480 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	481 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	483 Orthanc::PixelFormat_Grayscale16>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	484 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	485 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 else if (source.GetFormat() == Orthanc::PixelFormat_SignedGrayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	487 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	488 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	490 Orthanc::PixelFormat_SignedGrayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	492 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	493 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	494 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	495 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
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 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	498 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	500 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
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 else
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 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 }
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
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	511 void VolumeReslicer::CheckIterators(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	512 const CoordinateSystem3D& plane,
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	513 const OrientedVolumeBoundingBox& box) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 {
1624 59f95b9ea858 fix build Sebastien Jodogne <s.jodogne@gmail.com> parents: 1610 diff changeset	515 const unsigned int width = slice_->GetWidth();
59f95b9ea858 fix build Sebastien Jodogne <s.jodogne@gmail.com> parents: 1610 diff changeset	516 const unsigned int height = slice_->GetHeight();
59f95b9ea858 fix build Sebastien Jodogne <s.jodogne@gmail.com> parents: 1610 diff changeset	517
59f95b9ea858 fix build Sebastien Jodogne <s.jodogne@gmail.com> parents: 1610 diff changeset	518 for (unsigned int y = 0; y < height; y++)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	520 FastRowIterator fast(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	521 SlowRowIterator slow(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	522
1624 59f95b9ea858 fix build Sebastien Jodogne <s.jodogne@gmail.com> parents: 1610 diff changeset	523 for (unsigned int x = 0; x < width; x++)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	525 float px, py, pz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	526 fast.GetVolumeCoordinates(px, py, pz);
153 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 float qx, qy, qz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	529 slow.GetVolumeCoordinates(qx, qy, qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531 Vector d;
158 a053ca7fa5c6 LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 154 diff changeset	532 LinearAlgebra::AssignVector(d, px - qx, py - qy, pz - qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533 double norm = boost::numeric::ublas::norm_2(d);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534 if (norm > 0.0001)
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 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 fast.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	540 slow.Next();
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
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 void VolumeReslicer::Reset()
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 success_ = false;
1610 b7630b1a0253 ISceneLayer::GetBoundingBox() returns void Sebastien Jodogne <s.jodogne@gmail.com> parents: 1598 diff changeset	549 extent_.Clear();
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 slice_.reset(NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552
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 float VolumeReslicer::GetMinOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	555 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	556 switch (outputFormat_)
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 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 return 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562 break;
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 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	565 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	566 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	567 }
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	570 float VolumeReslicer::GetMaxOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	571 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	572 switch (outputFormat_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	573 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	574 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	575 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	576 return static_cast<float>(std::numeric_limits<uint8_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	577 break;
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 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 return static_cast<float>(std::numeric_limits<uint16_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	583 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	584 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585 }
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 VolumeReslicer::VolumeReslicer() :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590 outputFormat_(Orthanc::PixelFormat_Grayscale8),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 interpolation_(ImageInterpolation_Nearest),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 fastMode_(true),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	593 success_(false)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 ResetLinearFunction();
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
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 void VolumeReslicer::GetLinearFunction(float& scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	600 float& offset) const
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 if (hasLinearFunction_)
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 scaling = scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	605 offset = offset_;
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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	608 {
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	615 void VolumeReslicer::ResetLinearFunction()
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_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 scaling_ = 1.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 offset_ = 0.0f;
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::SetLinearFunction(float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	625 float offset)
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 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 hasLinearFunction_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 scaling_ = scaling;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630 offset_ = offset;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 }
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	634 void VolumeReslicer::SetWindow(float low,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	635 float high)
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 //printf("Range in pixel values: %f->%f\n", low, high);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638 float scaling = (GetMaxOutputValue() - GetMinOutputValue()) / (high - low);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 float offset = GetMinOutputValue() - scaling * low;
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 SetLinearFunction(scaling, offset);
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 /float x = scaling_ low + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	644 float y = scaling_ * high + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645 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	646 }
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	649 void VolumeReslicer::FitRange(const ImageBuffer3D& image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 float minInputValue, maxInputValue;
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 if (!image.GetRange(minInputValue, maxInputValue) \|\|
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	654 maxInputValue < 1)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	655 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	656 ResetLinearFunction();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 else
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 SetWindow(minInputValue, maxInputValue);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 }
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
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 void VolumeReslicer::SetWindowing(ImageWindowing windowing,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	666 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	667 float rescaleSlope,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	668 float rescaleIntercept)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	669 {
589 3080ec4ec6b9 removed enum value: ImageWindowing_Default Sebastien Jodogne <s.jodogne@gmail.com> parents: 439 diff changeset	670 if (windowing == ImageWindowing_Custom)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	671 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	672 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
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 float center, width;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	676 ComputeWindowing(center, width, windowing, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	677
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	678 float a = (center - width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	679 float b = (center + width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680 SetWindow(a, b);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 }
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 void VolumeReslicer::SetOutputFormat(Orthanc::PixelFormat format)
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 if (format != Orthanc::PixelFormat_Grayscale8 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 format != Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688 format != Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
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 if (hasLinearFunction_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695 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	696 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	697
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	698 outputFormat_ = format;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 }
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 void VolumeReslicer::SetInterpolation(ImageInterpolation interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705 if (interpolation != ImageInterpolation_Nearest &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706 interpolation != ImageInterpolation_Bilinear &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707 interpolation != ImageInterpolation_Trilinear)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	709 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710 }
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 interpolation_ = interpolation;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	713 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714 }
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	717 const Extent2D& VolumeReslicer::GetOutputExtent() const
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 if (success_)
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 return extent_;
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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 783 diff changeset	725 LOG(ERROR) << "VolumeReslicer::GetOutputExtent(): (!success_)";
a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 783 diff changeset	726
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 }
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
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 const Orthanc::ImageAccessor& VolumeReslicer::GetOutputSlice() const
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 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	737 return *slice_;
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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 783 diff changeset	741 LOG(ERROR) << "VolumeReslicer::GetOutputSlice(): (!success_)";
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 }
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
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 Orthanc::ImageAccessor* VolumeReslicer::ReleaseOutputSlice()
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 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	750 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	751 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 success_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753 return slice_.release();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 783 diff changeset	757 LOG(ERROR) << "VolumeReslicer::ReleaseOutputSlice(): (!success_)";
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
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 }
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 void VolumeReslicer::Apply(const ImageBuffer3D& source,
735 c3bbb130abc4 removing dependencies in ImageBuffer3D Sebastien Jodogne <s.jodogne@gmail.com> parents: 684 diff changeset	764 const VolumeImageGeometry& geometry,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765 const CoordinateSystem3D& plane)
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 // Choose the default voxel size as the finest voxel dimension
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 // of the source volumetric image
1640 52b8b96cb55f cleaning namespaces Sebastien Jodogne <s.jodogne@gmail.com> parents: 1624 diff changeset	769 const Vector dim = geometry.GetVoxelDimensions(VolumeProjection_Axial);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	770 double voxelSize = dim[0];
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 if (dim[1] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	773 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	774 voxelSize = dim[1];
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 if (dim[2] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779 voxelSize = dim[2];
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782 if (voxelSize <= 0)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	786
735 c3bbb130abc4 removing dependencies in ImageBuffer3D Sebastien Jodogne <s.jodogne@gmail.com> parents: 684 diff changeset	787 Apply(source, geometry, plane, voxelSize);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	788 }
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791 void VolumeReslicer::Apply(const ImageBuffer3D& source,
735 c3bbb130abc4 removing dependencies in ImageBuffer3D Sebastien Jodogne <s.jodogne@gmail.com> parents: 684 diff changeset	792 const VolumeImageGeometry& geometry,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794 double voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 Reset();
783 cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	797 pixelSpacing_ = voxelSize;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 // Firstly, compute the intersection of the source volumetric
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	800 // 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	801 // to 6 vertices. We compute the extent of the intersection
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802 // polygon, with respect to the coordinate system of the reslicing
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 // plane.
742 fa5febe0f0c2 moved OrientedBoundingBox in the Volumes folder Sebastien Jodogne <s.jodogne@gmail.com> parents: 735 diff changeset	804 OrientedVolumeBoundingBox box(geometry);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	806 if (!box.ComputeExtent(extent_, plane))
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	807 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 // 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	809 slice_.reset(new Orthanc::Image(outputFormat_, 0, 0, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 return;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	812 }
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 // Secondly, the extent together with the voxel size gives the
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 // size of the output image
172 316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	816 unsigned int width = boost::math::iround(extent_.GetWidth() / voxelSize);
316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	817 unsigned int height = boost::math::iround(extent_.GetHeight() / voxelSize);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	818
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 slice_.reset(new Orthanc::Image(outputFormat_, width, height, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	821 //CheckIterators(source, plane, box);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	822
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	823 if (fastMode_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	824 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	825 ProcessImage<FastRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	826 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	827 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	829 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	830 ProcessImage<SlowRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835 }
783 cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	836
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	837
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	838 double VolumeReslicer::GetPixelSpacing() const
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	839 {
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	840 if (success_)
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	841 {
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	842 return pixelSpacing_;
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	843 }
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	844 else
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	845 {
956 a7351ad54960 Made IsContextLost automatically set the flag by checking with the emscripten Benjamin Golinvaux <bgo@osimis.io> parents: 783 diff changeset	846 LOG(ERROR) << "VolumeReslicer::GetPixelSpacing(): (!success_)";
783 cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	847 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	848 }
cd13a062c9bd DicomVolumeImage Sebastien Jodogne <s.jodogne@gmail.com> parents: 742 diff changeset	849 }
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	850 }

Mercurial > hg > orthanc-stone

annotate OrthancStone/Sources/Volumes/VolumeReslicer.cpp @ 1772:6c246f862b00