VolumeMarchingCubes.h

// Simple Marching Cubes for VolumeViz
 
//
// This class is part of the Open Inventor Medical utility library.
//
// The medical utility classes are provided as a prebuilt library named
// "fei_inventor_medical", that can be used directly in an Open Inventor
// application. The classes in the prebuilt library are documented and
// supported by Thermo Fisher Scientific. These classes are also provided as source code.
//
// Please see $OIVHOME/include/Medical/InventorMedical.h for the full text.
//
 
#ifndef _VOLUME_MARCHING_CUBES_
#define _VOLUME_MARCHING_CUBES_
 
#include <Medical/InventorMedical.h>
 
#include <Inventor/nodes/SoTriangleSet.h>
#include <VolumeViz/nodes/SoVolumeData.h>
#include <LDM/SoLDMTileID.h>
#include <LDM/SoLDMDataAccess.h>
 
#include <vector>
 
class  SoLDMTopoOctree;
 
 
class INVENTORMEDICAL_API VolumeMarchingCubes
{
public:
  VolumeMarchingCubes();
 
  ~VolumeMarchingCubes();
 
  SoTriangleSet* getIsosurface( const SoVolumeData& volume, float isovalue );
 
  size_t getNumTriangles() const;
 
  float getIsovalue() const;
 
private:
  void marchVolume();
 
  void marchTile( SbVec3i32 tileDim, SbVec3f tileMin, SbVec3f tileMax, void* tileData );
 
  int marchCube( SbVec3f& corner, SbVec3f& size, float values[8], float isoval, SbVec3f vertices[15] );
 
  float getVoxelValue( int i, int j, int k);
 
  // Volume properties
  SbVec3i32 m_voldim;     // IJK dimensions
  SbBox3f   m_volext;     // XYZ extent
  SbVec3i32 m_tiledim;    // IJK dimensions
  SbVec3f   m_tileSize;   // XYZ size
  SbVec3f   m_voxelSize;  // XYZ size
  unsigned int        m_voxelBytes;
  SoDataSet::DataType m_voxelType;
 
  // Stuff to handle multiple tiles
  SoVolumeData*    m_volData;
  SoLDMDataAccess::DataInfo* m_dataInfo;  // Data for current locked tile (if any).
  SbVec3i32 m_tileMin;
  SbVec3i32 m_tileMax;
 
  // Only used during extraction.  Should be hidden in an impl class.
  std::vector<SbVec3f> m_vertices;
 
  // Other
  float m_isoValue;
  int   m_resolution;
  size_t m_numTriangles;
};
 
#endif