advancedusfusion.frag 5.47 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29
// ================================================================================================
// 
// This file is part of the CAMPVis Software Framework.
// 
// If not explicitly stated otherwise: Copyright (C) 2012, all rights reserved,
//      Christian Schulte zu Berge <christian.szb@in.tum.de>
//      Chair for Computer Aided Medical Procedures
//      Technische Universitt Mnchen
//      Boltzmannstr. 3, 85748 Garching b. Mnchen, Germany
// For a full list of authors and contributors, please refer to the file "AUTHORS.txt".
// 
// The licensing of this softare is not yet resolved. Until then, redistribution in source or
// binary forms outside the CAMP chair is not permitted, unless explicitly stated in legal form.
// However, the names of the original authors and the above copyright notice must retain in its
// original state in any case.
// 
// Legal disclaimer provided by the BSD license:
// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR
// IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY 
// AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR 
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR 
// CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR 
// SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY 
// THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR 
// OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
// POSSIBILITY OF SUCH DAMAGE.
// 
// ================================================================================================

30
#version 330
31

32 33
in vec3 ex_TexCoord;
out vec4 out_Color;
34

35
#include "tools/colorspace.frag"
36
#include "tools/texture2d.frag"
37 38
#include "tools/texture3d.frag"
#include "tools/transferfunction.frag"
39

40 41 42 43 44 45 46 47 48
#ifdef USE_3D_TEX
#define SAMPLER_TYPE sampler3D
#define TEXPARAMS_TYPE TextureParameters3D
#define TEXTURE_LOOKUP_FUNC getElement3DNormalized
#else
#define SAMPLER_TYPE sampler2D
#define TEXPARAMS_TYPE TextureParameters2D
#define TEXTURE_LOOKUP_FUNC getElement2DNormalized
#endif
49

50 51
uniform SAMPLER_TYPE _usImage;
uniform TEXPARAMS_TYPE _usTextureParams;
52

53 54 55 56 57
uniform SAMPLER_TYPE _blurredImage;
uniform TEXPARAMS_TYPE _blurredTextureParams;

uniform SAMPLER_TYPE _confidenceMap;
uniform TEXPARAMS_TYPE _confidenceTextureParams;
58 59 60

uniform sampler1D _transferFunction;
uniform TFParameters1D _transferFunctionParams;
61

62
uniform int _sliceNumber;
63 64
uniform int _viewIndex;
uniform float _blurredScaling;
65

66
void main() {
67
#ifdef USE_3D_TEX
68
    vec3 texCoord = vec3(ex_TexCoord.xy, _usTextureParams._sizeRCP.z * (_sliceNumber + 0.5));
69 70 71
#else
    vec2 texCoord = ex_TexCoord.xy;
#endif
72

73 74 75
    vec4 texel = TEXTURE_LOOKUP_FUNC(_usImage, _usTextureParams, texCoord);
    vec4 blurred = TEXTURE_LOOKUP_FUNC(_blurredImage, _blurredTextureParams, texCoord) * _blurredScaling;
    float confidence = TEXTURE_LOOKUP_FUNC(_confidenceMap, _confidenceTextureParams, texCoord).a;
76

77 78
    switch (_viewIndex) {
        case 0:
79
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
80 81
            break;
        case 1:
82
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, blurred.a);
83 84
            break;
        case 2:
85
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a + texel.a - blurred.a);
86 87
            break;
        case 3:
88
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
89
            vec3 hsv = rgb2hsv(out_Color.xyz);
90
            hsv.x = 0.15;
91 92 93 94
            hsv.y = 1.0 - confidence;
            out_Color.xyz = hsv2rgb(hsv);
            break;
        case 4:
95
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
96 97 98 99 100 101
            vec3 hsl = rgb2hsl(out_Color.xyz);
            hsl.x = 0.15;
            hsl.y = 1.0 - confidence;
            out_Color.xyz = hsl2rgb(hsl);
            break;
        case 5:
102 103 104 105 106 107 108 109 110 111 112 113 114 115
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
            vec3 tsl = rgb2tsl(out_Color.xyz);
            //tsl.x = 0.15;
            //tsl.y = 1.0 - confidence;
            out_Color.xyz = tsl2rgb(tsl);
            break;
        case 6:
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
            vec3 hcl = rgb2hcl(out_Color.xyz);
            hcl.x = 0.15;
            hcl.y = 1.0 - confidence;
            out_Color.xyz = hcl2rgb(hcl);
            break;
        case 7:
116
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
117 118 119 120 121
            vec3 hcy = rgb2hcy(out_Color.xyz);
            hcy.x = 0.15;
            hcy.y = 1.0 - confidence;
            out_Color.xyz = hcy2rgb(hcy);
            break;
122 123 124 125 126 127 128 129
        case 8:
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, texel.a);
            vec3 lab = rgb2lab(out_Color.xyz);
            lab.y = (1.0 - confidence);
            //lab.z = 50.0 * (1.0 - confidence);
            out_Color.xyz = lab2rgb(lab);
            break;
        case 9:
130
            float intensity = mix(blurred.a, 2.0 * texel.a - blurred.a, confidence);
131
            out_Color = lookupTF(_transferFunction, _transferFunctionParams, intensity);
132 133
            break;
    }
134
}