isosurface.h

Go to the documentation of this file.

/***************************************************************************
                         isosurface.h  -  description
                             -------------------
    begin                : Fri Mar 11 2005
    copyright            : (C) 2005-2006 by Ben Swerts
    email                : bswerts@users.sourceforge.net
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/


#ifndef ISOSURFACE_H
#define ISOSURFACE_H


// STL includes
#include <map>
#include <vector>
using std::map;
using std::vector;

// Xbrabo includes
#include <point3d.h>

class IsoSurface
{
  public:
        IsoSurface();                       // constructor
          ~IsoSurface();                      // destructor

          void setParameters(const std::vector<double>* values, const Point3D<unsigned int>& pointDimension, const Point3D<float>& pointDelta, const Point3D<float>& pointOrigin);         // set up the parameters for the surface 
          void addSurface(const double isoDensity); // calculates a new surface
    void changeSurface(const unsigned int surface, const double isoDensity);      // recalculates a surface

    bool densityPresent() const;          // returns whether a density has been loaded
    unsigned int numSurfaces() const;     // returns the number of calculated surfaces
    unsigned int numTriangles(const unsigned int surface) const;        // returns the number of triangles a certain surface consists of
          unsigned int numVertices(const unsigned int surface) const;         // returns the number of points a certain surface consists of
          void getTriangle(const unsigned int surface, const unsigned int index, Point3D<float>& point1, Point3D<float>& point2, Point3D<float>& point3, 
                     Point3D<float>& normal1, Point3D<float>& normal2, Point3D<float>& normal3) const;// return the data of a triangle of a surface    
    Point3D<float> getPoint(const unsigned int surface, const unsigned int index) const;    // returns the coordinates of a point on a surface
    void clearParameters();               // clear all data
    void clearSurfaces();                 // removes all existing surfaces
    void removeSurface(const unsigned int surface); // removes a certain surface
    Point3D<float> getOrigin() const;               // returns the set origin
    Point3D<float> getDelta() const;                // returns the set deltas
    Point3D<unsigned int> getNumPoints() const;     // returns the number of points in all directions

  private:
          struct Triangle
          {
      unsigned int pointID[3];
          };

    void calculateSurface(const double isoDensity); // does the basic surface calculation
    Point3D<float> intersection(const unsigned int x, const unsigned int y, const unsigned int z, const unsigned int edge);    // calculates the intersection
    Point3D<float> interpolate(const Point3D<float> point1, const Point3D<float> point2, const double var1, const double var2);   // linear interpolation
          unsigned int getEdgeID(const unsigned int x, const unsigned int y, const unsigned int z, const unsigned int edge);  // returns the ID of the edge 
          unsigned int getVertexID(const unsigned int x, const unsigned int y, const unsigned int z);     // returns the ID of the vertex 
          void renameVerticesAndTriangles(vector<Point3D<float> >* singleVerticesList, vector<unsigned int>* singleTriangleIndices);  // renames the vertices and triangles
    void calculateNormals(vector<float>* singleNormals, const unsigned int surface);        // calculates the normals
    unsigned int getArrayIndex(const unsigned int x, const unsigned int y, const unsigned int z) const;         // returns the index into the densityValues array

    vector<double> densityValues;         
    Point3D<unsigned int> numPoints;      
    Point3D<float> delta;                 
    Point3D<float> origin;                
    map<unsigned int, Point3D<float> > vertices;    
    vector<Triangle> triangles;           
    double currentIsoLevel;               
    vector<double> isoLevels;             
    vector< vector<Point3D<float> >* > verticesList;
    vector< vector<unsigned int>* > triangleIndices;
    vector< vector<float>* > normals;     

          static const unsigned int edgeTable[256];        
          static const int triTable[256][16];     

};

#endif

---