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

annotate Framework/Volumes/VolumeReslicer.cpp @ 177:83200c4d07ca wasm

fix interpolation

author	Sebastien Jodogne <s.jodogne@gmail.com>
date	Fri, 09 Mar 2018 12:32:01 +0100
parents	316324f42848
children	db21c1810c89

rev	line source
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1 #include "VolumeReslicer.h"
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	2
159 0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	3 #include "../Toolbox/GeometryToolbox.h"
0a73d76333db populating LinearAlgebra Sebastien Jodogne <s.jodogne@gmail.com> parents: 158 diff changeset	4
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	5 #include <Core/Logging.h>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	6 #include <Core/OrthancException.h>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	7
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	8 #include <boost/math/special_functions/round.hpp>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	9
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	10 #if defined(_MSC_VER)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	11 # define ORTHANC_STONE_FORCE_INLINE __forceinline
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	12 #elif defined(__GNUC__) \|\| defined(__clang__) \|\| defined(__EMSCRIPTEN__)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	13 # define ORTHANC_STONE_FORCE_INLINE inline __attribute((always_inline))
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	14 #else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	15 # error Please support your compiler here
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	16 #endif
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	17
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	18
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	19 namespace OrthancStone
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	20 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	21 // Anonymous namespace to avoid clashes between compilation modules
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	22 namespace
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	23 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	24 enum TransferFunction
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	25 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	26 TransferFunction_Copy,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	27 TransferFunction_Float,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	28 TransferFunction_Linear
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	29 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	30
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	31
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	32 template <Orthanc::PixelFormat Format>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	33 struct PixelTraits;
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 template <>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	36 struct PixelTraits<Orthanc::PixelFormat_Grayscale8>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	37 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	38 typedef uint8_t PixelType;
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 static void SetOutOfVolume(PixelType& target)
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 target = 0;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	43 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	44
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	45 static void GetVoxel(PixelType& target,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	46 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	47 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	48 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	49 unsigned int z)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	50 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	51 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	52 target = image.GetVoxelGrayscale8Unchecked(x, y, z);
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	55 static float GetFloatVoxel(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	56 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	57 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	58 unsigned int z)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	59 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	60 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	61 return static_cast<float>(image.GetVoxelGrayscale8Unchecked(x, y, z));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	62 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	63 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	64
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	65 template <>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	66 struct PixelTraits<Orthanc::PixelFormat_Grayscale16>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	67 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	68 typedef uint16_t PixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	69
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	70 static void SetOutOfVolume(PixelType& target)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	71 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	72 target = 0;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	73 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	74
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	75 static void GetVoxel(PixelType& target,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	76 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	77 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	78 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	79 unsigned int z)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	80 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	81 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	82 target = image.GetVoxelGrayscale16Unchecked(x, y, z);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	83 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	84
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	85 static float GetFloatVoxel(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	86 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	87 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	88 unsigned int z)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	89 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	90 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	91 return static_cast<float>(image.GetVoxelGrayscale16Unchecked(x, y, z));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	92 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	93 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	94
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	95
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	96 template <>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	97 struct PixelTraits<Orthanc::PixelFormat_SignedGrayscale16>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	98 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	99 typedef int16_t PixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	100
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	101 static void SetOutOfVolume(PixelType& target)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	102 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	103 target = std::numeric_limits<PixelType>::min();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	104 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	105
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	106 static void GetVoxel(PixelType& target,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	107 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	108 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	109 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	110 unsigned int z)
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 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	113 target = image.GetVoxelSignedGrayscale16Unchecked(x, y, z);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	114 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	115
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	116 static float GetFloatVoxel(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	117 unsigned int x,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	118 unsigned int y,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	119 unsigned int z)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	120 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	121 assert(x < image.GetWidth() && y < image.GetHeight() && z < image.GetDepth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	122 return static_cast<float>(image.GetVoxelSignedGrayscale16Unchecked(x, y, z));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	123 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	124 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	125
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	126
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	127 template <>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	128 struct PixelTraits<Orthanc::PixelFormat_BGRA32>
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 struct PixelType
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	131 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	132 uint8_t blue_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	133 uint8_t green_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	134 uint8_t red_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	135 uint8_t alpha_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	136 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	137 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	138
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	139
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	140
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	141 template <Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	142 Orthanc::PixelFormat OutputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	143 class PixelWriter
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	144 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	145 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	146 typedef typename PixelTraits<InputFormat>::PixelType InputPixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	147 typedef PixelTraits<OutputFormat> OutputPixelTraits;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	148 typedef typename PixelTraits<OutputFormat>::PixelType OutputPixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	149
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	150 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	151 template <typename T>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	152 static void SetValueInternal(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	153 const T& value)
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 if (value < std::numeric_limits<OutputPixelType>::min())
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	156 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	157 *pixel = std::numeric_limits<OutputPixelType>::min();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	158 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	159 else if (value > std::numeric_limits<OutputPixelType>::max())
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	160 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	161 *pixel = std::numeric_limits<OutputPixelType>::max();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	162 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	163 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	164 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	165 *pixel = static_cast<OutputPixelType>(value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	166 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	167 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	168
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	169 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	170 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	171 void SetFloatValue(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	172 float value) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	173 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	174 SetValueInternal<float>(pixel, value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	175 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	176
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	177 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	178 void SetValue(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	179 const InputPixelType& value) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	180 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	181 SetValueInternal<InputPixelType>(pixel, value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	182 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	183 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	184
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	185
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	186 template <Orthanc::PixelFormat InputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	187 class PixelWriter<InputFormat, Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	188 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	189 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	190 typedef typename PixelTraits<InputFormat>::PixelType InputPixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	191 typedef PixelTraits<Orthanc::PixelFormat_BGRA32> OutputPixelTraits;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	192 typedef PixelTraits<Orthanc::PixelFormat_BGRA32>::PixelType OutputPixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	193
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	194 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	195 template <typename T>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	196 static void SetValueInternal(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	197 const T& value)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	198 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	199 uint8_t v;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	200 if (value < 0)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	201 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	202 v = 0;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	203 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	204 else if (value >= 255.0f)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	205 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	206 v = 255;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	207 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	208 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	209 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	210 v = static_cast<uint8_t>(value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	211 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	212
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	213 pixel->blue_ = v;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	214 pixel->green_ = v;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	215 pixel->red_ = v;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	216 pixel->alpha_ = 255;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	217 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	218
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	219 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	220 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	221 void SetFloatValue(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	222 float value) const
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 SetValueInternal<float>(pixel, value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	225 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	226
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	227 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	228 void SetValue(OutputPixelType* pixel,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	229 const InputPixelType& value) const
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 SetValueInternal<InputPixelType>(pixel, value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	232 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	233 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	234
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	235
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	236
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	237 class VoxelReaderBase : public boost::noncopyable
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	238 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	239 private:
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	240 const ImageBuffer3D& image_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	241 unsigned int width_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	242 unsigned int height_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	243 unsigned int depth_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	244 float widthFloat_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	245 float heightFloat_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	246 float depthFloat_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	247
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	248 public:
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	249 VoxelReaderBase(const ImageBuffer3D& image) :
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	250 image_(image),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	251 width_(image.GetWidth()),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	252 height_(image.GetHeight()),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	253 depth_(image.GetDepth()),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	254 widthFloat_(static_cast<float>(image.GetWidth())),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	255 heightFloat_(static_cast<float>(image.GetHeight())),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	256 depthFloat_(static_cast<float>(image.GetDepth()))
153 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 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	259
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	260 const ImageBuffer3D& GetImage() const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	261 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	262 return image_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	263 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	264
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	265 const unsigned int GetImageWidth() const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	266 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	267 return width_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	268 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	269
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	270 const unsigned int GetImageHeight() const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	271 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	272 return height_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	273 }
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	274
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	275 const unsigned int GetImageDepth() const
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	276 {
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	277 return depth_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	278 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	279
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	280 bool GetNearestCoordinates(unsigned int& imageX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	281 unsigned int& imageY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	282 unsigned int& imageZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	283 float& fractionalX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	284 float& fractionalY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	285 float& fractionalZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	286 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	287 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	288 float volumeZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	289 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	290 if (volumeX >= 0 &&
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	291 volumeY >= 0 &&
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	292 volumeZ >= 0)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	293 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	294 const float x = volumeX * widthFloat_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	295 const float y = volumeY * heightFloat_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	296 const float z = volumeZ * depthFloat_;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	297
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	298 imageX = static_cast<unsigned int>(std::floor(x));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	299 imageY = static_cast<unsigned int>(std::floor(y));
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	300 imageZ = static_cast<unsigned int>(std::floor(z));
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	301
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	302 if (imageX < width_ &&
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	303 imageY < height_ &&
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	304 imageZ < depth_)
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	305 {
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	306 fractionalX = x - static_cast<float>(imageX);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	307 fractionalY = y - static_cast<float>(imageY);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	308 fractionalZ = z - static_cast<float>(imageZ);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	309 return true;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	310 }
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	311 else
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	312 {
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	313 return false;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	314 }
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	315 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	316 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	317 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	318 return false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	319 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	320 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	321 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	322
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	323
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	324 template <Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	325 ImageInterpolation Interpolation>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	326 class VoxelReader;
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	329 template <Orthanc::PixelFormat InputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	330 class VoxelReader<InputFormat, ImageInterpolation_Nearest> :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	331 public VoxelReaderBase
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	332 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	333 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	334 typedef typename PixelTraits<InputFormat>::PixelType InputPixelType;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	335
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	336 VoxelReader(const ImageBuffer3D& image) :
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	337 VoxelReaderBase(image)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	338 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	339 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	340
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	341 ORTHANC_STONE_FORCE_INLINE
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	342 float GetFloatValue(float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	343 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	344 float volumeZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	345 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	346 InputPixelType value;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	347 GetValue(value, volumeX, volumeY, volumeZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	348 return static_cast<float>(value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	349 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	350
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	351 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	352 void GetValue(InputPixelType& target,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	353 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	354 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	355 float volumeZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	356 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	357 unsigned int imageX, imageY, imageZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	358 float fractionalX, fractionalY, fractionalZ; // unused
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	359
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	360 if (GetNearestCoordinates(imageX, imageY, imageZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	361 fractionalX, fractionalY, fractionalZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	362 volumeX, volumeY, volumeZ))
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	363 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	364 PixelTraits<InputFormat>::GetVoxel(target, GetImage(), imageX, imageY, imageZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	365 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	366 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	367 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	368 PixelTraits<InputFormat>::SetOutOfVolume(target);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	369 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	370 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	371 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	372
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	373
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	374 template <Orthanc::PixelFormat InputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	375 class VoxelReader<InputFormat, ImageInterpolation_Bilinear> :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	376 public VoxelReaderBase
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	377 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	378 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	379 float outOfVolume_;
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 public:
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	382 VoxelReader(const ImageBuffer3D& image) :
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	383 VoxelReaderBase(image)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	384 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	385 typename PixelTraits<InputFormat>::PixelType value;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	386 PixelTraits<InputFormat>::SetOutOfVolume(value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	387 outOfVolume_ = static_cast<float>(value);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	388 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	389
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	390 void SampleVoxels(float& f00,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	391 float& f01,
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	392 float& f10,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	393 float& f11,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	394 unsigned int imageX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	395 unsigned int imageY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	396 unsigned int imageZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	397 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	398 f00 = PixelTraits<InputFormat>::GetFloatVoxel(GetImage(), imageX, imageY, imageZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	399
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	400 if (imageX + 1 < GetImageWidth())
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	401 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	402 f01 = PixelTraits<InputFormat>::GetFloatVoxel(GetImage(), imageX + 1, imageY, imageZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	403 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	404 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	405 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	406 f01 = f00;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	407 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	408
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	409 if (imageY + 1 < GetImageWidth())
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	410 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	411 f10 = PixelTraits<InputFormat>::GetFloatVoxel(GetImage(), imageX, imageY + 1, imageZ);
153 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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	414 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	415 f10 = f00;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	416 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	417
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	418 if (imageX + 1 < GetImageWidth() &&
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	419 imageY + 1 < GetImageHeight())
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	420 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	421 f11 = PixelTraits<InputFormat>::GetFloatVoxel(GetImage(), imageX + 1, imageY + 1, imageZ);
153 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 f11 = f00;
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 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	428
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	429 float GetOutOfVolume() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	430 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	431 return outOfVolume_;
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
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	434 float GetFloatValue(float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	435 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	436 float volumeZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	437 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	438 unsigned int imageX, imageY, imageZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	439 float fractionalX, fractionalY, fractionalZ;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	440
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	441 if (GetNearestCoordinates(imageX, imageY, imageZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	442 fractionalX, fractionalY, fractionalZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	443 volumeX, volumeY, volumeZ))
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	444 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	445 float f00, f01, f10, f11;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	446 SampleVoxels(f00, f01, f10, f11, imageX, imageY, imageZ);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	447 return GeometryToolbox::ComputeBilinearInterpolationUnitSquare
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	448 (fractionalX, fractionalY, f00, f01, f10, f11);
153 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 else
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 return outOfVolume_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	453 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	454 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	455 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	456
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	457
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	458 template <Orthanc::PixelFormat InputFormat>
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	459 class VoxelReader<InputFormat, ImageInterpolation_Trilinear> :
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	460 public VoxelReaderBase
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	461 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	462 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	463 typedef VoxelReader<InputFormat, ImageInterpolation_Bilinear> Bilinear;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	464
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	465 Bilinear bilinear_;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	466 unsigned int imageDepth_;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	467
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	468 public:
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	469 VoxelReader(const ImageBuffer3D& image) :
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	470 VoxelReaderBase(image),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	471 bilinear_(image),
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	472 imageDepth_(image.GetDepth())
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	473 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	474 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	475
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	476 float GetFloatValue(float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	477 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	478 float volumeZ) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	479 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	480 unsigned int imageX, imageY, imageZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	481 float fractionalX, fractionalY, fractionalZ;
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	482
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	483 if (GetNearestCoordinates(imageX, imageY, imageZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	484 fractionalX, fractionalY, fractionalZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	485 volumeX, volumeY, volumeZ))
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	486 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	487 float f000, f001, f010, f011;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	488 bilinear_.SampleVoxels(f000, f001, f010, f011, imageX, imageY, imageZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	489
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	490 if (imageZ + 1 < imageDepth_)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	491 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	492 float f100, f101, f110, f111;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	493 bilinear_.SampleVoxels(f100, f101, f110, f111, imageX, imageY, imageZ + 1);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	494 return GeometryToolbox::ComputeTrilinearInterpolationUnitSquare
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	495 (fractionalX, fractionalY, fractionalZ,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	496 f000, f001, f010, f011, f100, f101, f110, f111);
153 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 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	499 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	500 return GeometryToolbox::ComputeBilinearInterpolationUnitSquare
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	501 (fractionalX, fractionalY, f000, f001, f010, f011);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	502 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	503 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	504 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	505 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	506 return bilinear_.GetOutOfVolume();
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 template <typename VoxelReader,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	513 typename PixelWriter,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	514 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	515 class PixelShader;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	516
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	517
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	518 template <typename VoxelReader,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	519 typename PixelWriter>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	520 class PixelShader<VoxelReader, PixelWriter, TransferFunction_Copy>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	521 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	522 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	523 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	524 PixelWriter writer_;
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 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	527 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	528 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	529 float /* offset */) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	530 reader_(image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	531 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	532 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	533
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	534 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	535 void Apply(typename PixelWriter::OutputPixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	536 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	537 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	538 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	539 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	540 typename VoxelReader::InputPixelType image;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	541 reader_.GetValue(image, volumeX, volumeY, volumeZ);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	542 writer_.SetValue(pixel, image);
153 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 template <typename VoxelReader,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	548 typename PixelWriter>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	549 class PixelShader<VoxelReader, PixelWriter, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	550 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	551 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	552 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	553 PixelWriter writer_;
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 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	556 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	557 float /* scaling */,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	558 float /* offset */) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	559 reader_(image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	560 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	561 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	562
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	563 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	564 void Apply(typename PixelWriter::OutputPixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	565 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	566 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	567 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	568 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	569 writer_.SetFloatValue(pixel, reader_.GetFloatValue(volumeX, volumeY, volumeZ));
153 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 };
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	574 template <typename VoxelReader,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	575 typename PixelWriter>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	576 class PixelShader<VoxelReader, PixelWriter, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	577 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	578 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	579 VoxelReader reader_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	580 PixelWriter writer_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	581 float scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	582 float offset_;
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 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	585 PixelShader(const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	586 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	587 float offset) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	588 reader_(image),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	589 scaling_(scaling),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	590 offset_(offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	591 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	592 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	593
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	594 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	595 void Apply(typename PixelWriter::OutputPixelType* pixel,
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	596 float volumeX,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	597 float volumeY,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	598 float volumeZ)
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	599 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	600 writer_.SetFloatValue(pixel, scaling_ * reader_.GetFloatValue(volumeX, volumeY, volumeZ) + offset_);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	601 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	602 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	603
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	604
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	605
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	606 class FastRowIterator : public boost::noncopyable
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 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	609 float position_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	610 float offset_[3];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	611
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	612 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	613 FastRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	614 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	615 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	616 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	617 unsigned int y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	618 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	619 const double width = static_cast<double>(slice.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	620 const double height = static_cast<double>(slice.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	621 assert(y < height);
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 Vector q1 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	624 (extent.GetX1() + extent.GetWidth() * static_cast<double>(0) / static_cast<double>(width + 1),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	625 extent.GetY1() + extent.GetHeight() * static_cast<double>(y) / static_cast<double>(height + 1));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	626
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	627 Vector q2 = plane.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	628 (extent.GetX1() + extent.GetWidth() * static_cast<double>(width - 1) / static_cast<double>(width + 1),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	629 extent.GetY1() + extent.GetHeight() * static_cast<double>(y) / static_cast<double>(height + 1));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	630
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	631 Vector r1, r2;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	632 box.ToInternalCoordinates(r1, q1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	633 box.ToInternalCoordinates(r2, q2);
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 position_[0] = static_cast<float>(r1[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	636 position_[1] = static_cast<float>(r1[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	637 position_[2] = static_cast<float>(r1[2]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	638
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	639 Vector tmp = (r2 - r1) / static_cast<double>(width - 1);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	640 offset_[0] = static_cast<float>(tmp[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	641 offset_[1] = static_cast<float>(tmp[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	642 offset_[2] = static_cast<float>(tmp[2]);
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	645 ORTHANC_STONE_FORCE_INLINE
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	646 void Next()
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 position_[0] += offset_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	649 position_[1] += offset_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	650 position_[2] += offset_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	651 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	652
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	653 ORTHANC_STONE_FORCE_INLINE
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	654 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	655 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	656 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	657 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	658 x = position_[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	659 y = position_[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	660 z = position_[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	661 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	662 };
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	663
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	664
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	665 class SlowRowIterator : public boost::noncopyable
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	666 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	667 private:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	668 const Orthanc::ImageAccessor& slice_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	669 const Extent2D& extent_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	670 const CoordinateSystem3D& plane_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	671 const OrientedBoundingBox& box_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	672 unsigned int x_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	673 unsigned int y_;
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 public:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	676 SlowRowIterator(const Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	677 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	678 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	679 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	680 unsigned int y) :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	681 slice_(slice),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	682 extent_(extent),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	683 plane_(plane),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	684 box_(box),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	685 x_(0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	686 y_(y)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	687 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	688 assert(y_ < slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	689 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	690
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	691 void Next()
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 x_++;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	694 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	695
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	696 void GetVolumeCoordinates(float& x,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	697 float& y,
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	698 float& z) const
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	699 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	700 assert(x_ < slice_.GetWidth());
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 const double width = static_cast<double>(slice_.GetWidth());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	703 const double height = static_cast<double>(slice_.GetHeight());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	704
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	705 Vector q = plane_.MapSliceToWorldCoordinates
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	706 (extent_.GetX1() + extent_.GetWidth() * static_cast<double>(x_) / (width + 1.0),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	707 extent_.GetY1() + extent_.GetHeight() * static_cast<double>(y_) / (height + 1.0));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	708
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	709 Vector r;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	710 box_.ToInternalCoordinates(r, q);
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 x = static_cast<float>(r[0]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	713 y = static_cast<float>(r[1]);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	714 z = static_cast<float>(r[2]);
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
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	719 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	720 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	721 Orthanc::PixelFormat OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	722 ImageInterpolation Interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	723 TransferFunction Function>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	724 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	725 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	726 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	727 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	728 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	729 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	730 float offset)
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 typedef VoxelReader<InputFormat, Interpolation> Reader;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	733 typedef PixelWriter<InputFormat, OutputFormat> Writer;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	734 typedef PixelShader<Reader, Writer, Function> Shader;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	735
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	736 const unsigned int outputWidth = slice.GetWidth();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	737 const unsigned int outputHeight = slice.GetHeight();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	738
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	739 Shader shader(source, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	740
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	741 for (unsigned int y = 0; y < outputHeight; y++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	742 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	743 typename Writer::OutputPixelType* p = reinterpret_cast<typename Writer::OutputPixelType*>(slice.GetRow(y));
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 RowIterator it(slice, extent, plane, box, y);
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 for (unsigned int x = 0; x < outputWidth; x++, p++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	748 {
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	749 float volumeX, volumeY, volumeZ;
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	750 it.GetVolumeCoordinates(volumeX, volumeY, volumeZ);
83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	751 shader.Apply(p, volumeX, volumeY, volumeZ);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	752 it.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	753 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	754 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	755 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	756
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	757
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	758 template <typename RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	759 Orthanc::PixelFormat InputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	760 Orthanc::PixelFormat OutputFormat>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	761 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	762 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	763 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	764 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	765 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	766 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	767 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	768 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	769 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	770 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	771 if (hasLinearFunction)
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 switch (interpolation)
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 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	776 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	777 ImageInterpolation_Nearest, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	778 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	779 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	780
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	781 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	782 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	783 ImageInterpolation_Bilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	784 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	785 break;
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 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	788 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	789 ImageInterpolation_Trilinear, TransferFunction_Linear>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	790 (slice, extent, source, plane, box, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	791 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	792
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	793 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	794 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	795 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	796 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	797 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	798 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	799 switch (interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	800 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	801 case ImageInterpolation_Nearest:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	802 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	803 ImageInterpolation_Nearest, TransferFunction_Copy>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	804 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	805 break;
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 case ImageInterpolation_Bilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	808 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	809 ImageInterpolation_Bilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	810 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	811 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	812
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	813 case ImageInterpolation_Trilinear:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	814 ProcessImage<RowIterator, InputFormat, OutputFormat,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	815 ImageInterpolation_Trilinear, TransferFunction_Float>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	816 (slice, extent, source, plane, box, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	817 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	818
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	819 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	820 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
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 }
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
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 template <typename RowIterator>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	827 static void ProcessImage(Orthanc::ImageAccessor& slice,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	828 const Extent2D& extent,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	829 const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	830 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	831 const OrientedBoundingBox& box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	832 ImageInterpolation interpolation,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	833 bool hasLinearFunction,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	834 float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	835 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	836 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	837 if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	838 slice.GetFormat() == Orthanc::PixelFormat_Grayscale8)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	839 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	840 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	841 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	842 Orthanc::PixelFormat_Grayscale8>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	843 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	844 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	845 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	846 slice.GetFormat() == Orthanc::PixelFormat_Grayscale16)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	847 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	848 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	849 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	850 Orthanc::PixelFormat_Grayscale16>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	851 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	852 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	853 else if (source.GetFormat() == Orthanc::PixelFormat_SignedGrayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	854 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	855 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	856 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	857 Orthanc::PixelFormat_SignedGrayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	858 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	859 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	860 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	861 else if (source.GetFormat() == Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	862 slice.GetFormat() == Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	863 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	864 ProcessImage<RowIterator,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	865 Orthanc::PixelFormat_Grayscale16,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	866 Orthanc::PixelFormat_BGRA32>
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	867 (slice, extent, source, plane, box, interpolation, hasLinearFunction, scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	868 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	869 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	870 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	871 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	872 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	873 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	874 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	875
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	876
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	877
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	878 void VolumeReslicer::CheckIterators(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	879 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	880 const OrientedBoundingBox& box) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	881 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	882 for (unsigned int y = 0; y < slice_->GetHeight(); y++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	883 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	884 FastRowIterator fast(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	885 SlowRowIterator slow(*slice_, extent_, plane, box, y);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	886
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	887 for (unsigned int x = 0; x < slice_->GetWidth(); x++)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	888 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	889 float px, py, pz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	890 fast.GetVolumeCoordinates(px, py, pz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	891
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	892 float qx, qy, qz;
177 83200c4d07ca fix interpolation Sebastien Jodogne <s.jodogne@gmail.com> parents: 172 diff changeset	893 slow.GetVolumeCoordinates(qx, qy, qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	894
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	895 Vector d;
158 a053ca7fa5c6 LinearAlgebra toolbox Sebastien Jodogne <s.jodogne@gmail.com> parents: 154 diff changeset	896 LinearAlgebra::AssignVector(d, px - qx, py - qy, pz - qz);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	897 double norm = boost::numeric::ublas::norm_2(d);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	898 if (norm > 0.0001)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	899 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	900 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	901 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	902
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	903 fast.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	904 slow.Next();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	905 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	906 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	907 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	908
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	909
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	910 void VolumeReslicer::Reset()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	911 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	912 success_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	913 extent_.Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	914 slice_.reset(NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	915 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	916
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	917
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	918 float VolumeReslicer::GetMinOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	919 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	920 switch (outputFormat_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	921 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	922 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	923 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	924 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	925 return 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	926 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	927
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	928 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	929 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	930 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	931 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	932
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	933
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	934 float VolumeReslicer::GetMaxOutputValue() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	935 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	936 switch (outputFormat_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	937 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	938 case Orthanc::PixelFormat_Grayscale8:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	939 case Orthanc::PixelFormat_BGRA32:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	940 return static_cast<float>(std::numeric_limits<uint8_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	941 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	942
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	943 case Orthanc::PixelFormat_Grayscale16:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	944 return static_cast<float>(std::numeric_limits<uint16_t>::max());
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	945 break;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	946
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	947 default:
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	948 throw Orthanc::OrthancException(Orthanc::ErrorCode_NotImplemented);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	949 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	950 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	951
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	952
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	953 VolumeReslicer::VolumeReslicer() :
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	954 outputFormat_(Orthanc::PixelFormat_Grayscale8),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	955 interpolation_(ImageInterpolation_Nearest),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	956 fastMode_(true),
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	957 success_(false)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	958 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	959 ResetLinearFunction();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	960 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	961
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	962
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	963 void VolumeReslicer::GetLinearFunction(float& scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	964 float& offset) const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	965 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	966 if (hasLinearFunction_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	967 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	968 scaling = scaling_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	969 offset = offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	970 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	971 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	972 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	973 scaling = 1.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	974 offset = 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	975 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	976 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	977
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	978
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	979 void VolumeReslicer::ResetLinearFunction()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	980 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	981 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	982 hasLinearFunction_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	983 scaling_ = 1.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	984 offset_ = 0.0f;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	985 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	986
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	987
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	988 void VolumeReslicer::SetLinearFunction(float scaling,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	989 float offset)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	990 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	991 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	992 hasLinearFunction_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	993 scaling_ = scaling;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	994 offset_ = offset;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	995 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	996
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	997
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	998 void VolumeReslicer::SetWindow(float low,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	999 float high)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1000 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1001 //printf("Range in pixel values: %f->%f\n", low, high);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1002 float scaling = (GetMaxOutputValue() - GetMinOutputValue()) / (high - low);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1003 float offset = GetMinOutputValue() - scaling * low;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1004
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1005 SetLinearFunction(scaling, offset);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1006
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1007 /float x = scaling_ low + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1008 float y = scaling_ * high + offset_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1009 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	1010 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1011
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1012
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1013 void VolumeReslicer::FitRange(const ImageBuffer3D& image)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1014 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1015 float minInputValue, maxInputValue;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1016
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1017 if (!image.GetRange(minInputValue, maxInputValue) \|\|
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1018 maxInputValue < 1)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1019 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1020 ResetLinearFunction();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1021 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1022 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1023 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1024 SetWindow(minInputValue, maxInputValue);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1025 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1026 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1027
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1028
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1029 void VolumeReslicer::SetWindowing(ImageWindowing windowing,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1030 const ImageBuffer3D& image,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1031 float rescaleSlope,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1032 float rescaleIntercept)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1033 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1034 if (windowing == ImageWindowing_Custom \|\|
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1035 windowing == ImageWindowing_Default)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1036 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1037 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1038 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1039
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1040 float center, width;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1041 ComputeWindowing(center, width, windowing, 0, 0);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1042
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1043 float a = (center - width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1044 float b = (center + width / 2.0f - rescaleIntercept) / rescaleSlope;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1045 SetWindow(a, b);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1046 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1047
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1048
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1049 void VolumeReslicer::SetOutputFormat(Orthanc::PixelFormat format)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1050 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1051 if (format != Orthanc::PixelFormat_Grayscale8 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1052 format != Orthanc::PixelFormat_Grayscale16 &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1053 format != Orthanc::PixelFormat_BGRA32)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1054 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1055 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1056 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1057
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1058 if (hasLinearFunction_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1059 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1060 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	1061 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1062
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1063 outputFormat_ = format;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1064 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1065 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1066
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1067
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1068 void VolumeReslicer::SetInterpolation(ImageInterpolation interpolation)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1069 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1070 if (interpolation != ImageInterpolation_Nearest &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1071 interpolation != ImageInterpolation_Bilinear &&
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1072 interpolation != ImageInterpolation_Trilinear)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1073 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1074 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1075 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1076
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1077 interpolation_ = interpolation;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1078 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1079 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1080
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1081
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1082 const Extent2D& VolumeReslicer::GetOutputExtent() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1083 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1084 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1085 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1086 return extent_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1087 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1088 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1089 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1090 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1091 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1092 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1093
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1094
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1095 const Orthanc::ImageAccessor& VolumeReslicer::GetOutputSlice() const
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1096 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1097 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1098 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1099 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1100 return *slice_;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1101 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1102 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1103 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1104 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1105 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1106 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1107
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1108
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1109 Orthanc::ImageAccessor* VolumeReslicer::ReleaseOutputSlice()
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1110 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1111 if (success_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1112 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1113 assert(slice_.get() != NULL);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1114 success_ = false;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1115 return slice_.release();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1116 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1117 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1118 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1119 throw Orthanc::OrthancException(Orthanc::ErrorCode_BadSequenceOfCalls);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1120 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1121 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1122
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1123
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1124 void VolumeReslicer::Apply(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1125 const CoordinateSystem3D& plane)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1126 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1127 // Choose the default voxel size as the finest voxel dimension
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1128 // of the source volumetric image
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1129 const OrthancStone::Vector dim = source.GetVoxelDimensions(OrthancStone::VolumeProjection_Axial);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1130 double voxelSize = dim[0];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1131
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1132 if (dim[1] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1133 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1134 voxelSize = dim[1];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1135 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1136
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1137 if (dim[2] < voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1138 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1139 voxelSize = dim[2];
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1140 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1141
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1142 if (voxelSize <= 0)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1143 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1144 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1145 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1146
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1147 Apply(source, plane, voxelSize);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1148 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1149
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1150
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1151 void VolumeReslicer::Apply(const ImageBuffer3D& source,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1152 const CoordinateSystem3D& plane,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1153 double voxelSize)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1154 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1155 Reset();
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1156
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1157 // Firstly, compute the intersection of the source volumetric
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1158 // 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	1159 // to 6 vertices. We compute the extent of the intersection
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1160 // polygon, with respect to the coordinate system of the reslicing
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1161 // plane.
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1162 OrientedBoundingBox box(source);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1163
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1164 if (!box.ComputeExtent(extent_, plane))
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1165 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1166 // 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	1167 slice_.reset(new Orthanc::Image(outputFormat_, 0, 0, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1168 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1169 return;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1170 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1171
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1172 // Secondly, the extent together with the voxel size gives the
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1173 // size of the output image
172 316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	1174 unsigned int width = boost::math::iround(extent_.GetWidth() / voxelSize);
316324f42848 fix Sebastien Jodogne <s.jodogne@gmail.com> parents: 171 diff changeset	1175 unsigned int height = boost::math::iround(extent_.GetHeight() / voxelSize);
153 ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1176
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1177 slice_.reset(new Orthanc::Image(outputFormat_, width, height, false));
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1178
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1179 //CheckIterators(source, plane, box);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1180
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1181 if (fastMode_)
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1182 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1183 ProcessImage<FastRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1184 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1185 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1186 else
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1187 {
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1188 ProcessImage<SlowRowIterator>(*slice_, extent_, source, plane, box,
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1189 interpolation_, hasLinearFunction_, scaling_, offset_);
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1190 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1191
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1192 success_ = true;
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1193 }
ae531ab5dcd9 new class: VolumeReslicer Sebastien Jodogne <s.jodogne@gmail.com> parents: diff changeset	1194 }

Mercurial > hg > orthanc-stone

annotate Framework/Volumes/VolumeReslicer.cpp @ 177:83200c4d07ca wasm