Lévy Flight-Driven Simulated Annealing for B-spline Curve Fitting

Point cloud approximation by spline models, also called curve and surface reconstruction, is an active research field in computer-aided design and manufacturing (CAD/CAM). Due to the physical and mechanical processes used to obtain the data, the measurements are often affected by noise and other distortions. Obtaining a suitable spline model to reconstruct the underlying shape of the data while maintaining a low design complexity leads to a multivariate and highly non-linear optimization problem, also known to be non-convex and multi-modal. In this work, we propose a method to fit a given point cloud by means of a B-spline curve model. Our approach to solve the optimization problem is based on a powerful thermodynamics-driven metaheuristic known as the Simulated Annealing. We compute the model parameters by combining traditional SA techniques with Levy flights (random walks based on the Levy distribution). The ability to perform such a flight allows the algorithm to escape from local minima and energy plateaus, a strong requirement when dealing with highly multi-modal problems. The performance and robustness of our algorithm is tested against three illustrative examples. Our experimental results show that our method is able to reconstruct the underlying shape of the data, even in the presence of noise, with acceptable accuracy and in a completely automated way.

[1]  Steven G. Louie,et al.  A Monte carlo simulated annealing approach to optimization over continuous variables , 1984 .

[2]  Andrés Iglesias,et al.  Particle swarm optimization for non-uniform rational B-spline surface reconstruction from clouds of 3D data points , 2012, Inf. Sci..

[3]  Andrés Iglesias,et al.  Cuckoo Search with Lévy Flights for Weighted Bayesian Energy Functional Optimization in Global-Support Curve Data Fitting , 2014, TheScientificWorldJournal.

[4]  M. Locatelli Simulated Annealing Algorithms for Continuous Global Optimization , 2002 .

[5]  Toshinobu Harada,et al.  Automatic knot placement by a genetic algorithm for data fitting with a spline , 1999, Proceedings Shape Modeling International '99. International Conference on Shape Modeling and Applications.

[6]  Andrés Iglesias,et al.  Global-Support Rational Curve Method for Data Approximation with Bat Algorithm , 2015, AIAI.

[7]  Kathryn A. Dowsland,et al.  Simulated Annealing , 1989, Encyclopedia of GIS.

[8]  A. Iglesias,et al.  Firefly Algorithm for Explicit B-Spline Curve Fitting to Data Points , 2013 .

[9]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[10]  Two Simulated Annealing Optimization Schemas for Rational Bézier Curve Fitting in the Presence of Noise , 2016 .

[11]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[12]  Cecilia R. Aragon,et al.  Optimization by Simulated Annealing: An Experimental Evaluation; Part II, Graph Coloring and Number Partitioning , 1991, Oper. Res..

[13]  Andrés Iglesias,et al.  Elitist clonal selection algorithm for optimal choice of free knots in B-spline data fitting , 2015, Appl. Soft Comput..

[14]  Angel Cobo,et al.  Bézier Curve and Surface Fitting of 3D Point Clouds Through Genetic Algorithms, Functional Networks and Least-Squares Approximation , 2007, ICCSA.

[15]  H. Szu Fast simulated annealing , 1987 .

[16]  Andrés Iglesias,et al.  Simulated Annealing and Natural Neighbor for Rational Bézier Surface Reconstruction from Scattered Data Points , 2017, ICHSA.

[17]  Xin-She Yang,et al.  Flower Pollination Algorithm for Global Optimization , 2012, UCNC.

[18]  Hyungjun Park,et al.  An error-bounded approximate method for representing planar curves in B-splines , 2004, Comput. Aided Geom. Des..

[19]  Xin-She Yang,et al.  Cuckoo Search via Lévy flights , 2009, 2009 World Congress on Nature & Biologically Inspired Computing (NaBIC).

[20]  Andrés Iglesias,et al.  A new iterative mutually coupled hybrid GA-PSO approach for curve fitting in manufacturing , 2013, Appl. Soft Comput..

[21]  Chuyi Song,et al.  Parameter Optimization for Bezier Curve Fitting Based on Genetic Algorithm , 2013, ICSI.

[22]  Andrés Iglesias,et al.  Memetic firefly algorithm for data fitting with rational curves , 2015, 2015 IEEE Congress on Evolutionary Computation (CEC).

[23]  Erkan Ülker B-Spline Curve Approximation Using Pareto Envelope- Based Selection Algorithm - PESA , 2013 .

[24]  Hyungjun Park,et al.  B-spline curve fitting based on adaptive curve refinement using dominant points , 2007, Comput. Aided Des..

[25]  M. E. Johnson,et al.  Generalized simulated annealing for function optimization , 1986 .

[26]  R. Mantegna,et al.  Fast, accurate algorithm for numerical simulation of Lévy stable stochastic processes. , 1994, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  Andrés Iglesias,et al.  Firefly Algorithm for Polynomial Bézier Surface Parameterization , 2013, J. Appl. Math..

[28]  Andrés Iglesias,et al.  New memetic self-adaptive firefly algorithm for continuous optimisation , 2016 .

[29]  Andries Petrus Engelbrecht,et al.  Fundamentals of Computational Swarm Intelligence , 2005 .

[30]  Andrés Iglesias,et al.  Efficient particle swarm optimization approach for data fitting with free knot B-splines , 2011, Comput. Aided Des..

[31]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[32]  D. Jupp Approximation to Data by Splines with Free Knots , 1978 .

[33]  Andrés Iglesias,et al.  Simulated Annealing Algorithm for Bezier Curve Approximation , 2014, 2014 International Conference on Cyberworlds.

[34]  Ahmet Arslan,et al.  Automatic knot adjustment using an artificial immune system for B-spline curve approximation , 2009, Inf. Sci..

[35]  Andrés Iglesias,et al.  Iterative two-step genetic-algorithm-based method for efficient polynomial B-spline surface reconstruction , 2012, Inf. Sci..

[36]  Mounib Mekhilef,et al.  Optimization of a representation , 1993, Comput. Aided Des..

[37]  Cecilia R. Aragon,et al.  Optimization by Simulated Annealing: An Experimental Evaluation; Part I, Graph Partitioning , 1989, Oper. Res..

[38]  Xin-She Yang,et al.  Engineering Optimization: An Introduction with Metaheuristic Applications , 2010 .

[39]  Andrés Iglesias,et al.  Four Adaptive Memetic Bat Algorithm Schemes for Bézier Curve Parameterization , 2016, Trans. Comput. Sci..

[40]  Toshinobu Harada,et al.  Data fitting with a spline using a real-coded genetic algorithm , 2003, Comput. Aided Des..

[41]  Andrés Iglesias,et al.  Memetic Simulated Annealing for Data Approximation with Local-Support Curves , 2017, ICCS.

[42]  Andrés Iglesias,et al.  Hybridizing mesh adaptive search algorithm and artificial immune systems for discrete rational Bézier curve approximation , 2015, The Visual Computer.

[43]  M. Powell A Direct Search Optimization Method That Models the Objective and Constraint Functions by Linear Interpolation , 1994 .

[44]  Andrés Iglesias,et al.  Cuckoo search with Lévy flights for reconstruction of outline curves of computer fonts with rational Bézier curves , 2016, 2016 IEEE Congress on Evolutionary Computation (CEC).

[45]  Lester Ingber,et al.  Adaptive simulated annealing (ASA): Lessons learned , 2000, ArXiv.

[46]  G. Wiselin Jiji,et al.  An enhanced particle swarm optimization with levy flight for global optimization , 2016, Appl. Soft Comput..

[47]  Andrés Iglesias,et al.  Nature-Inspired Swarm Intelligence for Data Fitting in Reverse Engineering: Recent Advances and Future Trends , 2016, Nature-Inspired Computation in Engineering.

[48]  Ralph R. Martin,et al.  Reverse engineering of geometric models - an introduction , 1997, Comput. Aided Des..