An Octree Isosurface Codification Based On Discrete Planes
Imma Boada
Institut Informatics i Aplicacions
Universitat de Girona
Girona, Spain
imma@ima.udg.es
Isabel Navazo
Dept. Llenguatges i Sistemes Informatics
Universitat Politkcnica de Catalunya
Barcelona, Spain
isabel@lsi.upc.es
Abstract
ing on a single structure surface and volume data is still in
a developing stage.
The purpose of this work is to describe a method to code
a decimated model of an isosurjiace on an octree representation while maintaining volume data if it is needed. The proposed technique is based on grouping the Marching Cubes
(MC) patterns in five configurations according the topology and the number ofplanes of the sutface contained into
the cell. Moreover; a discrete number of planes, where the
surface lain on, is fixed. Starting from a complete volume
octree with the isosurface codijied at terminal nodes according to the new configurations, a bottom-up strategy is
taken in merging cells. Such strategy allows to implicitly
represent coplanar faces in upper octree levels without introducing error: At the end of this merging process, when
is required, a reconstruction strategy is applied to generate the surjiace contained in the octree intersected leaves.
Some examples with medical data demonstrate that reduction u p to 50 per cent on the the number of polygons could
be achieved.
The method we are going to present, is thejirst step on
the definition of an hybrid framework able to maintain multiresolution surface and volume data renderings. The model
allows to codify a surface on an octree data structure that,
different of previous octree based methods, maintains the
surface-volume data connection.
Based on the classification of original Marching Cubes
patterns in five classes, characterized by the number and
the position of the planes that define the surface contained
into an intersected cell, the surface is initially represented at
maximal division nodes of an octree. Then, a merging process is applied in a bottom-up traversal and sibling nodes,
which can be represented in the parent node, in one of the
five possible configurations, without introducing error, are
purged. At the end of this process, the surface is codified
in a set of nodes distributed at different levels of the octree.
To generate the surface, a reconstruction technique is applied at each terminal intersected node, achiving, in the best
situation, simplifications of almost 50% on the number of
polygons.
1. Introduction
The main advantages of this codification are that not
only reduces surface output fragmentation it also provides
a framework able to support multiresolution renderings, if
a n error is introduced in the merging process. Moreover,
for each isosurface terminal node its volume information is
maintained on its lower octree nodes.
Multiresolution has become one of the most promising optimization methods for resolving the conflict among
the increasing size of volume data sets and the capabilities of available hardware. Encoding in a single structure
a large number of different representations, multiresolution
schemes allow the extraction of variable resolution representations according to the user or the viewing parameters
specifications. Many different approaches have been proposed for the multiresolution management of surfaces (see
[6] for a survey). Multiresolution volume data management
has been also extendely analised, with methods based on
classical wavelets [26, 141, on multidimensional trees [29],
or more sophisticated methods that take advantage of texture hardware capabilities [ 12, 2, 301. However, despite the
improvements achieved in these areas the idea of combin-
0-7695-1215-l/Ol $lO.OO 0 2001 IEEE
The paper is organized as follows. After a brief description in Section 2 of previous work, the proposed surface
codification technique is introduced in Section 3. A complete description of the algorithm proposed to codify the
surface in the octree is given in Section 4. In Section 5, the
reconstruction algorithm applied to generate the surface is
described. Several results achieved with the proposed technique are shown in Section 6. Finally, concluding remarks
and future work are presented in Section 7.
130
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
2. Previous Work
codification of the MC patterns on the octree.
The Marching Cubes (MC) has become the best known
of the surface-based rendering techniques [ 1 I]. Since its
publication, in 1987, to overcome its main disadvantages
several improvements have been proposed. Algorithm's efficiency has been enhanced with strategies that avoid processing all volume data cells [3,4, 13,24,28]. Several techniques to solve ambiguities have been proposed [7,25, 161.
And a lot of strategies to reduce the large amount of generated polygons defining surfaces have been also presented
[8,9, 10, 221.
Despite the hight quality renderings generated with the
method, interactive renderings are difficult to be achieved.
A common solution to overcome this problem is the combination of the MC with surface simplification strategies. In
the literature, a variety of surface simplification techniques
have been proposed. These techniques can be classified into
two categories.
The first category, is characterized by the application of
a simplification process previously to the surface generation
phase. In general, strategies in this group exploit spatial coherence defining coarser representations of the original volume data, just averaging or subsampling neighboring vertices of a given grid into the vertices of a new coarser grid.
As the new grid has fewer cubes, the number of faces defining the surface decreases. However, the new surface may
show differences of the surface generated from the original
data. Such approaches, combined with adaptive isosurface
extraction techniques, are an efficient tool to achieve interactive exploration of volume models. Nevertheless, the surface reconstruction process becomes more complex, due to
the possible cracks that could appear on the generated surface [ 18, 23, 20, 271.
The second group, is characterized by the application
of the simplification process on the reconstructed surface,
commonly represented as a triangular mesh. Most of the
techniques in this group, simplify triangular meshes, either
by merging elements or by re-sampling vertices using different error criteria to measure the fitness of the approximated surfaces [22, 10, 61.
The surface codification we propose, in next section,
exploits the simplifications introduced in the Discretized
Marching Cubes (DiscMC). The DiscMC [15, 171, is a
derivation of the MC, that replaces the edge interpolation,
applied for the vertex polygon computation in the standard
MC, by a midpoint selection. This simplification reduces
the set of vertices that can be generated and implies that the
generated facets, that define the surface, can lie only on a
finite number of distinct plane incidences. The ability to
merge co-planar elements of this method achieves a high
simplification rate. In our approach this merging process is
based on the proposed configurations and the geometrical
3. The Octree Surface Codification
The hierarchical structure of octrees has lead this representation scheme a widespread used framework to maintain multiresolution on regular data sets. In the classical octree representation [21] nodes are labeled as black, white or
grey. Black nodes correspond to regions completely inside
the object. White nodes correspond to regions completely
outside the object. White and black nodes are leaves, i.e.
they are no further subdivided. Nodes containing part of the
surface are labeled as grey and are recursively subdivided.
Leaf grey nodes are called terminal grey (TG).
In the recent years, several extensions of the octree
model have been proposed. One of them is the Extended
Octree (EO) model [19]. The EO incorporates new terminal nodes that contain part of the surface object, yielding exact representations for polyhedra and reducing the required storage. In particular, the EO introduces as a terminal nodes: Face nodes (crossed by a single planar face
of the solid), Edge nodes (contain part of two neighboring
faces and a part of their common edge), Vertex nodes (contain a vertex of the polyhedron and a part of the faces and
edges converging to it) and Nearly Vertex nodes (contain
two or more faces that converge to a single vertex outside
the node).
Driven by the simplification introduced in the DiscMC
and by the original MC patterns, we propose a variation of
the EO to codify the surface. The specific terminal surface
nodes (TS) of the proposed octree codification are described
in next section.
3.1. Characterization of Terminal Surface Nodes
The selection of the terminal surface nodes has been
done according the number and the topology of the planes
that define the surface contained in an intersected cell (see
figure 2).
.
..
..
..
Figure 1. The surface contained in a cell can
codified as a Face, an Edge,a Band", a Band*
or a DoubleEdge.
Let consider P* = { p ~ , p l ., . , p , } the set of polygons
vertices obtained from the application of the DiscMC to an
intersected cell. According the position of these polygons
this cell can be classified as:
13 1
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
0
0
0
a Face, denoted as F , if exists a 7r plane such that, all
pi E P* vertices lye on it. Patterns 1,2,8 and 9 of the
original MC are classified as F cells (see figure 2).
an Edge, denoted as E , if exist two planes TO and 7r1
such that, given a point pi E P* one of the following
situations is given: (a) pi lies on the “0 plane; or (b)
pi lies on the 7r1 plane, or (c) pi lies on “ 0 and 7r1 (this
situation appears when the vertex lies on the T O ,7r1 intersection). Pattern 5 of the MC configurations corresponds to an Edge cell (see figure 2).
2
I
6
a Double Edge, denoted as D E , if exist three planes
T O ,7r1 and 7r2 such that, given a point pi E P* one of
the following situations is satisfied: (a)pz lies on the TO
plane, or (b)pi lies on the 7r1 plane, or (c) pi lies on the
7r2 plane, or (d) pi lies on TO and 7r1 (this situation appears when the vertex lies on the ~ 0 , 7 r 1intersection),
or (e) pi lies on 7r1 and 7r2 (this situation appears when
the vertex lies on the 7r1, 7r2 intersection). Patterns I 1
and 14 of figure 2 are examples of DoubleEdge cells.
in
12
0
a Band, denoted as B n ,if exist n plane equations, that
never intersect into the cell and given a point pi E P*
always lies in one and only one of these planes. The
number of planes, n, contained in a Bandn could be
2,3 or 4. Patterns 3, 4, 6, 7, 10 and 13 of figure 2
correspond to a Band cell configuration.
Figure 2. The patterns of the original MC.
1
I
a Special Band, denoted as B * ,if it can be represented
as the union of a Face cell and an Edge cell. In this
case, three plane equations are required to represent
the surface contained in the cell. Pattern 12 of figure 2
is an example of a Special Band.
ClassConfiguration
Face
Edge
DoubleEdge
Band
SpecialBand
Null
I MCpattern
I
I
1.2.8.9
5
11,14
I
3.4.6.7.10.
. . . . . 13
12
l o
Table 1. The correspondence between MC
patterns and the five possible configuration.
The correspondence between the MC patterns and the
five configurations is represented in Table 1. The non intersected node is denoted as Null.
Based on this classification we define as terminal surface
nodes those that can be codified as a Face, an Edge, a DoubleEdge, a Band or a SpecialBand (see figure I).
the number of plane incidences is limited, and thus, each Vi
could be represented as an integer pointer to a pre-computed
look-up-table that stores all possible plane incidences. This
integer pointer is also used to implicitly maintain the information of the inner and outer side of the object. The
codification of the surface is reduced to a set of integers.
3.2. Terminal Surface Nodes Codification
The surface of a terminal surface node 1zi is codified as
4. The Octree Construction Algorithm
The key point of the surface codification algorithm is
the application of a simplification strategy that reduces the
number of polygons defining the surface. This simplification is based on the detection and merging of coplanar polygons.
The algorithm starts with the surface (initially obtained
with the DiscMC), codified at maximal subdivision inter-
where K represents the class at which the node belongs to
(i.e. F, E , D E , Bn or B * )and (< Vo,do >, < Vl, d l >
, . . . , < V,, d, >) represents the set of plane equations that
di > corresponds
define the surface. Each couple < Vi,
to a 7ri plane where V , = ( A i ,Bi, Ci)is the normal vector and di is the plane distance to the node origin. Note
that, applying the simplification introduced by the DiscMC,
132
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
61
4567
5
0
I
I23
3
Figure 3. The reduction of a set of Face nodes to an Edge. The merging process starts grouping these
, ( n 2n3),
,
(n4,n5) and ( n 6n7).
,
At this first step the Plane Position test
descendents by pairs (n onl),
hasn't to be applied, the resultant configuration can be obtained from the MT table directly (see Table
2, where (N,
F ) -+ F ) . The resultant configurations { F ,F, F, F } are grouped again and evaluated. In
this case, the Plane Position test has to be applied to select if the pair ( F , F ) generates a Face, an
Edge or a Band". As planes intersected into the node, the generated codificationsare Edges. Finally,
(no123,n4567) is evaluated. In this case, the ( E ,E ) configurations generate an Edge.
sected nodes according to its configuration. On a bottomup traversal of this octree, coplanar polygons are detected
and merged taking into account the proposed terminal surface configurations. Starting at level n - 1 (i.e. parents
of maximal subdivion nodes) the eight descendants of each
intersected node ni are evaluated:
of resultant codificable configurations that can be generated
(see Table 2).
Given a pair of codificable nodes the resultant codification not always can be obtained directly from the table. See
for example the ( F ,F ) pair. In this case, three possible configurations can be obtained a F , an E or a B". To select
which of these possible codifications correspond to the analyzed pair, a Plane Position test has to be applied.
Ifni : Aix+Biy+Ciz = di and7r: : A:x+B:y+C,'z = d:,
are the plane equations defining the surface contained in
F nodes, respectively, the resultant codification after the
merging is: (i) F , if the planes are Coplanar; (ii) B2 if
the planes are Parallel; (iii) E if the planes are Incident and
its intersection falls within the parent node; (iv) B2 if the
planes are Incident and the intersection does not fall within
the parent node. Note that, as all the plane equations are
codified as pointers to a table of planes, the Plane Position test is reduced to a set of integer comparisons.
In figure 3, we have illustrated the application of the
..,n7}
merging process to a node ni, where {no,n1,n~,
and {N,F, NIF, N,F, NIF } represent its codificable
child nodes and their configurations.
if child nodes can be codified in the parent node, ni, as
a Face, an Edge, a DoubleEdge, a Band" or a Band',
descendants information is represented in ni, and these
child nodes are pruned. To represent the surface in the
parent node, the planes equations defining the surface
have to be updated.
conversely, if this codification is not possible the process stops and no more ancestors of the node are processed. In this case, the child nodes correspond to a TS
node.
Summarizing, if we define a codificable node as a node
that corresponds to one of the five defined TSN configurations, the simplification of child nodes in the parent, is possible if they can be reduced to a codificable node. Note that
the merging process always starts with a set of eight cod@cable nodes and the simplification is possible if they can be
reduced to a codificable node. All possible configurations
can be pre-computed and stored in a table to facilitate the
application of this process. This look-up-table is described
in next section.
5. Surface Reconstruction Algorithm
Once the octree has been codified if an explicit representation of the surface is needed, for rendering purposes
for example, a surface reconstruction algorithm has to be
applied.This algorithm could be decomposed in two steps.
The first, traverses in a top-down manner the octree and
selects all terminal intersected nodes. The set of selected
nodes is represented as S .
The second step applies a reconstruction process to each
node ni E S. Before the description of this process some
4.1. The Merging Table
The Merging Table (MT) determines when the merging
is possible and how has to be done. To generate the M T
we have considered all possible combinations of two codificable nodes and for each pair of configurations the set
133
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
@
;.
:
...
e0
(C)
B*
B*
B*
B*
NC
NC
NC
(b)
NC
Figure 4. The surface information of the node
is s(ni) = (E,7ro,7rl). (a) In a first step the
intersection of each plane with the node is
and
computed obtaining: Po = {po.pl,p~,p3}
PI = {pb.p:,pb,pi}. The edge vertices (eo,el)
are identified. (b) Point selection is done according the plane signs: T > 0 + p b , p ; are
eliminated ;T’ < 0 -+ p l , p 2 are eliminated. (c)
{po~p3relr
and {P;I&I e l , e o } polygons define the surface.
Table 2. Given a pair of nodes the MT Table
indicates the set of codificable configurations
that can be generated. NC represents a no
codificable configuration.
considerations have to be taken into account: (i) each plane
equation represented in ~ ( n ,has
) to be represented by one
polygon; (ii) this polygon is defined from the connection
o f all or some of the intersection points between the plane
and the edges of the node. (iii) the number of polygons required to define the surface varies with the configuration of
the node. In particular, the configuration-polygon relation
is F: I , E:2, DE:3, B”:n and B*:3.
To define the surface of an intersected node ni we apply
the three following steps:
follow the order established for the intersection point
computation. To guarantee the correct connection possible configurations of 4,s and 6 intersection points
have been evaluated and stored in a look-up-table (see
figure 4(c)).
The error involved in all the surface extraction process is
the error introduced by the DiscMC, i.e. an error of 1 / 2 of
the cell edge size.
This reconstruction algorithm can be extended to allow
multiresolution renderings. The first step, is substituted by
a selection process that selects nodes according to the distance to a focus of interest. The reconstruction step is also
modified to reduce the number of polygons that define the
surface of non-interesting nodes [ I ] .
Ititersection Point Coniputution. On a first step s(rii)
is evaluated and for each plane equation < Vi,rli >E
s ( n , ) (see Subsection 3.2), the intersection points between the plane and the edges of the node are computed. The limited number of plane incidences allows
the definition of a look-up-table in which the intersection points generated for each plane are stored. Thus,
the set of intersection points P
i of each plane in s ( n i )
will be obtained directly from this table just applying a translation respect the origin of the node. When
the node corresponds to an Edge, a DoubleEdge or a
Band* the points defining the generated edges are also
computed. (see figure 4(a))
6. Experimental Results
Four data models have been used for the tests: a CTvertebra of 128x128~80;a CT-head of 9 6 x 9 6 ~ 6 9 ;a CTscanned jaw of 1 2 8 x 1 2 8 ~ 4 0and a sphere of 643. The
method has been implemented on a Pentium I11 at 450MHz.
In our tests we have evaluated the degree of simplification that can be achieved with respect to standard MC and
the quality of the renderings. In table 3 the comparison of
the performance of the proposed simplification algorithm
and the original MC in terms o f number of polygons, number of intersection points defining the surface and surface
Point Selection. Once the intersection points are
known, a selection process is applied to determine
which of these points are members of the surface. This
selection depends of the node’s configuration and i t is
done evaluating the normal vector direction. Selected
vertices are represented in Pi” (see figure 4(b)).
Points Connection. Finally, P,’ vertices are connected.
To proper counter clock-wise ordering these points we
134
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
Figure 5. The CT-vertebra wire-mesh model before and after the simplification process.
extraction time, are presented. Best results are obtained on
the CT-vertebra (see figure 8). In this case the reduction on
the number of polygons is of a 50%. Despite the simplification on the number of polygons (see figure 5 ) the quality of
the renderings is comparable with the MC renderings. The
simplification on the CT-head only achieves a reduction of
15 (3% (see figure 7). While in the CT-jaw the reduction is
approximately of a 38% (see figure 6). The sphere achieves
a reduction of a 35%.
~
Data set
CT-vert. (MC)
CT-vert.*
CT-head (MC)
CT-head*
CT-jaw (MC)
CT-jaw*
sphere (MC)
sDhere*
The execution time is reasonable longer than that of the
MC, however, the reduction on the number of polygons will
reduce the complexity of posterior operations.
num.
polygons
13.597
7.501
33.480
28.990
38.491
23.656
24.464
15.720
num. inter.
points
52.755
28.637
120.629
104.033
144.470
86.357
78.936
47.732
time
339
33
68
8,4
89
9,83
2,7
2,9
Table 3. Experimental results. All times are in
seconds
We are conscious of the similarity with the DiscMC algorithm and that the degree of simplification of our method
is lower than the one obtained with the DiscMC. Despite
this similarity, note that in the DiscMC the face merging is
performed on all adjacent faces which lie on the same plane
while we are constrained to the eight descendants of an OCtree node. On the other hand, note that in [17] to give in
output a triangle-cased representation, it needs to triangulate these merged polygons; but these polygons might be
non convex, or can contain holes. Therefore, it adopts an
efficient solution for triangulating them on the fly (and thus
times are fast, in most cases slightly better than standard
MC). A disadvantage of this approach is that it usually produces in output very thin triangles, which can give problems
both in rendering and in further geometrical processing of
the mesh. Our approach, in contrast is based on a hierarchical approach so, by definition, the merged region should
be much more regular (even if in some cases of lower extension than the DiscMC); therefore also the triangulated
representation should posses a more regular shape.
Moreover, we have to keep in mind that our main objective is twofold: to reduce the number of polygons defining
the surface and to maintain the connection between surface
and volume data information.
Figure 6. a rendering of the CT-jaw.
135
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
7. Conclusions and Future Work
We have presented a method to codify a surface in an
octree data structure. The proposed surface coditication is
based on the classification o f the original MC patterns in
five configurations according the position and the number
of planes that define the surface contained in an intersected
node. Starting with the surface codified at maximal subdivision nodes, we have defined a simplification strategy
that detects coplanar polygons and represents them in upper
levels of the octree. Such simplification process achieves
reductions of almost a 50% on the number of facets with
respect to the standard MC in real medical data sets. The
algorithm to reconstruct the surfacc o f an intersected node
has been also described.
The main advantages of the proposed codification are
that not only reduces the output fragmentation it also maintains the connection between surface and volume data. The
method is the first step on the definition of an hybrid framework able to maintain multiresolution surtace and volume
data renderings. Our future work will be centered on niultiresolution surface fitting and volume hybrid renderings,
Figure 7. A rendering of the CT-head.
References
[ I ] 1.Boada and LNavazo: Multiresolution Isosurface Fit-
ting using an Octree based Surface Hieruchy. Research
Report IIiA 01-02-RR, Institut Infornihtica i Aplicacions, University of Girona, February 2001.
[ 2 ] I.B oada. I. Navazo and R .Sc o pi g n o . Multi re so I u t i on
Volume Visualization with Texture-based Octree. The
Visual Computer, Springer International. 200 1 (in
press).
[3] C.L. Bajaj, V. Pascucci and D.R. Schikore: Fast Isocontoruing tor improved interactivity. In 1996 ACM Symposium on Volume Visualization, pages 39-46. IEEE
Computer Society Press, Los Alamitos, CA, October
1996.
[4] P. Cignoni,P. Marno, C.Montani, E.Puppo and
R.Scopigno: Speeding up isosurface extraction using
interval trees. IEEE Transactions on Visualization and
Computer Graphics, 3(2): 158-170, 1997.
[5] P. Cignoni, C.Montani, E.Puppo and R.Scopigno: Optimal Isosurface Extraction from Irregular Volume Data.
In 1996 Symposium on Volume Visualization, pages
3 1-49. IEEE Computer Society, IEEE Computer Society Press. Los Alamitos, CA, October 1996.
[6] P.Cignoni, C. Montani and R.Scopigno. A comparison of mesh simplification algorithms. Computer and
Graphics 22 ( l ) , pp.37-54, 1998
Figlure 8. Different views of the CT-vertebra.
136
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
[7] M.J. Durst: Additional Reference to Marching Cubes
(letter). Computer Graphics 22(4), pages 72-74, 1988.
[21] H.Samet. Applications of Spatial Data Structures. Addison Wesley, Reading, MA, 1990.
[8] M. Deering: Geometry Compression. In Computer
Graphics (SIGGRAPH '95 Proceedings), pages 13-25,
1995.
[22] WSchroeder, J.A.Zarge and W.E. Lorensen: Decimation of Triangle Meshes. Proceedings SIGGRAPH 92,
65-70, 1992.
[9] F. Evans, S. Skiens and A. Varshney: Optimizing triangle strips for fast rendering. In IEEE Visualization
1996, pages 319-326, 1996.
[23] R.Shekhar, E.Fayyad, R.Yage1 and J.Cornhil1. Octree
based Decimation of Marching Cubes surfaces. Visualization 96, 335-342, 1996.
[IO] Hugues Hoppe: Progressive Meshes. Computer
Graphics (SIGGRAPH 96 Proceedings), 99- 108, 1996.
[24] H.W. Shen, C.Hansen, Y.Livnat and C.R.Johson: Isosurfacing in Span Space with Utmost Efficiency (ISSUE). Proceedings IEEE Visualization 96, 287-294,
1996.
[ l l ] W.Lorensen and H.Cline: Marching cubes a high
resolution 3D surface construction algorithm, ACM
Computer Graphics (Proceedings of SIGGRAPH '87),
v01.21, n 4, pp 163-170, 1987.
[25] A.Van Gelder and J. Wilhems: Topological Considerations on Isosurface Generation. ACM Transactions on
Graphics, 13(4). 337-375, October 1994.
[ 121 E.LaMar, B.Hamman and K.Joy. Multiresolution
Techniques for interactive texture-based volume visualization. In IEEE Visualization 99. IEEE CS Press,
pp.355-361 October 1999.
[26] R.Westemann. A Multiresolution Framework for Volume Rendering. In Proceedings of 1994 Sysmposium
on Volume Visualization, pp. 51-58. ACM Press October 17-18 1994
[I31 Y.Livnat, H.W. Shen and C.R. Johsons. A Near Optimal Isosurface Extraction Algorithm using the Span
Space. IEEE Transactions on Visualization and Computer Graphics, 2(1):73-84, 1996.
[27] R.Westemann, L.Kobbalt and T.Erl: Real-time exploration of Regular Volume Data by Adaptative Reconstruction of Isosurfaces. The Visual Computing 1998
[ 141 [Mur95] S.Muraki. Multiscale Volume Representation
[28] J. Wilhems and A. Van Gelder: Octrees for Faster
Isosurface generation. ACM Transactions on Graphics,
1 l(3). 201-227, July 1992.
by a dog wavelet. IEEE Trans. on Visualization and
Computer Graphics, 1(2):109:116, June 1995.
[ 151 C.Montani, R.Scateni and R.Scopigno: Discretized
[29] J.Wilhems and A. Van Gelder. Multi-dimensional
Trees for Controlled Volume Rendering and Compression. In proceedings of 1994 Symposium on Volume
Visualization, pp 27-34. ACM Press, October 17- 18
1994.
Marching Cubes, in Visualization '94 Proceedings,
R.D. Bergeron and A.E.Kaufman, Eds. 1994, pp.281287, IEEE Computer Society Press.
[ 161 C.Montani, R.Scateni and RScopigno: A Modified
Look-Up-Table for Implicit Disambiguation of Marching Cubes. The Visual Computer 1994, (10) pp.353355, IEEE Computer Society Press.
[30] M. Weiler, R. Westermann, C. Hansen, K. Zimmerman, T.Ertl: Level-of-Detail Volume Rendering via 3D
Textures. In Volume Visualization and Graphics Symposium 2000, pp.7- 13,2000.
[ 171 C.Montani, R.Scateni and RScopigno: Decreasing
Isosurface Complexity via Discrete Fitting, Computer
Aided Geometric Design, 17 (2000) 207-232.
[ 181 H.Mueller and M.Stark: Adaptive Generation of Sur-
face in Volume Data. The Visual Computer 9 (4), 182199.
[ 191 LNavazo: Extended Octree Representation of General
Solids with Plane Faces: Model Structure and Algorithms, Computer and Graphics, vol 13, l , 1989, pp.516.
[20] M.Ohlberger and M. Rumpf Hierarchical and Adaptative Visualization on Nested Grids. Computing ,
59(4), pages 269-285,1997
137
Authorized licensed use limited to: UNIVERSITAT DE GIRONA. Downloaded on April 27,2010 at 07:20:03 UTC from IEEE Xplore. Restrictions apply.
View publication stats