基于差分进化算法的B 样条曲线曲面拟合

应用B 样条曲线曲面拟合内在形状带有间断或者尖点的数据时,最小二乘法得到的 拟合结果往往在间断和尖点处误差较大,原因在于最小二乘法将拟合函数B 样条的节点固定。本 文在利用3 次B 样条曲线和曲面拟合数据时,应用差分进化算法设计出一种能够自适应地设置B 样条节点的方法,同时对节点的数量和位置进行优化,使得B 样条拟合曲线曲面在间断和尖点处 产生拟多重节点,实现高精度地拟合采样于带有间断或尖点的曲线和曲面数据。     

Elitist clonal selection algorithm for optimal choice of free knots in B-spline data fitting

Data fitting with B-splines is a challenging problem in reverse engineering for CAD/CAM, virtual reality, data visualization, and many other fields. It is well-known that the fitting improves greatly if knots are considered as free variables. This leads, however, to a very difficult multimodal and multivariate continuous nonlinear optimization problem, the so-called knot adjustment problem. In this context, the present paper introduces an adapted elitist clonal selection algorithm for automatic knot adjustment of B-spline curves. Given a set of noisy data points, our method determines the number and location of knots automatically in order to obtain an extremely accurate fitting of data. In addition, our method minimizes the number of parameters required for this task. Our approach performs very well and in a fully automatic way even for the cases of underlying functions requiring identical multiple knots, such as functions with discontinuities and cusps. To evaluate its performance, it has been applied to three challenging test functions, and results have been compared with those from other alternative methods based on AIS and genetic algorithms. Our experimental results show that our proposal outperforms previous approaches in terms of accuracy and flexibility. Some other issues such as the parameter tuning, the complexity of the algorithm, and the CPU runtime are also discussed.     

Knot calculation for spline fitting via sparse optimization

Curve fitting with splines is a fundamental problem in computer-aided design and engineering. However, how to choose the number of knots and how to place the knots in spline fitting remain a difficult issue. This paper presents a framework for computing knots (including the number and positions) in curve fitting based on a sparse optimization model. The framework consists of two steps: first, from a dense initial knot vector, a set of active knots is selected at which certain order derivative of the spline is discontinuous by solving a sparse optimization problem; second, we further remove redundant knots and adjust the positions of active knots to obtain the final knot vector. Our experiments show that the approximation spline curve obtained by our approach has less number of knots compared to existing methods. Particularly, when the data points are sampled dense enough from a spline, our algorithm can recover the ground truth knot vector and reproduce the spline.     

基于遗传算法的B样条曲线和Bézier曲线的最小二乘拟合

考虑用B样条曲线拟合平面有序数据使得最小二乘拟合误差最小.一般有两种考虑,一种是保持B样条基函数的节点不变,选择参数使得拟合较优.参数的选择方法包括均匀取值、累加弦长法、centripetal model、Gauss-Newton迭代法等.另一种则是先确定好参数值(一般用累加弦长法),然后再用.某一算法计算出节点,使得拟合较优.同时把两者统一考虑,用遗传算法同时求出参数、节点使得拟合在最小二乘误差意义下最优.与Gauss-Newton迭代法、Piegl算法相比,本方法具有较好的鲁棒性(拟合曲线与初始值无关)、较高的精度及控制顶点少等优点.实验结果说明采用遗传算法得到的曲线逼近效果更好.用遗传算法对Bezier曲线拟合平面有序数据也进行了研究.     

Automatic knot adjustment using an artificial immune system for B-spline curve approximation

Reverse engineering transforms real parts into engineering concepts or models. First, sampled points are mapped from the object’s surface by using tools such as laser scanners or cameras. Then, the sampled points are fitted to a free-form surface or a standard shape by using one of the geometric modeling techniques. The curves on the surface have to be modeled before surface modeling. In order to obtain a good B-spline curve model from large data, the knots are usually respected as variables. A curve is then modeled as a continuous, nonlinear and multivariate optimization problem with many local optima. For this reason it is very difficult to reach a global optimum. In this paper, we convert the original problem into a discrete combinatorial optimization problem like in Yoshimoto et al. [F. Yoshimoto, M. Moriyama, T. Harada, Automatic knot placement by a genetic algorithm for data fitting with a spline, in: Proceedings of the International Conference on Shape Modeling and Applications, IEEE Computer Society Press, 1999, pp. 162-169] and Sarfraz et al. [M. Sarfraz, S.A. Raza, Capturing outline of fonts using genetic algorithm and splines, in: Fifth International Conference on Information Visualisation (IV’01), 2001, pp. 738-743]. Then, we suggest a new method that solves the converted problem by artificial immune systems. We think the candidates of the locations of knots as antibodies. We define the affinity measure benefit from Akaike’s Information Criterion (AIC). The proposed method determines the appropriate location of knots automatically and simultaneously. Furthermore, we do not need any subjective parameter or good population of initial location of knots for a good iterative search. Some examples are also given to demonstrate the efficiency and effectiveness of our method.     

Algorithms for smoothing data on the sphere with tensor product splines

Algorithms are presented for fitting data on the sphere by using tensor product splines which satisfy certain boundary constraints. First we consider the least-squares problem when the knots are given. Then we discuss the construction of smoothing splines on the sphere. Here the knots are located automatically. A Fortran IV implementation of these two algorithms is described.     

A study of equivalent electrical circuit fitting to electrochemical impedance using a stochastic method

Modeling electrochemical impedance spectroscopy is usually done using equivalent electrical circuits. These circuits have parameters that need to be estimated properly in order to make possible the simulation of impedance data. Despite the fitting procedure is an optimization problem solved recurrently in the literature, rarely statistical significance of the estimated parameters is evaluated. In this work, the optimization process for the equivalent electrical circuit fitting to the impedance data is detailed. First, a mathematical development regarding the minimization of residual least squares is presented in order to obtain a statistically valid objective function of the complex nonlinear regression problem. Then, the optimization method used in this work is presented, the Differential Evolution, a global search stochastic method. Furthermore, it is shown how a population-based stochastic method like this can be used directly to obtain confidence regions to the estimated parameters. A sensitivity analysis was also conducted. Finally, the equivalent circuit fitting is done to model synthetic experimental data, in order to demonstrate the adopted procedure.     

Fitting optimal piecewise linear functions using genetic algorithms

Constructing a model for data in R² is a common problem in many scientific fields, including pattern recognition, computer vision, and applied mathematics. Often little is known about the process which generated the data or its statistical properties. For example, in fitting a piecewise linear model, the number of pieces, as well as the knot locations, may be unknown. Hence, the method used to build the statistical model should have few assumptions, yet, still provide a model that is optimal in some sense. Such methods can be designed through the use of genetic algorithms. We examine the use of genetic algorithms to fit piecewise linear functions to data in R². The number of pieces, the location of the knots, and the underlying distribution of the data are assumed to be unknown. We discuss existing methods which attempt to solve this problem and introduce a new method which employs genetic algorithms to optimize the number and location of the pieces. Experimental results are presented which demonstrate the performance of our method and compare it to the performance of several existing methods, We conclude that our method represents a valuable tool for fitting both robust and nonrobust piecewise linear functions     

Algorithm/algorithmus 42 an algorithm for cubic spline fitting with convexity constraints

In this paper an algorithm is presented for fitting a cubic spline satisfying certain local concavity and convexity constraints, to a given set of data points. When using theL ₂ norm, this problem results in a quadratic programming problem which is solved by means of the Theil-Van de Panne procedure. The algorithm makes use of the well-conditioned B-splines to represent the cubic splines. The knots are located automatically, as a function of a given upper limit for the sum of squared residuals. A Fortran IV implementation is given.     

基于相邻体素选择的脑白质纤维交叉分叉问题解决方法

为了解决脑白质纤维交叉分叉问题,在传统算法的启发下,提出一种基于相邻体素选择的盘状张量分解算法。首先,选择合适的起点进行非分叉纤维的追踪,建立拟合函数数据集,得到拟合函数;其次,在纤维追踪出现交叉分叉问题时,建立该交叉分叉点及周围区域体素所对应的棋盘图;然后,计算以交叉点为中心的相邻张量的夹角,结合夹角的大小并利用得到的拟合函数进行纤维整体走行方向的估计,实现盘状张量的分解。算法既保证了局部信息的合理适用,又考虑了整体信息的影响,能够更加精确完整地跟踪纤维路径,解决纤维交叉分叉问题。与传统方法相比,本文算法可以更有效解决纤维分叉及交叉处的跟踪问题,从而使得到的纤维路径更加真实。     

连续等距区间上积分值的二次样条插值

目的 在现实中,某些插值问题结点处的函数值往往是未知的,而仅仅已知一些区间上的积分值。为此提出一种给定已知函数在连续等距区间上的积分值构造二次样条插值函数的方法。方法 首先,利用二次B样条基函数的线性组合去满足给定的积分值和两个端点插值条件,该插值问题等价于求解n+2个方程带宽为3的线性方程组。然后,运用算子理论给出二次样条插值函数的误差估计,继而得到二次样条函数逼近结点处的函数值时具有超收敛性。最后,通过等距区间上积分值的线性组合逼近两个端点的函数值方法实现了不带任何边界条件的积分型二次样条插值问题。结果 选取低频率函数,对积分型二次样条插值方法和改进方法分别进行数值测试,发现这两种方法逼近效果都是良好的。同样,选取高频率函数对积分型二次样条插值方法进行数值实验,得到数值收敛阶与理论值相一致。结论 实验结果表明,本文算法相比已有的方法更简单有效,对改进前后的二次样条插值函数在逼近结点处的函数值时的超收敛性得到了验证。该方法对连续等距区间上积分值的函数重构具有普适性。     

Surface Reconstruction Based on Hierarchical Floating Radial Basis Functions

In this paper we address the problem of optimal centre placement for scattered data approximation using radial basis functions (RBFs) by introducing the concept of floating centres. Given an initial least‐squares solution, we optimize the positions and the weights of the RBF centres by minimizing a non‐linear error function. By optimizing the centre positions, we obtain better approximations with a lower number of centres, which improves the numerical stability of the fitting procedure. We combine the non‐linear RBF fitting with a hierarchical domain decomposition technique. This provides a powerful tool for surface reconstruction from oriented point samples. By directly incorporating point normal vectors into the optimization process, we avoid the use of off‐surface points which results in less computational overhead and reduces undesired surface artefacts. We demonstrate that the proposed surface reconstruction technique is as robust as recent methods, which compute the indicator function of the solid described by the point samples. In contrast to indicator function based methods, our method computes a global distance field that can directly be used for shape registration.     

基于能量优化的数据参数化方法

构造一条通过一组给定数据点的参数样条曲线的关键之一是选取节点.本文提出了一种基于离散能量模型的确定节点参数的新方法,该方法首先通过极小化离散能量模型的目标函数确定参数样条曲线在节点处的最佳切向角,然后以最佳切向角为参变量来计算节点参数.该方法所得到的方程为线性方程,便于求解.本文最后通过实例对新方法与累加弦长方法、向心参数化方法、修正弦长参数化方法以及ZCM方法进行了比较.     

A fitting method with generalized Erlang distributions

We present a fitting technique that fits trace data into a generalized Erlang distribution class using an EM method. A generalized Erlang (GEr) distribution can be made by convolution of the third order ME distributions similar to the formulation of an Erlang distribution with exponential distributions. We give a sufficient condition for the representation to make a probability density function and we implement a fitting algorithm into a GEr distribution set by solving a nonlinear optimization problem with the EM algorithm. The effectiveness of the proposed fitting algorithm is presented by applying fitting methods to sets of synthetic data and measurement data. We present comparative numerical simulation results of our approach and other methods.     

带法向约束的隐式T样条曲线重构

目的 隐式曲线能够描述复杂的几何形状和拓扑结构,而传统的隐式B样条曲线的控制网格需要大量多余的控制点满足拓扑约束。有些情况下,获取的数据点不仅包含坐标信息,还包含相应的法向约束条件。针对这个问题,提出了一种带法向约束的隐式T样条曲线重建算法。方法 结合曲率自适应地调整采样点的疏密,利用二叉树及其细分过程从散乱数据点集构造2维T网格;基于隐式T样条函数提出了一种有效的曲线拟合模型。通过加入偏移数据点和光滑项消除额外零水平集,同时加入法向项减小曲线的法向误差,并依据最优化原理将问题转化为线性方程组求解得到控制系数,从而实现隐式曲线的重构。在误差较大的区域进行T网格局部细分,提高重建隐式曲线的精度。结果 实验在3个数据集上与两种方法进行比较,实验结果表明,本文算法的法向误差显著减小,法向平均误差由10^-3数量级缩小为10^-4数量级,法向最大误差由10^-2数量级缩小为10^-3数量级。在重构曲线质量上,消除了额外零水平集。与隐式B样条控制网格相比,3个数据集的T网格的控制点数量只有B样条网格的55.88%、39.80%和47.06%。结论 本文算法能在保证数据点精度的前提下,有效降低法向误差,消除了额外的零水平集。与隐式B样条曲线相比,本文方法减少了控制系数的数量,提高了运算速度。     

Approximation and estimation bounds for free knot splines

The fitting to data by splines has long been known to improve dramatically if the knots can be adjusted adaptively. To demonstrate the quality of the obtained free knot spline, it is essential to characterize its generalization ability. By utilizing the powerful techniques of the empirical process and approximation theory to address the estimation and approximation error bounds, respectively, the generalization ability of the free knot spline learning strategy is successfully characterized. We show that the Pseudo-dimension of free knot splines is essentially a linear function of the number of knots. A class of rather general loss functions is considered here and the squared loss is specially treated for its excellent property. We also provide some numerical results to demonstrate the utility of these theoretical results in guiding the process of choosing the appropriate knot numbers through the training data to avoid the overfitting/underfitting problem.     

Suboptimal Robot Joint Interpolation Within User-Specified Knot Tolerances

Approximation of a desired robot path can be accomplished by interpolating a curve through a sequence of joint-space knots. A smooth interpolated trajectory can be realized by using trigonometric splines. But, sometimes the joint trajectory is not required to exactly pass through the given knots. The knots may rather be centers of tolerances near which the trajectory is required to pass. In this article, we optimize trigonometric splines through a given set of knots subject to user-specified knot tolerances. The contribution of this article is the straightforward way in which intermediate constraints (i.e., knot angles) are incorporated into the parameter optimization problem. Another contribution is the exploitation of the decoupled nature of trigonometric splines to reduce the computational expense of the problem. The additional freedom of varying the knot angles results in a lower objective function and a higher computational expense compared to the case in which the knot angles are constrained to exact values. The specific objective functions considered are minimum jerk and minimum torque. In the minimum jerk case, the optimization problem reduces to a quadratic programming problem. Simulation results for a two-link manipulator are presented to support the results of this article.     

基于极小化二阶导矢确定节点

构造参数拟合曲线的关键问题之一是为每个数据点指定一个参数值（节点）。提出了一种确定节点的新方法。对于每个数据点,新方法由相邻的三个数据点构造一条二次多项式曲线,二次曲线的节点通过极小化其二阶导矢的平方确定。两个相邻数据点间的节点区间由两条二次曲线确定。为使节点计算公式能有效反映出相邻数据点的变化情况,新方法改进了修正弦长方法并应用于节点计算。新方法是一个局部化方法,因此适合于曲线曲面的交互设计。实验结果说明,新方法比其他节点计算方法有效。     

Surface reconstruction using bivariate simplex splines on Delaunay configurations

Recently, a new bivariate simplex spline scheme based on Delaunay configuration has been introduced into the geometric computing community, and it defines a complete spline space that retains many attractive theoretic and computational properties. In this paper, we develop a novel shape modeling framework to reconstruct a closed surface of arbitrary topology based on this new spline scheme. Our framework takes a triangulated set of points, and by solving a linear least-square problem and iteratively refining parameter domains with newly added knots, we can finally obtain a continuous spline surface satisfying the requirement of a user-specified error tolerance. Unlike existing surface reconstruction methods based on triangular B-splines (or DMS splines), in which auxiliary knots must be explicitly added in advance to form a knot sequence for construction of each basis function, our new algorithm completely avoids this less-intuitive and labor-intensive knot generating procedure. We demonstrate the efficacy and effectiveness of our algorithm on real-world, scattered datasets for shape representation and computing.     

Numerische Lösung gewöhnlicher Differentialgleichungen mit Splinefunktionen

In this paper a general procedure to obtain spline approximations for the solutions of initial value problems for ordinary differential equations is presented. Several well-known spline approximation methods are included as special cases. It is common practice to partition the interval for which the initial value problem is defined into equidistant subintervals and to construct successively the spline approximation; thereby the spline function has to satisfy certain conditions at the knots. In the general procedure presented here additional knots are admitted in every subinterval. At these points which need not be equally spaced the spline approximation has to fulfill analogous conditions as at the original knots. Convergence and divergence theorems are proved; especially the influence of the additional knots on convergence and divergence of the method is investigated.