0
|
1 /**
|
|
2 * Stone of Orthanc
|
|
3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics
|
|
4 * Department, University Hospital of Liege, Belgium
|
|
5 *
|
|
6 * This program is free software: you can redistribute it and/or
|
|
7 * modify it under the terms of the GNU General Public License as
|
|
8 * published by the Free Software Foundation, either version 3 of the
|
|
9 * License, or (at your option) any later version.
|
|
10 *
|
|
11 * In addition, as a special exception, the copyright holders of this
|
|
12 * program give permission to link the code of its release with the
|
|
13 * OpenSSL project's "OpenSSL" library (or with modified versions of it
|
|
14 * that use the same license as the "OpenSSL" library), and distribute
|
|
15 * the linked executables. You must obey the GNU General Public License
|
|
16 * in all respects for all of the code used other than "OpenSSL". If you
|
|
17 * modify file(s) with this exception, you may extend this exception to
|
|
18 * your version of the file(s), but you are not obligated to do so. If
|
|
19 * you do not wish to do so, delete this exception statement from your
|
|
20 * version. If you delete this exception statement from all source files
|
|
21 * in the program, then also delete it here.
|
|
22 *
|
|
23 * This program is distributed in the hope that it will be useful, but
|
|
24 * WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
25 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
26 * General Public License for more details.
|
|
27 *
|
|
28 * You should have received a copy of the GNU General Public License
|
|
29 * along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
30 **/
|
|
31
|
|
32
|
|
33 #include "GrayscaleFrameRenderer.h"
|
|
34
|
|
35 #include "../Orthanc/Core/OrthancException.h"
|
|
36
|
|
37 namespace OrthancStone
|
|
38 {
|
|
39 CairoSurface* GrayscaleFrameRenderer::GenerateDisplay(const RenderStyle& style)
|
|
40 {
|
|
41 std::auto_ptr<CairoSurface> result;
|
|
42
|
|
43 float windowCenter, windowWidth;
|
|
44 style.ComputeWindowing(windowCenter, windowWidth,
|
|
45 defaultWindowCenter_, defaultWindowWidth_);
|
|
46
|
|
47 float x0 = windowCenter - windowWidth / 2.0f;
|
|
48 float x1 = windowCenter + windowWidth / 2.0f;
|
|
49
|
|
50 //LOG(INFO) << "Window: " << x0 << " => " << x1;
|
|
51
|
|
52 result.reset(new CairoSurface(frame_->GetWidth(), frame_->GetHeight()));
|
|
53
|
|
54 const uint8_t* lut = NULL;
|
|
55 if (style.applyLut_)
|
|
56 {
|
|
57 if (Orthanc::EmbeddedResources::GetFileResourceSize(style.lut_) != 3 * 256)
|
|
58 {
|
|
59 // Invalid colormap
|
|
60 throw Orthanc::OrthancException(Orthanc::ErrorCode_InternalError);
|
|
61 }
|
|
62
|
|
63 lut = reinterpret_cast<const uint8_t*>(Orthanc::EmbeddedResources::GetFileResourceBuffer(style.lut_));
|
|
64 }
|
|
65
|
|
66 Orthanc::ImageAccessor target = result->GetAccessor();
|
|
67 for (unsigned int y = 0; y < target.GetHeight(); y++)
|
|
68 {
|
|
69 const float* p = reinterpret_cast<const float*>(frame_->GetConstRow(y));
|
|
70 uint8_t* q = reinterpret_cast<uint8_t*>(target.GetRow(y));
|
|
71
|
|
72 for (unsigned int x = 0; x < target.GetWidth(); x++, p++, q += 4)
|
|
73 {
|
|
74 uint8_t v = 0;
|
|
75 if (windowWidth >= 0.001f) // Avoid division by zero
|
|
76 {
|
|
77 if (*p >= x1)
|
|
78 {
|
|
79 v = 255;
|
|
80 }
|
|
81 else if (*p <= x0)
|
|
82 {
|
|
83 v = 0;
|
|
84 }
|
|
85 else
|
|
86 {
|
|
87 // https://en.wikipedia.org/wiki/Linear_interpolation
|
|
88 v = static_cast<uint8_t>(255.0f * (*p - x0) / (x1 - x0));
|
|
89 }
|
|
90
|
|
91 if (style.reverse_)
|
|
92 {
|
|
93 v = 255 - v;
|
|
94 }
|
|
95 }
|
|
96
|
|
97 if (style.applyLut_)
|
|
98 {
|
|
99 assert(lut != NULL);
|
|
100 q[3] = 255;
|
|
101 q[2] = lut[3 * v];
|
|
102 q[1] = lut[3 * v + 1];
|
|
103 q[0] = lut[3 * v + 2];
|
|
104 }
|
|
105 else
|
|
106 {
|
|
107 q[3] = 255;
|
|
108 q[2] = v;
|
|
109 q[1] = v;
|
|
110 q[0] = v;
|
|
111 }
|
|
112 }
|
|
113 }
|
|
114
|
|
115 return result.release();
|
|
116 }
|
|
117
|
|
118
|
|
119 GrayscaleFrameRenderer::GrayscaleFrameRenderer(Orthanc::ImageAccessor* frame,
|
|
120 const DicomFrameConverter& converter,
|
|
121 const SliceGeometry& viewportSlice,
|
|
122 const SliceGeometry& frameSlice,
|
|
123 double pixelSpacingX,
|
|
124 double pixelSpacingY,
|
|
125 bool isFullQuality) :
|
|
126 FrameRenderer(viewportSlice, frameSlice, pixelSpacingX, pixelSpacingY, isFullQuality),
|
|
127 frame_(frame),
|
|
128 defaultWindowCenter_(converter.GetDefaultWindowCenter()),
|
|
129 defaultWindowWidth_(converter.GetDefaultWindowWidth())
|
|
130 {
|
|
131 if (frame == NULL)
|
|
132 {
|
|
133 throw Orthanc::OrthancException(Orthanc::ErrorCode_ParameterOutOfRange);
|
|
134 }
|
|
135
|
|
136 converter.ConvertFrame(frame_);
|
|
137 assert(frame_.get() != NULL);
|
|
138
|
|
139 if (frame_->GetFormat() != Orthanc::PixelFormat_Float32)
|
|
140 {
|
|
141 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat);
|
|
142 }
|
|
143 }
|
|
144 }
|