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

annotate OrthancStone/Sources/Volumes/VolumeReslicer.cpp @ 1870:3889ae96d2e9

added copyright UCLouvain

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Tue, 11 Jan 2022 11:04:09 +0100
parents	95ece40bb298
children	7053b8a0aaec

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

Mercurial > hg > orthanc-stone

annotate OrthancStone/Sources/Volumes/VolumeReslicer.cpp @ 1870:3889ae96d2e9