Priority Program 1253
|
|
- Rosamund Ball
- 5 years ago
- Views:
Transcription
1 Deutsche Forschungsgemeinschaft Priority Program 1253 Optimization with Partial Differential Equations N. D. Botkin, K.-H. Hoffmann, N. Mayer, L. Kovacs, M. Eder, S. Raith and A. Volf Application Oriented Methods for Computing Shortest Paths and Principal Curvatures on Triangle Surfaces April 2014 Preprint-Number SPP
2
3 Application oriented methods for computing shortest paths and principal curvatures on triangle surfaces N. D. Botkin 1, K.-H. Hoffmann 1, N. Mayer 2, L. Kovacs 3,4, M. Eder 3,4, S. Raith 3, and A. Volf 3 1 TUM*, Mathematics M6, Garching, Germany 2 TUM, Mechanical engineering LCC, Garching, Germany 3 TUM, Computer Aided Plastic Surgery (CAPS), Garching, Germany 4 TUM, Klinikum rechts der Isar, Klinik und Poliklinik für Plastische Chirurgie und Handchirurgie, München, Germany * TUM = Technische Universität München Abstract In recent years, imaging techniques and three-dimensional (3-D) surface scans in plastic and aesthetic surgery become more and more popular. The data generated by these methods can be used for planning of operations and clinical evaluation of the surgical success. Conventional aesthetic assessments are purely empirical and based on the subjective evaluation through the surgeon or the patient itself, which may be very individual. Commonly accepted objective and quantifiable evaluation approaches do not exist. Preliminary objective 3-D evaluation methods are based on the simple rightto-left side comparison of distances, volumes, or surface areas. It is desirable to develop an objective evaluation index that takes into account a set of relevant measurements and collect them with weights in the so-called Symmetry Index. This would allow to use the objective analysis of surgical results, which would make the doctor-patient communication more constructive. It should be noticed that the Symmetry Index should contain the intrinsic characteristics of 3-D surfaces coming from the topology and differential geometry such as geodesic distances and curvatures. Another important problem is setting characteristic anatomical points (called landmarks) on 3-D face and body scans to estimate the surgical results. The problem here botkin@ma.tum.de hoffmann@ma.tum.de mayer@lcc.mw.tum.de l.kovacs@lrz.tum.de m.eder@lrz.tum.de stefan.raith@caps.me.tum.de alexander.volf@caps.me.tum.de
4 is that the process of landmark positioning is manual, which is very cumbersome and time-consuming, depends on the particular evaluator, and therefore is only partially reproducible. Thus, an objective, standardized, and reproducible method for the fully automatic anthropometric analysis of three-dimensional data sets is required. There is the believe that a stable algorithm for the computation of curvature tensors can solve this challenging problem. Finally, the problem of sparse representation of aesthetic regions by semantic surface models containing additional meta-information about the anatomical peculiarities of these regions is important and challenging. The development of such semantic surface models requires a semi-analytical description of triangle surface regions. Here, local maps and parametric shape functions can be used. Fitting these objects to the specific regions is then carried out by the variation of the corresponding parameter. The decomposition of the surface into the aesthetic regions and construction of the shape functions should be based on the analysis of curvature. Roughly speaking, the areas with the main curvature signatures (+ +) - convex, (- -) - concave, and (- +) - saddle-type, have to be distinguished. The calculation of semantic basis functions can be done using the kriging method that analytically approximates large data sets with low numerical effort. This paper considers two problems: computation of shortest paths on triangle surfaces and calculation of principal curvatures at vertices of triangle surfaces. Existing methods compute shortest paths along edges so that the result may be completely unacceptable. We propose a technique that allows us to compute smoother and shorter paths using intermediate admissible directions and linear interpolation of potential functions. Conventional methods that compute the curvature in nodes of triangle surfaces are based on processing only adjacent vertices. They are non stable because triangle meshes contain frequently oscillating errors in positions of the vertices. We propose a simple practical method that accounts for larger neighborhoods of vertices. The kriging method aimed to the sparse representation of specific surface parts (aesthetic regions) is sketched. Preliminary numerical results are presented. 1 Introduction In the last decade, discrete differential geometry became very popular because of its application to visualization of technical and medical objects. The basic aim of the discrete differential geometry is to extend geometrical characteristics inherent to smooth manifolds to the case of polyhedron surfaces. We consider two problems: computation of shortest paths on triangle surfaces and calculation of principal curvatures at vertices of triangle surfaces. There are many algorithms for computing shortest paths between two points on graphs (see e.g. a very exhaustive survey [1]). Nevertheless, these methods are not perfectly applicable to practical problems because shortest paths along edges may be very jagged. Moreover, such paths may lie far away from geodesic lines impressed by the observer. We propose a technique that allows us to compute smoother and shorter paths using intermediate admissible directions and linear interpolation of the potential function. This approach is related to the numerical treatment of time independent Hamilton-Jcobi equations arising from optimal 2
5 control theory. As for computing the curvature in nodes of a triangle surface, there are several fundamental principles of extension of the notion of curvature to discrete objects. For example, the Gaussian total curvature at a vertex is given by K = 2π j α j, where α j are the angles between the edges adjacent to the vertex. The mean curvature vector at a vertex p i is given by the cotan formula H i = 1/2 j (cot α ij + cot α ji )(p i p j ) (Pinkall and Polthier, see e.g. [4]), where α ij and α ji are angles opposite to the edge p i p j. There are works that propose formulaes for computing the principal curvatures at a vertex using the edges adjacent to this vertex (see e.g. [2, 3]). Nevertheless, these results are not perfectly applicable to practical problems where triangle surfaces arise from scans, say, of human faces. Such scans yield clouds of points that are processed to obtain triangle surfaces representing human faces. It is clear that such triangle meshes contain frequently oscillating errors in positions of the vertices, which restricts application of local methods that use only adjacent vertices. We propose a simple practical method that accounts for larger neighborhoods of vertices. 2 Computing a shortest path between two points 2.1 Shortest path along edges In this subsection, basic ideas of computing shortest paths are outlined. Consider a triangle surface and denote the set of all vertices (resp. all points) by S (resp. by Σ). The majority of existing algorithms are based on the direct or indirect construction of the potential function d(s), s S, that express the distance between any vertex s and a fixed vertex s f along edges (see Fig. 1). Thus, d(s) = dist(s, s f ) along edges. It is easily to prove that the potential function satisfies the following optimality condition: d(s f ) = 0, d(s) = min l L(s) {d(s + l) + l }, s s f. (1) Here, L(s) denotes the set of all edges adjacent to the vertex s, i.e. the set of all admissible directions (see Fig. 1), l denotes the length of the edge l. Such a nonlinear equation can be solved using the following iterative relaxsation procedure: d 0 (s) = 0 d n (s f ) = 0 d n+1 (s) = min l L(s) {dn (s + l) + l }, s s f. (2) s, n, 3
6 L(s) = {l 1, l 2,..., l 5 } Figure 1: A triangle surface with a selected vertex s f. The set L(s) consists of edges adjacent to s. It easily to prove the following properties of the sequence d n (s): i) d n (s) d(s) s, ii) d n+1 (s) d n (s) s. These properties follow from the monotonicity of the operator q Π[q] defined as Π[q](s) := min l L(s) {q(s + l) + l }, i.e. if q 1 (s) q 2 (s) for all s, then Π[q 1 ](s) Π[q 2 ](s) for all s. Thus, conditions i) and ii) together with relations (1) and (2) imply that d n (s) monotonically converges to d(s) for all s. Construction of a shortest path between s and s f is shown in Fig. 2. Assume that the path is already constructed to the vertex s i. To obtain the next vertex, chose an optimal direction l min as a minimizer in min {d(s i +l)+ l } and set s i+1 = s i +l min. If the minimizer l L(s i ) is not unique, choose anyone. Remark 1 Note that the potential function d( ) is used only in a small neighborhood of a shortest path. Therefore, the construction of the potential function can be restricted to a subregion containing all vertices of a shortest path, which essentially speedups the calculation. Moreover, the computations at different vertices in (2) are completely independent, and therefore very efficient parallelization is possible. 4
7 Figure 2: Construction of a shortest path along edges. 2.2 Shortest path in the case of large sets of admissible directions Now, large sets of admissible directions will be used to compute a modified potential function at the vertices. Fig. 3 shows different ways of choosing admissible directions. On the picture to the left, the ends of admissible directions lie on the edges of the cluster about s, whereas, on the picture to the right, the ends of admissible directions lie inside the cluster. In both cases, if the end point of an admissible direction does not meet a vertex, the linear interpolation is used to compute a value in such a point. Figure 3: Choice of admissible directions. The ends of admissible directions may lie on the edges of the cluster about s (see the picture to the left), or they may lie inside the cluster (see the picture to the right). 5
8 Consider the following iterative procedure: d 0 (s) = 0 d n (s f ) = 0 d n+1 (s) = min l L(s) { d n (s + l) + l }, s s f, (3) where L(s), s S, are large sets of admissible directions, and denotes the operation of linear interpolation of a grid function. It is easily to prove that the sequence d n (s) satisfies properties i) and ii), and therefore d n (s) monotonically converges to some ˆd(s) for all s. Construction of a shortest path between s and s f is explained in Fig. 4. Assume that the path is already constructed to a point s i (not necessary a vertex). To obtain the next point, chose an optimal direction l min as a minimizer in min l L(s i ) { ˆd(si + l) + l } and set s i+1 = s i + l min. Here, L(s), s Σ, are large sets of admissible directions when constructing the path. Remember that Σ is the set of all points of the triangle surface. s, n, Figure 4: Construction of a shortest path between s and s f using intermediate directions and linear interpolation of the potential function ˆd. 3 Computing principal curvatures Assume that the principal curvatures are to be computed at a vertex s labeled by a large green circle in Fig. 5. Consider the first order cluster consisting of triangles adjacent to s (the bounding vertices of this cluster are shown in red). The second order cluster consists 6
9 of all triangles adjacent to the first order cluster (the bounding vertices are shown in blue). Moreover, it is easily to define the third order cluster, and so on. Compute a normal vector n s at s as a weighted sum of the normal vectors to the triangles adjacent to s and draw a cutting plane passing through n s. It is clear that such a plane is uniquely defined by the rotation angle about n s. Let a value of the rotation angle be fixed. To compute the curvature corresponding to this angle, find all intersection points of the cutting plane with the edges of the first, second, and third order clusters (such points are shown in yellow in Fig. 5). Consider the cutting plane and the intersection points on it (see Fig. 6), chose a local reference system in the cutting plane, and construct a polynomial (reasonable degrees are 2,3,4, and 5) that provides the best approximation of the intersection points in the sense of weighted minimum squares. That is, P (x) = a 0 + a 1 x + a 2 x a m x m, and the coefficients are obtained from the condition min a 0,a 1,...a m N ) 2, w j (P (x j ) y j j=1 where w j are some weights, and (x j, y j ) are the coordinates of the intersection points in the local reference system chosen on the cutting plane. The curvature of the polynomial graph is then given by k = 2a 2 /(1 + a 2 1). This value approximates the curvature of the triangle surface at s for a fixed rotation angle of the cutting plane. Rotating the cutting plane about the normal n s, we compute the principal (maximal and minimal) curvatures. 7
10 n s s s Figure 5: The first and second order clusters (red and blue vertices, respectively) about the vertex s, a normal vector n s at s, and the intersection points (yellow) of the cutting plane and the edges of the clusters. Figure 6: Two cutting planes, the corresponding intersection points, and the graphs of approximating polynomials are shown. The picture to the left (resp. to the right) corresponds to the maximal (resp. minimal) value of the curvature. 4 Preliminary numerical results Figure 7 shows results of computing shortest paths. The red line represents a shortest path along edges. The picture at the upper right explains the appearance of jags in the path. It is seen that there are not good admissible direction in the shown area. The green line is constructed using a global refinement of the triangle grid. In this case, the computation time essentially increases. The blue line is constructed using larger sets of admissible directions. 8
11 The picture at the down right shows that the path can cross edges in this case. Note that this path is 15% shorter than the path along edges and really looks like the geodesic line. The method works also well for relatively rough triangle meshes (see Fig. 8). Figure 7: Computation of shortest paths on a triangle surface: along edges (the red line), along edges but using a global refinement of the triangle mesh (the green line), and using larger sets of admissible directions (the blue line). The numbers of vertices and triangles of the surface are about and , respectively. 9
12 Figure 8: Computation of shortest paths on a relatively rough triangle surface. A shortest path along edges is given by the red line, and the green line shows a shortest path computed using larger sets of admissible directions. The number of vertices of the surface is about Figure 9 presents numerical experiments with rough, inhomogeneous, and very distorted triangle meshes. The picture to the left shows two shortest paths along edges. The corresponding paths computed using larger sets of admissible directions are shown in the picture in the middle. Note that the red line is essentially better in the case of larger sets of admissible directions, whereas the green line avoids to cross large distorted triangles. This effect is caused by the interpolation error that is large in the case of big distorted triangles. The way to improve that is related to a modification of the optimality principal when the constructed path meets a big distorted triangle. In the case of relatively homogenous meshes, the construction of shortest paths using larger sets of admissible directions gives good results even in the case of very distorted triangles (see the picture to the right). 10
13 Figure 9: Numerical experiments with rough, inhomogeneous, and very distorted triangle meshes. Shortest paths along edges (to the left) and the corresponding paths computed using larger sets of admissible directions (in the middle) are presented. The computation of a shortest path for a relatively homogeneous but very distorted triangle mesh using larger sets of admissible directions is shown on the picture to the right. Figures 10, and 11 show results of computing the curvature. Figure 10: Computing the maximal (to the left) and mean (to the right) curvatures in the case of an asymmetric face. 11
14 saddle type concavity convexity Figure 11: Computing the Gaussian curvature (to the left) and the classification of vertices according to the signs of the principal curvatures (to the right). Figure 12: Landmarks (see [5]) describing indicative points of a face. 12
15 Figure 13: Distance measurement between landmarks as described by Farkas [5] with the new shortest path method. Figure 14: Localization of the left mouth landmark with the minimum of the maximal curvature. 13
16 Figure 15: Localization of the right mouth landmark with the minimum of the maximal curvature. Figure 16: The result of the localization of the mouth landmarks. 14
17 Surface representation using kriging models Reconstructing free-form surfaces from 3-D measurements is a central problem in many applications. A vary appropriate and commonly used tool here are kriging models. It should be noticed that these models have been developed as estimators of randomly disturbed data. Nevertheless, they work well with deterministic computer models (see e.g. [6]). Kriging models show very good flexibility to approximate many different and complex response functions. They are also a good choice for analytic representation of grid surfaces using few parameters. The mathematical form of a kriging model with, e.g. scalar response, see (4), has two parts. The first part is a linear regression of the data with k regressors modeling the trend over the domain. Most engineering applications use a constant trend model and rely upon the second part of the model to pull the response surface through the observed data by quantifying the correlation of nearby points. ŷ( X) = k β i f i ( X) + Z( X). (4) i=1 Here X is a p-dimensional input vector, e.g. a set of parameters of a computer PDE-based model, f i known functions (polynomials), β i regression coefficients, and Z( X) a model of a Gaussian and stationary random process with zero mean and covariance: cov(x 1, x 2 ) = σ 2 R(x 1, x 2 ). (5) The process variance, σ 2, scales the spatial correlation function R(x 1, x 2 ) that controls the smoothness of the resulting kriging model, the influence of other nearby points, and the differentiability of the response surface by quantifying the correlation between the two observations x 1 and x 2. There are many potential functions that can be used to quantify the correlation between the observations. For example, R(x 1, x 2, θ) = exp ( θ (x 1 x 2 ) 2), θ > 0. (6) The range parameter, 1/ θ, indicates the distance at which the influence equals e 1 = or approximately 37%. To construct a multivariate correlation function, a univariate correlation function is used for each of the p input dimensions, and a product correlation rule is used: p R( X 1, X 2, θ) = R(x 1,i, x 2,i, θ i ). (7) Therefore, a kriging model contains the unknown parameters i=1 γ = { β, σ, θ} to be determined from the n-vector, Y, of observations. The most popular method of estimating kriging model parameters is Maximum Likelihood Estimation (see [6]). The logarithm of 15
18 the multivariate Gaussian likelihood function is given by: l( γ Y ) = n 2 ln[2πσ2 ] 1 2 ln[ R ] 1 2σ 2 ( Y F β) T R 1 ( Y F β), (8) where R is a matrix quantifying the correlation of all of the observations, R ij = R( X i, X j, θ), i, j = 1,...n, and F a matrix of f( X) evaluated at each observation, F ij = f j (X i ), i = 1...n, j = 1...k. Minimization of this function yields relations β = (F T R 1 F ) 1 F T R 1 Y, σ 2 = 1 n ( Y F β) T ( Y F β). (9) Notice that the matrix R in (9) depends on the parameter set γ. To obtain a closed solution, the relations (9) can be substituted into (8), and the last can be minimized over θ. This procedure is being iteratively repeated (see [6]). There is a wide list of software for constructing kriging models (see e.g. [7] for the toolbox DACE). This toolbox provides regression models with polynomials of orders 0, 1, and 2 and several correlation models: Gaussian, see (6), exponential, spherical, cubic, etc. The next figure shows the usage of a kriging estimator for the prediction of outputs of a physical model. The toolbox DACE (see [7]) was utilized for the construction of a kriging estimator. Figure 17: Three surfaces representing outputs of a physical model. Each surface is computed using a kriging estimator whose regression part utilizes polynomials of the first order. The kriging estimator itself is constructed on the basis of 20 observations, simulations of the physical model with different parameter sets. The black dots represent the values obtained from the model simulations but not involved into the construction of the kriging estimator. One can see a very good prediction. 16
19 Conclusion The results presented in this paper are targeted for medical applications. First, the computation of smooth shortest paths on triangle surfaces is supposed to be used for the development of a Symmetry Index providing the objective analysis of surgical results. The measurement of the geodesic distances between characteristic anatomical points, landmarks, on 3-D face and body scans provides intrinsic characteristics of aesthetic regions. Second, the computation of curvatures for sampled directions at each node allows us to develop an objective, standardized, and reproducible method for fully automatic positioning of land marks. And finally, the problem of sparse representation of aesthetic regions by semantic surface models can be solved using the decomposition of the surface into curvature specific regions (convex, concave, and saddle-type) and applying kriging models to represent them analytically. References [1] Giorgio Gallo. Shortest Path Algorithms. Annals of Operations Research 13(1988)3-79 [2] Gabriel Taubin. Estimating the tensor of curvature of a surface from a polyhedral approximation. In Proceedings of Fifth International Conference on Computer Vision (ICCV 95), [3] Meyer, M., Desbrun, M., Schröder, P. and Barr A.H. (2003), Discrete differentialgeometry operators for triangulated 2-manifolds. In H.-C. Hege and K. Polthier, editors, Visualization and Mathematics III, pp , Springer-Verlag, Heidelberg. [4] Max Wardetzky. Convergence of the Cotangent Formula: An Overview (p ). In. Discrete Differential Geometry (ebook). Alexander I. Bobenko, Peter Schröder, John M. Sullivan, Günter M. Ziegler (Eds.) Birkhäuser [5] L. Kovacs, A. Zimmermann, G. Brockmann, H. Baurecht, K. Schwenzer-Zimmerer, N. A. Papadopulos, M. A. Papadopoulos, R. Sader, E. Biemer, and H. F. Zeilhofer, Accuracy and Precision of the Three-Dimensional Assessment of the Facial Surface Using a 3-D Laser Scanner, IEEE Transactions on medical imaging, 25 (2006) No. 6. [6] J. D. Martin and T. W. Simpson. Use of Kriging Models to Approximate Deterministic Computer Models. AIAA Journal, 43(4): , [7] Søren N Lophaven, Hans Bruun Nielsen, Jacob Søndergaard. DACE a Matlab Kriging Toolbox. Technical report, Technical University of Denmark, Informatics and Mathematical Modelling,
Shape 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 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 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 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 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 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 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 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 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 informationMULTIVARIATE TEXTURE DISCRIMINATION USING A PRINCIPAL GEODESIC CLASSIFIER
MULTIVARIATE TEXTURE DISCRIMINATION USING A PRINCIPAL GEODESIC CLASSIFIER A.Shabbir 1, 2 and G.Verdoolaege 1, 3 1 Department of Applied Physics, Ghent University, B-9000 Ghent, Belgium 2 Max Planck Institute
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 informationGeodesic Paths on Triangular Meshes
Geodesic Paths on Triangular Meshes Dimas Martínez Luiz Velho Paulo Cezar Carvalho IMPA Instituto Nacional de Matemática Pura e Aplicada Estrada Dona Castorina, 110, 22460-320 Rio de Janeiro, RJ, Brasil
More 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 informationFree-Form Shape Optimization using CAD Models
Free-Form Shape Optimization using CAD Models D. Baumgärtner 1, M. Breitenberger 1, K.-U. Bletzinger 1 1 Lehrstuhl für Statik, Technische Universität München (TUM), Arcisstraße 21, D-80333 München 1 Motivation
More informationRecent advances in Metamodel of Optimal Prognosis. Lectures. Thomas Most & Johannes Will
Lectures Recent advances in Metamodel of Optimal Prognosis Thomas Most & Johannes Will presented at the Weimar Optimization and Stochastic Days 2010 Source: www.dynardo.de/en/library Recent advances in
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 informationGeodesic Paths on Triangular Meshes
Geodesic Paths on Triangular Meshes DIMAS MARTÍNEZ, LUIZ VELHO, PAULO CEZAR CARVALHO IMPA Instituto Nacional de Matemática Pura e Aplicada - Estrada Dona Castorina, 110, 22460-320 Rio de Janeiro, RJ, Brasil
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 Tobler s Law All places are related, but nearby places are related more than distant places Corollary:
More informationSUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS
SUBDIVISION ALGORITHMS FOR MOTION DESIGN BASED ON HOMOLOGOUS POINTS M. Hofer and H. Pottmann Institute of Geometry Vienna University of Technology, Vienna, Austria hofer@geometrie.tuwien.ac.at, pottmann@geometrie.tuwien.ac.at
More informationSpatial Interpolation - Geostatistics 4/3/2018
Spatial Interpolation - Geostatistics 4/3/201 (Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Distance between pairs of points Lag Mean Tobler s Law All places are related, but nearby places
More informationComparing different interpolation methods on two-dimensional test functions
Comparing different interpolation methods on two-dimensional test functions Thomas Mühlenstädt, Sonja Kuhnt May 28, 2009 Keywords: Interpolation, computer experiment, Kriging, Kernel interpolation, Thin
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
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 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 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 informationDijkstra's Algorithm
Shortest Path Algorithm Dijkstra's Algorithm To find the shortest path from the origin node to the destination node No matrix calculation Floyd s Algorithm To find all the shortest paths from the nodes
More informationA GENTLE INTRODUCTION TO THE BASIC CONCEPTS OF SHAPE SPACE AND SHAPE STATISTICS
A GENTLE INTRODUCTION TO THE BASIC CONCEPTS OF SHAPE SPACE AND SHAPE STATISTICS HEMANT D. TAGARE. Introduction. Shape is a prominent visual feature in many images. Unfortunately, the mathematical theory
More informationMATH3016: OPTIMIZATION
MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is
More informationDiscrete Optimization. Lecture Notes 2
Discrete Optimization. Lecture Notes 2 Disjunctive Constraints Defining variables and formulating linear constraints can be straightforward or more sophisticated, depending on the problem structure. The
More informationAccurate Fitting of Noisy Irregular Beam Data for the Planck Space Telescope
Proceedings of the 5th European Conference on Antennas and Propagation (EUCAP) Accurate Fitting of Noisy Irregular Beam Data for the Planck Space Telescope Oscar Borries, Per Heighwood Nielsen, Jan Tauber,
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationElastic registration of medical images using finite element meshes
Elastic registration of medical images using finite element meshes Hartwig Grabowski Institute of Real-Time Computer Systems & Robotics, University of Karlsruhe, D-76128 Karlsruhe, Germany. Email: grabow@ira.uka.de
More informationA Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
Scientific Papers, University of Latvia, 2010. Vol. 756 Computer Science and Information Technologies 207 220 P. A Modified Spline Interpolation Method for Function Reconstruction from Its Zero-Crossings
More informationNormals of subdivision surfaces and their control polyhedra
Computer Aided Geometric Design 24 (27 112 116 www.elsevier.com/locate/cagd Normals of subdivision surfaces and their control polyhedra I. Ginkel a,j.peters b,,g.umlauf a a University of Kaiserslautern,
More informationCurvature Berkeley Math Circle January 08, 2013
Curvature Berkeley Math Circle January 08, 2013 Linda Green linda@marinmathcircle.org Parts of this handout are taken from Geometry and the Imagination by John Conway, Peter Doyle, Jane Gilman, and Bill
More informationRecovering Primitives in 3D CAD meshes
Recovering Primitives in 3D CAD / R.Bénière 1/17 Recovering Primitives in 3D CAD meshes Roseline Bénière G. Subsol, G. Gesquière, F. Le Breton and W. Puech LIRMM, Montpellier, France C4W, Montpellier,
More informationCSE 554 Lecture 6: Fairing and Simplification
CSE 554 Lecture 6: Fairing and Simplification Fall 2012 CSE554 Fairing and simplification Slide 1 Review Iso-contours in grayscale images and volumes Piece-wise linear representations Polylines (2D) and
More informationLocally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling
Locally Weighted Least Squares Regression for Image Denoising, Reconstruction and Up-sampling Moritz Baecher May 15, 29 1 Introduction Edge-preserving smoothing and super-resolution are classic and important
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 informationConjectures concerning the geometry of 2-point Centroidal Voronoi Tessellations
Conjectures concerning the geometry of 2-point Centroidal Voronoi Tessellations Emma Twersky May 2017 Abstract This paper is an exploration into centroidal Voronoi tessellations, or CVTs. A centroidal
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 informationA Geometric Approach to the Bisection Method
Geometric pproach to the isection Method laudio Gutierrez 1, Flavio Gutierrez 2, and Maria-ecilia Rivara 1 1 epartment of omputer Science, Universidad de hile lanco Encalada 2120, Santiago, hile {cgutierr,mcrivara}@dcc.uchile.cl
More informationApproximation of 3D-Parametric Functions by Bicubic B-spline Functions
International Journal of Mathematical Modelling & Computations Vol. 02, No. 03, 2012, 211-220 Approximation of 3D-Parametric Functions by Bicubic B-spline Functions M. Amirfakhrian a, a Department of Mathematics,
More informationDiscrete differential geometry: Surfaces made from Circles
Discrete differential geometry: Alexander Bobenko (TU Berlin) with help of Tim Hoffmann, Boris Springborn, Ulrich Pinkall, Ulrike Scheerer, Daniel Matthes, Yuri Suris, Kevin Bauer Papers A.I. Bobenko,
More informationBayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis
Bayesian Spherical Wavelet Shrinkage: Applications to Shape Analysis Xavier Le Faucheur a, Brani Vidakovic b and Allen Tannenbaum a a School of Electrical and Computer Engineering, b Department of Biomedical
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 informationNote Set 4: Finite Mixture Models and the EM Algorithm
Note Set 4: Finite Mixture Models and the EM Algorithm Padhraic Smyth, Department of Computer Science University of California, Irvine Finite Mixture Models A finite mixture model with K components, for
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationModular network SOM : Theory, algorithm and applications
Modular network SOM : Theory, algorithm and applications Kazuhiro Tokunaga and Tetsuo Furukawa Kyushu Institute of Technology, Kitakyushu 88-96, Japan {tokunaga, furukawa}@brain.kyutech.ac.jp Abstract.
More informationGeneralized Principal Component Analysis CVPR 2007
Generalized Principal Component Analysis Tutorial @ CVPR 2007 Yi Ma ECE Department University of Illinois Urbana Champaign René Vidal Center for Imaging Science Institute for Computational Medicine Johns
More informationA Data Dependent Triangulation for Vector Fields
A Data Dependent Triangulation for Vector Fields Gerik Scheuermann Hans Hagen Institut for Computer Graphics and CAGD Department of Computer Science University of Kaiserslautern, Postfach 3049, D-67653
More informationNumerical Treatment of Geodesic Differential. Equations on a Surface in
International Mathematical Forum, Vol. 8, 2013, no. 1, 15-29 Numerical Treatment of Geodesic Differential Equations on a Surface in Nassar H. Abdel-All Department of Mathematics, Faculty of Science Assiut
More informationGrade 9 Math Terminology
Unit 1 Basic Skills Review BEDMAS a way of remembering order of operations: Brackets, Exponents, Division, Multiplication, Addition, Subtraction Collect like terms gather all like terms and simplify as
More information6.2 Classification of Closed Surfaces
Table 6.1: A polygon diagram 6.1.2 Second Proof: Compactifying Teichmuller Space 6.2 Classification of Closed Surfaces We saw that each surface has a triangulation. Compact surfaces have finite triangulations.
More informationHyperbolic structures and triangulations
CHAPTER Hyperbolic structures and triangulations In chapter 3, we learned that hyperbolic structures lead to developing maps and holonomy, and that the developing map is a covering map if and only if the
More informationLAPLACIAN MESH SMOOTHING FOR TETRAHEDRA BASED VOLUME VISUALIZATION 1. INTRODUCTION
JOURNAL OF MEDICAL INFORMATICS & TECHNOLOGIES Vol.4/2002, ISSN 642-6037 Rafał STĘGIERSKI *, Paweł MIKOŁAJCZAK * volume data,triangle mesh generation, mesh smoothing, marching tetrahedra LAPLACIAN MESH
More informationComparison and affine combination of generalized barycentric coordinates for convex polygons
Annales Mathematicae et Informaticae 47 (2017) pp. 185 200 http://ami.uni-eszterhazy.hu Comparison and affine combination of generalized barycentric coordinates for convex polygons Ákos Tóth Department
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 information7. The Gauss-Bonnet theorem
7. The Gauss-Bonnet theorem 7.1 Hyperbolic polygons In Euclidean geometry, an n-sided polygon is a subset of the Euclidean plane bounded by n straight lines. Thus the edges of a Euclidean polygon are formed
More informationTensor Based Approaches for LVA Field Inference
Tensor Based Approaches for LVA Field Inference Maksuda Lillah and Jeff Boisvert The importance of locally varying anisotropy (LVA) in model construction can be significant; however, it is often ignored
More informationAuto-Extraction of Modelica Code from Finite Element Analysis or Measurement Data
Auto-Extraction of Modelica Code from Finite Element Analysis or Measurement Data The-Quan Pham, Alfred Kamusella, Holger Neubert OptiY ek, Aschaffenburg Germany, e-mail: pham@optiyeu Technische Universität
More informationPlanar Symmetry Detection by Random Sampling and Voting Process
Planar Symmetry Detection by Random Sampling and Voting Process Atsushi Imiya, Tomoki Ueno, and Iris Fermin Dept. of IIS, Chiba University, 1-33, Yayo-cho, Inage-ku, Chiba, 263-8522, Japan imiya@ics.tj.chiba-u.ac.jp
More informationDIGITAL TERRAIN MODELS
DIGITAL TERRAIN MODELS 1 Digital Terrain Models Dr. Mohsen Mostafa Hassan Badawy Remote Sensing Center GENERAL: A Digital Terrain Models (DTM) is defined as the digital representation of the spatial distribution
More informationShape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow
Shape-based Diffeomorphic Registration on Hippocampal Surfaces Using Beltrami Holomorphic Flow Abstract. Finding meaningful 1-1 correspondences between hippocampal (HP) surfaces is an important but difficult
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 informationCIS 520, Machine Learning, Fall 2015: Assignment 7 Due: Mon, Nov 16, :59pm, PDF to Canvas [100 points]
CIS 520, Machine Learning, Fall 2015: Assignment 7 Due: Mon, Nov 16, 2015. 11:59pm, PDF to Canvas [100 points] Instructions. Please write up your responses to the following problems clearly and concisely.
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 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 informationPlane Based Free Stationing for Building Models
Christian MANTHE, Germany Key words: plane based building model, plane detection, plane based transformation SUMMARY 3D Building models are used to construct, manage and rebuild buildings. Thus, associated
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 informationData Representation in Visualisation
Data Representation in Visualisation Visualisation Lecture 4 Taku Komura Institute for Perception, Action & Behaviour School of Informatics Taku Komura Data Representation 1 Data Representation We have
More informationClustering. Robert M. Haralick. Computer Science, Graduate Center City University of New York
Clustering Robert M. Haralick Computer Science, Graduate Center City University of New York Outline K-means 1 K-means 2 3 4 5 Clustering K-means The purpose of clustering is to determine the similarity
More informationWeighted Neighborhood Sequences in Non-Standard Three-Dimensional Grids Parameter Optimization
Weighted Neighborhood Sequences in Non-Standard Three-Dimensional Grids Parameter Optimization Robin Strand and Benedek Nagy Centre for Image Analysis, Uppsala University, Box 337, SE-7505 Uppsala, Sweden
More informationTwo-Dimensional Fitting of Brightness Profiles in Galaxy Images with a Hybrid Algorithm
Two-Dimensional Fitting of Brightness Profiles in Galaxy Images with a Hybrid Algorithm Juan Carlos Gomez, Olac Fuentes, and Ivanio Puerari Instituto Nacional de Astrofísica Óptica y Electrónica Luis Enrique
More informationDesigning Cylinders with Constant Negative Curvature
Designing Cylinders with Constant Negative Curvature Ulrich Pinkall Abstract. We describe algorithms that can be used to interactively construct ( design ) surfaces with constant negative curvature, in
More information3.3 Interpolation of rotations represented by quaternions
3 S 3.3 Interpolation of rotations represented by quaternions : set of unit quaternions ; S :set of unit 3D vectors Interpolation will be performed on 3 S (direct linear interpolation produces nonlinear
More informationNorbert Schuff VA Medical Center and UCSF
Norbert Schuff Medical Center and UCSF Norbert.schuff@ucsf.edu Medical Imaging Informatics N.Schuff Course # 170.03 Slide 1/67 Objective Learn the principle segmentation techniques Understand the role
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationClipping Polygons. A routine for clipping polygons has a variety of graphics applications.
The Mathematica Journal Clipping Polygons Garry Helzer Department of Mathematics University of Maryland College Park, MD 20742 gah@math.umd.edu A routine for clipping polygons has a variety of graphics
More informationParameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia
Parameterization II Some slides from the Mesh Parameterization Course from Siggraph Asia 2008 1 Non-Convex Non Convex Boundary Convex boundary creates significant distortion Free boundary is better 2 Fixed
More informationTwo Counterexamples of Global Differential Geometry for Polyhedra
Two Counterexamples of Global Differential Geometry for Polyhedra Abdênago Barros, Esdras Medeiros and Romildo Silva Departatamento de Matemática Universidade Federal do Ceará 60450-520 Fortaleza-Ceará,
More informationChapter 4. Clustering Core Atoms by Location
Chapter 4. Clustering Core Atoms by Location In this chapter, a process for sampling core atoms in space is developed, so that the analytic techniques in section 3C can be applied to local collections
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 informationing and enhancing operations for shapes represented by level sets (isosurfaces) and applied unsharp masking to the level set normals. Their approach t
Shape Deblurring with Unsharp Masking Applied to Mesh Normals Hirokazu Yagou Λ Alexander Belyaev y Daming Wei z Λ y z ; ; Shape Modeling Laboratory, University of Aizu, Aizu-Wakamatsu 965-8580 Japan fm50534,
More 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 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 informationAlgorithm research of 3D point cloud registration based on iterative closest point 1
Acta Technica 62, No. 3B/2017, 189 196 c 2017 Institute of Thermomechanics CAS, v.v.i. Algorithm research of 3D point cloud registration based on iterative closest point 1 Qian Gao 2, Yujian Wang 2,3,
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 informationAnna Porazilová THE SHORTEST PATH. 1 Introduction 25. KONFERENCE O GEOMETRII A POČÍTAČOVÉ GRAFICE
25. KONFERENCE O GEOMETRII A POČÍTAČOVÉ GRAFICE Anna Porazilová THE SHORTEST PATH Abstract This paper describes the shortest path problem, its classification and the best known algorithms. A new algorithm
More informationGeostatistics 2D GMS 7.0 TUTORIALS. 1 Introduction. 1.1 Contents
GMS 7.0 TUTORIALS 1 Introduction Two-dimensional geostatistics (interpolation) can be performed in GMS using the 2D Scatter Point module. The module is used to interpolate from sets of 2D scatter points
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 informationFeature Extraction for Illustrating 3D Stone Tools from Unorganized Point Clouds
Feature Extraction for Illustrating 3D Stone Tools from Unorganized Point Clouds Enkhbayar Altantsetseg 1) Yuta Muraki 2) Katsutsugu Matsuyama 2) Fumito Chiba 3) Kouichi Konno 2) 1) Graduate School of
More informationOn Smooth Bicubic Surfaces from Quad Meshes
On Smooth Bicubic Surfaces from Quad Meshes Jianhua Fan and Jörg Peters Dept CISE, University of Florida Abstract. Determining the least m such that one m m bi-cubic macropatch per quadrilateral offers
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 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(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 informationStudy of Panelization Techniques to Inform Freeform Architecture
Study of Panelization Techniques to Inform Freeform Architecture Daniel Hambleton, Crispin Howes, Jonathan Hendricks, John Kooymans Halcrow Yolles Keywords 1 = Freeform geometry 2 = Planar quadrilateral
More informationSeparation of Surface Roughness Profile from Raw Contour based on Empirical Mode Decomposition Shoubin LIU 1, a*, Hui ZHANG 2, b
International Conference on Advances in Mechanical Engineering and Industrial Informatics (AMEII 2015) Separation of Surface Roughness Profile from Raw Contour based on Empirical Mode Decomposition Shoubin
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 information