Geodesics in heat: A new approach to computing distance
|
|
- Mae Rice
- 5 years ago
- Views:
Transcription
1 Geodesics in heat: A new approach to computing distance based on heat flow Diana Papyan Faculty of Informatics - Technische Universität München Abstract In this report we are going to introduce new method of computing geodesic distance from one point to the specified subset of a given domain, using heat flow. The heat method is simple to implement in comparison with other methods, as it is based on solving standard linear elliptic problems. The other great advantage of this method is that it is possible to apply on any type of domain (point clouds, triangular meshes, polygonal meshes, etc.) since it is using only standard differential operators. At the end of the report we will present a graph which will provide you with information about efficiency of two different methods (Heat method, Fast Marching) of computing geodesic distance. There it is clearly visible that out of those two, heat method is faster then the others and in the same place mean error in most of cases is better then in fast marching. 1 Introduction There exists already several methods to compute geodesics distance, but they are either too slow, or hard to implement or not efficient enough. The heat method can be solution for all of this problems. When applying heat on one point of shape it will spread all over the surface after some time. To describe that heat flow lets take a hot needle and attach it to a point on bunny s surface. The function which will show the amount of heat reaching from source point (hot needle, x) to destination point (y) after time (t) is called heat kernel and is noted as k t,x (y). One of the most famous methods to count geodesic distance using heat kernel is Varadhan s formula 1. Basically the formula uses only pointwise transformation of heat kernel to get distance between any two points. φ(x, y) = lim 4t log k t,x (y) (1) t 0 In this equation (x) represents the source Figure 1: Geodesic distance from heat source to point in Riemannian manifold, (y) - destination the surface point in Riemannian manifold, (t) - the time passed after which we are measuring the amount of heat reaching destination point. The reason why we can t use Varadhan s formula in every case is that it is highly dependent on heat kernel. Any approximation of it will bring to quite big error in distance. You can see that in the Figure 1. Heat method is avoiding this problem by using wider range of inputs. The only requirement for using function as input in heat method is that functions gradient should be parallel to geodesics. Except this, heat method has two major advantages to already existing algorithms. First one is that you can use it on any shape (point cloud, triangular mesh, etc). The other advantage is that it is possible to prefactor parts of solution and use it later, as heat method includes only sparse linear systems. 1
2 Figure 2: Top left is exact reconstruction of heat kernel and under it is counted geodesic distance using that heat kernel. Given an exact reconstruction of the heat kernel (top left) Varadhan s formula can be used to recover geodesic distance (bottom left) but fails in the presence of approximation or numerical error (middle, right), as shown here for a point source in 1D. The robustness of the heat method stems from the fact that it depends only on the direction of the gradient. 2 Existing methods The base approach for solving distance problem is using eikonal equation 2 with boundary condition φ γ = 0 where γ is the subset of the domain. φ(x) = 1 (2) φ is showing the distance from the point where heat was applied to the x point. denotes the gradient of φ function. But this equation is non-linear partially differential hyperbolic equation. Solving it can be really challenging. Two best known methods to solve this equation are Fast Marching and Fast Sweeping. Though Fast Sweeping is easier to implement but it is slower in comparison with Fast Marching. [2] But still both of them are slow and not efficient as they are not reusing any computation, so for each subset of domain you need to compute everything from scratch. Also this method is almost not possible to make parallel as it is using priority queue. So even if you have a big cluster the program still will run on only one processor, which is a huge problem. The heat method is mostly close to Rangarajan and Gurumoorthy method 3, which is using Schrödinger wave equation formalism to solve Eikonal equation. φ = log ψ (3) In equation is stated for Planck s constant. But the problem of this method is that you need to use Fast Fourier Transformation to solve it, which is automatically restricting you to use regular grids only. And also there is no any information on how should be chosen for getting best results. 3 The Heat Method From here on will be stated as negative, semi definite Laplace-Beltrami operator acting on differentiable real valued functions over the Riemannian manifold(m). 2
3 Figure 3: Heat method brief outline: u heat is spreading over domain. u is the gradient of heat function which is point to the the heat source. X is a unit vector field whose direction is opposite to gradients direction. Function φ is final recovered distance. Now lets discuss heat method. Heat Method consists of three main steps which are 1) Solve heat equation u = u 2) Normalize gradient X = u u 3) Solve Poisson equation φ = X φ approximates geodesic distance. It is most precise when time (t) is approaching 0 in heat equation. Also when solving the Poisson equation you re not getting unique answer. All the values should be later shifted so that the smallest number will become 0. So basically in heat method we don t have problem of not having precise value of heat kernel as we did in Varadhans formula, as we only need correct direction of gradient. Also we can see from Eikonal equation that the length of true distance function has always length of 1, so no need to have precise value of gradients magnitude at all. So we need to compute the normalized gradient field and find the closest scalar potential by minimizing M φ X 2. Step descriptions can be seen in Figure Time Discretization For being able to solve first step (heat equation) we need to discretize it in time. For doing that we can use single backward Euler step for some fixed time t. So we need to solve 4 linear equation over the whole domain M. (id t )u t = u 0 (4) Here id is stated for identity operator. Now lets consider elliptic boundary value problem (id t )υ t = 0 on M\γ υ t = 1 on γ υ t is equal to u t up to a multiplicative constant for a point source. υ t also has a close relationship with distance, which is 6 t lim t 0 2 log υ t = φ (6) outside of cut locus. Above mentioned relationship shows that transformation applied to υ t preserves the direction of the gradient, and this means that prerequisites are satisfied for step 2 and 3. (5) 3
4 3.2 Spatial Discretization Heat method can be applied to any type of shape with discrete gradient ( ), divergence ( ) and Laplace operator ( ). Underneath is it shown how it can be done for several types of shapes, but it is not limited to those types, you can apply it to other meshes too. The computational cost of the heat method mostly depends on the intrinsic dimension n of M. The other methods like fast marching requires a grid of the same dimension m as the ambient space. This is a great advantage of heat method when using it in machine learning projects where m can be significantly larger than n Simplicial Meshes Let u R V specify a piecewise linear function on a triangulated surface. Standart discretization of Laplacian at the vertex i is defined by 7 (Lu) i = 1 (cot α i,j + cot β i,j )(u j u i ) (7) 2A i where A i is defined by one third the area of all triangles incident on vertex i, the sum is counted on all neighboring vertices j, and α i,j, β i,j are the angles opposing the corresponding edge. To get feeling of what was said you can see figure It can also be rewritten as matrix operation 8 j L = A 1 L C (8) where A, L C R V V. A is a diagonal matrix containing the vertex areas, L C is named as cotan operator and is the sum part from the equation 7. Heat flow can be computed using equation 9 (A tl C )u = δ γ (9) where δ is a Kronecker delta over γ (Kronecker delta gives the integrated value of a Dirac delta). The gradient in following triangle is expressed as equation 10 u = 1 u i (N e i ) (10) 2A f i where A f is the area of the face, N is its unit normal, e i is the i th edge vector (oriented counterclockwise), and u i is the value of u at the opposing vertex. The divergence associated with i vertex of triangle can be computed with equation X = 1 cot θ 1 (e 1 X j ) + cot θ 2 (e 2 X j ) (11) 2 j where the sum is taken over incident triangles j each with a vector X j. e 1 and e 2 are the two edge vectors of triangle j containing i vertex, and θ 1, θ 2 are the opposing angles. Now lets take b R V, and lets assume that it is a vector of (integrated) divergences of the normalized vector field X. In that case we can compute final distance by solving Poisson equation which is the last step of heat method 12 L C φ = b (12) These solution can be generalized to higher dimensional meshes(tetrahedral meshes) using discrete operators. 4
5 Diana Papyan Figure 4: In the examples it is visible that heat method is working quite good on point clouds. Even if there is some missing points like in left image. And both images don t have any points connecting information Point Clouds Now we need to solve the heat equation for a discrete point sample P Rn of M with no connection information. For a discrete point sample P Rn of M with no connectivity information, there is a recent method for computing Laplacian which is A 1 ν LP C, where Aν is a diagonal matrix of Voronoi areas and LP C is symmetric positive semidefinite. To solve step 2 of heat method (computing vector field) we need to declare new function u : p R which would be height function over approximate tangent planes Tp at each point of point cloud and then count the gradient of a weighted least squares (WLS) approximation of u. Computing tangent planes can be done in different ways, but it is suggested to use moving least squares (MLS) approximation as it is simple to implement. Step 3 of heat method (getting real distance) is solved by positive semi-definite Poisson equation 13. LP C φ = DT Aν X 3.3 (13) Choosing Time Step Accuracy of heat method is highly dependent on time step t which is chosen by person. To go back to 6 equation, we can see that in smooth setting, choosing smaller value of time step would bring to better accuracy in geodesic distance. But the same time lets look into equation 4. The solution for this equation is a function of the combinatorial distance, independent of how we discretize the Laplacian. So if we have a fixed mesh making the t value smaller not always brings to best results. There is an option to over-go that problem by refining the mesh and then decreasing t accordingly. But there is a reason for not doing that, because for getting smoothed approximation of geodesic distance we need to pick large values of t. That is why we will try to find an optimal t which is not too big and not to small. In that case we will get smooth approximation and good accuracy together. Coming up with optimal expression for t is quite difficult because of complexity of analyses involving cut lotus. But practice shows that t = mh2 where h is the mean spacing between adjacent nodes and m is just a positive constant, works quite well in most cases. This estimate is motivated by the fact that h2 is invariant to scaling and refinement. Some experiments were 5
6 made to check which is the smallest possible m which can recover real distance. You can see that in 3.3 and 3.3 figures. Figure 5: The experiment is done with h = 1. The solution approaches the combinatorial distance as t goes to zero So the smallest m which recovers distance with good accuracy is when m = 1. So in all upcoming examples we use 14 expression for time step. t = h 2 (14) This equation is not bringing to good results only if we need to get really smooth approximation of geodesic distance. 3.4 Smoothed Distance Geodesic distance is not smooth at cut locuss points, which are for example points from where there is no unique shortest path to the source. Non-smoothness can bring to numerical difficulty in cases when there will be need to compute derivatives of the distance function, or may simply be undesirable aesthetically. Figure 6: Horizontal axis shows different values of m, vertical axis shows mean error. Best results are reaching when m = 1 There are some other methods to get really smooth distances, which are for example diffusion distance or biharmonic distance, but they are far away from being close to real geodesic distance as indicated by uneven spacing of isolines. You can check 3.4 figures middle part for proof. Figure 7: When the heat is applied from front side of shape it results in non-smooth edges on the opposite side (left). Going to right side those are images with longer t duration. The best smoothed approximations of geodesic distance is reached when t is biggest (right). 6
7 Diana Papyan Figure 8: Top: heat method smoothed geodesic approximation (left) and biharmonic distance (center) both mitigate sharp edges, yet isoline spacing of the biharmonic distance can vary dramatically. Bottom: biharmonic distance tends to exhibit elliptical level lines near the heat source, while heat method smoothed distance maintains isotropic circular profiles as seen in the exact distance (right). But we can reach smooth approximation with heat method too. The trick is to choose big value for t. The results can be seen in 3.4 figure. The computational cost remains the same. As can be seen from figure isolines are evenly spaced for any value of t. That comes from normalization step of heat method (step 2). When using equation t = mh2 meaningful values of m is varying from 1 till 106. When choosing bigger values little visible change can be noticed. 4 Accuracy Tests were done to compare accuracy of fast marching and heat method. In figure 3.4 you can see the visual difference between exact geodesic distance, distance computed with heat method and distance computed by fast marching method. More information about accuracy can 7 Figure 9: Left: exact geodesic distance. Using default parameters, the heat method (middle) and fast marching (right) both produce results of comparable accuracy, here within less than 1% of the exact distance
8 be viewed in 4 figure. References [1] Keenan Crane, Clarisse Weischedel, and Max Wardetzky. Geodesics in Heat: A New Approach to Computing Distance Based on Heat Flow. [2] Pierre A. Gremaud and Christopher M. Kuster. Computational study of fast methods for the eikonal equation. 8
A Primer on Laplacians. Max Wardetzky. Institute for Numerical and Applied Mathematics Georg-August Universität Göttingen, Germany
A Primer on Laplacians Max Wardetzky Institute for Numerical and Applied Mathematics Georg-August Universität Göttingen, Germany Warm-up: Euclidean case Warm-up The Euclidean case Chladni s vibrating plates
More information05 - Surfaces. Acknowledgements: Olga Sorkine-Hornung. CSCI-GA Geometric Modeling - Daniele Panozzo
05 - Surfaces Acknowledgements: Olga Sorkine-Hornung Reminder Curves Turning Number Theorem Continuous world Discrete world k: Curvature is scale dependent is scale-independent Discrete Curvature Integrated
More informationShape Modeling and Geometry Processing
252-0538-00L, Spring 2018 Shape Modeling and Geometry Processing Discrete Differential Geometry Differential Geometry Motivation Formalize geometric properties of shapes Roi Poranne # 2 Differential Geometry
More informationSHORTEST PATHS ON SURFACES GEODESICS IN HEAT
SHORTEST PATHS ON SURFACES GEODESICS IN HEAT INF555 Digital Representation and Analysis of Shapes 28 novembre 2015 Ruoqi He & Chia-Man Hung 1 INTRODUCTION In this project we present the algorithm of a
More informationCS 523: Computer Graphics, Spring Shape Modeling. Differential Geometry of Surfaces
CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry of Surfaces Andrew Nealen, Rutgers, 2011 2/22/2011 Differential Geometry of Surfaces Continuous and Discrete Motivation Smoothness
More informationSpectral Surface Reconstruction from Noisy Point Clouds
Spectral Surface Reconstruction from Noisy Point Clouds 1. Briefly summarize the paper s contributions. Does it address a new problem? Does it present a new approach? Does it show new types of results?
More informationPreparation Meeting. Recent Advances in the Analysis of 3D Shapes. Emanuele Rodolà Matthias Vestner Thomas Windheuser Daniel Cremers
Preparation Meeting Recent Advances in the Analysis of 3D Shapes Emanuele Rodolà Matthias Vestner Thomas Windheuser Daniel Cremers What You Will Learn in the Seminar Get an overview on state of the art
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 informationGeometric Modeling Assignment 3: Discrete Differential Quantities
Geometric Modeling Assignment : Discrete Differential Quantities Acknowledgements: Julian Panetta, Olga Diamanti Assignment (Optional) Topic: Discrete Differential Quantities with libigl Vertex Normals,
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 informationSOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS
SOLVING PARTIAL DIFFERENTIAL EQUATIONS ON POINT CLOUDS JIAN LIANG AND HONGKAI ZHAO Abstract. In this paper we present a general framework for solving partial differential equations on manifolds represented
More informationCS 468 (Spring 2013) Discrete Differential Geometry
CS 468 (Spring 2013) Discrete Differential Geometry 1 Math Review Lecture 14 15 May 2013 Discrete Exterior Calculus Lecturer: Justin Solomon Scribe: Cassidy Saenz Before we dive into Discrete Exterior
More informationarxiv: v1 [cs.gr] 23 May 2018
The Vector Heat Method NICHOLAS SHARP, YOUSUF SOLIMAN, and KEENAN CRANE, Carnegie Mellon University arxiv:1805.09170v1 [cs.gr] 23 May 2018 This paper describes a method for efficiently computing parallel
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 informationLaplacian Meshes. COS 526 Fall 2016 Slides from Olga Sorkine and Yaron Lipman
Laplacian Meshes COS 526 Fall 2016 Slides from Olga Sorkine and Yaron Lipman Outline Differential surface representation Ideas and applications Compact shape representation Mesh editing and manipulation
More informationDef De orma f tion orma Disney/Pixar
Deformation Disney/Pixar Deformation 2 Motivation Easy modeling generate new shapes by deforming existing ones 3 Motivation Easy modeling generate new shapes by deforming existing ones 4 Motivation Character
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More informationLectures in Discrete Differential Geometry 3 Discrete Surfaces
Lectures in Discrete Differential Geometry 3 Discrete Surfaces Etienne Vouga March 19, 2014 1 Triangle Meshes We will now study discrete surfaces and build up a parallel theory of curvature that mimics
More informationScanning Real World Objects without Worries 3D Reconstruction
Scanning Real World Objects without Worries 3D Reconstruction 1. Overview Feng Li 308262 Kuan Tian 308263 This document is written for the 3D reconstruction part in the course Scanning real world objects
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 informationEdge and local feature detection - 2. Importance of edge detection in computer vision
Edge and local feature detection Gradient based edge detection Edge detection by function fitting Second derivative edge detectors Edge linking and the construction of the chain graph Edge and local feature
More informationGeometric Modeling and Processing
Geometric Modeling and Processing Tutorial of 3DIM&PVT 2011 (Hangzhou, China) May 16, 2011 6. Mesh Simplification Problems High resolution meshes becoming increasingly available 3D active scanners Computer
More informationParallel Computation of Spherical Parameterizations for Mesh Analysis. Th. Athanasiadis and I. Fudos University of Ioannina, Greece
Parallel Computation of Spherical Parameterizations for Mesh Analysis Th. Athanasiadis and I. Fudos, Greece Introduction Mesh parameterization is a powerful geometry processing tool Applications Remeshing
More informationAssignment 4: Mesh Parametrization
CSCI-GA.3033-018 - Geometric Modeling Assignment 4: Mesh Parametrization In this exercise you will Familiarize yourself with vector field design on surfaces. Create scalar fields whose gradients align
More information(Discrete) Differential Geometry
(Discrete) Differential Geometry Motivation Understand the structure of the surface Properties: smoothness, curviness, important directions How to modify the surface to change these properties What properties
More informationLarge Mesh Deformation Using the Volumetric Graph Laplacian
Large Mesh Deformation Using the Volumetric Graph Laplacian Kun Zhou1 Jin Huang2 John Snyder3 Xinguo Liu1 Hujun Bao2 Baining Guo1 Heung-Yeung Shum1 1 Microsoft Research Asia 2 Zhejiang University 3 Microsoft
More information1.7.1 Laplacian Smoothing
1.7.1 Laplacian Smoothing 320491: Advanced Graphics - Chapter 1 434 Theory Minimize energy functional total curvature estimate by polynomial-fitting non-linear (very slow!) 320491: Advanced Graphics -
More informationComputing and Processing Correspondences with Functional Maps
Computing and Processing Correspondences with Functional Maps SIGGRAPH 2017 course Maks Ovsjanikov, Etienne Corman, Michael Bronstein, Emanuele Rodolà, Mirela Ben-Chen, Leonidas Guibas, Frederic Chazal,
More informationGreedy Routing with Guaranteed Delivery Using Ricci Flow
Greedy Routing with Guaranteed Delivery Using Ricci Flow Jie Gao Stony Brook University Joint work with Rik Sarkar, Xiaotian Yin, Wei Zeng, Feng Luo, Xianfeng David Gu Greedy Routing Assign coordinatesto
More informationDigital Geometry Processing
Digital Geometry Processing Spring 2011 physical model acquired point cloud reconstructed model 2 Digital Michelangelo Project Range Scanning Systems Passive: Stereo Matching Find and match features in
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 informationParameterization of Meshes
2-Manifold Parameterization of Meshes What makes for a smooth manifold? locally looks like Euclidian space collection of charts mutually compatible on their overlaps form an atlas Parameterizations are
More informationFast marching methods
1 Fast marching methods Lecture 3 Alexander & Michael Bronstein tosca.cs.technion.ac.il/book Numerical geometry of non-rigid shapes Stanford University, Winter 2009 Metric discretization 2 Approach I:
More informationSurfaces, meshes, and topology
Surfaces from Point Samples Surfaces, meshes, and topology A surface is a 2-manifold embedded in 3- dimensional Euclidean space Such surfaces are often approximated by triangle meshes 2 1 Triangle mesh
More informationDigital Geometry Processing Parameterization I
Problem Definition Given a surface (mesh) S in R 3 and a domain find a bective F: S Typical Domains Cutting to a Disk disk = genus zero + boundary sphere = closed genus zero Creates artificial boundary
More informationCorrespondence. CS 468 Geometry Processing Algorithms. Maks Ovsjanikov
Shape Matching & Correspondence CS 468 Geometry Processing Algorithms Maks Ovsjanikov Wednesday, October 27 th 2010 Overall Goal Given two shapes, find correspondences between them. Overall Goal Given
More informationJustin Solomon MIT, Spring Numerical Geometry of Nonrigid Shapes
Justin Solomon MIT, Spring 2017 Numerical Geometry of Nonrigid Shapes Intrinsically far Extrinsically close Geodesic distance [jee-uh-des-ik dis-tuh-ns]: Length of the shortest path, constrained not to
More information: Mesh Processing. Chapter 8
600.657: Mesh Processing Chapter 8 Handling Mesh Degeneracies [Botsch et al., Polygon Mesh Processing] Class of Approaches Geometric: Preserve the mesh where it s good. Volumetric: Can guarantee no self-intersection.
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationDigital Image Processing. Prof. P. K. Biswas. Department of Electronic & Electrical Communication Engineering
Digital Image Processing Prof. P. K. Biswas Department of Electronic & Electrical Communication Engineering Indian Institute of Technology, Kharagpur Lecture - 21 Image Enhancement Frequency Domain Processing
More informationSurface Parameterization
Surface Parameterization A Tutorial and Survey Michael Floater and Kai Hormann Presented by Afra Zomorodian CS 468 10/19/5 1 Problem 1-1 mapping from domain to surface Original application: Texture mapping
More informationDiffusion Wavelets for Natural Image Analysis
Diffusion Wavelets for Natural Image Analysis Tyrus Berry December 16, 2011 Contents 1 Project Description 2 2 Introduction to Diffusion Wavelets 2 2.1 Diffusion Multiresolution............................
More informationCS 177 Homework 1. Julian Panetta. October 22, We want to show for any polygonal disk consisting of vertex set V, edge set E, and face set F:
CS 177 Homework 1 Julian Panetta October, 009 1 Euler Characteristic 1.1 Polyhedral Formula We want to show for any polygonal disk consisting of vertex set V, edge set E, and face set F: V E + F = 1 First,
More informationImage Sampling and Quantisation
Image Sampling and Quantisation Introduction to Signal and Image Processing Prof. Dr. Philippe Cattin MIAC, University of Basel 1 of 46 22.02.2016 09:17 Contents Contents 1 Motivation 2 Sampling Introduction
More informationImage Sampling & Quantisation
Image Sampling & Quantisation Biomedical Image Analysis Prof. Dr. Philippe Cattin MIAC, University of Basel Contents 1 Motivation 2 Sampling Introduction and Motivation Sampling Example Quantisation Example
More informationDifferential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder]
Differential Geometry: Circle Patterns (Part 1) [Discrete Conformal Mappinngs via Circle Patterns. Kharevych, Springborn and Schröder] Preliminaries Recall: Given a smooth function f:r R, the function
More informationEdge detection. Convert a 2D image into a set of curves. Extracts salient features of the scene More compact than pixels
Edge Detection Edge detection Convert a 2D image into a set of curves Extracts salient features of the scene More compact than pixels Origin of Edges surface normal discontinuity depth discontinuity surface
More informationNumerical Geometry of Nonrigid Shapes. CS 468, Spring 2013 Differential Geometry for Computer Science Justin Solomon and Adrian Butscher
Numerical Geometry of Nonrigid Shapes CS 468, Spring 2013 Differential Geometry for Computer Science Justin Solomon and Adrian Butscher Intrinsically far Extrinsically close Straightest Geodesics on Polyhedral
More information04 - Normal Estimation, Curves
04 - Normal Estimation, Curves Acknowledgements: Olga Sorkine-Hornung Normal Estimation Implicit Surface Reconstruction Implicit function from point clouds Need consistently oriented normals < 0 0 > 0
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 informationHomework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in class hard-copy please)
Virginia Tech. Computer Science CS 5614 (Big) Data Management Systems Fall 2014, Prakash Homework 4: Clustering, Recommenders, Dim. Reduction, ML and Graph Mining (due November 19 th, 2014, 2:30pm, in
More informationRobotic Motion Planning: Review C-Space and Start Potential Functions
Robotic Motion Planning: Review C-Space and Start Potential Functions Robotics Institute 16-735 http://www.cs.cmu.edu/~motionplanning Howie Choset http://www.cs.cmu.edu/~choset What if the robot is not
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationDelaunay Triangulations
Delaunay Triangulations (slides mostly by Glenn Eguchi) Motivation: Terrains Set of data points A R 2 Height ƒ(p) defined at each point p in A How can we most naturally approximate height of points not
More informationTutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass
Tutorial 3 Comparing Biological Shapes Patrice Koehl and Joel Hass University of California, Davis, USA http://www.cs.ucdavis.edu/~koehl/ims2017/ What is a shape? A shape is a 2-manifold with a Riemannian
More informationReal-Time Shape Editing using Radial Basis Functions
Real-Time Shape Editing using Radial Basis Functions, Leif Kobbelt RWTH Aachen Boundary Constraint Modeling Prescribe irregular constraints Vertex positions Constrained energy minimization Optimal fairness
More informationA TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3
A TESSELLATION FOR ALGEBRAIC SURFACES IN CP 3 ANDREW J. HANSON AND JI-PING SHA In this paper we present a systematic and explicit algorithm for tessellating the algebraic surfaces (real 4-manifolds) F
More informationINTRODUCTION TO THE HOMOLOGY GROUPS OF COMPLEXES
INTRODUCTION TO THE HOMOLOGY GROUPS OF COMPLEXES RACHEL CARANDANG Abstract. This paper provides an overview of the homology groups of a 2- dimensional complex. It then demonstrates a proof of the Invariance
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 informationSurfaces: notes on Geometry & Topology
Surfaces: notes on Geometry & Topology 1 Surfaces A 2-dimensional region of 3D space A portion of space having length and breadth but no thickness 2 Defining Surfaces Analytically... Parametric surfaces
More informationAssignment 5: Shape Deformation
CSCI-GA.3033-018 - Geometric Modeling Assignment 5: Shape Deformation Goal of this exercise In this exercise, you will implement an algorithm to interactively deform 3D models. You will construct a two-level
More informationManipulating the Boundary Mesh
Chapter 7. Manipulating the Boundary Mesh The first step in producing an unstructured grid is to define the shape of the domain boundaries. Using a preprocessor (GAMBIT or a third-party CAD package) you
More informationDr. Allen Back. Nov. 21, 2014
Dr. Allen Back of Nov. 21, 2014 The most important thing you should know (e.g. for exams and homework) is how to setup (and perhaps compute if not too hard) surface integrals, triple integrals, etc. But
More information274 Curves on Surfaces, Lecture 5
274 Curves on Surfaces, Lecture 5 Dylan Thurston Notes by Qiaochu Yuan Fall 2012 5 Ideal polygons Previously we discussed three models of the hyperbolic plane: the Poincaré disk, the upper half-plane,
More informationMesh data structures. Abstract. 1 Introduction. 2 Method. 2.1 Implement the half edge mesh LIU TNM
Mesh data structures Fredrik Salomonsson fresa517@student.liu.se LIU TNM079 2010-05-28 Abstract When dealing with meshes containing large amount of triangles or quadrilaterals one need to have some kind
More informationChapter 13. Boundary Value Problems for Partial Differential Equations* Linz 2002/ page
Chapter 13 Boundary Value Problems for Partial Differential Equations* E lliptic equations constitute the third category of partial differential equations. As a prototype, we take the Poisson equation
More informationSurface Reconstruction with MLS
Surface Reconstruction with MLS Tobias Martin CS7960, Spring 2006, Feb 23 Literature An Adaptive MLS Surface for Reconstruction with Guarantees, T. K. Dey and J. Sun A Sampling Theorem for MLS Surfaces,
More informationGeometric Computations for Simulation
1 Geometric Computations for Simulation David E. Johnson I. INTRODUCTION A static virtual world would be boring and unlikely to draw in a user enough to create a sense of immersion. Simulation allows things
More informationParameterization. Michael S. Floater. November 10, 2011
Parameterization Michael S. Floater November 10, 2011 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to generate from point
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 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 informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationTutorial 9 - Fast Marching Methods
236861 Numerical Geometry of Images Tutorial 9 Fast Marching Methods c 2012 Why measure distances? Shape analysis Image analysis Registration Morphology Navigation And many other uses.. There are many
More informationParallel Lines Investigation
Year 9 - The Maths Knowledge Autumn 1 (x, y) Along the corridor, up the stairs (3,1) x = 3 Gradient (-5,-2) (0,0) y-intercept Vertical lines are always x = y = 6 Horizontal lines are always y = Parallel
More informationHow Much Geometry Lies in The Laplacian?
How Much Geometry Lies in The Laplacian? Encoding and recovering the discrete metric on triangle meshes Distance Geometry Workshop in Bad Honnef, November 23, 2017 Maks Ovsjanikov Joint with: E. Corman,
More informationIntersection of an Oriented Box and a Cone
Intersection of an Oriented Box and a Cone David Eberly, Geometric Tools, Redmond WA 98052 https://www.geometrictools.com/ This work is licensed under the Creative Commons Attribution 4.0 International
More informationMultimedia Computing: Algorithms, Systems, and Applications: Edge Detection
Multimedia Computing: Algorithms, Systems, and Applications: Edge Detection By Dr. Yu Cao Department of Computer Science The University of Massachusetts Lowell Lowell, MA 01854, USA Part of the slides
More informationGraphs. The ultimate data structure. graphs 1
Graphs The ultimate data structure graphs 1 Definition of graph Non-linear data structure consisting of nodes & links between them (like trees in this sense) Unlike trees, graph nodes may be completely
More informationAdaptive and Smooth Surface Construction by Triangular A-Patches
Adaptive and Smooth Surface Construction by Triangular A-Patches Guoliang Xu Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, China Abstract
More informationSpace Filling Curves and Hierarchical Basis. Klaus Speer
Space Filling Curves and Hierarchical Basis Klaus Speer Abstract Real world phenomena can be best described using differential equations. After linearisation we have to deal with huge linear systems of
More informationMultiresolution Computation of Conformal Structures of Surfaces
Multiresolution Computation of Conformal Structures of Surfaces Xianfeng Gu Yalin Wang Shing-Tung Yau Division of Engineering and Applied Science, Harvard University, Cambridge, MA 0138 Mathematics Department,
More informationAs a consequence of the operation, there are new incidences between edges and triangles that did not exist in K; see Figure II.9.
II.4 Surface Simplification 37 II.4 Surface Simplification In applications it is often necessary to simplify the data or its representation. One reason is measurement noise, which we would like to eliminate,
More informationAnalysis of high dimensional data via Topology. Louis Xiang. Oak Ridge National Laboratory. Oak Ridge, Tennessee
Analysis of high dimensional data via Topology Louis Xiang Oak Ridge National Laboratory Oak Ridge, Tennessee Contents Abstract iii 1 Overview 1 2 Data Set 1 3 Simplicial Complex 5 4 Computation of homology
More informationDISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU /858B Fall 2017
DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU 15-458/858B Fall 2017 LECTURE 10: DISCRETE CURVATURE DISCRETE DIFFERENTIAL GEOMETRY: AN APPLIED INTRODUCTION Keenan Crane CMU 15-458/858B
More informationEstimating normal vectors and curvatures by centroid weights
Computer Aided Geometric Design 21 (2004) 447 458 www.elsevier.com/locate/cagd Estimating normal vectors and curvatures by centroid weights Sheng-Gwo Chen, Jyh-Yang Wu Department of Mathematics, National
More informationLecture 15: Segmentation (Edge Based, Hough Transform)
Lecture 15: Segmentation (Edge Based, Hough Transform) c Bryan S. Morse, Brigham Young University, 1998 000 Last modified on February 3, 000 at :00 PM Contents 15.1 Introduction..............................................
More informationTopic 6 Representation and Description
Topic 6 Representation and Description Background Segmentation divides the image into regions Each region should be represented and described in a form suitable for further processing/decision-making Representation
More informationAn introduction to mesh generation Part IV : elliptic meshing
Elliptic An introduction to mesh generation Part IV : elliptic meshing Department of Civil Engineering, Université catholique de Louvain, Belgium Elliptic Curvilinear Meshes Basic concept A curvilinear
More informationAccurate and Dense Wide-Baseline Stereo Matching Using SW-POC
Accurate and Dense Wide-Baseline Stereo Matching Using SW-POC Shuji Sakai, Koichi Ito, Takafumi Aoki Graduate School of Information Sciences, Tohoku University, Sendai, 980 8579, Japan Email: sakai@aoki.ecei.tohoku.ac.jp
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 informationTHIS paper presents the recent advances in mesh deformation
1 On Linear Variational Surface Deformation Methods Mario Botsch Computer Graphics Laboratory ETH Zurich Olga Sorkine Computer Graphics Group TU Berlin Abstract This survey reviews the recent advances
More informationHow shapes are represented in 3D Graphics. Aims and objectives By the end of the lecture you will be able to describe
Today s lecture Today we will learn about The mathematics of 3D space vectors How shapes are represented in 3D Graphics Modelling shapes as polygons Aims and objectives By the end of the lecture you will
More information1 Appendix to notes 2, on Hyperbolic geometry:
1230, notes 3 1 Appendix to notes 2, on Hyperbolic geometry: The axioms of hyperbolic geometry are axioms 1-4 of Euclid, plus an alternative to axiom 5: Axiom 5-h: Given a line l and a point p not on l,
More informationDiFi: Distance Fields - Fast Computation Using Graphics Hardware
DiFi: Distance Fields - Fast Computation Using Graphics Hardware Avneesh Sud Dinesh Manocha UNC-Chapel Hill http://gamma.cs.unc.edu/difi Distance Fields Distance Function For a site a scalar function f:r
More informationMotion Planning. Howie CHoset
Motion Planning Howie CHoset Questions Where are we? Where do we go? Which is more important? Encoders Encoders Incremental Photodetector Encoder disk LED Photoemitter Encoders - Incremental Encoders -
More informationHi everyone. I hope everyone had a good Fourth of July. Today we're going to be covering graph search. Now, whenever we bring up graph algorithms, we
Hi everyone. I hope everyone had a good Fourth of July. Today we're going to be covering graph search. Now, whenever we bring up graph algorithms, we have to talk about the way in which we represent the
More informationData Analysis 3. Support Vector Machines. Jan Platoš October 30, 2017
Data Analysis 3 Support Vector Machines Jan Platoš October 30, 2017 Department of Computer Science Faculty of Electrical Engineering and Computer Science VŠB - Technical University of Ostrava Table of
More informationCS 532: 3D Computer Vision 14 th Set of Notes
1 CS 532: 3D Computer Vision 14 th Set of Notes Instructor: Philippos Mordohai Webpage: www.cs.stevens.edu/~mordohai E-mail: Philippos.Mordohai@stevens.edu Office: Lieb 215 Lecture Outline Triangulating
More informationSIMPLICIAL ENERGY AND SIMPLICIAL HARMONIC MAPS
SIMPLICIAL ENERGY AND SIMPLICIAL HARMONIC MAPS JOEL HASS AND PETER SCOTT Abstract. We introduce a combinatorial energy for maps of triangulated surfaces with simplicial metrics and analyze the existence
More information