Mercurial > hg > orthanc-stone

comparison Framework/Layers/GrayscaleFrameRenderer.cpp @ 0:351ab0da0150

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
initial commit
author Sebastien Jodogne <s.jodogne@gmail.com>
date Fri, 14 Oct 2016 15:34:11 +0200
parents
children ff1e935768e7
comparison
equal deleted inserted replaced
-1:000000000000 0:351ab0da0150
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 }