mitkBSPTree.h

/*=========================================================================

  Program:   3DMed
  Date:      $Date: 2014-02-25 18:30:00 +0800 $
  Version:   $Version: 4.6.0 $
  Copyright: MIPG, Institute of Automation, Chinese Academy of Sciences

=========================================================================*/


#ifndef __mitkBSPTree_h
#define __mitkBSPTree_h

template <class BSPTreeNode> 
class mitkBSPTree
{
public:
    float plane[4];
    BSPTreeNode *node;
    mitkBSPTree *front;
    mitkBSPTree *back;  

    mitkBSPTree() : node(NULL), front(NULL), back(NULL) {}
    ~mitkBSPTree() { delete node; node=NULL; delete front; front=NULL; delete back; back=NULL; }
};

template <class BSPTreeNode, class NodeProc>
inline void mitkBSPTreeTravelFrontToBack(mitkBSPTree<BSPTreeNode> *tree, float eyepos[3])
{
    if (tree->node)
    {
        NodeProc(tree->node);
    }
    else
    {
        float const *p = tree->plane;
        if (p[0]*eyepos[0]+p[1]*eyepos[1]+p[2]*eyepos[2]+p[3] <= 0)
        {
            if (tree->back) mitkBSPTreeTravelFrontToBack(tree->back, eyepos);
            if (tree->front) mitkBSPTreeTravelFrontToBack(tree->front, eyepos);
        }
        else
        {
            if (tree->front) mitkBSPTreeTravelFrontToBack(tree->front, eyepos);
            if (tree->back) mitkBSPTreeTravelFrontToBack(tree->back, eyepos);
        }
    }
}

template <class BSPTreeNode, class NodeProc>
inline void mitkBSPTreeTravelBackToFront(mitkBSPTree<BSPTreeNode> *tree, float eyepos[3])
{
    if (tree->node)
    {
        NodeProc(tree->node);
    }
    else
    {
        float const *p = tree->plane;
        if (p[0]*eyepos[0]+p[1]*eyepos[1]+p[2]*eyepos[2]+p[3] <= 0)
        {
            if (tree->back) mitkBSPTreeTravelBackToFront(tree->back, eyepos);
            if (tree->front) mitkBSPTreeTravelBackToFront(tree->front, eyepos);
        }
        else
        {
            if (tree->front) mitkBSPTreeTravelBackToFront(tree->front, eyepos);
            if (tree->back) mitkBSPTreeTravelBackToFront(tree->back, eyepos);
        }
    }
}

#endif

---