THE GEOMETRIC HEAT EQUATION AND SURFACE FAIRING
|
|
- Justina Curtis
- 6 years ago
- Views:
Transcription
1 THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING ANREW WILLIS BROWN UNIVERSITY, IVISION OF ENGINEERING, PROVIENCE, RI 02912, USA 1. INTROUCTION This paper concentrates on analysis and discussion of the heat equation as it pertains to smoothing of geometric shapes and its relationship to the problem of surface fairing. The geometric heat equation distorts a given shape in order to obtain scale-space representation of a shape. This scale-space provides a complete description of the original shape in terms of small to large scale structures. These structures may then be interrogated by recognition algorithms in an attempt to classify the shape. In computer science, researchers have been working on methods for eliminating noise from mesh data obtained via 3 measurements. In this case one wishes to eliminate the noise present in the surface measurements in order to obtain an improved approximation of the true surface. This process is called surface fairing. These two seemingly different problems actually have many common goals. This report will focus on a generic implementation of the geometric heat equation and a particular implementation of surface fairing suggested by Gabriel Taubin [8]. Each method will be explained and their inter-relationships will be made clear. These relationships will be supported via experimental results obtained from 3 measurements on a set of human faces. The input data consisted of 3 mesh data sets obtained via a 3 laser range scanner, specifically the ShapeGrabber product of Vitana Corporation [1]. 2. THE GEOMETRIC HEAT EQUATION In 1986, Mokhtarian et. al. described a way to modify the standard heat equation for heat conductivity as a method of smoothing shapes in order to extract a multi-scale shape representation of a given object [6]. In doing so, one modifies the original boundary of the shape,, by moving the contour in the direction of the normal by some positive factor. multiplied by the signed curvature of the contour. (2.1) Here controls the rate of the geometric diffusion, or equivalently, curvature flow. Both of these phrases describe in different words or terminologies what is happening to the boundary. The phrase geometric diffusion is derived from the terminology of heat conduction and indicates that the curvature at each point along the boundary is being dissipated among it s neighbors in the same sense that heat dissipates through a conductive substrate by exchanging energy among its neighbors. The phrase curvature flow is derived from the terminology of differential geometry. In this case, the curvature function is defining the extent of motion or flow which the contour experiences in the direction of the normal vector field. In practice, the contour is represented approximately by a set of linear segments specified by ordered points. This provides a polygonal approximation of the actual 2 shape. From this we compute and via standard numerical methods as applied to equations provided in textbooks on Classical ifferential Geometry such as Struik s [7]. Of note here is the fact that we are using methods of differential geometry to represent the boundary. This alternative representation of the shape makes computation somewhat more difficult. Since each value is differential with respect to the boundary moving frame, they depend upon local regions of the selected parametrization of the curve. This makes the numerical solutions of the PE more difficult to compute since it implies that the grid is non-uniform for each iteration and also may have varying domains for any two sequential timesteps. Hence, constant care must be taken to preserve the spatial structure of the curve. As time passes, the geometric heat equation causes the curvature signal of the boundary, parametrized according to arc length,, to converge to a constant. This implies that the boundary will converge to a circle, and for large enough, this circle eventually reduces to a point 1. In the same way, the standard 1-dimensional PE for heat conduction (2.2) decreases the curvature at each point of the signal with increasing time. (2.2) "!#! $ % " & Hence we observe the same phenomenon with as observed with. However, since the actual values of the heat signal correspond to the curvature signature of the boundary, the net effect of this heat& conduction results in modification of the shape boundary. Hence, by using the heat equation on a differential 1 representation of shape, we have combined classical physics and differential geometry to generate a method for evolving shape contours. 1 In 3 this phenomenon is referred to as shrinkage and is considered to be an undesirable effect. 1
2 !! % $ THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 2 I = 0 I = I = I = I = I = FIGURE 2.1. The Geometric Heat Equation : parabolic curves are superimposed on evolved surfaces for dataset Agnes.mat. Here denotes the number of integration steps necessary to compute the solution. This method also extends to 3 surfaces in Monge form, i.e. of the form %. Since this is a function graph, we may apply the standard methods of elementary differential geometry to obtain representations for curvature. For the face data, we had the benefit of a uniform spatial grid in, which allowed for direct computation of the corresponding form of the geometric heat equation: % Here, the height function at a point % %! % %! % %! % % is modified proportional to the mean curvature. Since the RHS of the equation above differs from the classical definition of mean curvature by only % % [2]. This has the implication that the above equation is equivalent to the exact geometric heat equation for a Monge patch with a non-constant dissipation term! "#. Consequently, application of this PE will smooth the surface according to the value of the mean curvature using a non-uniform timestep scheme. This PE was implemented and executed with the initial condition % being several examples of face surfaces. The results of these experiments over time is provided in Fig. (2.1). As demonstrated in [5], visualization of surface evolution can be enhanced by tracking the evolution of specific lines of curvature on the face. In particular, parabolic lines of maximum curvature and minimum curvature on the face surface are indicated on each of the plots. These correspond to surface regions where one of the two principal directions %'&! % %!#! " %
3 THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 3 I = 6000 I = I = I = I = I = FIGURE 2.2. The Geometric Heat Equation : parabolic curves are superimposed on evolved surfaces for dataset Agnes.mat. Here denotes the number of integration steps necessary to compute the solution. or % is close to zero. Initially, we attempted to compute parabolic lines with respect to the function graph % This may be accomplished by looking for points on the surface which satisfy the relation!#!!. However, in many instances the face was not oriented such that the parabolic lines would lie along areas of importance. Hence, the only meaningful parabolic curves could be extracted via coordinate-free methods using %'&! or %. One can see that due to the necessarily small timestep required by using the forward Euler integration method [4], the geometric heat equation evolves the surface very slowly. For this reason, additional results are provided in Fig. (2.2) showing the extension of the surface evolution up to integration iterations. One could attempt to align the face such that the face is more orthogonal to the plane which should improve the results obtained via the function graph method of computing parabolic lines. Additionally, there are two possible ways of enforcing boundary conditions. The experiments below show results obtained by enforcing vonneumann style boundary conditions, i.e.! for all points on the boundary of the domain. It may be interesting to experiment with irichlet style boundary conditions, i.e. % at all boundary points for some constant. 3. SURFACE FAIRING The problem of surface fairing is somewhat different than that discussed in 2. Here, the object is to eliminate unwanted noise or geometry from a surface which may be a-priori known or unknown. In the case of a known surface, one may wish to amplify specific aspects of the surface while attenuating others for the various purposes. In general, surface fairing is used as a pre-processing step prior to computation or extraction of higher order shape statistics. In these cases, the fairing step allows for more accurate, and consequently more stable, numerical computations based on the surface mesh. Examples of such computations are extraction of surface normals, principal curvatures, and geodesic curves along the surface.
4 ', % 2 = % THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 4 The critical difference between the concepts outlined in 2 and those of surface fairing is that we are attempting to extract a more accurate estimate of the true underlying geometry given a surface mesh. The surface mesh is an approximate description of the true surface via a connected set of 3 points,. Each 3 point or vertex of the mesh corresponds to an actual measurement of the true 3 surface and the connections between the 3 points form planar surface patch approximations of the true surface between vertices. These points may be distributed in a non-uniform manner about the surface and the only requirement on the graph is that it be fully connected (i..e each point must have at least 1 connecting edge to another neighboring point ). The problem of surface fairing has a compelling relationship with the geometric heat equation from 2. For surface fairing methods, the de-facto standard method is called Laplacian smoothing. This method is discussed in Taubin s paper [9] and is considered insufficient due to the shrinkage effect. In order to understand the term Laplacian smoothing we must first define the Laplacian operator and understand some of it s properties and specifically how this operator is related to the solution of the heat equation The Laplacian Operator and the Heat Equation in 1. We start by examining the standard heat equation as discussed in 2. The RHS of this equation may be re-written in terms of the Laplacian operator: "!#!. The corresponding finite difference expression for the Laplacian operator applied to the function at point, time is given by!!!!. This can be computed simply by defining the discrete Laplacian kernel filter,! " $#, and performing convolution: By observing that the filter! "!!! "& is clearly a high-pass filter, it is initially difficult to understand why the Laplacian operator would be associated with any smoothing effect. Yet, by examining the numerical solution of the heat equation we are finally able to establish a connection between the Laplacian operator and smoothing. In order to solve the 1 heat equation, we may integrate the PE using a forward Euler finite difference equation, where all known terms are collected to the RHS of the equation (here ' denotes the timestep). (3.1) )( '* Using the identity matrix and the! " operators, we may re-express (3.1) as (3.2) below. )( (3.2) '+! " & Finally we see that the solution to our partial differential equation now appears numerically the same as filtering. Further, one may simply add the two filters '+! " 10 #-, ' /.. This kernel corresponds to a scaled version of the well known Gaussian kernel 0 for some values of. This is expected since the Gaussian kernel is a Green s function of the heat equation. As ' the diffusion process is accelerated, or equivalently, solutions are computed with larger time steps. Consequently, the solution to this system is obtained by smoothing the signal via convolution with some specific realization of a scaled Gaussian filter. For this reason, the Laplacian operator, which is essentially a second derivative operator has the unlikely and rather inappropriate association with the smoothing effect of Gaussian filtering. This is commonly referred to as Laplacian smoothing, or equivalently, Gaussian smoothing The Shrinkage Effect. The shrinkage effect is analogous to the phenomenon observed in the geometric heat equation where the boundary eventually converges to the shape s center of mass. In the case of explicit integration, at each iteration we applying an averaging filter. which causes each of the grid points move in the direction of their statistical average. It is interesting to note that without external forces applied (i.e. boundary conditions are not forced) 687:9<; this system will eventually converge to the average of all the points as time goes to is given by (3.3). infinity, i.e. for 32, the solution on the domain 4 5 (3.3) 5 For the curvature signature of the boundary, this causes the shape to converge approximately to a circle which then shrinks in size until it finally converges to to the shape s center of mass. This analogy continues to hold true in 3 and the resulting effect of monotonically decreasing perimeter (2) or area (3) is referred to as the shrinkage effect The Laplacian Operator and the Geometric Heat Equation in 2. For the 2 heat equation, this method may be trivially extended via the kernel (3.4) which is often used in image processing applications [10]. (3.4)! " BAC However, when we attempt to apply the Laplacian operator in the case of the 2 geometric heat equation, we encounter a problem which entirely relates to the form of the data. As mentioned in 2 we simply wish to apply the 1 heat equation to the curvature signature of the 6 => &? %
5 0 THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 5 boundary,. This incurrs several numerical nuisances for computation which add to the complexity of computing a stable solution quickly. Specifically, at every time step one must re-compute the normal field and curvature for every point on the boundary. However, if we apply an approximate geometric Gaussian filter to the spatial location of the points, we may compute an approximate solution without the need of additional estimates of curvature and arc length. This is a huge computational advantage for both speed and stability of the computed solution while sacrificing accuracy. Taubin s method of surface fairing suggests one possibility for this filter :. '... The surface may then be smoothed by applying this filter to the boundary (2) or surface (3). For the sake of simplicity, the smoothing algorithm given applies to the smoothing of 2 contours. This algorithm is then modified trivially for direct use on 3 surfaces. # enote the vertices of the boundary curve,, as the ordered set of 2 points. For any given point, the neighborhood of the point is defined as all points which share a common edge with. In this case, since the points are ordered, the neighborhood of ". In 2, the fairing algorithm works as follows: (1) For each point compute the midpoint of the segment connecting the neighborhood as shown in (3.5). (3.5) " (2) Compute the weighting function for the point. This is an arbitrary function which defines the magnitude of the curvature flow at the point. However, the weighting function should have unit gain, i.e., and operate solely upon those points within the neighborhood of. In his paper, Taubin suggests using of. (3) Move in the direction of by a factor of '. (4) Move 0 backwards along the same trajectory by a factor of., where is the number of points in the neighborhood Steps 1-3 above describe the process of standard Laplacian smoothing or, given an appropriate weighting function, the geometric heat equation from 2. Specifically, we have equivalence when the weighting function from Step 2 is an estimate of the curvature such as, where ( ( is the change in the exterior angle of the tangent and is the change in arc length. Step 4 is an attempt to reduce shrinkage. In the problem of surface fairing, one wishes to preserve the volume of the object. Hence, an additional step is necessary after applying the approximate geometric heat equation. This step essentially expands the smoothed surface in the direction of the original surface The 3 Extension : Laplacian Smoothing. The extension of the above algorithm to 3 meshes is trivial. In the case of a surface mesh, instead of computing the midpoint of Step 1, one computes the centroid of all the points in the neighborhood. Step 2 is identical, save for the fact that the neighborhood may now consist of many 3 points, all of which share a common edge with the point of interest. Steps 3 and 4 also generalize directly, where the vectors, and are now 3-dimensional vectors. Since the emergence of this method in 1995, there have been many researchers which have suggested a wide range of possible weighting functions for step 2 of the algorithm. Many of these functions attempt to approximate surface curvature in some way using various combinations of features obtained in the neighborhood of the vertex such as edge lengths, triangle areas or interior / exterior angles of triangles. A list of several of the most important functions is provided in [9]. Unfortunately, this method as defined in Taubin s paper will always either converge to a point, the centroid of the data, or diverge to infinite size depending upon the relative values of ' and. This effect is entirely analogous to the effects of the heat equation summarized by equation (3.3). Later work by esbrun et. al. at Caltech [3] demonstrated an implicit integration method (backward Euler), for solving (3.1) which is unconditionally stable. However, if one wishes to maintain a shape of constant volume the method becomes non-linear and leads to some considerable computational complexity. However, the fact that Taubin s algorithm may be computed quickly and easily for general 3 surfaces of arbitrary topology makes his method extremely powerful. Fig. (3.1) shows results obtained by applying Taubin s surface fairing method to a 3 face dataset. The six figures show how the surface evolves for the weighting function mentioned in Step 2. It is clear that the approximate smoothing evolves the surface in a distinctly different manner. However, the smoothing operation appears to proceed at a much faster rate than observed with the geometric heat equation. It is interesting to note that this averaging process allows for more uniform dissipation due to the definition of the weighting function. For example, in Figs (2.1) and (2.2) the eyes of the surface model are surrounded by very well defined parabolic curves which are almost circular in shape. These shapes remain almost unchanged for the entire duration of the surface evolution. This may be explained by the fact that by definition, the mean curvature at a parabolic point is equal to one of the principal curvatures alone. If this curvature is small, the surface evolution in these regions will proceed extremely slowly. Contrarywise, the surface fairing algorithm does not respect the mean curvature in computing the magnitude of the normal flow. Here, the magnitude of the curvature flow is given roughly by average exterior angle of the connected edges to the given vertex. Hence, one can see that the corresponding parabolic curves quickly evolve as the eyes are smoothed. The lines of curvature seem to combine and split as the surface evolves in both cases destroying some lines of curvature and creating new ones in other instances as saddle points are formed between regions which were previously well separated by non-parabolic points. One example of this can be seen in Fig (3.2) for the surface evolutions and. Fig. (3.2) shows the continuation of Taubin s
6 THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 6 I = I = I = I = I = I = FIGURE 3.1. Smoothing parameters are ' Surface Fairing : parabolic curves are superimposed on the evolved surfaces for the dataset Agnes.mat., denotes the number of integration steps. and surface method after many iterations. For these results, stability was maintained via a surface area normalization step which followed the standard algorithm. 4. CONCLUSIONS Although the objectives of the geometric heat equation and the surface fairing algorithms are somewhat different, it is clear that the underlying theory is very much the same. In essence, the two methods described here are exactly the same given a specific conditions on the surface fairing algorithm of 3. Additionally, results were provided which demonstrated the differences between the two methods. Here, we demonstrated that the approximate method of surface fairing does not as accurately compute curvature flows over the sample mesh. This results in a different surface evolution which can be clearly tracked by the observation of parabolic ridges on the surface. REFERENCES [1] Shapegrabber inc. [2] Wolfram research web site, 2. [3] M. esbrun, M. Meyer, P. Schroder, and A. Barr. Implicit fairing of irregular meshes using diffusion and curvature flow. In SIGGRAPH, [4] B. Gustafsson, H. O. Kriess, and J. Oliger. Time ependent Problems and ifference Methods. John Wiley and Sons Inc., Canada, [5] Peter Hallinan, Gaile Gordon, A.L. Yuille, P. Giblin, and. Mumford. Two-and-Three imensional Patterns of the Face. A.K. Peters Ltd., Natick, MA., USA, [6] F. Mokhtarian and A. K. Mackworth. Scale-based description and recognition of planar curves and two-dimensional shapes. IEEE Trans. on Pattern Anal. Machine Intell., 8(1):34 43, [7] irk J. Struik. Lectures on Classical ifferential Geometry, chapter Elementary Theory of Surfaces. over, NY, 19. [8] Gabriel Taubin. A signal processing approach to fair surface design. In SIGGRAPH, [9] Gabriel Taubin. Geometric signal processing on polygonal meshes. In EUROGRAPHICS, 0. [10] E. Trucco and A. Verri. Introductory Techniques for 3- Computer Vision, chapter Appendix A.2. Prentice Hall, New Jersey, USA, 1998.
7 THE GEOMETRIC HEAT EQUATION AN SURFACE FAIRING 7 I = 400 I = I = I = I = I = FIGURE 3.2. Surface Fairing : parabolic curves are superimposed on the evolved surfaces for the dataset Agnes.mat., denotes the number of integration steps. Smoothing parameters are ' and
Multi-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 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 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 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 informationTHE HALF-EDGE DATA STRUCTURE MODELING AND ANIMATION
THE HALF-EDGE DATA STRUCTURE MODELING AND ANIMATION Dan Englesson danen344@student.liu.se Sunday 12th April, 2011 Abstract In this lab assignment which was done in the course TNM079, Modeling and animation,
More informationCurvature Estimation on Smoothed 3-D Meshes
Curvature Estimation on Smoothed 3-D Meshes Peter Yuen, Nasser Khalili and Farzin Mokhtarian Centre for Vision, Speech and Signal Processing School of Electronic Engineering, Information Technology and
More informationIEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 17, NO. 5, MAY
IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 17, NO. 5, MAY 2008 645 A Real-Time Algorithm for the Approximation of Level-Set-Based Curve Evolution Yonggang Shi, Member, IEEE, and William Clem Karl, Senior
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 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 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 informationMesh Processing Pipeline
Mesh Smoothing 1 Mesh Processing Pipeline... Scan Reconstruct Clean Remesh 2 Mesh Quality Visual inspection of sensitive attributes Specular shading Flat Shading Gouraud Shading Phong Shading 3 Mesh Quality
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 informationA NON-TRIGONOMETRIC, PSEUDO AREA PRESERVING, POLYLINE SMOOTHING ALGORITHM
A NON-TRIGONOMETRIC, PSEUDO AREA PRESERVING, POLYLINE SMOOTHING ALGORITHM Wayne Brown and Leemon Baird Department of Computer Science The United States Air Force Academy 2354 Fairchild Dr., Suite 6G- USAF
More informationSmoothing & Fairing. Mario Botsch
Smoothing & Fairing Mario Botsch Motivation Filter out high frequency noise Desbrun, Meyer, Schroeder, Barr: Implicit Fairing of Irregular Meshes using Diffusion and Curvature Flow, SIGGRAPH 99 2 Motivation
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 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 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 informationA New Approach to Computation of Curvature Scale Space Image for Shape Similarity Retrieval
A New Approach to Computation of Curvature Scale Space Image for Shape Similarity Retrieval Farzin Mokhtarian, Sadegh Abbasi and Josef Kittler Centre for Vision Speech and Signal Processing Department
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 informationPTE 519 Lecture Note Finite Difference Approximation (Model)
PTE 519 Lecture Note 3 3.0 Finite Difference Approximation (Model) In this section of the lecture material, the focus is to define the terminology and to summarize the basic facts. The basic idea of any
More informationExtracting Axially Symmetric Geometry From Limited 3D Range Data
Extracting Axially Symmetric Geometry From Limited 3D Range Data Andrew Willis, Xavier Orriols, Senem Velipasalar, Xavier Binefa, David B. Cooper Brown University, Providence, RI 02912 Computer Vision
More informationProcessing 3D Surface Data
Processing 3D Surface Data Computer Animation and Visualisation Lecture 17 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing
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 informationOutline. Visualization Discretization Sampling Quantization Representation Continuous Discrete. Noise
Fundamentals Data Outline Visualization Discretization Sampling Quantization Representation Continuous Discrete Noise 2 Data Data : Function dependent on one or more variables. Example Audio (1D) - depends
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 informationIn this lecture we introduce the Gauss-Bonnet theorem. The required section is The optional sections are
Math 348 Fall 2017 Lectures 20: The Gauss-Bonnet Theorem II Disclaimer. As we have a textbook, this lecture note is for guidance and supplement only. It should not be relied on when preparing for exams.
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 informationLevel-set MCMC Curve Sampling and Geometric Conditional Simulation
Level-set MCMC Curve Sampling and Geometric Conditional Simulation Ayres Fan John W. Fisher III Alan S. Willsky February 16, 2007 Outline 1. Overview 2. Curve evolution 3. Markov chain Monte Carlo 4. Curve
More informationA Singular Example for the Averaged Mean Curvature Flow
To appear in Experimental Mathematics Preprint Vol. No. () pp. 3 7 February 9, A Singular Example for the Averaged Mean Curvature Flow Uwe F. Mayer Abstract An embedded curve is presented which under numerical
More informationProcessing 3D Surface Data
Processing 3D Surface Data Computer Animation and Visualisation Lecture 15 Institute for Perception, Action & Behaviour School of Informatics 3D Surfaces 1 3D surface data... where from? Iso-surfacing
More informationVariational Methods II
Mathematical Foundations of Computer Graphics and Vision Variational Methods II Luca Ballan Institute of Visual Computing Last Lecture If we have a topological vector space with an inner product and functionals
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 informationThe goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a
The goal is the definition of points with numbers and primitives with equations or functions. The definition of points with numbers requires a coordinate system and then the measuring of the point with
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationGAUSS-BONNET FOR DISCRETE SURFACES
GAUSS-BONNET FOR DISCRETE SURFACES SOHINI UPADHYAY Abstract. Gauss-Bonnet is a deep result in differential geometry that illustrates a fundamental relationship between the curvature of a surface and its
More informationMöbius Transformations in Scientific Computing. David Eppstein
Möbius Transformations in Scientific Computing David Eppstein Univ. of California, Irvine School of Information and Computer Science (including joint work with Marshall Bern from WADS 01 and SODA 03) Outline
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 informationAn Algorithm to Generate Repeating Hyperbolic Patterns
An Algorithm to Generate Repeating Hyperbolic Patterns Douglas Dunham Department of Computer Science University of innesota, Duluth Duluth, N 5581-3036, USA E-mail: ddunham@d.umn.edu Web Site: http://www.d.umn.edu/
More informationShape Modeling. Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell. CS 523: Computer Graphics, Spring 2011
CS 523: Computer Graphics, Spring 2011 Shape Modeling Differential Geometry Primer Smooth Definitions Discrete Theory in a Nutshell 2/15/2011 1 Motivation Geometry processing: understand geometric characteristics,
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 18: Putting Lines Together: Polylines and Polygons Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book
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 informationAn Efficient, Geometric Multigrid Solver for the Anisotropic Diffusion Equation in Two and Three Dimensions
1 n Efficient, Geometric Multigrid Solver for the nisotropic Diffusion Equation in Two and Three Dimensions Tolga Tasdizen, Ross Whitaker UUSCI-2004-002 Scientific Computing and Imaging Institute University
More informationLecture Notes (Reflection & Mirrors)
Lecture Notes (Reflection & Mirrors) Intro: - plane mirrors are flat, smooth surfaces from which light is reflected by regular reflection - light rays are reflected with equal angles of incidence and reflection
More informationChapter 11 Representation & Description
Chain Codes Chain codes are used to represent a boundary by a connected sequence of straight-line segments of specified length and direction. The direction of each segment is coded by using a numbering
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 informationPrinciples of Architectural and Environmental Design EARC 2417 Lecture 2 Forms
Islamic University-Gaza Faculty of Engineering Architecture Department Principles of Architectural and Environmental Design EARC 2417 Lecture 2 Forms Instructor: Dr. Suheir Ammar 2016 1 FORMS ELEMENTS
More informationSelf-formation, Development and Reproduction of the Artificial System
Solid State Phenomena Vols. 97-98 (4) pp 77-84 (4) Trans Tech Publications, Switzerland Journal doi:.48/www.scientific.net/ssp.97-98.77 Citation (to be inserted by the publisher) Copyright by Trans Tech
More informationMedical Image Segmentation using Level Sets
Medical Image Segmentation using Level Sets Technical Report #CS-8-1 Tenn Francis Chen Abstract Segmentation is a vital aspect of medical imaging. It aids in the visualization of medical data and diagnostics
More informationCOMPUTER AND ROBOT VISION
VOLUME COMPUTER AND ROBOT VISION Robert M. Haralick University of Washington Linda G. Shapiro University of Washington A^ ADDISON-WESLEY PUBLISHING COMPANY Reading, Massachusetts Menlo Park, California
More informationLecture notes: Object modeling
Lecture notes: Object modeling One of the classic problems in computer vision is to construct a model of an object from an image of the object. An object model has the following general principles: Compact
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationGeometric Modeling and Processing
Geometric Modeling and Processing Tutorial of 3DIM&PVT 2011 (Hangzhou, China) May 16, 2011 4. Geometric Registration 4.1 Rigid Registration Range Scanning: Reconstruction Set of raw scans Reconstructed
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationEXAM SOLUTIONS. Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006,
School of Computer Science and Communication, KTH Danica Kragic EXAM SOLUTIONS Image Processing and Computer Vision Course 2D1421 Monday, 13 th of March 2006, 14.00 19.00 Grade table 0-25 U 26-35 3 36-45
More informationGeometric Transformations
Geometric Transformations CS 4620 Lecture 9 2017 Steve Marschner 1 A little quick math background Notation for sets, functions, mappings Linear and affine transformations Matrices Matrix-vector multiplication
More informationCOMPUTER VISION > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE
COMPUTER VISION 2017-2018 > OPTICAL FLOW UTRECHT UNIVERSITY RONALD POPPE OUTLINE Optical flow Lucas-Kanade Horn-Schunck Applications of optical flow Optical flow tracking Histograms of oriented flow Assignment
More informationComputational QC Geometry: A tool for Medical Morphometry, Computer Graphics & Vision
Computational QC Geometry: A tool for Medical Morphometry, Computer Graphics & Vision Part II of the sequel of 2 talks. Computation C/QC geometry was presented by Tony F. Chan Ronald Lok Ming Lui Department
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationLecture 16: Voronoi Diagrams and Fortune s Algorithm
contains q changes as a result of the ith insertion. Let P i denote this probability (where the probability is taken over random insertion orders, irrespective of the choice of q). Since q could fall through
More informationLecture 25: Bezier Subdivision. And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10
Lecture 25: Bezier Subdivision And he took unto him all these, and divided them in the midst, and laid each piece one against another: Genesis 15:10 1. Divide and Conquer If we are going to build useful
More informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationSimulation Details for 2D
Appendix B Simulation Details for 2D In this appendix we add some details two-dimensional simulation method. The details provided here describe the method used to obtain results reported in Chapters 3
More informationMass-Spring Systems. Last Time?
Mass-Spring Systems Last Time? Implicit Surfaces & Marching Cubes/Tetras Collision Detection & Conservative Bounding Regions Spatial Acceleration Data Structures Octree, k-d tree, BSF tree 1 Today Particle
More informationGeometry and Gravitation
Chapter 15 Geometry and Gravitation 15.1 Introduction to Geometry Geometry is one of the oldest branches of mathematics, competing with number theory for historical primacy. Like all good science, its
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 informationBiometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong)
Biometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong) References: [1] http://homepages.inf.ed.ac.uk/rbf/hipr2/index.htm [2] http://www.cs.wisc.edu/~dyer/cs540/notes/vision.html
More informationTracking Minimum Distances between Curved Objects with Parametric Surfaces in Real Time
Tracking Minimum Distances between Curved Objects with Parametric Surfaces in Real Time Zhihua Zou, Jing Xiao Department of Computer Science University of North Carolina Charlotte zzou28@yahoo.com, xiao@uncc.edu
More informationSubdivision. Outline. Key Questions. Subdivision Surfaces. Advanced Computer Graphics (Spring 2013) Video: Geri s Game (outside link)
Advanced Computer Graphics (Spring 03) CS 83, Lecture 7: Subdivision Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs83/sp3 Slides courtesy of Szymon Rusinkiewicz, James O Brien with material from Denis
More informationMath background. 2D Geometric Transformations. Implicit representations. Explicit representations. Read: CS 4620 Lecture 6
Math background 2D Geometric Transformations CS 4620 Lecture 6 Read: Chapter 2: Miscellaneous Math Chapter 5: Linear Algebra Notation for sets, functions, mappings Linear transformations Matrices Matrix-vector
More informationSTATISTICS AND ANALYSIS OF SHAPE
Control and Cybernetics vol. 36 (2007) No. 2 Book review: STATISTICS AND ANALYSIS OF SHAPE by H. Krim, A. Yezzi, Jr., eds. There are numerous definitions of a notion of shape of an object. These definitions
More informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationIsophote-Based Interpolation
Isophote-Based Interpolation Bryan S. Morse and Duane Schwartzwald Department of Computer Science, Brigham Young University 3361 TMCB, Provo, UT 84602 {morse,duane}@cs.byu.edu Abstract Standard methods
More informationGeodesics in heat: A new approach to computing distance
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
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 informationSIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE
78 Proceedings of the 4 th International Conference on Informatics and Information Technology SIMULATION OF ARTIFICIAL SYSTEMS BEHAVIOR IN PARAMETRIC EIGHT-DIMENSIONAL SPACE D. Ulbikiene, J. Ulbikas, K.
More informationA Fast and Stable Approach for Restoration of Warped Document Images
A Fast and Stable Approach for Restoration of Warped Document Images Kok Beng Chua, Li Zhang, Yu Zhang and Chew Lim Tan School of Computing, National University of Singapore 3 Science Drive 2, Singapore
More informationDATA MODELS IN GIS. Prachi Misra Sahoo I.A.S.R.I., New Delhi
DATA MODELS IN GIS Prachi Misra Sahoo I.A.S.R.I., New Delhi -110012 1. Introduction GIS depicts the real world through models involving geometry, attributes, relations, and data quality. Here the realization
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 information2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into
2D rendering takes a photo of the 2D scene with a virtual camera that selects an axis aligned rectangle from the scene. The photograph is placed into the viewport of the current application window. A pixel
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
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 informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationEECS150 - Digital Design Lecture 14 FIFO 2 and SIFT. Recap and Outline
EECS150 - Digital Design Lecture 14 FIFO 2 and SIFT Oct. 15, 2013 Prof. Ronald Fearing Electrical Engineering and Computer Sciences University of California, Berkeley (slides courtesy of Prof. John Wawrzynek)
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 informationGaussian and Mean Curvature Planar points: Zero Gaussian curvature and zero mean curvature Tangent plane intersects surface at infinity points Gauss C
Outline Shape Analysis Basics COS 526, Fall 21 Curvature PCA Distance Features Some slides from Rusinkiewicz Curvature Curvature Curvature κof a curve is reciprocal of radius of circle that best approximates
More informationu 0+u 2 new boundary vertex
Combined Subdivision Schemes for the design of surfaces satisfying boundary conditions Adi Levin School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel. Email:fadilev@math.tau.ac.ilg
More informationPrentice Hall Mathematics: Geometry 2007 Correlated to: Arizona Academic Standards for Mathematics (Grades 9-12)
Strand 1: Number Sense and Operations Every student should understand and use all concepts and skills from the previous grade levels. The standards are designed so that new learning builds on preceding
More informationLocal Image Registration: An Adaptive Filtering Framework
Local Image Registration: An Adaptive Filtering Framework Gulcin Caner a,a.murattekalp a,b, Gaurav Sharma a and Wendi Heinzelman a a Electrical and Computer Engineering Dept.,University of Rochester, Rochester,
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 informationPrime Time (Factors and Multiples)
CONFIDENCE LEVEL: Prime Time Knowledge Map for 6 th Grade Math Prime Time (Factors and Multiples). A factor is a whole numbers that is multiplied by another whole number to get a product. (Ex: x 5 = ;
More informationImproved Curvature Estimation on Triangular Meshes
Washington University in St. Louis Washington University Open Scholarship All Computer Science and Engineering Research Computer Science and Engineering Report Number: WUCSE-2004-9 2004-09-01 Improved
More informationAFFINE INVARIANT CURVE MATCHING USING NORMALIZATION AND CURVATURE SCALE-SPACE. V. Giannekou, P. Tzouveli, Y. Avrithis and S.
AFFINE INVARIANT CURVE MATCHING USING NORMALIZATION AND CURVATURE SCALE-SPACE V. Giannekou, P. Tzouveli, Y. Avrithis and S.Kollias Image Video and Multimedia Systems Laboratory, School of Electrical and
More informationTHE preceding chapters were all devoted to the analysis of images and signals which
Chapter 5 Segmentation of Color, Texture, and Orientation Images THE preceding chapters were all devoted to the analysis of images and signals which take values in IR. It is often necessary, however, to
More informationCURVES OF CONSTANT WIDTH AND THEIR SHADOWS. Have you ever wondered why a manhole cover is in the shape of a circle? This
CURVES OF CONSTANT WIDTH AND THEIR SHADOWS LUCIE PACIOTTI Abstract. In this paper we will investigate curves of constant width and the shadows that they cast. We will compute shadow functions for the circle,
More informationLevel set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field
Level Set Methods Overview Level set methods Formulation of Interface Propagation Boundary Value PDE Initial Value PDE Motion in an externally generated velocity field Convection Upwind ddifferencingi
More informationShape fitting and non convex data analysis
Shape fitting and non convex data analysis Petra Surynková, Zbyněk Šír Faculty of Mathematics and Physics, Charles University in Prague Sokolovská 83, 186 7 Praha 8, Czech Republic email: petra.surynkova@mff.cuni.cz,
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 informationOutline 7/2/201011/6/
Outline Pattern recognition in computer vision Background on the development of SIFT SIFT algorithm and some of its variations Computational considerations (SURF) Potential improvement Summary 01 2 Pattern
More informationComparison between Various Edge Detection Methods on Satellite Image
Comparison between Various Edge Detection Methods on Satellite Image H.S. Bhadauria 1, Annapurna Singh 2, Anuj Kumar 3 Govind Ballabh Pant Engineering College ( Pauri garhwal),computer Science and Engineering
More informationCURVILINEAR MESH GENERATION IN 3D
CURVILINEAR MESH GENERATION IN 3D Saikat Dey, Robert M. O'Bara 2 and Mark S. Shephard 2 SFA Inc. / Naval Research Laboratory, Largo, MD., U.S.A., dey@cosmic.nrl.navy.mil 2 Scientific Computation Research
More information