source: trunk/packages/vizservers/nanovis/Volume.cpp @ 2852

/* -*- 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-2006  Purdue Research Foundation
 *
 *  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 "Volume.h"
#include "Trace.h"

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

Volume::Volume(float x, float y, float z,
               int w, int h, int d, float s, 
               int n, float *data,
               double v0, double v1, double nz_min) :
    aspect_ratio_width(1),
    aspect_ratio_height(1),
    aspect_ratio_depth(1),
    id(0),
    width(w),
    height(h),
    depth(d),
    size(s),
    pointsetIndex(-1),
    _tfPtr(NULL),
    _specular(6.),              // default value
    _diffuse(3.),               // default value
    _opacity_scale(10.),        // default value
    _name(NULL),
    _data(NULL),
    _n_components(n),
    _nonzero_min(nz_min),
    _tex(NULL),
    _location(x, y, z),
    _n_slices(512),             // default value
    _enabled(true),
    _data_enabled(true),        // default value
    _outline_enabled(true),     // default value
    _outline_color(1., 1., 1.), // default value
    _volume_type(CUBIC),        // default is a cubic volume
    _iso_surface(0)
{
    _tex = new Texture3D(w, h, d, GL_FLOAT, GL_LINEAR, n);
    int fcount = width * height * depth * _n_components;
    _data = new float[fcount];
    if (data != NULL) {
        TRACE("data is copied\n");
        memcpy(_data, data, fcount * sizeof(float));
        _tex->initialize(_data);
    } else {
        TRACE("data is null\n");
        memset(_data, 0, sizeof(width * height * depth * _n_components * 
                                sizeof(float)));
        _tex->initialize(_data);
    }

    id = _tex->id;

    wAxis.SetRange(v0, v1);

    // VOLUME
    aspect_ratio_width  = s * _tex->aspect_ratio_width;
    aspect_ratio_height = s * _tex->aspect_ratio_height;
    aspect_ratio_depth =  s * _tex->aspect_ratio_depth;

    //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();
    add_cutplane(1, 0.5f);
    add_cutplane(2, 0.5f);
    add_cutplane(3, 0.5f);

    //initialize the labels  
    label[0] = "X Label";
    label[1] = "Y Label";
    label[2] = "Z Label";

    TRACE("End -- Volume constructor\n");
}

Volume::~Volume()
{ 
    if (pointsetIndex != -1) {
        // TBD...
    }

    delete [] _data;
    delete _tex;
}