glreductiontest.cpp 4.97 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
30
31
32
33
34
35
36
37
38
39
40
41
42
// ================================================================================================
// 
// This file is part of the CAMPVis Software Framework.
// 
// If not explicitly stated otherwise: Copyright (C) 2012-2013, all rights reserved,
//      Christian Schulte zu Berge <christian.szb@in.tum.de>
//      Chair for Computer Aided Medical Procedures
//      Technische Universität München
//      Boltzmannstr. 3, 85748 Garching b. München, Germany
// 
// For a full list of authors and contributors, please refer to the file "AUTHORS.txt".
// 
// Licensed under the Apache License, Version 2.0 (the "License"); you may not use this file 
// except in compliance with the License. You may obtain a copy of the License at
// 
// http://www.apache.org/licenses/LICENSE-2.0
// 
// Unless required by applicable law or agreed to in writing, software distributed under the 
// License is distributed on an "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, 
// either express or implied. See the License for the specific language governing permissions 
// and limitations under the License.
// 
// ================================================================================================

#include "glreductiontest.h"

#include "tgt/event/keyevent.h"
#include "core/datastructures/imagedata.h"
#include "core/datastructures/renderdata.h"

#include "core/classification/geometry1dtransferfunction.h"
#include "core/classification/tfgeometry1d.h"
#include "core/tools/glreduction.h"


namespace campvis {

    GlReductionTest::GlReductionTest(DataContainer* dc)
        : AutoEvaluationPipeline(dc)
        , _imageReader()
        , _resampler(&_canvasSize)
        , _ve(&_canvasSize)
43
44
        , _glr(nullptr)
        , _glr2(nullptr)
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
    {
        addProcessor(&_imageReader);
        addProcessor(&_resampler);
        addProcessor(&_ve);

        addEventListenerToBack(&_ve);
    }

    GlReductionTest::~GlReductionTest() {
    }

    void GlReductionTest::init() {
        AutoEvaluationPipeline::init();
    
        _glr = new GlReduction(GlReduction::MAX);
        _glr2 = new GlReduction(GlReduction::MIN);
        _imageReader.s_validated.connect(this, &GlReductionTest::onProcessorValidated);

        _ve.p_outputImage.setValue("result");
        _renderTargetID.setValue("result");

        _imageReader.p_url.setValue(CAMPVIS_SOURCE_DIR "/modules/vis/sampledata/smallHeart.mhd");
        _imageReader.p_targetImageID.setValue("reader.output");
        _imageReader.p_targetImageID.addSharedProperty(&_resampler.p_inputImage);

        _resampler.p_outputImage.setValue("resampled");
        _resampler.p_outputImage.addSharedProperty(&_ve.p_inputVolume);
        _resampler.s_validated.connect(this, &GlReductionTest::onProcessorValidated);

        Geometry1DTransferFunction* dvrTF = new Geometry1DTransferFunction(128, tgt::vec2(0.f, .05f));
        dvrTF->addGeometry(TFGeometry1D::createQuad(tgt::vec2(.1f, .125f), tgt::col4(255, 0, 0, 32), tgt::col4(255, 0, 0, 32)));
        dvrTF->addGeometry(TFGeometry1D::createQuad(tgt::vec2(.4f, .5f), tgt::col4(0, 255, 0, 128), tgt::col4(0, 255, 0, 128)));
77
        static_cast<TransferFunctionProperty*>(_ve.getNestedProperty("VolumeRendererProperties::RaycasterProps::TransferFunction"))->replaceTF(dvrTF);
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120

    }

    void GlReductionTest::deinit() {
        AutoEvaluationPipeline::deinit();

        delete _glr;
        delete _glr2;
    }


    void GlReductionTest::onProcessorValidated(AbstractProcessor* processor) {
        if (processor == &_resampler) {
            ImageRepresentationGL::ScopedRepresentation img(getDataContainer(), _resampler.p_outputImage.getValue());
            if (img != 0) {
                float* foo = reinterpret_cast<float*>(img->getTexture()->downloadTextureToBuffer(GL_RED, GL_FLOAT));
                std::vector<float> v(foo, foo + tgt::hmul(img->getTexture()->getDimensions()));
                float mmm = *std::max_element(v.begin(), v.end());

                float max = _glr->reduce(img->getTexture()).front();
                LDEBUG("3D: " << mmm << ", " << max << ", DIFF: " << (mmm - max));
            }

            ScopedTypedData<RenderData> result(getDataContainer(), _ve.p_outputImage.getValue() + ".raycaster");
            if (result != 0) {
                const ImageData* depthImg = result->getDepthTexture();
                if (depthImg != 0) {
                    const ImageRepresentationGL* depthRep = depthImg->getRepresentation<ImageRepresentationGL>();
                    if (depthRep != 0) {
                        float* foo = reinterpret_cast<float*>(depthRep->getTexture()->downloadTextureToBuffer());
                        std::vector<float> v(foo, foo + depthRep->getNumElements());
                        float min1 = *std::min_element(v.begin(), v.end());

                        float min2 = _glr2->reduce(depthRep->getTexture()).front();
                        LDEBUG("2D: " << min1 << ", " << min2 << ", DIFF: " << std::abs(min2 - min1));
                    }
                }
            }
        }
    }


}