Manifold T-spline. Ying He 1 Kexiang Wang 2 Hongyu Wang 2 Xianfeng David Gu 2 Hong Qin 2. Geometric Modeling and Processing 2006

Similar documents
Manifold splines. Xianfeng Gu, Ying He *, Hong Qin

Manifold Splines. Abstract. 1 Introduction and Motivation. Ying He. Hong Qin. Xianfeng Gu Stony Brook University

CONFORMAL SPHERICAL PARAMETRIZATION FOR HIGH GENUS SURFACES

Polycube Splines. Abstract. 1 Introduction and Motivation. Hong Qin Stony Brook. Ying He NTU. Xianfeng Gu Stony Brook. Hongyu Wang Stony Brook

Conformal Spherical Parametrization for High Genus Surfaces

Computational Conformal Geometry and Its Applications

Polycube splines. Hongyu Wang a, Ying He b,, Xin Li a, Xianfeng Gu a, Hong Qin a. 1. Introduction and motivation

Geometric structures on 2-orbifolds

Manifold Splines with Single Extraordinary Point

Topic: Orientation, Surfaces, and Euler characteristic

Polycube Splines. Hongyu Wang a, Ying He b,, Xin Li a, Xianfeng Gu a,hong Qin a, 1. Introduction and Motivation

Meshless Modeling, Animating, and Simulating Point-Based Geometry

Möbius Transformations in Scientific Computing. David Eppstein

Manifold Splines with Single Extraordinary Point

Geometric structures on manifolds

COMPUTING SURFACE UNIFORMIZATION USING DISCRETE BELTRAMI FLOW

05 - Surfaces. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Daniele Panozzo

CS Object Representation. Aditi Majumder, CS 112 Slide 1

Parameterization of triangular meshes

Introduction to geometry

Topological Data Analysis - I. Afra Zomorodian Department of Computer Science Dartmouth College

Human Brain Mapping with Conformal Geometry and Multivariate Tensor-based Morphometry

Geometric structures on manifolds

DISCRETE DIFFERENTIAL GEOMETRY

3D Modeling Parametric Curves & Surfaces

274 Curves on Surfaces, Lecture 5

Surfaces Beyond Classification

Geometric Modeling Mortenson Chapter 11. Complex Model Construction

Discrete Surface Ricci Flow

COMPUTING CONFORMAL INVARIANTS: PERIOD MATRICES

Computing Geodesic Spectra of Surfaces

Greedy Routing with Guaranteed Delivery Using Ricci Flow

Flat Surfaces, Teichmueller Discs, Veech Groups, and the Veech Tessellation

3D Modeling Parametric Curves & Surfaces. Shandong University Spring 2013

Three Points Make a Triangle Or a Circle

Constrained Willmore Tori in the 4 Sphere

Lectures in Discrete Differential Geometry 3 Discrete Surfaces

Generalized PolyCube Trivariate Splines

Surfaces for CAGD. FSP Tutorial. FSP-Seminar, Graz, November

Digital Geometry Processing Parameterization I

Freeform Curves on Spheres of Arbitrary Dimension

Discrete Ricci Flow and Its Applications

Globally Optimal Surface Mapping for Surfaces with Arbitrary Topology

Shape Modeling and Geometry Processing

E. Vinberg, V. Kac, Quasi-homogeneous cones. (Russian) Mat. Zametki 1 (1967)

Dynamic Points: When Geometry Meets Physics. Xiaohu Guo Stony Brook

2D Spline Curves. CS 4620 Lecture 18

From isothermic triangulated surfaces to discrete holomorphicity

Automated Surface Matching using Mutual Information Applied to Riemann Surface Structures

Geometric Methods in Engineering Applications

Introduction to Computational Manifolds and Applications

Globally Optimal Surface Mapping for Surfaces with Arbitrary Topology

Rational Bezier Curves

CS354 Computer Graphics Rotations and Quaternions

Symmetric Conformal Mapping for Surface Matching and Registration

Greedy Routing in Wireless Networks. Jie Gao Stony Brook University

Parameterization. Michael S. Floater. November 10, 2011

Tutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass

Curves and Surfaces for Computer-Aided Geometric Design

THE UNIFORMIZATION THEOREM AND UNIVERSAL COVERS

Multiresolution Computation of Conformal Structures of Surfaces

Introduction p. 1 What Is Geometric Modeling? p. 1 Computer-aided geometric design Solid modeling Algebraic geometry Computational geometry

Brain Surface Conformal Spherical Mapping

HOMEOMORPHISMS OF 3-MANIFOLDS AND THE REALIZATION OF NIELSEN NUMBER

Two Connections between Combinatorial and Differential Geometry

2D Spline Curves. CS 4620 Lecture 13

Conformal Surface Parameterization Using Euclidean Ricci Flow

Shape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow

Curves-on-Surface: A General Shape Comparison Framework

Embedded graphs. Sasha Patotski. Cornell University November 24, 2014

Dynamics on some Z 2 -covers of half-translation surfaces

Implicit Matrix Representations of Rational Bézier Curves and Surfaces

Cannon s conjecture, subdivision rules, and expansion complexes

The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a

Pythagorean - Hodograph Curves: Algebra and Geometry Inseparable

Lecture notes for Topology MMA100

Bezier Curves, B-Splines, NURBS

Don t just read it; fight it! Ask your own questions, look for your own examples, discover your own proofs. Is the hypothesis necessary?

Brain Surface Conformal Parameterization with Algebraic Functions

Parameterization of Star-Shaped Volumes Using Green s Functions

A Construction of Rational Manifold Surfaces of Arbitrary Topology and Smoothness from Triangular Meshes

Mathematical Research Letters 1, (1994) MÖBIUS CONE STRUCTURES ON 3-DIMENSIONAL MANIFOLDS. Feng Luo

Fathi El-Yafi Project and Software Development Manager Engineering Simulation

Barycentric Coordinates and Parameterization

Introduction to Riemann surfaces, moduli spaces and its applications

Surface Parameterization

The Convex Hull of a Space Curve. Bernd Sturmfels, UC Berkeley (two papers with Kristian Ranestad)

Optimal Möbius Transformation for Information Visualization and Meshing

Tripod Configurations

Geometric Queries for Ray Tracing

Translation surfaces: saddle connections, triangles and covering constructions

GEOMETRIC TOOLS FOR COMPUTER GRAPHICS

Curve and Surface Basics

Algebra II: Review Exercises

Lecture IV Bézier Curves

GEOMETRIC LIBRARY. Maharavo Randrianarivony

Spline Thin-Shell Simulation of Manifold Surfaces

Lecture 0: Reivew of some basic material

CS-184: Computer Graphics. Today

Parameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia

Transcription:

Ying He 1 Kexiang Wang 2 Hongyu Wang 2 Xianfeng David Gu 2 Hong Qin 2 1 School of Computer Engineering Nanyang Technological University, Singapore 2 Center for Visual Computing (CVC) Stony Brook University, USA Geometric Modeling and Processing 2006

Geometry Erlangen Program - F. Klein 1872 Different geometries study the invariants under different transformation groups. Euclidean Geometry : Rigid motion on R 2. Distances between arbitrary two points are the invariants. Affine Geometry: Affine transformations. Parallelism and barry centric coordinates are the invariants. Real Projective Geometry: Real projective transformations. Collinearity and cross ratios are the invariants.

Geometry Erlangen Program - F. Klein 1872 Different geometries study the invariants under different transformation groups. Euclidean Geometry : Rigid motion on R 2. Distances between arbitrary two points are the invariants. Affine Geometry: Affine transformations. Parallelism and barry centric coordinates are the invariants. Real Projective Geometry: Real projective transformations. Collinearity and cross ratios are the invariants.

Algorithms vs. Geometries Central Problem Can different geometries be defined on general surfaces? Can different planar algorithms be generalized to surface domains directly? The answers are yes and yes. The major theoretic tool is the Geometric Structure.

Algorithms vs. Geometries Central Problem Can different geometries be defined on general surfaces? Can different planar algorithms be generalized to surface domains directly? The answers are yes and yes. The major theoretic tool is the Geometric Structure.

Main Ideas Geometry Structure A surface is covered by local coordinate charts. Geometric construction is invariant during the transition from one local coordinate to another.

Manifold U α U β g replacements φ α S φ β φ αβ φ α (U α ) φ β (U β )

Manifold Definition (Manifold) A manifold is a topological space Σ covered by a set of open sets {U α }. A homeomorphism φ α : U α R n maps U α to the Euclidean space R n. (U α, φ α ) is called a chart of Σ, the set of all charts {(U α, φ α )} form the atlas of Σ. Suppose U α U β, then φ αβ = φ β φ α : φ α (U α U β ) φ β (U α U β ) is a transition map. Transition maps satisfy cocycle condition, suppose U α U β U γ, then φ βγ φ αβ = φ αγ.

(X, G) structure Definition ((X,G) Atlas) Suppose X is a topological space, G is the transformation group of X. A manifold Σ with an atlas A = {(U α, φ α )} is an (X, G) atlas if 1 φ α (U α ) X, for all charts (U α, φ α ). 2 Transition maps φ αβ G.

(X,G) structure Definition (Equivalent (X, G) atlases) Two (X, G) atlases A 1 and A 2 of Σ are equivalent, if their union is still an (X, G) atlas of Σ. Definition ((X,G) structure) An (X, G) structure of a manifold Σ is an equivalent class of its (X, G) atlases.

(X,G) structure Definition (Equivalent (X, G) atlases) Two (X, G) atlases A 1 and A 2 of Σ are equivalent, if their union is still an (X, G) atlas of Σ. Definition ((X,G) structure) An (X, G) structure of a manifold Σ is an equivalent class of its (X, G) atlases.

Geometries M.C.Esher s art works: Angels and Devils Regular divisin Sphere with Angels Circle limit IV of the plane and Devils Heaven and Hell

Geometries defined on surfaces

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Spherical Structure X: Unit sphere S 2. G: Rotation group. Surfaces: Genus zero closed surfaces; any open surfaces. Harmonic maps.

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Spherical Structure X: Unit sphere S 2. G: Rotation group. Surfaces: Genus zero closed surfaces; any open surfaces. Harmonic maps.

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Affine Structure X: Real plane R 2. G: Affine transformation group. Surfaces: Genus one closed surface and open surfaces. Holomorphic 1-forms.

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Affine Structure X: Real plane R 2. G: Affine transformation group. Surfaces: Genus one closed surface and open surfaces. Holomorphic 1-forms.

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Hyperbolic Structure X: Hyperbolic plane H 2. G: Möbius transformation group. Surfaces: with negative Euler number. Hyperbolic Ricci flow

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Hyperbolic Structure X: Hyperbolic plane H 2. G: Möbius transformation group. Surfaces: with negative Euler number. Hyperbolic Ricci flow

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Real Projective Structure X: Real projective plane RP 2. G: Real projective transformation group. Surfaces: any surface. Hyperbolic Ricci flow.

Common (X,G) structure lacements g T(P) T(g(P)) lacements g T(P) T(g(P)) Real Projective Structure X: Real projective plane RP 2. G: Real projective transformation group. Surfaces: any surface. Hyperbolic Ricci flow.

Pseudo (X,G) structure lacements g T(P) T(g(P)) Conformal Structure X: Complex plane C. G: Biholomorphic maps. Surfaces: any surface. Holomorphic 1-forms lacements g T(P) T(g(P))

Pseudo (X,G) structure lacements g T(P) T(g(P)) Conformal Structure X: Complex plane C. G: Biholomorphic maps. Surfaces: any surface. Holomorphic 1-forms lacements g T(P) T(g(P))

Conformal Structure cementspsfrag replacementspsfrag replacements g g g T(P) T(g(P)) T(P) T(g(P)) T(P) T(g(P)) cementspsfrag replacementspsfrag replacements g g g T(P) T(g(P)) T(P) T(g(P)) T(P) T(g(P))

(X,G) Invariant Algorithms Definition ((X,G) invariant Algorithm) Suppose X is a topological space, G is the transformation group on X. A geometric operator Ω defined on X is (X, G) invariant, if and only if Examples: Ω g = g Ω, g G. Convex Hull: Projective invariant. Voronoi Diagram: Rigid motion invariant. Polar form : Affine invariant.

(X,G) structure Theorem Suppose a manifold with an (X, G) structure, then any (X, G) invariant algorithms can be generalized on the manifold. Corollary (Manifold Splines - Gu,He,Qin 2005) Spline schemes based on polar forms can be defined on a manifold, if and only if the manifold has an affine structure.

Planar Splines cements g T(P) T(g(P)) PSfrag replacements g T(P) T(g(P)) Parametric Affine Invariant The spline is invariants under the affine transformations of the knots and the parameters.

Manifold SPlines F ents C α, F α φ α (U α ) φ β (U β ) φ α U α U β M Z φ β C β, F β φ α,β = φ β φ 1 α Idea: Geometry Structure A mesh is covered by local coordinate charts. Geometric construction is invariant during the transition from one local coordinate to another. Global Parameterization Find atlas with special transition functions.

Manifold SPlines F ents C α, F α φ α (U α ) φ β (U β ) φ α U α U β M Z φ β C β, F β φ α,β = φ β φ 1 α Idea: Geometry Structure A mesh is covered by local coordinate charts. Geometric construction is invariant during the transition from one local coordinate to another. Global Parameterization Find atlas with special transition functions.

Topological Obstructions Theorem (Benzécri 1959) If a closed surface admits an affine structure, it has zero Euler class. Real projective structure Real projective structure is general, it exists for all surfaces. Real projective structure is simple, all transitions are linear rational functions. Real projective structure is suitable for designing manifold spline schemes.

Topological Obstructions Theorem (Benzécri 1959) If a closed surface admits an affine structure, it has zero Euler class. Real projective structure Real projective structure is general, it exists for all surfaces. Real projective structure is simple, all transitions are linear rational functions. Real projective structure is suitable for designing manifold spline schemes.

Conformal Structure Global Tensor Product Structure

Holomorphic 1-forms Definition (Holomorphic 1-form) Suppose Σ is a Riemann surface, {z α } is a local complex parameter, a holomorphic 1-form ω has a local representation as ω = f (z α )dz α, where f (z α ) is a holomorphic function. PSfrag replacements g Locally, ω is the derivative of a T(P) holomorphic function. Globally, T(g(P)) it is not.

Holomorphic 1-forms Original Surface

Holomorphic 1-forms One basis holomorphic 1- form

Holomorphic 1-forms Another one basis holomorphic 1-form

Holomorphic 1-forms Summation of ω 1 and ω 2

Holomorphic 1-forms Difference between ω 1 and ω 2

Holomorphic 1-forms Holomorphic 1-form induces a conformal parameterization.

Holomorphic 1-forms Holomorphic 1-form induces a conformal parameterization.

Holomorphic 1-forms Theorem (Holomorphic 1-forms) All holomorphic 1-forms form a linear space Ω(Σ) which is isomorphic to the first cohomology group H 1 (Σ, R).

Affine Structure Theorem Holomorphic 1-form induces affine structure By integrating a holomorphic 1-form, local coordinate charts can be established. The charts covers the surface without the singularises, the transition maps are translations.

Manifold Splines SPM2005

Manifold Powell-Sabin Spline

Manifold Powell-Sabin Spline

T-spline T-spline is the superset of tensor-product B-spline and industry-standard NURBS Allow T-junction in parametric domain and control net Natural hierarchical structure Much more flexible than NURBS!

Critical Graph

Critical Graph and Local Charts

Critical Graph

Given the domain manifold M with conformal structure φ : M R 2, the manifold T-spline can be formulated as follows: F(u) = n C i B i (φ(u)), u M, (1) i=1 where B i s are basis functions and C i = (x i, y i, z i, w i ) are control points in P 4 whose weights are w i, and whose Cartesian coordinates are 1 w i (x i, y i, z i ). The cartesian coordinates of points on the surface are given by n i=1 (x i, y i, z i )B i (φ(u)) n i=1 w. (2) ib i (φ(u)) Given a parameter u M, the evaluation can be carried out on arbitrary charts covering u.

Hierarchical Surface Reconstruction Minimize a linear combination of interpolation and fairness functionals, min E = E dist + λe fair, (3) where E dist = m F(u i ) p i 2 i=1 and E fair in (3) is a smoothing term.

Hierarchical Surface Reconstruction P, N v = 200K Conformal structure N 1 c = 105 L 1 = 9.6% N 2 c = 295 L 2 = 5.7%

Hierarchical Surface Reconstruction N 3 c = 950 L 3 = 3.8% N 4 c = 2130 L 4 = 2.4% N 5 c = 5087 L 5 = 1.3% N 6 c = 7706 L 6 = 0.74%

Hierarchical Surface Reconstruction Figure: Close-up view of the reconstructed details

Examples

Examples

Examples

Examples

Examples

Statistics Table: Statistics of test cases. N p, # of points in the polygonal mesh; N c, # of control points; rms, root-mean-square error; L, maximal error. The execution time measures in minutes. Object N p N c rms L Time David 200, 000 7, 706 0.08% 0.74% 39m Bunny 34, 000 1, 304 0.09% 0.81% 18m Iphegenia 150, 000 9, 907 0.06% 0.46% 53m Rocker Arm 50, 000 2, 121 0.04% 0.36% 26m Kitten 40, 000 765 0.05% 0.44% 12m

Future Directions Manifold Splines with single singularity. Manifold Splines which are polynomials everywhere with C k continuity. Planar splines based on projective invariants.

Thanks For more information, please email to gu@cs.sunysb.edu. Thank you!