Bounding sphere
![](../I/m/Smallest_circle_problem.svg.png)
In mathematics, given a non-empty set of objects of finite extension in -dimensional space, for example a set of points, a bounding sphere, enclosing sphere or enclosing ball for that set is an -dimensional solid sphere containing all of these objects.
In the plane the terms bounding or enclosing circle are used.
Used in computer graphics and computational geometry, a bounding sphere is a special type of bounding volume. There are several fast and simple bounding sphere construction algorithms with a high practical value in real-time computer graphics applications.[1]
In statistics and operations research, the objects are typically points, and generally the sphere of interest is the minimal bounding sphere, that is, the sphere with minimal radius among all bounding spheres. It may be proven that such a sphere is unique: If there are two of them, then the objects in question lie within their intersection. But an intersection of two non-coinciding spheres of equal radius is contained in a sphere of smaller radius.
The problem of computing the center of a minimal bounding sphere is also known as the "unweighted Euclidean 1-center problem".
Applications
Clustering
Such spheres are useful in clustering, where groups of similar data points are classified together.
In statistical analysis the scattering of data points within a sphere may be attributed to measurement error or natural (usually thermal) processes, in which case the cluster represents a perturbation of an ideal point. In some circumstances this ideal point may be used as a substitute for the points in the cluster, advantageous in reducing calculation time.
In operations research the clustering of values to an ideal point may also be used to reduce the number of inputs in order to obtain approximate values for NP-hard problems in a reasonable time. The point chosen is not usually the center of the sphere, as this can be biased by outliers, but instead some form of average location such as a least squares point is computed to represent the cluster.
Algorithms
There are exact and approximate algorithms for solving bounding sphere problem.
Ritter's bounding sphere
In 1990, Jack Ritter proposed a simple algorithm to find a non-minimal bounding sphere.[2] It is widely used in various applications for its simplicity. The algorithm works in this way:
- Pick a point from , search a point in , which has the largest distance from ;
- Search a point in , which has the largest distance from . Set up an initial ball , with its centre as the midpoint of and , the radius as half of the distance between and ;
- If all points in are within ball , then we get a bounding sphere. Otherwise, let be a point outside the ball, construct a new ball covering both point and previous ball. Repeat this step until all points are covered.
Ritter's algorithm runs in time on inputs consisting of points in -dimensional space, which makes it very efficient. However it gives only a coarse result which is usually 5% to 20% larger than the optimum.
Linear programming
In 1983, Nimrod Megiddo proposed a "prune and search" algorithm which finds the optimum bounding sphere and runs in linear time if the dimension is fixed as a constant. When dimension is taken into account, the execution time complexity is , impractical for high-dimensional applications. Megiddo used this approach to solve linear programming in linear time when dimension is fixed.
In 1991, Emo Welzl proposed a much simpler randomized algorithm based in the extension of a randomized linear programming algorithm by Raimund Seidel. It runs in expected linear time and provided experimental results demonstrating its practicality in higher dimensions.[3]
The open-source Computational Geometry Algorithms Library (CGAL) contains an implementation of this algorithm.[4]
Core-set based approximation
Bădoiu et al. presented a approximation to the bounding sphere problem,[5] where a approximation means that the constructed sphere has radius at most , where is the smallest possible radius of a bounding sphere.
A coreset is a small subset, that a expansion of the solution on the subset is a bounding sphere of the whole set. The coreset is constructed incrementally by adding the farthest point into the set in each iteration.
Kumar et al. improved this approximation algorithm[6] so that it runs in time .
Fischer's exact solver
Fischer et al. (2003) proposed an exact solver to this minimal enclosing ball problem, though the algorithm does not have a polynomial running time in the worst case.[7] The algorithm is purely combinatorial and implements a pivoting scheme similar to the simplex method for linear programming, used earlier in some heuristics. It starts with a big ball that covers all points and gradually shrink it until it can not be shrunk further. Its major contributions are correct termination rules in cases of degeneracies, overlooked by other authors, and efficient handling of partial solutions, which produced a major speed-up. The authors verified that the algorithm is efficient in practice in low and moderately low (up to 10,000) dimensions and claim it does not exhibit numerical stability problems in its floating-point operations.[7] A C++ implementation of the algorithm is available as an open-source project.[8]
Extremal points optimal sphere
Larsson (2008) proposed the "extremal points optimal sphere" method with controllable speed to accuracy approximation to solve the bounding sphere problem. This method works by taking a suitable subset of the points, apply an exact solver over them, iterates over the remaining points if any, and grows the sphere if necessary. This method is orders of magnitude faster than exact methods while giving comparable results, with a worse case time of , where is an accuracy variable representing the number of direction considered in the algorithm.[1]
See also
References
- 1 2 Larsson, Thomas (2008), "Fast and tight fitting bounding spheres", SIGRAD 2008: The Annual SIGRAD Conference, Special Theme: Interaction, November 27-28, 2008, Stockholm, Sweden, Linköping Electronic Conference Proceedings, 34, Linköping, Sweden: Linköping University
- ↑ Ritter, Jack (1990), "An efficient bounding sphere", in Glassner, Andrew S., Graphics Gems, San Diego, CA, US: Academic Press Professional, Inc., pp. 301–303, ISBN 0-12-286166-3
- ↑ Welzl, Emo (1991), "Smallest enclosing disks (balls and ellipsoids)", in Maurer, Hermann, New Results and New Trends in Computer Science: Graz, Austria, June 20–21, 1991, Proceedings, Lecture Notes in Computer Science, 555, Berlin, Germany: Springer, pp. 359–370, doi:10.1007/BFb0038202, MR 1254721
- ↑ CGAL 4.3 - Bounding Volumes - Min_sphere_of_spheres_d, retrieved 2014-03-30.
- ↑ Bādoiu, Mihai; Har-Peled, Sariel; Indyk, Piotr (2002), "Approximate clustering via core-sets" (PDF), Proceedings of the Thirty-Fourth Annual ACM Symposium on Theory of Computing, New York, NY, US: ACM, pp. 250–257, doi:10.1145/509907.509947, MR 2121149
- ↑ Kumar, Piyush; Mitchell, Joseph S. B.; Yıldırım, E. Alper (2003), "Computing core-sets and approximate smallest enclosing hyperspheres in high dimensions", in Ladner, Richard E., Proceedings of the Fifth Workshop on Algorithm Engineering and Experiments, Baltimore, MD, USA, January 11, 2003, Philadelphia, PA, US: SIAM, pp. 45–55
- 1 2 Fischer, Kaspar; Gärtner, Bernd; Kutz, Martin (2003), "Fast smallest-enclosing-ball computation in high dimensions" (PDF), in Battista, Giuseppe Di; Zwick, Uri, Algorithms: ESA 2003, 11th Annual European Symposium, Budapest, Hungary, September 16-19, 2003, Proceedings, Lecture Notes in Computer Science, 2832, Springer, Berlin, pp. 630–641, doi:10.1007/978-3-540-39658-1_57
- ↑ miniball open-source project
External links
- Smallest Enclosing Circle Problem – describes several algorithms for enclosing a point set, including Megiddo's linear-time algorithm