Marching Cubes

最新推荐文章于 2017-06-29 09:44:00 发布

iteye_15675

最新推荐文章于 2017-06-29 09:44:00 发布

阅读量145

点赞数

Marching Cubes
created: May 2000, updated: November 2001

Using the Open Inventor libraries I created a subclass ofSoEnginefor computing an isosurface usingMarching Cubes: A High Resolution 3D Surface Construction Algorithm, byWilliam Lorensenand Harvey E. Cline. This implementation isVERYdifferent from a previous version posted on this site. This version properly subclasses Open Inventor and efficiently creates an isosurface. The tables provided in the code found atPaul Bourke's pagePolygonising a scalar fieldwere used in this implementation.

The engine takes as input a scalar field and an isovalue. Then as output are the points, normals, and indexes into the points array used to create the triangles. The engine can easily be connected to other Inventor nodes to create anSoIndexedFaceSet. Below is an example of connecting the engine into a scene graph.

// create MarchingCubes object and connect fields
MarchingCubes *mcubes = new MarchingCubes();
mcubes->ref();
mcubes->data.setValue(dims, data);
mcubes->isoValue = 1.0;

// MarchingCubes gives the triangles with the
// vertices ordered clockwise
SoShapeHints *hints = new SoShapeHints();
hints->vertexOrdering.
  setValue(SoShapeHints::CLOCKWISE);
root->addChild(hints);

// connect the MarchingCubes points and
// normals to an SoVertexProperty node
SoVertexProperty *vprop = new SoVertexProperty();
vprop->ref();
vprop->vertex.connectFrom(&mcubes->points);
vprop->normal.connectFrom(&mcubes->normals);

// connect the indexes from MarchingCubes
SoIndexedFaceSet *faceSet = new SoIndexedFaceSet();
faceSet->vertexProperty = vprop;
faceSet->coordIndex.connectFrom(&mcubes->indexes);
root->addChild(faceSet);

Thanks to Gokhan Kisacikoglu, a vertex blending technique was added. The blending is nothing but a simple test to the closest voxel corner based on a percentage from the corner. Below is an example with the blending and without. Also a gif animation showing the blending can be foundhere.

Without Blending

With Blending

The blending reduces the number of triangles of a surface by 40% in some cases. When tested with the MRI data to produce the brain surface, the triangles were reduced from 1,065,244 to 636,702 without loss to the mesh structure. An interactive example of the blending technique can be found by downloading the3D GrapherOpen Inventor application.

TheMarchingCubesengine was used to create all the images on this page. Below are links to the MarchingCubes source code as well as an example which creates a sphere dataset and displays it with the engine node.(TheSFScalarFieldclass is needed by MarchingCubes.)

MarchingCubes.h
MarchingCubes.cpp
main.cpp
main.cpp.in
All code needed to compile example (including a configure script):mcubesdemo-1.0.tar.gz
All code needed to compile example on Irix (no configure script):mcubes.tar.gz

http://www.wheatchex.com/projects/mcubes/