Mercurial > hg > orthanc-stone
annotate Framework/Scene2D/Internals/OpenGLFloatTextureProgram.cpp @ 770:e6850c866469
fix windowing
| author | Sebastien Jodogne <s.jodogne@gmail.com> |
|---|---|
| date | Fri, 24 May 2019 11:02:20 +0200 |
| parents | e3f21a265be5 |
| children | ed783cbc92d7 |
| rev | line source |
|---|---|
| 591 | 1 /** |
| 2 * Stone of Orthanc | |
| 3 * Copyright (C) 2012-2016 Sebastien Jodogne, Medical Physics | |
| 4 * Department, University Hospital of Liege, Belgium | |
| 5 * Copyright (C) 2017-2019 Osimis S.A., Belgium | |
| 6 * | |
| 7 * This program is free software: you can redistribute it and/or | |
| 8 * modify it under the terms of the GNU Affero General Public License | |
| 9 * as published by the Free Software Foundation, either version 3 of | |
| 10 * the License, or (at your option) any later version. | |
| 11 * | |
| 12 * This program is distributed in the hope that it will be useful, but | |
| 13 * WITHOUT ANY WARRANTY; without even the implied warranty of | |
| 14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU | |
| 15 * Affero General Public License for more details. | |
| 16 * | |
| 17 * You should have received a copy of the GNU Affero General Public License | |
| 18 * along with this program. If not, see <http://www.gnu.org/licenses/>. | |
| 19 **/ | |
| 20 | |
| 21 | |
| 22 #include "OpenGLFloatTextureProgram.h" | |
|
611
e3f21a265be5
Added version directive to GLSL shader code + glew init function in sample code
Benjamin Golinvaux <bgo@osimis.io>
parents:
592
diff
changeset
|
23 #include "OpenGLShaderVersionDirective.h" |
| 591 | 24 |
| 25 #include <Core/OrthancException.h> | |
| 26 #include <Core/Images/Image.h> | |
| 27 #include <Core/Images/ImageProcessing.h> | |
| 28 | |
| 29 | |
| 30 static const char* FRAGMENT_SHADER = | |
|
611
e3f21a265be5
Added version directive to GLSL shader code + glew init function in sample code
Benjamin Golinvaux <bgo@osimis.io>
parents:
592
diff
changeset
|
31 ORTHANC_STONE_OPENGL_SHADER_VERSION_DIRECTIVE |
| 591 | 32 "uniform float u_offset; \n" |
| 33 "uniform float u_slope; \n" | |
| 34 "uniform float u_windowCenter; \n" | |
| 35 "uniform float u_windowWidth; \n" | |
| 36 "uniform sampler2D u_texture; \n" | |
| 37 "varying vec2 v_texcoord; \n" | |
| 38 "void main() \n" | |
| 39 "{ \n" | |
| 40 " vec4 t = texture2D(u_texture, v_texcoord); \n" | |
| 41 " float v = (t.r * 256.0 + t.g) * 256.0; \n" | |
| 42 " v = v * u_slope + u_offset; \n" // (*) | |
| 43 " float a = u_windowCenter - u_windowWidth; \n" | |
| 770 | 44 " float dy = 1.0 / u_windowWidth; \n" |
| 591 | 45 " if (v <= a) \n" |
| 46 " v = 0.0; \n" | |
| 47 " else \n" | |
| 48 " { \n" | |
| 49 " v = (v - a) * dy; \n" | |
| 50 " if (v >= 1.0) \n" | |
| 51 " v = 1.0; \n" | |
| 52 " } \n" | |
| 53 " gl_FragColor = vec4(v, v, v, 1); \n" | |
| 54 "}"; | |
| 55 | |
| 56 | |
| 57 namespace OrthancStone | |
| 58 { | |
| 59 namespace Internals | |
| 60 { | |
| 61 OpenGLFloatTextureProgram::Data::Data(const Orthanc::ImageAccessor& texture, | |
| 62 bool isLinearInterpolation) | |
| 63 { | |
| 64 if (texture.GetFormat() != Orthanc::PixelFormat_Float32) | |
| 65 { | |
| 66 throw Orthanc::OrthancException(Orthanc::ErrorCode_IncompatibleImageFormat); | |
| 67 } | |
| 68 | |
| 69 float minValue, maxValue; | |
| 70 Orthanc::ImageProcessing::GetMinMaxFloatValue(minValue, maxValue, texture); | |
| 71 | |
| 72 offset_ = minValue; | |
| 73 | |
| 74 if (LinearAlgebra::IsCloseToZero(maxValue - minValue)) | |
| 75 { | |
| 76 slope_ = 1; | |
| 77 } | |
| 78 else | |
| 79 { | |
| 80 slope_ = (maxValue - minValue) / 65536.0f; | |
| 81 assert(!LinearAlgebra::IsCloseToZero(slope_)); | |
| 82 } | |
| 83 | |
| 84 const unsigned int width = texture.GetWidth(); | |
| 85 const unsigned int height = texture.GetHeight(); | |
| 86 | |
| 87 Orthanc::Image converted(Orthanc::PixelFormat_RGB24, width, height, true); | |
| 88 | |
| 89 for (unsigned int y = 0; y < height; y++) | |
| 90 { | |
| 91 const float *p = reinterpret_cast<const float*>(texture.GetConstRow(y)); | |
| 92 uint8_t *q = reinterpret_cast<uint8_t*>(converted.GetRow(y)); | |
| 93 | |
| 94 for (unsigned int x = 0; x < width; x++) | |
| 95 { | |
| 96 /** | |
| 97 * At (*), the floating-point "value" is reconstructed as | |
| 98 * "value = texture * slope + offset". | |
| 99 * <=> texture = (value - offset) / slope | |
| 100 **/ | |
| 101 | |
| 102 float texture = (*p - offset_) / slope_; | |
| 103 if (texture < 0) | |
| 104 { | |
| 105 texture = 0; | |
| 106 } | |
| 107 else if (texture >= 65535.0f) | |
| 108 { | |
| 109 texture = 65535.0f; | |
| 110 } | |
| 111 | |
| 112 uint16_t t = static_cast<uint16_t>(texture); | |
| 113 | |
| 114 q[0] = t / 256; // red | |
| 115 q[1] = t % 256; // green | |
| 116 q[2] = 0; // blue is unused | |
| 117 | |
| 118 p++; | |
| 119 q += 3; | |
| 120 } | |
| 121 } | |
| 122 | |
| 123 texture_.Load(converted, isLinearInterpolation); | |
| 124 } | |
| 125 | |
| 126 | |
| 592 | 127 OpenGLFloatTextureProgram::OpenGLFloatTextureProgram(OpenGL::IOpenGLContext& context) : |
| 128 program_(context, FRAGMENT_SHADER) | |
| 129 { | |
| 130 } | |
| 131 | |
| 132 | |
| 591 | 133 void OpenGLFloatTextureProgram::Apply(Data& data, |
| 134 const AffineTransform2D& transform, | |
| 135 float windowCenter, | |
| 136 float windowWidth) | |
| 137 { | |
| 138 OpenGLTextureProgram::Execution execution(program_, data.GetTexture(), transform); | |
| 139 | |
| 140 glUniform1f(execution.GetUniformLocation("u_slope"), data.GetSlope()); | |
| 141 glUniform1f(execution.GetUniformLocation("u_offset"), data.GetOffset()); | |
| 142 glUniform1f(execution.GetUniformLocation("u_windowCenter"), windowCenter); | |
| 143 glUniform1f(execution.GetUniformLocation("u_windowWidth"), windowWidth); | |
| 144 | |
| 145 execution.DrawTriangles(); | |
| 146 } | |
| 147 } | |
| 148 } |