Solving the initial value problem of discrete geodesics. (January 2016)

• Record Type: Journal Article
• Authors: Cheng, Peng; Miao, Chunyan; Liu, Yong-Jin; Tu, Changhe; He, Ying
• Is Part Of: Computer aided design. Volume 70(2016)
• Page Start: 144
• Page End: 152
• Publication Date: 2016-01
• Subjects: Discrete geodesics -- Initial value problem -- Geodesic polar coordinates; Computer-aided design -- Periodicals; Engineering design -- Data processing -- Periodicals; Computer graphics -- Periodicals; Conception technique -- Informatique -- Périodiques; Infographie -- Périodiques; Computer graphics; Engineering design -- Data processing; Periodicals; Electronic journals; 620.00420285
• DOI: 10.1016/j.cad.2015.07.012
• Languages: English
• ISSNs: 0010-4485

Abstract

Abstract: Computing geodesic paths and distances is a common operation in computer graphics and computer-aided geometric design. The existing discrete geodesic algorithms are mainly designed to solve the boundary value problem, i.e., to find the shortest path between two given points. In this paper, we focus on the initial value problem, i.e., finding a uniquely determined geodesic path from a given point in any direction. Since the shortest paths do not provide the unique solution on triangle meshes, we solve the initial value problem in an indirect manner: given a fixed point and an initial tangent direction on a triangle mesh M, we first compute a geodesic curve γ ̂ on a piecewise smooth surface M ̂, which well approximates the input mesh M and can be constructed at little cost. Then, we solve a first-order ODE of the tangent vector using the fourth-order Runge–Kutta method, and parallel transport it along γ ̂ . When the geodesic curve reaches the boundary of the current patch, its tangent can be directly transported to the neighboring patch, thanks to the G 1 -continuity along the common boundary of two adjacent patches. Finally, once the geodesic curve γ ̂ is available, we project it onto the underlying mesh M, producing the discrete geodesic path γ, which is guaranteed to be unique on M . It is worth noting that our method is different from the conventional methods of directly solving the geodesic equation (i.e., a second-order ODE of the position) on piecewise smooth surfaces, which are difficult to implement due to the complicated representation of the geodesic equation involving Christoffel symbols. The proposed method, based on the first-order ODE of the tangent vector, is intuitive and easy for implementation. Our method is particularly useful for computing geodesic paths on low-resolution meshes which may have large and/or skinny triangles, since the conventional straightest geodesic paths are usually far from the ground truth. Highlights: Shortest geodesic is not able to solve the initial value problem of discrete geodesic. Geodesic equation are second-order ODEs. We solve the initial value problem on triangle meshes by solving a first-order ODE The computed discrete geodesic path converges to the one on the smooth surface.