source: nanovis/tags/1.2.3/Volume.cpp @ 6369

/* -*- mode: c++; c-basic-offset: 4; indent-tabs-mode: nil -*- */
/*
 * ----------------------------------------------------------------------
 * Volume.cpp: 3d volume class
 *
 * ======================================================================
 *  AUTHOR:  Wei Qiao <qiaow@purdue.edu>
 *           Purdue Rendering and Perceptualization Lab (PURPL)
 *
 *  Copyright (c) 2004-2013  HUBzero Foundation, LLC
 *
 *  See the file "license.terms" for information on usage and
 *  redistribution of this file, and for a DISCLAIMER OF ALL WARRANTIES.
 * ======================================================================
 */
#include <memory.h>
#include <assert.h>
#include <float.h>

#include <vrmath/Vector4f.h>
#include <vrmath/Matrix4x4d.h>

#include "Volume.h"
#include "Trace.h"

using namespace nv;
using namespace vrmath;

bool Volume::updatePending = false;
double Volume::valueMin = 0.0;
double Volume::valueMax = 1.0;

Volume::Volume(int width, int height, int depth,
               int numComponents, float *data,
               double vmin, double vmax, double nonZeroMin) :
    _id(0),
    _width(width),
    _height(height),
    _depth(depth),
    _transferFunc(NULL),
    _ambient(0.6f),
    _diffuse(0.4f),
    _specular(0.3f),
    _specularExp(90.0f),
    _lightTwoSide(false),
    _opacityScale(0.5f),
    _data(NULL),
    _numComponents(numComponents),
    _nonZeroMin(nonZeroMin),
    _cutplanesVisible(true),
    _tex(NULL),
    _position(0, 0, 0),
    _scale(1, 1, 1),
    _numSlices(512),
    _enabled(true),
    _dataEnabled(true),
    _outlineEnabled(true),
    _volumeType(CUBIC),
    _isosurface(0)
{
    TRACE("Enter: %dx%dx%d", _width, _height, _depth);

    _outlineColor[0] = _outlineColor[1] = _outlineColor[2] = 1.0f;

    _tex = new Texture3D(_width, _height, _depth, GL_FLOAT, GL_LINEAR, _numComponents);
    int fcount = _width * _height * _depth * _numComponents;
    _data = new float[fcount];
    memcpy(_data, data, fcount * sizeof(float));
    _tex->initialize(_data);

    _id = _tex->id();

    wAxis.setRange(vmin, vmax);

    //Add cut planes. We create 3 default cut planes facing x, y, z directions.
    //The default location of cut plane is in the middle of the data.
    _plane.clear();
    addCutplane(CutPlane::X_AXIS, 0.5f);
    addCutplane(CutPlane::Y_AXIS, 0.5f);
    addCutplane(CutPlane::Z_AXIS, 0.5f);

    TRACE("Leave");
}

Volume::~Volume()
{
    TRACE("Enter");

    delete [] _data;
    delete _tex;
}

void Volume::setData(float *data, double vmin, double vmax, double nonZeroMin)
{
    int fcount = _width * _height * _depth * _numComponents;
    memcpy(_data, data, fcount * sizeof(float));
    _tex->update(_data);
    wAxis.setRange(vmin, vmax);
    _nonZeroMin = nonZeroMin;
    updatePending = true;
}

void Volume::getWorldSpaceBounds(Vector3f& bboxMin, Vector3f& bboxMax) const
{
    Vector3f scale = getPhysicalScaling();

    Matrix4x4d mat;
    mat.makeTranslation(_position);
    Matrix4x4d mat2;
    mat2.makeScale(scale);

    mat.multiply(mat2);

    bboxMin.set(FLT_MAX, FLT_MAX, FLT_MAX);
    bboxMax.set(-FLT_MAX, -FLT_MAX, -FLT_MAX);

    Vector3f modelMin(0, 0, 0);
    Vector3f modelMax(1, 1, 1);

    Vector4f bvert[8];
    bvert[0] = Vector4f(modelMin.x, modelMin.y, modelMin.z, 1);
    bvert[1] = Vector4f(modelMax.x, modelMin.y, modelMin.z, 1);
    bvert[2] = Vector4f(modelMin.x, modelMax.y, modelMin.z, 1);
    bvert[3] = Vector4f(modelMin.x, modelMin.y, modelMax.z, 1);
    bvert[4] = Vector4f(modelMax.x, modelMax.y, modelMin.z, 1);
    bvert[5] = Vector4f(modelMax.x, modelMin.y, modelMax.z, 1);
    bvert[6] = Vector4f(modelMin.x, modelMax.y, modelMax.z, 1);
    bvert[7] = Vector4f(modelMax.x, modelMax.y, modelMax.z, 1);

    for (int i = 0; i < 8; i++) {
        Vector4f worldVert = mat.transform(bvert[i]);
        if (worldVert.x < bboxMin.x) bboxMin.x = worldVert.x;
        if (worldVert.x > bboxMax.x) bboxMax.x = worldVert.x;
        if (worldVert.y < bboxMin.y) bboxMin.y = worldVert.y;
        if (worldVert.y > bboxMax.y) bboxMax.y = worldVert.y;
        if (worldVert.z < bboxMin.z) bboxMin.z = worldVert.z;
        if (worldVert.z > bboxMax.z) bboxMax.z = worldVert.z;
    }
}