Adaptive Mesh Extraction using Simplification and Refinement
|
|
- Abner Barrett
- 6 years ago
- Views:
Transcription
1 Adaptive Mesh Extraction using Simplification and Refinement Adelailson Peixoto and Luiz Velho Abstract. This work presents a method for multiresolution mesh extraction with important mathematical properties. The generated mesh represents a regular surface from 3D Euclidian space. The input surface may be specified either implicitly or directly as a volumetric object. The method applies simplification operations to obtain a low resolution initial mesh and applies refinement operations to obtain a multiresolution representation from initial mesh. These combined operations provide nice properties to the mesh. 1. Introduction In several areas like medical visualization, geosciences and others, mesh extraction is crucial for many applications. In these areas, initial data may be represented as a volume data that encodes samples of a function over a grid. Each function value defines an iso-surface inside the volume that may be extracted as a mesh. Most existing algorithms for mesh extraction are classified as voxel based and as cross section based algorithms. The former generates the mesh by creating triangles inside voxels that intersect the iso-surface [14]. The latter generates triangles between consecutive pairs of parallel slices, where each slice contains a set of contours that represent the iso-surface [1,2,3,5,7,8,11]. In general, both methods generate meshes with constant resolution and uniform sampling. Because large volume are very common, these methods may result in meshes with millions of triangles. This may encumber several applications like rendering, denoising, network transmission and others. Multiresolution mesh representation helps to overcome constant resolution drawbacks [6]. It is defined as a set of meshes {M 0,..., M n } such that the number of vertices, edges and triangles of mesh M i increases Curve and Surface Fitting: Saint-Malo XXX (eds.), pp Copyright oc 2002 by Nashboro Press, Nashville, TN. ISBN XXX. All rights of reproduction in any form reserved.
2 2 A. Peixoto and L. Velho monotonically with i and there is a dependency relation between triangles at two subsequent levels i and i + 1. In general multiresolution mesh is obtained by simplification or by refinement methods. Simplification methods act over a supersampled initial mesh and coarsify it, creating lower resolution meshes [9,12,13,17]. Initial mesh may arise from 3D scanning techniques or from iso-surface extraction, like Marching Cubes (MC) [14]. One way to obtain the multiresolution is from a repeated decimation process applied to initial mesh, where vertices are iteratively removed. The problem is that most decimation algorithms produce highly irregular meshes, that are undesirable to many applications. An alternative to building multiresolution through decimation is remeshing, where the problem of an oversampled mesh is repaired through remeshing techniques. Although these methods are compelling, in many applications they seem to be inefficient and inelegant in the settings of iso-surface extraction from volumes. If one desires to extract an adaptively sampled hierarchical mesh, e.g., it is more natural that samples are adaptively added to the mesh in accord with the specific features of the underlying iso-surface. Refinement is applied over an initial coarse mesh and progressively computes higher resolution meshes [10,20]. In many applications it is more suitable than simplification: e.g. in network transmission it is more natural transmit first low resolution information, followed by details. So, combining simplification and refinement seems to be a good strategy to generate meshes with nice properties and available to more aplications. This work presents a combined simplification and refinement method such that the generated mesh presents nice properties (Section 2). 2. Algorithm Overview In our approach, the mesh generation is performed in two main steps: simplification and refinement. The simplification step generates a coarse mesh, called base mesh, containing the correct topology of the iso-surface. First, we generate a supersampled and structured set of points G, called the connectivity graph, which contains the geometry and topology of the iso-surface. Now we construct a set of disks over G such that the centers of the disks define vertices of the base mesh. From this set of disks we generate a set, called the Voronoi covering, that will be used to triangulate the base mesh. The simplification step is described in details in Section 3. The refinement step is applied to the base mesh to generate a mesh hierarchy. It consists of two operators that are applied simultaneously to the mesh: subdivision and adaptation of triangles. The former subdivides triangles, adding new vertices, edges and triangles to the mesh. The latter computes positions of new vertices and adjust them. Refinement is described in Section 4. Figure 1 shows a general scheme of our method. The proposed method presents a set of important properties: gener-
3 Mesh Extraction using Simplification and Refinement 3 Fig. 1. General Scheme. ated mesh exhibits correct topology; triangles have a good aspect ratios; meshes have a semi-regular structure; the method applies to surfaces with boundary and with multiple connected components; the mesh has a multiresolutionrepresentation; the mesh converges correctly to the iso-surface; refinement is adaptive; triangles subdivision is consistent. 3. Simplification This step is based on [10], where extracted multiresolution meshes produce triangles with a good aspect ratio. We employ a similiar process to extracted a corse mesh with constant resolution, instead of a multiresolution representation. Connectivity Graph A surfel is a surface patch formed by the intersection of the iso-surface with a voxel [10,20]. A node is the intersection of a surfel with a voxel edge. A c-vertex (candidate vertex) v is associated to each surfel such that the coordinates of v are the average of the node coordinates of the surfel associated to v (Figure 2a). Two c-vertices are neighbors if their associated surfels share at least one node. The connectivity graph G is defined by all c-vertices, and its edges are defined by pairs of neighbors c-vertices. The weight of each edge is the Euclidian distance between the c-vertices it connects. G represents the highest resolution of the iso-surface. Disk Covering Given a point p S, wheres is a regular surface, there exist an open neighborhood of p homeomorphic to a planar disk [4]. This neighorhood D S is called a disk over S. If for all p S there is a disk D over S with p D, we say that there is a disk covering C S associated to S. The goal is to construct a disk covering C G over the graph G. Toachievethe
4 4 A. Peixoto and L. Velho Fig. 2. a) Surfel and c-vertex. b) Disk constraints. c) Region classification. mesh properties, two topological constraints are imposed to C G (Figure 2b): first, C G is an independent disk covering, what means that no disk may contain the center of another disk. The second is that the radii of any two intersecting disks don t differ by more than a factor of two. The C G construction is made from discrete definitions of geodesic and disk. Ageodesicgeod u,v, u, v G, isacurveoverg composed by all c- vertices that minimize the distance between u and v. The length of geod u,v is denoted by leng u,v. A disk Dv, r with center v and radius r, isasetof c-vertices w G such that leng w,v r. The computation of a disk Dv r begins at the c-vertex v, and goes on, by a propagation process, until reaching all c-vertices w, withleng w,v r. The implementation of disk generation may be found in [10,16]. In order to assume that disk Dv r is homeomorphic to a planar disk, it is necessary to verify the occurrence of topological events during disk generation. A topological event occurs if the disk boundary self-intersects, or if the disk boundary vanishes. Let B G be the boundary of disk Dv r,andletu B be the closest c-vertex to v. Consider the sets V u = {x Dv x r is neighbor of u} and B u = {x B x is neighbor of u}. There is a boundary self-intersection if the set V u is not connected, and there is a boundary vanishing if B u = B. Inorderto maintain radius constraint when a topological event occurs, the radius of the disk must be halved. Let r max be the maximum radius allowed in disk covering C G. C G is constructed by classifying c-vertices of G in two types: forbidden or F, indicates c-vertices that cannot be taken as center of any disk, and allowable i or A i, indicates c-vertices that can be taken as center of a disk of radius r max /2 i (Figure 2c). Initially all c-vertices from G are classified as A 0. Suppose v is a A i c-vertex (v may be the center of a disk with radius r = r max /2 i ). But during generation of Dv r, topological event may occur, and the resulting disk Dv r will be halved until the event disappears, which means that the disk radius r r max /2 i. In order to maintain radius constraints, the region that surrounds a disk of radius r max /2 i must be classified as A i. An implementation of C G is shown in [10,16].
5 Mesh Extraction using Simplification and Refinement 5 Fig. 3. a) Incorrect Cylinder. b) Disk covering and Voronoi covering subdivision. Voronoi Covering Given a disk Dv Rv C G,aVoronoi cell VC G (v) is the set of c-vertices w G, such that for any disk Du Ru C G, leng w,v <leng w,u.avoronoi Covering VC G is the set of all Voronoi cells defined in C G. Two cells VC G(v) andv C G (u) areneighbors if there exist c-vertices x VC G(v) andy V C G (u) such that x and y are neighbors in G. The set of c-vertices x VC G(v) with the property that there exist y V C G(u) such that x and y are neighbors in G is called the u-boundary of VC G(v). Ac-vertexx VC G (v) is called a Voronoi vertex if it belongs to at least two separated boundaries of VC G (v). Voronoi vertices are very important to generating dual of Voronoi covering VC G. Sometimes, because of the amount of cells in VC G, it is not possible to reconstruct the base mesh correctly. In Figure 3a, the amount of disks (and cells) aren t sufficient to reconstruct a cylinder base mesh correctly. This problem didn t appear in [10] because C G was sufficiently refined to generate a final adapted mesh. To solve this problem here, we have not permitted the occurrence of violated cells in the Voronoi covering, i.e., cells with disconnected boundaries and cells with less than three boundaries. If they occur, we subdivide them: first we subdivide the disk covering, by halving disks that correspond to the violated cells and to their neighbor cells. When disks are halved, some forbidden c-vertices become A i and new disks must be created. After updating C G,weupdateVC G (Figure 3b). All this refinement process is applied locally to a neighborhood of violated cells, it means we don t need to recompute all cells of VC G. The subdivision proceeds until there are no more violated cells. [16] shows an implementation of VC G. Triangulation From the definition of a regular surface [4], it follows that the VC G has the same structure of a planar Voronoi diagram, so its dual is equivalent to the Delaunay triangulation. A triangulation may be computed by traversing all neighbors of each cell VC G (v) and generating all triangles sharing c-vertex v. To do this, neighbor cells of VC G (v) must be sorted. The simple projection of centers at a tangent plane is not a good way to sort cells: it may not reflect the correct order.
6 6 A. Peixoto and L. Velho Given a cell C cur with center c, we want to construct all triangles sharing c. For this, we must take a cell N cur with center n 1 and neighbor of C cur. We must compute a cell N nxt,neighborofc cur and of N cur simultaneously, with center n 2. The first triangle is c, n 1,n 2. The second triangle is computed by making N cur equals to N nxt and then computing anewn nxt, as before. We repeat this process until all triangles sharing c are created. It is important to take into account in the computation of N nxt the cases where C cur is on the boundary. Another important observation concerns the computation of holes: besides requiring N nxt to be simultaneous neighbor of C cur and N cur, we need the existence of neighboring Voronoi vertices c C cur,n 1 N cur and n 2 N nxt ; if there are no such Voronoi vertices, then C cur,n cur and N nxt form a hole. [16] shows an implementation of the triangulation. 4. Refinement Refinement is applied to a triangle T if its variance V T (d) >δ[20]. V T is computed over a discrete set of distances between T and the iso-surface: V T (d) = A(d 2 ) A(d), where A denotes the average of its argument. Refinement is applied simultaneously in two steps: subdivison (changes the mesh connectivity), and adaptation (changes the mesh geometry). Subdivision Triangle subdivision is based on variable resolution 4-8 meshes [19], where mesh connectivity is equivalent to [4.8 2 ] Laves tilings [18,19]. These tilings are defined by right isoscele triangles with vertices of valence 4 and 8, alternately. Its basic structure is a pair of such triangles that form a square as the basic block. The crucial advantage of this scheme is that it uses bisection as the basic refinement operation, rather than the more complex and commonly used face and vertex split. A bisection operation is applied to the internal edge of a basic block, subdividing the block into 4 triangles. The structure of 4-8 meshes is suitable for adaptive refinement because, even for non-uniform refinement, the mesh is still guaranteed to be conforming. A general framework for 4-8 meshes may be found in [18,19]. This refinement process permits the creation of semi-regular meshes. Adaptation When an edge e =(v 1,v 2 ) is subdivided by bisection, a new vertex v is added to the mesh. The coordinates of v must be computed such that triangles adapt to the iso-surface. The operations we ve applied to the mesh geometry are: moving vertices toward iso-surface, and reparameterizing vertices to improve mesh smoothness.
7 Mesh Extraction using Simplification and Refinement 7 Movement of vertex v (with coordinates p) towards the iso-surface S is performed by appling a local operator defined by F S (p) = T S (p) T S (p). This operator was inspired by [15], where T S is an implicit function. In our definition of F S, T S (p) is the signed distance between p and the isosurface, and T S (p) is the gradient that points in the normal direction of the iso-surface. If the subdivision edge e intersects the medial axis of the iso-surface, it is possible that the new vertex v will move in an incorrect direction. In the general case, this F S formulation is not able to decide the correct direction in which v moves. In [15], F S works well because it is guaranteed that vertices are in a neighborhood of the iso-surface, and subdivision edges don t intersect the medial axis. In [15] F S is applied to a constant resolution mesh. To move new vertices correctly we ve applied the geodesic geod v1,v2, instead of the F S operator: when the edge (v 1,v 2 ) is refined, the new vertex v is moved to the middle point of geod v1,v2. Points of geod v1,v2 may be computed like a disk generation process, where propagation begins at v 1 and proceeds until reaching v 2. If v 1 (or v 2 ) is on the boundary of the base mesh, it must be moved towards the boundary of the isosurface (connectivity graph). If v 1 and v 2 are both on the boundary, besides moving them towards the boundary of the iso-surface, we must compute geod v1,v2 on points restricted to the boundary. This process is implemented in [16]. After moving vertex v towards the iso-surface, a reparameterization is applied to it. This intends to maintain the mesh smoothness during the refinement process. Reparameterization may be defined as a tangential component F T of the Laplacian operator U, defined as [15,20]: U(p) = 1 n n i=1 (p i p), where p represents v coordinates and p i is are v neighbor is coordinates. The tangential component F T is defined as F T (p) =τ[u(p) (U(p). n) n], where τ is a constant and n is the normal to v. Application of F T to v optimizes triangles that share v, but may move v away from the iso-surface. As v is close to the iso-surface S, itmaybemovedback toward S by using F S operator. 5. Results, Conclusions and Future Work We have employed two error estimation functions. The first is an edge 1 deviation function, defined as ɛ e = U(e), where A(T )isthearea T A(T ) of each triangle T, U(e) =D(p) 2 (A(T 1 )+A(T 2 ))/2 is the error computed on edge e, D(p) is the distance from medial point p (of e) to the iso-surface and, T 1 and T 2 are triangles that share e. The second function estimates the deviation of the triangle normals 1 from T S. It is defined as [15] ɛ n = [A(T )(1 n(t ).m(t ) )], T e A(T ) T
8 8 A. Peixoto and L. Velho Fig. 4. Multiresolution mesh: base mesh (left), refined mesh (right). Fig. 5. Object Top: Multiresolution Method (left), Marching Cubes (right). Fig. 6. Normal and edge estimation and aspect ratio histogram. where m(t ) is the gradient computed at the centroid of T, and n(t ) is the unit normal to T. Figure 4 shows a base mesh generated from a synthetic volume and the multiresolution mesh. Figure 5 shows details of the top of the mesh, comparing our method with MC. Figure 6 shows the graphs of e and n applied to the above refined meshes in different resolutions: the graphs show that the mesh converges to the iso-surface during the refinement. Figure 6 (right) shows the aspect ratio histogram comparing our method with MC. We ve presented a method based on combined simplification and refinement operations for generating adaptive multiresolution meshes, with good properties. The method is mainly based on [10,20], and used important definitions and results of [15,19]. In [10] the main structures (DG, VCG and dual) are computed for
9 Mesh Extraction using Simplification and Refinement 9 every mesh extracted in a given resolution. Besides this disadvantage, there is no dependency relation between two generated meshes (there is no mesh hierarchy), and the mesh isn t semi-regular. Our method creates the main structures only once (base mesh), and creates a mesh hierarchy with a semi-regular structure. Some advantages of our work in relation to [20] are the capabilities of extracting meshes with boundary and with multiple connected components. Furthermore, we maintain efficiently and directly the consistence of meshes during the subdivision, due to bisection [19], while in [20] it is necessary to take into account some operations to maintain the mesh conformity. The combined simplification/refinement operations used by our method provide a set of important properties to the mesh (Section 2). This set of properties isn t provided by any other method. Future work will include the following: Definition of smooth operators that keep sharp features during refinement. Definition of criteria to decide when vertex displacement need not be computed by geodesics. After some refinement steps, displacement would be computed by the F S operator. Definition of criteria to choose disk centers. References 1. Bajaj, C., E. Coyle, and K. Link, Arbitrary topology shape reconstruction from planar cross sections, Graph. Models Image Process 58 (1996), Barequet, G., and M. Sharir, Piecewise-linear interpolation between polygonal slices,computer Vision and Image Undestanding 6 (1996), Boissonnat, J. D., Shape reconstruction from planar cross sections, Computer Vision, Graphics and Image Processing 44 (1988), Carmo,M.P.,Differential Geometry of Curves and Surfaces, Englewood Cliffs, N. J.: Prentice-Hall, Cong, G., and B. Parvin, An algebraic solution to surface recovery from cross-sectional contours, Graphical Models and Image Processing 61 (1999), Eck, M., T. DeRose, T. Duchamp, H. Hoppe, M. Lounsbery, and W. Stuertzle, Multiresolution analysis of arbitrary meshes, Proceedings of SIGGRAPH, 1995, Ekoule, A. B., F. C. Peyrin, and C. L. Odet, A triangulation algorithm from arbitrary shaped multiple planar contours, ACM Trans. on Graph 10 (1991),
10 10 A. Peixoto and L. Velho 8. Galin, E. and S. Akkouche, Fast surface reconstruction from contours using implicit surfaces, Implicit Surface 98. Eurographics and ACM SIGGRAPH, 1998, Garland, M. and P. S. Heckbert, Surface simplification using quadric error metrics, Proceedings of SIGGRAPH, 1996, Gavriliu, M., J. Carranza, D. Breen, and A. Barr, Fast extraction of adaptive multireslution meshes with guaranteed properties from volumetric data, Proceedings of IEEE Visualization, 2001, Gitlin, C., J. O Rourke, and V. Subramanian, On reconstructing polyhedral from parallel slices, International Journal of Comp. Geometry and Applications 6 (1996), Hoppe, H., Progressive meshes, Proceedings of SIGGRAPH, 1997, Lee, A. W. F., W. Sweldens, P. Schroder, L. Cowsar, and D. Dobkin, Maps: multiresolution adaptive parameterization of surfaces, Proceedings of SIGGRAPH, 1998, Lorensen, W. E., and H. E. Cline, Marching cubes: a high resolution 3D surface construction algorithm, Computer Graphics (Proceedings of SIGGRAPH) 21 (1987), Ohtake, Y., A. Belyaev, and A. Pasko, Dynamic meshes for accurate polygonization of implicit surfaces with sharp features, Shape Modeling International, 2001, Peixoto, A., Multiresolution adaptive mesh extraction from volumes, using simplification and refinement, Ph.D. Thesis, Department of Informatic, PUC-Rio- Catholic University of Rio de Janeiro, Shaffer, E., and M. Garland, Efficient adaptive simplification of massive meshes, IEEE Visualization, 2001, Velho, L., 4-8 subdivision, Computer-Aided Geometric Design 5 (2001), Velho, L., and J. Gomes, Variable resolution 4-k meshes: concepts and applications, Computer-Graphics Forum 3 (2000), Wood, Z. J., Semi-regular mesh extraction from volumes, Proceedings of Visualization, 2000, Adelailson Peixoto IMPA - Estrada D. Castorina, Rio de Janeiro - Brazil apeixoto@visgraf.impa.br Luiz Velho IMPA - Estrada D. Castorina, Rio de Janeiro - Brazil lvelho@visgraf.impa.br
Generalizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that g
Generalizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that generalizes the four-directional box spline of class
More informationA Robust Procedure to Eliminate Degenerate Faces from Triangle Meshes
A Robust Procedure to Eliminate Degenerate Faces from Triangle Meshes Mario Botsch, Leif P. Kobbelt Computer Graphics Group, RWTH Aachen, kobbelt,botsch @cs.rwth-aachen.de Abstract When using triangle
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationMesh Decimation. Mark Pauly
Mesh Decimation Mark Pauly Applications Oversampled 3D scan data ~150k triangles ~80k triangles Mark Pauly - ETH Zurich 280 Applications Overtessellation: E.g. iso-surface extraction Mark Pauly - ETH Zurich
More informationMesh Decimation Using VTK
Mesh Decimation Using VTK Michael Knapp knapp@cg.tuwien.ac.at Institute of Computer Graphics and Algorithms Vienna University of Technology Abstract This paper describes general mesh decimation methods
More information3-Dimensional Object Modeling with Mesh Simplification Based Resolution Adjustment
3-Dimensional Object Modeling with Mesh Simplification Based Resolution Adjustment Özgür ULUCAY Sarp ERTÜRK University of Kocaeli Electronics & Communication Engineering Department 41040 Izmit, Kocaeli
More informationA Developer s Survey of Polygonal Simplification algorithms. CS 563 Advanced Topics in Computer Graphics Fan Wu Mar. 31, 2005
A Developer s Survey of Polygonal Simplification algorithms CS 563 Advanced Topics in Computer Graphics Fan Wu Mar. 31, 2005 Some questions to ask Why simplification? What are my models like? What matters
More informationAdaptive Semi-Regular Remeshing: A Voronoi-Based Approach
Adaptive Semi-Regular Remeshing: A Voronoi-Based Approach Aymen Kammoun 1, Frédéric Payan 2, Marc Antonini 3 Laboratory I3S, University of Nice-Sophia Antipolis/ CNRS (UMR 6070) - France 1 kammoun@i3s.unice.fr
More informationMesh Repairing and Simplification. Gianpaolo Palma
Mesh Repairing and Simplification Gianpaolo Palma Mesh Repairing Removal of artifacts from geometric model such that it becomes suitable for further processing Input: a generic 3D model Output: (hopefully)a
More informationcoding of various parts showing different features, the possibility of rotation or of hiding covering parts of the object's surface to gain an insight
Three-Dimensional Object Reconstruction from Layered Spatial Data Michael Dangl and Robert Sablatnig Vienna University of Technology, Institute of Computer Aided Automation, Pattern Recognition and Image
More informationLevel Set Extraction from Gridded 2D and 3D Data
Level Set Extraction from Gridded 2D and 3D Data David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationApplications. Oversampled 3D scan data. ~150k triangles ~80k triangles
Mesh Simplification Applications Oversampled 3D scan data ~150k triangles ~80k triangles 2 Applications Overtessellation: E.g. iso-surface extraction 3 Applications Multi-resolution hierarchies for efficient
More informationDual/Primal Mesh Optimization for Polygonized Implicit Surfaces
Dual/Primal Mesh Optimization for Polygonized Implicit Surfaces Yutaka Ohtake y and Alexander G. Belyaev y;z y Computer Graphics Group, Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany z
More informationSurface Reconstruction. Gianpaolo Palma
Surface Reconstruction Gianpaolo Palma Surface reconstruction Input Point cloud With or without normals Examples: multi-view stereo, union of range scan vertices Range scans Each scan is a triangular mesh
More informationDual-Primal Mesh Optimization for Polygonized Implicit Surfaces with Sharp Features
Dual-Primal Mesh Optimization for Polygonized Implicit Surfaces with Sharp Features Yutaka Ohtake Ý and Alexander G. Belyaev Ý Þ Ý Computer Graphics Group, Max-Planck-Institut für Informatik, 66123 Saarbrücken,
More informationScalar Algorithms: Contouring
Scalar Algorithms: Contouring Computer Animation and Visualisation Lecture tkomura@inf.ed.ac.uk Institute for Perception, Action & Behaviour School of Informatics Contouring Scaler Data Last Lecture...
More informationVoronoi Diagram. Xiao-Ming Fu
Voronoi Diagram Xiao-Ming Fu Outlines Introduction Post Office Problem Voronoi Diagram Duality: Delaunay triangulation Centroidal Voronoi tessellations (CVT) Definition Applications Algorithms Outlines
More informationOutline. Reconstruction of 3D Meshes from Point Clouds. Motivation. Problem Statement. Applications. Challenges
Reconstruction of 3D Meshes from Point Clouds Ming Zhang Patrick Min cs598b, Geometric Modeling for Computer Graphics Feb. 17, 2000 Outline - problem statement - motivation - applications - challenges
More informationGeometry Processing & Geometric Queries. Computer Graphics CMU /15-662
Geometry Processing & Geometric Queries Computer Graphics CMU 15-462/15-662 Last time: Meshes & Manifolds Mathematical description of geometry - simplifying assumption: manifold - for polygon meshes: fans,
More informationA Short Survey of Mesh Simplification Algorithms
A Short Survey of Mesh Simplification Algorithms Jerry O. Talton III University of Illinois at Urbana-Champaign 1. INTRODUCTION The problem of approximating a given input mesh with a less complex but geometrically
More informationGeometric Modeling in Graphics
Geometric Modeling in Graphics Part 10: Surface reconstruction Martin Samuelčík www.sccg.sk/~samuelcik samuelcik@sccg.sk Curve, surface reconstruction Finding compact connected orientable 2-manifold surface
More informationMultiresolution Remeshing Using Weighted Centroidal Voronoi Diagram
Multiresolution Remeshing Using Weighted Centroidal Voronoi Diagram Chao-Hung Lin 1, Chung-Ren Yan 2, Ji-Hsen Hsu 2, and Tong-Yee Lee 2 1 Dept. of Geomatics, National Cheng Kung University, Taiwan 2 Dept.
More information2 Coarse Mesh Extraction
Semi-Regular Mesh Extraction from Volumes Zoë J. Wood Mathieu Desbrun Peter Schröder David Breen Abstract We present a novel method to extract iso-surfaces from distance volumes. It generates high quality
More informationGeometric Modeling. Mesh Decimation. Mesh Decimation. Applications. Copyright 2010 Gotsman, Pauly Page 1. Oversampled 3D scan data
Applications Oversampled 3D scan data ~150k triangles ~80k triangles 2 Copyright 2010 Gotsman, Pauly Page 1 Applications Overtessellation: E.g. iso-surface extraction 3 Applications Multi-resolution hierarchies
More informationVolume Illumination and Segmentation
Volume Illumination and Segmentation Computer Animation and Visualisation Lecture 13 Institute for Perception, Action & Behaviour School of Informatics Overview Volume illumination Segmentation Volume
More informationPolygonization of Implicit Surfaces
Polygonization of Implicit Surfaces Hongxin Zhang and Jieqing Feng 2007-01-11 State Key Lab of CAD&CG Zhejiang University Contents Polygonization of Implicit Surfaces Other Methods for Displaying Implicit
More informationTechnical Report. Mesh Coarsening based upon Multiplicatively Weighted Centroidal Voronoi Diagrams. Iurie Chiosa and Andreas Kolb
Computer Graphics and Multimedia Systems Technical Report Mesh Coarsening based upon Multiplicatively Weighted Centroidal Voronoi Diagrams Iurie Chiosa and Andreas Kolb Version: 1.0 (24.04.2006) Faculty
More informationParameterization of Triangular Meshes with Virtual Boundaries
Parameterization of Triangular Meshes with Virtual Boundaries Yunjin Lee 1;Λ Hyoung Seok Kim 2;y Seungyong Lee 1;z 1 Department of Computer Science and Engineering Pohang University of Science and Technology
More informationHierarchical Iso-Surface Extraction
Hierarchical Iso-Surface Extraction Ulf Labsik Kai Hormann Martin Meister Günther Greiner Computer Graphics Group University of Erlangen Multi-Res Modeling Group Caltech Computer Graphics Group University
More information2) For any triangle edge not on the boundary, there is exactly one neighboring
Triangulating Trimmed NURBS Surfaces Chang Shu and Pierre Boulanger Abstract. This paper describes techniques for the piecewise linear approximation of trimmed NURBS surfaces. The problem, called surface
More informationMoving Least Squares Multiresolution Surface Approximation
Moving Least Squares Multiresolution Surface Approximation BORIS MEDEROS LUIZ VELHO LUIZ HENRIQUE DE FIGUEIREDO IMPA Instituto de Matemática Pura e Aplicada Estrada Dona Castorina 110, 22461-320 Rio de
More informationSurface Simplification Using Quadric Error Metrics
Surface Simplification Using Quadric Error Metrics Authors: Michael Garland & Paul Heckbert Presented by: Niu Xiaozhen Disclaimer: Some slides are modified from original slides, which were designed by
More informationSurvey of Surface Reconstruction Algorithms
Journal of Signal and Information Processing, 2014, 5, 63-79 Published Online August 2014 in SciRes. http://www.scirp.org/journal/jsip http://dx.doi.org/10.4236/jsip.2014.53009 Survey of Surface Reconstruction
More informationMultiresolution Meshes. COS 526 Tom Funkhouser, Fall 2016 Slides by Guskov, Praun, Sweldens, etc.
Multiresolution Meshes COS 526 Tom Funkhouser, Fall 2016 Slides by Guskov, Praun, Sweldens, etc. Motivation Huge meshes are difficult to render store transmit edit Multiresolution Meshes! [Guskov et al.]
More informationFinal Project, Digital Geometry Processing
Final Project, Digital Geometry Processing Shayan Hoshyari Student #: 81382153 December 2016 Introduction In this project an adaptive surface remesher has been developed based on the paper [1]. An algorithm
More informationDigital Geometry Processing. Computer Graphics CMU /15-662
Digital Geometry Processing Computer Graphics CMU 15-462/15-662 Last time: Meshes & Manifolds Mathematical description of geometry - simplifying assumption: manifold - for polygon meshes: fans, not fins
More informationVolume Illumination, Contouring
Volume Illumination, Contouring Computer Animation and Visualisation Lecture 0 tkomura@inf.ed.ac.uk Institute for Perception, Action & Behaviour School of Informatics Contouring Scaler Data Overview -
More informationCS354 Computer Graphics Surface Representation IV. Qixing Huang March 7th 2018
CS354 Computer Graphics Surface Representation IV Qixing Huang March 7th 2018 Today s Topic Subdivision surfaces Implicit surface representation Subdivision Surfaces Building complex models We can extend
More informationFreeform Shape Representations for Efficient Geometry Processing
Freeform Shape Representations for Efficient Geometry Processing Leif Kobbelt Mario Botsch Computer Graphics Group RWTH Aachen Abstract The most important concepts for the handling and storage of freeform
More informationSemi-Regular 4 8 Refinement and Box Spline Surfaces
Semi-Regular 4 8 Refinement and Box Spline Surfaces LUIZ VELHO IMPA Instituto de Matemática Pura e Aplicada, lvelho@visgraf.impa.br Abstract. In this paper we introduce a new mesh refinement method for
More informationMulti-View Matching & Mesh Generation. Qixing Huang Feb. 13 th 2017
Multi-View Matching & Mesh Generation Qixing Huang Feb. 13 th 2017 Geometry Reconstruction Pipeline RANSAC --- facts Sampling Feature point detection [Gelfand et al. 05, Huang et al. 06] Correspondences
More informationGeometric Modeling. Bing-Yu Chen National Taiwan University The University of Tokyo
Geometric Modeling Bing-Yu Chen National Taiwan University The University of Tokyo What are 3D Objects? 3D Object Representations What are 3D objects? The Graphics Process 3D Object Representations Raw
More informationSurface Curvature Estimation for Edge Spinning Algorithm *
Surface Curvature Estimation for Edge Spinning Algorithm * Martin Cermak and Vaclav Skala University of West Bohemia in Pilsen Department of Computer Science and Engineering Czech Republic {cermakm skala}@kiv.zcu.cz
More informationSpace Subdivision for the Adaptive Edge Spinning Polygonization
Space Subdivision for the Adaptive Edge Spinning Polygonization MARTIN CERMAK 1, VACLAV SKALA Department of Computer Science and Engineering University of West Bohemia in Pilsen Univerzitni 8, 306 14 Plzen
More informationCloning Skeleton-driven Animation to Other Models
Cloning Skeleton-driven Animation to Other Models Wan-Chi Luo Jian-Bin Huang Bing-Yu Chen Pin-Chou Liu National Taiwan University {maggie, azar, toby}@cmlab.csie.ntu.edu.tw robin@ntu.edu.tw Abstract-3D
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationMedial Scaffolds for 3D data modelling: status and challenges. Frederic Fol Leymarie
Medial Scaffolds for 3D data modelling: status and challenges Frederic Fol Leymarie Outline Background Method and some algorithmic details Applications Shape representation: From the Medial Axis to the
More informationSpecial Topics in Visualization
Special Topics in Visualization Final Project Report Dual contouring of Hermite Data Submitted By S M Shahed Nejhum 8589-1199 May 19, 2008 Introduction Iso-surface extraction from 3D volumetric data is
More informationEfficient Representation and Extraction of 2-Manifold Isosurfaces Using kd-trees
Efficient Representation and Extraction of 2-Manifold Isosurfaces Using kd-trees Alexander Greß and Reinhard Klein University of Bonn Institute of Computer Science II Römerstraße 164, 53117 Bonn, Germany
More informationCurves & Surfaces. Last Time? Progressive Meshes. Selective Refinement. Adjacency Data Structures. Mesh Simplification. Mesh Simplification
Last Time? Adjacency Data Structures Curves & Surfaces Geometric & topologic information Dynamic allocation Efficiency of access Mesh Simplification edge collapse/vertex split geomorphs progressive transmission
More informationGeometric Representations. Stelian Coros
Geometric Representations Stelian Coros Geometric Representations Languages for describing shape Boundary representations Polygonal meshes Subdivision surfaces Implicit surfaces Volumetric models Parametric
More informationIntroduction to Computer Graphics. Modeling (3) April 27, 2017 Kenshi Takayama
Introduction to Computer Graphics Modeling (3) April 27, 2017 Kenshi Takayama Solid modeling 2 Solid models Thin shapes represented by single polygons Unorientable Clear definition of inside & outside
More informationScalar Visualization
Scalar Visualization Visualizing scalar data Popular scalar visualization techniques Color mapping Contouring Height plots outline Recap of Chap 4: Visualization Pipeline 1. Data Importing 2. Data Filtering
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationing and enhancing operations for shapes represented by level sets (isosurfaces) and applied unsharp masking to the level set normals. Their approach t
Shape Deblurring with Unsharp Masking Applied to Mesh Normals Hirokazu Yagou Λ Alexander Belyaev y Daming Wei z Λ y z ; ; Shape Modeling Laboratory, University of Aizu, Aizu-Wakamatsu 965-8580 Japan fm50534,
More informationx1_nearest x2_nearest Implicit Surface:[ f(x)=0,b(x)= false ]
Marching Triangles: Delaunay Implicit Surface Triangulation A.Hilton 1 and J.Illingworth Vision, Speech and Signal Processing Group, Department of Electronic and Electrical Engineering University of Surrey,
More informationProcessing 3D Surface Data
Processing 3D Surface Data Computer Animation and Visualisation Lecture 12 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing
More informationLAPLACIAN MESH SMOOTHING FOR TETRAHEDRA BASED VOLUME VISUALIZATION 1. INTRODUCTION
JOURNAL OF MEDICAL INFORMATICS & TECHNOLOGIES Vol.4/2002, ISSN 642-6037 Rafał STĘGIERSKI *, Paweł MIKOŁAJCZAK * volume data,triangle mesh generation, mesh smoothing, marching tetrahedra LAPLACIAN MESH
More informationShape from LIDAR Data. Univ. of Florida
Shape from LIDAR Data Yusuf Sahillioğlu Alper Üngör Univ. of Florida 1. Introduction LIght Detection And Ranging systems, LIDAR, are capable of acquiring data to produce accurate digital elevation models
More informationMotivation. Freeform Shape Representations for Efficient Geometry Processing. Operations on Geometric Objects. Functional Representations
Motivation Freeform Shape Representations for Efficient Geometry Processing Eurographics 23 Granada, Spain Geometry Processing (points, wireframes, patches, volumes) Efficient algorithms always have to
More informationExtracting consistent and manifold interfaces from multi-valued volume data sets
Extracting consistent and manifold interfaces from multi-valued volume data sets Stephan Bischoff, Leif Kobbelt Lehrstuhl für Informatik 8, RWTH Aachen, 52056 Aachen Email: {bischoff,kobbelt}@informatik.rwth-aachen.de
More informationNormal Mesh Compression
Normal Mesh Compression Andrei Khodakovsky Caltech 549B (e:54, p:45db) 1225B (e:20, p:54db) Igor Guskov Caltech 3037B (e:8.1, p:62db) 18111B (e:1.77, p:75db) original Figure 1: Partial reconstructions
More informationTexture Mapping with Hard Constraints
EUROGRAPHICS 2001 / A.Chalmers and T.-M.Rhyne Volume 20 (2001), Number 3 (Guest Editors) Texture Mapping with Hard Constraints Ilya Eckstein Vitaly Surazhsky Craig Gotsman Computer Science Department,
More informationFast Stellar Mesh Simplification
Fast Stellar Mesh Simplification ANTÔNIO W. VIEIRA 1,2, LUIZ VELHO 3, HÉLIO LOPES 1, GEOVAN TAVARES 1 AND THOMAS LEWINER 1,4 1 PUC Rio Laboratório Matmídia Rio de Janeiro 2 UNIMONTES CCET Montes Claros
More informationPoint-based Simplification Algorithm
Point-based Simplification Algorithm Pai-Feng Lee 1, Bin-Shyan Jong 2 Department of Information Management, Hsing Wu College 1 Dept. of Information and Computer Engineering Engineering, Chung Yuan Christian
More informationLecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More informationA General Simplification Algorithm
A General Simplification Algorithm Boštjan Pivec, Vid Domiter Abstract In this article a new general algorithm for triangular mesh simplification is proposed. The algorithm extends Krivograd's work from
More informationContours & Implicit Modelling 1
Contouring & Implicit Modelling Visualisation Lecture 8 Institute for Perception, Action & Behaviour School of Informatics Contours & Implicit Modelling 1 Brief Recap Contouring Implicit Functions lecture
More informationSubdivision Surfaces. Course Syllabus. Course Syllabus. Modeling. Equivalence of Representations. 3D Object Representations
Subdivision Surfaces Adam Finkelstein Princeton University COS 426, Spring 2003 Course Syllabus I. Image processing II. Rendering III. Modeling IV. Animation Image Processing (Rusty Coleman, CS426, Fall99)
More informationGeodesic Paths on Triangular Meshes
Geodesic Paths on Triangular Meshes DIMAS MARTÍNEZ, LUIZ VELHO, PAULO CEZAR CARVALHO IMPA Instituto Nacional de Matemática Pura e Aplicada - Estrada Dona Castorina, 110, 22460-320 Rio de Janeiro, RJ, Brasil
More informationContours & Implicit Modelling 4
Brief Recap Contouring & Implicit Modelling Contouring Implicit Functions Visualisation Lecture 8 lecture 6 Marching Cubes lecture 3 visualisation of a Quadric toby.breckon@ed.ac.uk Computer Vision Lab.
More informationA Global Laplacian Smoothing Approach with Feature Preservation
A Global Laplacian Smoothing Approach with Feature Preservation hongping Ji Ligang Liu Guojin Wang Department of Mathematics State Key Lab of CAD&CG hejiang University Hangzhou, 310027 P.R. China jzpboy@yahoo.com.cn,
More informationAn Adaptive Subdivision Scheme On Composite Subdivision Meshes
An Adaptive Subdivision Scheme On Composite Subdivision Meshes Anh-Cang PHAN VinhLong College of Economics and Finance VinhLong, Vietnam pacang@vcef.edu.vn Romain RAFFIN Aix-Marseille University CNRS,
More informationManifold Parameterization
Manifold Parameterization Lei Zhang 1,2, Ligang Liu 1,2, Zhongping Ji 1,2, and Guojin Wang 1,2 1 Department of Mathematics, Zhejiang University, Hangzhou 310027, China 2 State Key Lab of CAD&CG, Zhejiang
More information03 - Reconstruction. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Spring 17 - Daniele Panozzo
3 - Reconstruction Acknowledgements: Olga Sorkine-Hornung Geometry Acquisition Pipeline Scanning: results in range images Registration: bring all range images to one coordinate system Stitching/ reconstruction:
More informationIndirect Volume Rendering
Indirect Volume Rendering Visualization Torsten Möller Weiskopf/Machiraju/Möller Overview Contour tracing Marching cubes Marching tetrahedra Optimization octree-based range query Weiskopf/Machiraju/Möller
More informationA General Algorithm for Triangular Meshes Simplification
Proceedings of the 11th WSEAS International Conference on COMPUTERS, Agios Nikolaos, Crete Island, Greece, July 26-28, 2007 613 A General Algorithm for Triangular Meshes Simplification BOŠTJAN PIVEC University
More informationInterpolating and Approximating Implicit Surfaces from Polygon Soup
Interpolating and Approimating Implicit Surfaces from Polygon Soup Chen Shen, James F. O Brien, Jonathan R. Shewchuk University of California, Berkeley Geometric Algorithms Seminar CS 468 Fall 2005 Overview
More information3D Modeling techniques
3D Modeling techniques 0. Reconstruction From real data (not covered) 1. Procedural modeling Automatic modeling of a self-similar objects or scenes 2. Interactive modeling Provide tools to computer artists
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationSubdivision Surfaces. Homework 1: Questions on Homework? Last Time? Today. Tensor Product. What s an illegal edge collapse?
Homework 1: Questions/Comments? Subdivision Surfaces Questions on Homework? Last Time? What s an illegal edge collapse? Curves & Surfaces Continuity Definitions 2 3 C0, G1, C1, C 1 a b 4 Interpolation
More informationTopology-Free Cut-and-Paste Editing over Meshes
Topology-Free Cut-and-Paste Editing over Meshes Hongbo Fu, Chiew-Lan Tai, Hongxin Zhang Department of Computer Science Hong Kong University of Science and Technology Abstract Existing cut-and-paste editing
More informationIntrinsic Morphing of Compatible Triangulations. VITALY SURAZHSKY CRAIG GOTSMAN
International Journal of Shape Modeling Vol. 9, No. 2 (2003) 191 201 c World Scientific Publishing Company Intrinsic Morphing of Compatible Triangulations VITALY SURAZHSKY vitus@cs.technion.ac.il CRAIG
More informationGeodesic Paths on Triangular Meshes
Geodesic Paths on Triangular Meshes Dimas Martínez Luiz Velho Paulo Cezar Carvalho IMPA Instituto Nacional de Matemática Pura e Aplicada Estrada Dona Castorina, 110, 22460-320 Rio de Janeiro, RJ, Brasil
More informationCS 523: Computer Graphics, Spring Differential Geometry of Surfaces
CS 523: Computer Graphics, Spring 2009 Shape Modeling Differential Geometry of Surfaces Andrew Nealen, Rutgers, 2009 3/4/2009 Recap Differential Geometry of Curves Andrew Nealen, Rutgers, 2009 3/4/2009
More informationSubdivision Surfaces. Homework 1: Questions/Comments?
Subdivision Surfaces Homework 1: Questions/Comments? 1 Questions on Homework? What s an illegal edge collapse? 1 2 3 a b 4 7 To be legal, the ring of vertex neighbors must be unique (have no duplicates)!
More informationOverview. Efficient Simplification of Point-sampled Surfaces. Introduction. Introduction. Neighborhood. Local Surface Analysis
Overview Efficient Simplification of Pointsampled Surfaces Introduction Local surface analysis Simplification methods Error measurement Comparison PointBased Computer Graphics Mark Pauly PointBased Computer
More informationFree-Form Deformation and Other Deformation Techniques
Free-Form Deformation and Other Deformation Techniques Deformation Deformation Basic Definition Deformation: A transformation/mapping of the positions of every particle in the original object to those
More informationPostprocessing of Compressed 3D Graphic Data
Journal of Visual Communication and Image Representation 11, 80 92 (2000) doi:10.1006/jvci.1999.0430, available online at http://www.idealibrary.com on Postprocessing of Compressed 3D Graphic Data Ka Man
More informationJoe Warren, Scott Schaefer Rice University
Joe Warren, Scott Schaefer Rice University Polygons are a ubiquitous modeling primitive in computer graphics. Their popularity is such that special purpose graphics hardware designed to render polygons
More informationImplicit Surfaces & Solid Representations COS 426
Implicit Surfaces & Solid Representations COS 426 3D Object Representations Desirable properties of an object representation Easy to acquire Accurate Concise Intuitive editing Efficient editing Efficient
More informationScientific Visualization Example exam questions with commented answers
Scientific Visualization Example exam questions with commented answers The theoretical part of this course is evaluated by means of a multiple- choice exam. The questions cover the material mentioned during
More informationMulti-Scale Free-Form Surface Description
Multi-Scale Free-Form Surface Description Farzin Mokhtarian, Nasser Khalili and Peter Yuen Centre for Vision Speech and Signal Processing Dept. of Electronic and Electrical Engineering University of Surrey,
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationA Sub-Atomic Subdivision Approach
A Sub-Atomic Subdivision Approach S. Seeger 1 K. Hormann 2 G. Häusler 1 G. Greiner 2 1 University of Erlangen, Chair for Optics Staudtstr. 7/B2, 91058 Erlangen, Germany Email: {sseeger,ghaeusler}@optik.uni-erlangen.de
More informationGeometric Modeling. Bing-Yu Chen National Taiwan University The University of Tokyo
Geometric Modeling Bing-Yu Chen National Taiwan University The University of Tokyo Surface Simplification Motivation Basic Idea of LOD Discrete LOD Continuous LOD Simplification Problem Characteristics
More informationEfficient Representation and Extraction of 2-Manifold Isosurfaces Using kd-trees
Efficient Representation and Extraction of 2-Manifold Isosurfaces Using kd-trees Alexander Greß gress@cs.uni-bonn.de Reinhard Klein rk@cs.uni-bonn.de University of Bonn Institute of Computer Science II
More informationSimplification. Stolen from various places
Simplification Stolen from various places The Problem of Detail Graphics systems are awash in model data: highly detailed CAD models high-precision surface scans surface reconstruction algorithms Available
More informationHexahedral Mesh Generation for Volumetric Image Data
Hexahedral Mesh Generation for Volumetric Image Data Jason Shepherd University of Utah March 27, 2006 Outline Hexahedral Constraints Topology Boundary Quality Zhang et al. papers Smoothing/Quality Existing
More information