Cylinders through Five Points: Computational Algebra and Geometry

We address the following question: Given five points in R 3 , determine a right circular cylinder containing those points. We obtain algebraic equations for the axial line and radius parameters and show that these give six solutions in the generic case. An even number (0, 2, 4, or 6) will be real valued and hence correspond to actual cylinders in R 3 . We will investigate computational and theoretical matters related to this problem. In particular we will show how exact and numeric Gr¨ obner bases, equation solving, and related symbolic-numeric methods may be used to advantage. We will also discuss some applications. obner bases, nonlinear systems, symbolic-numeric computation 1. Outline of the Problem and Related Work Given five points in R 3 , we are to determine all right circular cylinders containing those points. We do this by solving equations for the axial line and radius parameters. We will show that generically one obtains six solutions to these equations. Of these an even number are real valued, as the complex valued ones appear in conjugate pairs (an immediate consequence is that there is no "unique" real cylinder through five given points unless it a solution with multiplicity). Moreover there are open regions in the real configuration space that give each of these possibilities so none are disallowed. The basic problem of determining cylinders from five points may be recast in a computational geometry setting: Given five points in R 3 , find the smallest positive r and orientation parameters such that the cylinder of radius 2r with those parameters encloses tangentially the balls of radius r centered at the points. Here are some questions we will consider. The first three are classical; we address them here to illustrate the utility of symbolic computation in such investigations. The last ones are related to more recent work in computational and integral geometry.

[1]  Marek Teichmann,et al.  Smallest enclosing cylinders , 1996, SCG '96.

[2]  G. Alistair Watson,et al.  Fitting enclosing cylinders to data in Rn , 2006, Numerical Algorithms.

[3]  François Goulette,et al.  Extracting Cylinders in Full 3D Data Using a Random Sampling Method and the Gaussian Image , 2001, VMV.

[4]  François Goulette,et al.  A note on the construction of right circular cylinders through five 3D points , 2003 .

[5]  Michael Kalkbrenner,et al.  Implicitization of Rational Parametric Curves and Surfaces , 1990, AAECC.

[6]  Ilias S. Kotsireas Homotopies and polynomial system solving I: basic principles , 2001, SIGS.

[7]  Patrizia M. Gianni,et al.  Algebraic Solution of Systems of Polynomial Equations Using Groebner Bases , 1987, AAECC.

[8]  Robert M. Corless,et al.  Gröbner bases and matrix eigenproblems , 1996, SIGS.

[9]  Daniel Lichtblau,et al.  Solving finite algebraic systems using numeric Gröbner bases and eigenvalues , 2009 .

[10]  Bruno Buchberger,et al.  Applications of Gro¨bner bases in non-linear computational geometry , 1988 .

[11]  Bernd Sturmfels,et al.  Bernstein’s theorem in affine space , 1997, Discret. Comput. Geom..

[12]  Maria Grazia Marinari,et al.  The shape of the Shape Lemma , 1994, ISSAC '94.

[13]  Robert M Corless Editor's Corner: Gr?bner Bases and Matrix Eigenproblems , 1996 .

[14]  Michel Petitjean,et al.  About the algebraic solutions of smallest enclosing cylinders problems , 2010, Applicable Algebra in Engineering, Communication and Computing.

[15]  R. Janssen,et al.  Trends in Computer Algebra , 1988, Lecture Notes in Computer Science.

[16]  Franco P. Preparata,et al.  Evaluating the cylindricity of a nominally cylindrical point set , 2000, SODA '00.

[17]  M. Levine,et al.  Extracting geometric primitives , 1993 .

[18]  W. W. Adams,et al.  An Introduction to Gröbner Bases , 2012 .

[19]  J. Verschelde,et al.  Homotopies exploiting Newton polytopes for solving sparse polynomial systems , 1994 .

[20]  B. Buchberger Gröbner Bases and Applications: Introduction to Gröbner Bases , 1998 .

[21]  Bruno Buchberger,et al.  Applications of Gröbner Bases in Non-linear Computational Geometry , 1987, Trends in Computer Algebra.

[22]  Donal O'Shea,et al.  Ideals, varieties, and algorithms - an introduction to computational algebraic geometry and commutative algebra (2. ed.) , 1997, Undergraduate texts in mathematics.

[23]  Heinz Kredel,et al.  Gröbner Bases: A Computational Approach to Commutative Algebra , 1993 .

[24]  C. Hoffmann,et al.  Symbolic and numerical techniques for constraint solving , 1998 .

[25]  Marc Moreno Maza,et al.  On Solving Parametric Polynomial Systems , 2012, Mathematics in Computer Science.

[26]  John Fitch,et al.  Symbolic Computation and the Finite Element Method , 1989 .

[27]  Ralf Fröberg,et al.  An introduction to Gröbner bases , 1997, Pure and applied mathematics.

[28]  Daniel Lichtblau,et al.  Cylinders Through Five Points: Complex and Real Enumerative Geometry , 2006, Automated Deduction in Geometry.

[29]  János Pach,et al.  Common Tangents to Four Unit Balls in R3 , 2001, Discret. Comput. Geom..

[30]  N. S. Barnett,et al.  Private communication , 1969 .

[31]  Christoph M. Hoffmann,et al.  Geometric and Solid Modeling: An Introduction , 1989 .

[32]  N. Bose Gröbner Bases: An Algorithmic Method in Polynomial Ideal Theory , 1995 .

[33]  David A. Cox,et al.  Ideals, Varieties, and Algorithms: An Introduction to Computational Algebraic Geometry and Commutative Algebra, 3/e (Undergraduate Texts in Mathematics) , 2007 .

[34]  Patrizia M. Gianni,et al.  Properties of Gröbner bases under specializations , 1987, EUROCAL.

[35]  Hans J. Stetter,et al.  Stabilization of polynomial systems solving with Groebner bases , 1997, ISSAC.

[36]  Ralf Westphal,et al.  Pose Estimation of Cylindrical Fragments for Semi-automatic Bone Fracture Reduction , 2003, DAGM-Symposium.

[37]  Alessandro Zinani,et al.  The Expected Volume of a Tetrahedron whose Vertices are Chosen at Random in the Interior of a Cube , 2003 .

[38]  Frank Sottile,et al.  An Excursion From Enumerative Geometry to Solving Systems of Polynomial Equations with Macaulay 2 , 2000, math/0007142.

[39]  Stephen Wolfram,et al.  The Mathematica book, 5th Edition , 2003 .

[40]  O. Bottema,et al.  On the lines in space with equal distances to n given points , 1977 .

[41]  Michael Kalkbrener Solving systems of algebraic equations by using Gröbner bases , 1987, EUROCAL.

[42]  Franco P. Preparata,et al.  On circular Cylinders by Four or Five Points in Space , 2001 .

[43]  L. O'carroll AN INTRODUCTION TO GRÖBNER BASES (Graduate Studies in Mathematics 3) , 1996 .

[44]  J. Yorke,et al.  The cheater's homotopy: an efficient procedure for solving systems of polynomial equations , 1989 .

[45]  Thorsten Theobald,et al.  Algebraic Methods for Computing Smallest Enclosing and Circumscribing Cylinders of Simplices , 2002, Applicable Algebra in Engineering, Communication and Computing.

[46]  Kenneth Falconer,et al.  Unsolved Problems In Geometry , 1991 .

[47]  Bo Yuan,et al.  On Spatial Constraint Solving Approaches , 2000, Automated Deduction in Geometry.

[48]  Shojiro Sakata Applied Algebra, Algebraic Algorithms and Error-Correcting Codes, 9th International Symposium, AAECC-9, New Orleans, LA, USA, October 7-11, 1991, Proceedings , 1991, AAECC.