CVolvoxels.h

#ifndef __CVOLVOXELS_H__
#define __CVOLVOXELS_H__

#include <windows.h>
#include <math.h>
#include <string>
#include <iostream>
#include <fstream>
#include <vector>
#include <sstream>
#include "utilityfuncs.h"

class CVoxelValue
{
public:
    CVoxelValue() { density = 0; diffuseVal = 0; }
    CVoxelValue(float i, float j)
    {
        density    = i;
        diffuseVal = j;
    }

    float density;
    float diffuseVal;
};


class CVolVoxels
{
public:

    CVolVoxels()
    {
        _dimX = 0;
        _dimY = 0;
        _dimZ = 0;

        _pVoxelData = 0;
        _ptmpVD = 0;    
        _grad   = 0;
    }

    ~CVolVoxels()
    {
        Shutdown();
    }

    void Shutdown()
    {
        if ( _pVoxelData )
            delete[] _pVoxelData;

        if ( _grad )
            delete[] _grad;

        _pVoxelData = 0;
        _grad = 0;
        _dimX = 0;
        _dimY = 0;
        _dimZ = 0;
    }

    void Init( unsigned short * pVoxelData, unsigned short dimX, unsigned short dimY, unsigned short dimZ );

    unsigned long GetDimX() { return _dimX; }
    unsigned long GetDimY() { return _dimY; }
    unsigned long GetDimZ() { return _dimZ; }
    unsigned long GetBoundingDim() { return _boundingDim; }

    void ResetVolumePtr() { _ptmpVD = _pVoxelData; }

    void GotoSliceX(int iSliceNum) { _ptmpVD = _pVoxelData + iSliceNum; }
    void GotoSliceY(int iSliceNum) { _ptmpVD = _pVoxelData + _dimX * iSliceNum; }
    void GotoSliceZ(int iSliceNum) { _ptmpVD = _pVoxelData + _dimX * _dimY * iSliceNum; }

    void UseNextXSlice(int i)   { _ptmpVD = _pVoxelData + i; }   // i is the slice number
    void UseNextYSlice(int i)   { _ptmpVD = _pVoxelData + i * _dimX; }   // i is the slice number
    void UseNextZSlice(int i)   { _ptmpVD = _pVoxelData + i * _dimX * _dimY; }   // i is the slice number
    void UseNextZSlice(int i, int offset)   { _ptmpVD = _pVoxelData + i * _dimX;
                                              _ptmpVD += offset * _dimX * _dimY; }   // i is the slice number

    unsigned short NextX() { unsigned short val = (*_ptmpVD); _ptmpVD++; return val; }
    unsigned short NextY() { unsigned short val = (*_ptmpVD); _ptmpVD += _dimX; return val; }
        unsigned short NextZ() { unsigned short val = (*_ptmpVD); _ptmpVD += _dimX * _dimY; return val; }
    unsigned short GetVoxelValue(unsigned i, unsigned j, unsigned k)
    {
        //cout << endl;
        // cout << i << "  " << j << "  " << k << "  ";

        return (* (_pVoxelData +  _dimX * _dimY * k + _dimX * j + i) );
    }

    float CVolVoxels::GetDensityTriLinear(sPoint3f pos);

    unsigned short GetMaxValue();
    void ComputeGradient();
    void SetGradient(unsigned i, unsigned j, unsigned k, sVector3f grad)
    {
        (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).x = grad.x;
        (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).y = grad.y;
        (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).z = grad.z;
    }

    void GetGradient(unsigned i, unsigned j, unsigned k, sGrad * grad)
    {
        grad->x = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).x;
        grad->y = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).y;
        grad->z = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).z;
    }

    sGrad GetGradient(unsigned i, unsigned j, unsigned k)
    {
        sGrad grad;

        grad.x = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).x;
        grad.y = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).y;
        grad.z = (* (_grad +  _dimX * _dimY * k + _dimX * j + i) ).z;

        return grad;
    }

    void GetGradientTriLinear(float i, float j, float k, sGrad * pGrad);

private:
    unsigned long       _dimX, _dimY, _dimZ;        // dimensions of data read
    unsigned short  *   _pVoxelData;                        // pointer to memory block of data read, allocated by rawreader
    unsigned short  *   _ptmpVD;                    // volatile pointer into the 3D volume data

    unsigned long       _boundingDim;
    unsigned long       _MaxValue;
    unsigned long       _size;
    sGrad           *   _grad;


};

#endif /* __CVOLVOXELS_H__ */

---