二次均匀B样条曲线的双圆弧逼近方法*

提出了一种用双圆弧对二次均匀B样条曲线的分段逼近方法。首先,对一条具有n 1个控制顶点的二次均匀B样条曲线按照相邻两节点界定的区间分成n-1段只有三个控制顶点的二次均匀B样条曲线段;然后对每一曲线段构造一条双圆弧进行逼近。所构造的双圆弧满足端点及端点切向量条件,即双圆弧的两个端点分别是所逼近的曲线段的端点,而且双圆弧在两个端点处的切向量是所逼近的曲线段在端点处的单位切向量。同时,双圆弧的连接点是双圆弧连接点轨迹圆与其所逼近的曲线段的交点。这些新构造出来的双圆弧连接在一起构成了一条圆弧样条曲线,即二次均匀B样条曲线的逼近曲线。另外给出了逼近误差分析和实例说明。     

Efficient circular arc interpolation based on active tolerance control

In this paper, we present an efficient sub-optimal algorithm for fitting smooth planar parametric curves by G¹ arc splines. To fit a parametric curve by an arc spline within a prescribed tolerance, we first sample a set of points and tangents on the curve adaptively as well as with enough density, so that an interpolation biarc spline curve can be with any desired high accuracy. Then, we construct new biarc curves interpolating local triarc spirals explicitly based on the control of permitted tolerances. To reduce the segment number of fitting arc spline as much as possible, we replace the corresponding parts of the spline by the new biarc curves and compute active tolerances for new interpolation steps. By applying the local biarc curve interpolation procedure recursively and sequentially, the result circular arcs with no radius extreme are minimax-like approximation to the original curve while the arcs with radius extreme approximate the curve parts with curvature extreme well too, and we obtain a near optimal fitting arc spline in the end. Even more, the fitting arc spline has the same end points and end tangents with the original curve, and the arcs will be jointed smoothly if the original curve is composed of several smooth connected pieces. The algorithm is easy to be implemented and generally applicable to circular arc interpolation problem of all kinds of smooth parametric curves. The method can be used in wide fields such as geometric modeling, tool path generation for NC machining and robot path planning, etc. Several numerical examples are given to show the effectiveness and efficiency of the method.     

Biarc approximation of polygons within asymmetric tolerance bands

We present an algorithm for approximating a simple planar polygon by a tangent-continuous approximation curve that consists of biarcs. Our algorithm guarantees that the approximation curve lies within a user-specified tolerance from the original polygon. If requested, the algorithm can also guarantee that the original polygon lies within a user-specified distance from the approximation curve. Both symmetric and asymmetric tolerances can be handled. In either case, the approximation curve is guaranteed to be simple. Simplicity of the approximation curve is achieved by restricting it to a ‘tolerance band’ which represents the user-specified tolerance and which takes into account bottlenecks of the input polygon. The tolerance band itself is computed by means of a regular grid and so-called k-dops. The basic algorithm is readily extended to compute biarc approximations of collections of polygonal curves simultaneously. Experimental results demonstrate that this algorithm computes biarc approximations of an n-vertex polygon with a close-to-minimum number of biarcs in roughly time.     

单圆弧样条保形插值算法

该文以插值具有偶数个点的闭多边形为例提出了一种新的圆弧样条插值算法。这种算法具有以下3个特点：（1）生成的圆弧样条曲线具有保形的特点；（2）圆弧样条中圆弧的段数与型值点个数相同。（3）圆弧段之间的连接点不一定在插值的型值点上，这样就能用更多的自由度来控制拟合曲线的形状。同此文中还提出了一个优化的算法来得到光顺的插值曲线，同时还给出了几个例子加以说明。     

A practicable approach to G biarc approximations for making accurate, smooth and non-gouged profile features in CNC contouring

Extensive research on G¹ biarc approximations to free-form curves has been conducted for the production of accurate, smooth and non-gouged profile features in CNC contouring. However, all the published work has only focused on improving the fitting accuracy between the biarc curve and the nominal free-form curve of a profile and minimizing the biarc number; as a result, the radii of the concave arcs of some biarcs could be less than the pre-determined tool radius, and the tool would overcut these arcs in machining, eventually gouging the profile. In this work, a new, practicable approach is proposed to completely solve this problem. The main feature of this approach is to find the gouging-free parameter interval of a biarc family, among which the radii of all the concave arcs are larger than the tool radius, and then to search in this interval for a best fitting biarc so that its approximation accuracy is within the tolerance. This approach is robust and easy to implement and can substantially promote the use of G¹ biarc curves for CNC machining.     

平面列表点曲线的最优双圆弧拟合

本文利用相切的双圆弧拟合平面列表点曲线，最优原则取“应变能”与弧长加权之和为最小。这个方法克服了单圆弧样条拟合及其它优化原则方法的缺点，在数控加工中将得到很好的应用。     

An approximation of circular arcs by quartic Bézier curves

In this paper we propose an approximation method for circular arcs by quartic Bézier curves. Using an alternative error function, we give the closed form of the Hausdorff distance between the circular arc and the quartic Bézier curve. We also show that the approximation order is eight. By subdivision of circular arcs with equi-length, our method yields the curvature continuous spline approximation of the circular arc. We confirm that the approximation proposed in this paper has a smaller error than previous quartic Bézier approximations.     

Approximation by fat arcs and fat biarcs

A general discussion of the including approximation of a curve by a fat arc is given followed by an algorithm for constructing an including fat arc for a parametric Bézier curve. An example of applying the algorithm is given. The results for a fat arc are then used to develop an including approximation for a curve segment using a fat biarc. An algorithm for a fat biarc including approximation is provided followed by examples of Bézier curves being included by a fat biarc.     

A Recursive Subdivision Algorithm for Piecewise Circular Spline

We present an algorithm for generating a piecewise G ¹ circular spline curve from an arbitrary given control polygon. For every corner, a circular biarc is generated with each piece being parameterized by its arc length. This is the first subdivision scheme that produces a piecewise biarc curve that can interpolate an arbitrary set of points. It is easily adopted in a recursive subdivision surface scheme to generate surfaces with circular boundaries with pieces parameterized by arc length, a property not previously available. As an application, a modified version of Doo–Sabin subdivision algorithm is outlined making it possible to blend a subdivision surface with other surfaces having circular boundaries such as cylinders.     

An optimization approach for biarc curve-fitting of B-spline curves

We present an approach to the optimal fitting of a biarc-spline to a given B-spline curve. The objective is to minimize the area between the original B-spline curve and the fitted curve. Such an objective has obvious practical implications. This approach differs from conventional biarc curve-fitting techniques in two main aspects and has some desirable features. Firstly, it exploits the inherent freedom in the choice of the biarc that can be fitted to a given pair of end-points and their tangents. The conventional approach to biarc curve-fitting introduces additional constraints, such as the minimal difference in curvature or others to uniquely determine successive biarcs. In this approach, such constraints are not imposed. Instead, the freedom is exploited in the problem formulation to achieve a better fit. Secondly, the end-points do not lie on the curve so that appropriate tolerance control can be imposed through the use of additional constraints. Almost all previous biarc-fitting methods consider end-points that are on the original curve. As a result of these two aspects, the resulting biarc curve fits closely to the original curve with relatively fewer segments. This has a desirable effect on the surface finish, verification of CNC codes and memory requirement. Numerical results of the application of this approach to several examples are presented.     

Data Approximation Using Biarcs

. An algorithm for data approximation with biarcs is presented. The method uses a specific formulation of biarcs appropriate for parametric curves in Bézier or NURBS formulation. A base curve is applied to obtain tangents and anchor points for the individual arcs joining in G 1 continuity. Data sampled from circular arcs or straight line segments is represented precisely by one biarc. The method is most useful in numerical control to drive the cutter along straight line or circular paths.     

Curvature continuous offset approximation based on circle approximation using quadratic Bézier biarcs

We present a method for G² end-point interpolation of offset curves using rational Bézier curves. The method is based on a G² end-point interpolation of circular arcs using quadratic Bézier biarcs. We also prove the invariance of the Hausdorff distance between two compatible curves under convolution. Using this result, we obtain the exact Hausdorff distance between an offset curve and its approximation by our method. We present the approximation algorithm and give numerical examples.     

Polyline approach for approximating Hausdorff distance between planar free-form curves

This paper presents a practical polyline approach for approximating the Hausdorff distance between planar free-form curves. After the input curves are approximated with polylines using the recursively splitting method, the precise Hausdorff distance between polylines is computed as the approximation of the Hausdorff distance between free-form curves, and the error of the approximation is controllable. The computation of the Hausdorff distance between polylines is based on an incremental algorithm that computes the directed Hausdorff distance from a line segment to a polyline. Furthermore, not every segment on polylines contributes to the final Hausdorff distance. Based on the bound properties of the Hausdorff distance and the continuity of polylines, two pruning strategies are applied in order to prune useless segments. The R-Tree structure is employed as well to accelerate the pruning process. We experimented on Bezier curves, B-Spline curves and NURBS curves respectively with our algorithm, and there are 95% segments pruned on approximating polylines in average. Two comparisons are also presented: One is with an algorithm computing the directed Hausdorff distance on polylines by building Voronoi diagram of segments. The other comparison is with equation solving and pruning methods for computing the Hausdorff distance between free-form curves.     

平面NURBS曲线及其Offset的双圆弧逼近

除直线、圆弧、速端曲线等少数几种曲线外,平面参数曲线的offset曲线通常不能表示成有 理参数形式,因此在实际应用中,为了方便造型系统中数据结构和几何算法的统一表示,offse t曲线通常用低次曲线逼近来表示.通过用双圆弧逼近表示NURBS(non-uniform rational B -spline)曲线及其offset,并利用双圆弧逼近的特有性质,把offset的双圆弧逼近转化为原 曲线的双圆弧逼近,简化了问题的求解.同时考虑了双圆弧逼近算法中分割点的选取、公切点 的确定以及误差估计等主要问题.具体算     

Biarc approximation of NURBS curves

An algorithm for approximating arbitrary NURBS curves with biarcs is presented. The main idea is to approximate the NURBS curve with a polygon, and then to approximate the polygon with biarcs to within the required tolerance. The method uses a parametric formulation of biarcs appropriate in geometric design using parametric curves. The method is most useful in numerical control to drive the cutter along straight line or circular paths.     

Total curvature variation fairing for medial axis regularization

We present a new fairing method for planar curves, which is particularly well suited for the regularization of the medial axis of a planar domain. It is based on the concept of total variation regularization. The original boundary (given as a closed B-spline curve or several such curves for multiply connected domains) is approximated by another curve that possesses a smaller number of curvature extrema. Consequently, the modified curve leads to a smaller number of branches of the medial axis. In order to compute the medial axis, we use the state-of-the-art algorithm from [1] which is based on arc spline approximation and a domain decomposition approach. We improve this algorithm by using a different decomposition strategy that allows to reduce the number of base cases from 13 to only 5. Moreover, the algorithm reduces the number of conic arcs in the output by approx. 50%.     

基于最佳圆弧样条逼近的快速距离曲面计算

距离曲面是一种常用的隐式曲面,它在几何造型和计算机动画中具有重要的应用价值,但以往往在对距离曲面进行多边形化时速较慢,为了提高点到曲线最近距离计算的效率,提出了一种基于最佳圆弧样条逼近的快速线骨架距离曲面计算方法,该算法对于一条任意的二维NURBS曲线,在用户给定的误差范围内,先用最少量的圆弧样条来逼近给定的曲线,从而把点到NURBS曲线最近距离的计算问题转化为点到圆弧样条最近距离的计算问题,由于在对曲面进行多边形化时,需要大量的点到曲线最近距离的计算,而该处可以将点到圆弧样条最近距离很少的计算量来解析求得,故该算法效率很高,该实验表明,算法简单实用,具有很大的应用价值。     

Precise Hausdorff distance computation for planar freeform curves using biarcs and depth buffer

We present a real-time algorithm for computing the precise Hausdorff Distance (HD) between two planar freeform curves. The algorithm is based on an effective technique that approximates each curve with a sequence of G ¹ biarcs within an arbitrary error bound. The distance map for the union of arcs is then given as the lower envelope of trimmed truncated circular cones, which can be rendered efficiently to the graphics hardware depth buffer. By sampling the distance map along the other curve, we can estimate a lower bound for the HD and eliminate many redundant curve segments using the lower bound. For the remaining curve segments, we read the distance map and detect the pixel(s) with the maximum distance. Checking a small neighborhood of the maximum-distance pixel, we can reduce the computation to considerably smaller subproblems, where we employ a multivariate equation solver for an accurate solution to the original problem. We demonstrate the effectiveness of the proposed approach using several experimental results.     

Computing a compact spline representation of the medial axis transform of a 2D shape

We present a full pipeline for computing the medial axis transform of an arbitrary 2D shape. The instability of the medial axis transform is overcome by a pruning algorithm guided by a user-defined Hausdorff distance threshold. The stable medial axis transform is then approximated by spline curves in 3D to produce a smooth and compact representation. These spline curves are computed by minimizing the approximation error between the input shape and the shape represented by the medial axis transform. Our results on various 2D shapes suggest that our method is practical and effective, and yields faithful and compact representations of medial axis transforms of 2D shapes.     

一种参数曲线间Hausdorff 距离的计算方法

针对一般的连续参数曲线,提出一种快速计算曲线间Hausdorff 距离的方法。由 于曲线的近似折线能很好的表示曲线,所以,许多软件中,采用曲线的近似折线绘制曲线。为 此,证明了在任意给定误差范围下,可以将曲线间的Hausdorff 距离转化为折线间的Hausdorff 距离,进一步转化为点到线段间的距离进行计算,并辅之必要的剪枝策略和增量式算法以提高 计算效率。该方法计算速度快,逼近度高,基本解决了参数曲线间Hausdorff 距离的计算问题, 在几何设计、图像匹配、图像识别等领域有广泛应用。