Subdivision Schemes for Thin Plate Splines

Thin plate splines are a well known entity of geometric design. They are defined as the minimizer of a variational problem whose differential operators approximate a simple notion of bending energy. Therefore, thin plate splines approximate surfaces with minimal bending energy and they are widely considered as the standard "fair" surface model. Such surfaces are desired for many modeling and design applications.
Traditionally, the way to construct such surfaces is to solve the associated variational problem using finite elements or by using analytic solutions based on radial basis functions. This paper presents a novel approach for defining and computing thin plate splines using subdivision methods. We present two methods for the construction of thin plate splines based on subdivision: A globally supported subdivision scheme which exactly minimizes the energy functional as well as a family of strictly local subdivision schemes which only utilize a small, finite number of distinct subdivision rules and approximately solve the variational problem. A tradeoff between the accuracy of the approximation and the locality of the subdivision scheme is used to pick a particular member of this family of subdivision schemes.
Later, we show applications of these approximating subdivision schemes to scattered data interpolation and the design of fair surfaces. In particular we suggest an efficient methodology for finding control points for the local subdivision scheme that will lead to an interpolating limit surface and demonstrate how the schemes can be used for the effective and efficient design of fair surfaces.     

曲面三参数binary细分法

在经典四点细分法的基础上,通过在曲线细分过程中引入三个参数,给出一种改进的细分曲线构造的算法,利用生成多项式等方法对细分法的一致收敛性、Ck连续性进行了分析。并把该方法扩展到曲面上,进而提出了曲面三参数binary细分法。在给定初始控制数据的条件下,可以通过对形状参数的适当选择来实现对细分极限曲面形状的调控。数值实验表明该算法较容易控制曲面形状,可方便地应用于工程实际,解决曲线、曲面位置调整和控制问题。     

任意拓扑网的细分改进

在改进任意拓扑网构造光滑表面时,初始控制网格确定的情况下,生成的曲面形状惟一确定,最终的物体造型也随之确定,不具有可调性,因而在曲面细分过程中引入了控制参数和摄动。通过引入控制参数,调节一个参数值,使得所得的细分曲面的表达度可控,可以得到一系列的细分曲面。引入摄动是为了改进了空间位置,允许局部地调控约束曲面的形状。最后给出了曲面设计的实例,表明这种算法简单、有效。     

细分曲面拟合的局部渐进插值方法

逼近型细分方法生成的细分曲面其品质要优于插值型细分方法生成的细分曲面.然而,逼近型细分方法生成的细分曲面不能插值于初始控制网格顶点.为使逼近型细分曲面具有插值能力,一般通过求解全局线性方程组,使其插值于网格顶点.当网格顶点较多时,求解线性方程组的计算量很大,因此,难以处理稠密网格.与此不同,在不直接求解线性方程组的情况下,渐进插值方法通过迭代调整控制网格顶点,最终达到插值的效果.渐进插值方法可以处理稠密的任意拓扑网格,生成插值于初始网格顶点的光滑细分曲面.并且经证明,逼近型细分曲面渐进插值具有局部性质,也就是迭代调整初始网格的若干控制顶点,且保持剩余顶点不变,最终生成的极限细分曲面仍插值于初始网格中被调整的那些顶点.这种局部渐进插值性质给形状控制带来了更多的灵活性,并且使得自适应拟合成为可能.实验结果验证了局部渐进插值的形状控制以及自适应拟合能力.     

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.     

Doo-Sabin细分算法在动态模式下的推广

提出一种基于均匀三角多项式B样条的动态保凸细分算法，它可以看作Doo-Sabin细分算法在动态模式下的一个推广．其细分规则基于张量积曲面细分模式的几何意义，不仅可以生成旋转曲面等特殊曲面，而且可以根据参数来控制细分曲面的形状．最后运用传统的离散傅里叶技术和特征根方法证明了该细分算法的收敛性．     

Artifact analysis on triangular box-splines and subdivision surfaces defined by triangular polyhedra

Surface artifacts are features in a surface which cannot be avoided by movement of control points. They are present in B-splines, box splines and subdivision surfaces. We showed how the subdivision process can be used as a tool to analyse artifacts in surfaces defined by quadrilateral polyhedra ( [Sabin et al., 2005] and [Augsd?rfer et al., 2011]).In this paper we are utilising the subdivision process to develop a generic expression which can be employed to determine the magnitude of artifacts in surfaces defined by any regular triangular polyhedra. We demonstrate the method by analysing box-splines and regular regions of subdivision surfaces based on triangular meshes: Loop subdivision, Butterfly subdivision and a novel interpolating scheme with two smoothing stages. We compare our results for surfaces defined by triangular polyhedra to those for surfaces defined by quadrilateral polyhedra.     

基于四边形网格的可调细分曲面造型方法

提出了一种基于四边形网格的可调细分曲面造型方法。该方法不仅适合闭域拓扑结构，且对初始网格是开域的也能进行处理。细分算法中引入了可调参数，增加了曲面造型的灵活性。在给定初始数据的条件下，曲面造型时可以通过调节参数来控制极限曲面的形状。该方法可以生成C1连续的细分曲面。试验表明该方法生成光滑曲面是有效的。     

Integration of generalized B-spline functions on Catmull–Clark surfaces at singularities

Subdivision surfaces are a common tool in geometric modelling, especially in computer graphics and computer animation. Nowadays, this concept has become established in engineering too. The focus here is on quadrilateral control grids and generalized B-spline surfaces of Catmull–Clark subdivision type. In the classical theory, a subdivision surface is defined as the limit of the repetitive application of subdivision rules to the control grid. Based on Stam’s idea, the labour-intensive process can be avoided by using a natural parameterization of the limit surface. However, the simplification is not free of defects. At singularities, the smoothness of the classically defined limit surface has been lost. This paper describes how to rescue the parameterization by using a subdivision basis function that is consistent with the classical definition, but is expensive to compute. Based on this, we introduce a characteristic subdivision finite element and use it to discretize integrals on subdivision surfaces. We show that in the integral representation the complicated parameterization reduces to a decisive factor. We compare the natural and the characteristic subdivision finite element approach solving PDEs on surfaces. As model problem we consider the mean curvature flow, whereby the computation is done on the step-by-step changing geometry.     

Multiresolution Surfaces having Arbitrary Topologies by a Reverse Doo Subdivision Method

We have shown how to construct multiresolution structures for reversing subdivision rules using global least squares models (Samavati and Bartels, Computer Graphics Forum, 18(2):97–119, June 1999). As a result, semiorthogonal wavelet systems have also been generated. To construct a multiresolution surface of an arbitrary topology, however, biorthogonal wavelets are needed. In Bartels and Samavati (Journal of Computational and Applied Mathematics, 119:29–67, 2000) we introduced local least squares models for reversing subdivision rules to construct multiresolution curves and tensor product surfaces, noticing that the resulting wavelets were biorthogonal (under an induced inner product). Here, we construct multiresolution surfaces of arbitrary topologies by locally reversing the Doo subdivision scheme. In a Doo subdivision, a coarse surface is converted into a fine one by the contraction of coarse faces and the addition of new adjoining faces. We propose a novel reversing process to convert a fine surface into a coarse one plus an error. The conversion has the property that the subdivision of the resulting coarse surface is locally closest to the original fine surface, in the least squares sense, for two important face geometries. In this process, we first find those faces of the fine surface which might have been produced by the contraction of a coarse face in a Doo subdivision scheme. Then, we expand these faces. Since the expanded faces are not necessarily joined properly, several candidates are usually at hand for a single vertex of the coarse surface. To identify the set of candidates corresponding to a vertex, we construct a graph in such a way that any set of candidates corresponds to a connected component. The connected components can easily be identified by a depth first search traversal of the graph. Finally, vertices of the coarse surface are set to be the average of their corresponding candidates, and this is shown to be equivalent to local least squares approximation for regular arrangements of triangular and quadrilateral faces.     

形状可调的Loop 细分曲面渐进插值方法

针对Loop 细分无法调整形状与不能插值的问题,提出了一种形状可调的Loop 细分 曲面渐进插值方法。首先给出了一个既能对细分网格顶点统一调整又便于引入权因子实现细分曲 面形状可调的等价Loop 细分模板。其次,通过渐进迭代调整初始控制网格顶点生成新网格,运 用本文的两步Loop 细分方法对新网格进行细分,得到插值于初始控制顶点的形状可调的Loop 细分曲面。最后,证明了该方法的收敛性,并给出实例验证了该方法的有效性。     

Point-augmented biquadratic C1 subdivision surfaces

Shape artifacts, especially for convex input polyhedra, make Doo and Sabin’s generalization of bi-quadratic (bi-2) subdivision surfaces unattractive for general design. Rather than tuning the eigenstructure of the subdivision matrix, we improve shape by adding a point and enriching the refinement rules. Adding a guiding point can also yield a polar bi-2 subdivision algorithm. Both the augmented and the polar bi-2 subdivision are complemented by a new Primal Bi-2 Subdivision scheme. All surfaces are $C^1$ and can be combined.     

An adjustable subdivision surface modeling scheme

Based on triangle and quadrilateral meshes, this paper presents an adjustable subdivision surface scheme. The scheme can produce subdivision surface of Cl continuity of limit surface Since an adjustable parameter is introduced to the scheme, the surface modeling is flexible. Depended on given initial data, the limited surface shape can be adjusted and controlled through selecting appropriate parameters. The method is effective in generating smooth surfaces.     

基于双参数四点细分法的曲面造型

将双参数四点细分曲线方法进行推广，提出了基于双参数四点细分法的曲面造型方法，并对其收敛性进行了分析。该方法通过对两个参数的适当调节能够较容易地控制极限曲面的形状，极限曲面能够达到C4连续，可以应用到对曲面的连续性要求较高的曲面造型中去。在给定初始数据的条件下，可通过对形状参数的适当选择来实现对极限曲面的形状调整和控制，试验表明该算法生成光滑曲面是有效的。     

用逼近型√3细分方法构造闭三角网格的插值曲面

为了避免用逼近型3~(1/2)细分方法构造插值曲面过程中出现的烦琐运算,利用3细分方法极限点计算公式,提出一种用逼近型3~(1/2)细分方法构造闭三角网格插值曲面的方法.给定待插值的闭三角网格,先用一个新的几何规则与原3~(1/2)细分方法的拓扑规则细分一次得到一个初始网格,用3~(1/2)细分方法细分该初始网格得到插值曲面;新几何规则根据极限点公式确定,保证了初始网格的极限曲面插值待插值的三角网格.由于初始网格的顶点仅与待插值顶点2邻域内的点相关,所以插值曲面具有良好的局部性,即改变一个待插值点的位置时,只影响插值曲面在其附近的形状.该方法中只有确定初始网格顶点的几何规则与原3细分方法不同,故易于整合到原有的细分系统中.实验结果表明,该方法具有计算简单、有充分的自由度调整插值曲面的形状等特点,使得利用3~(1/2)细分方法构造三角网格的插值曲面变得极其简单.     

B样条的p-nary细分

有关B样条曲线曲面的binary细分技巧及其应用的研究已经获得了许多成果,建立在B样条binary细分基础上的binary细分法收敛性连续性分析的生成多项式法就是其中之一。该文研究了B样条曲线的p-nary细分问题,给出并证明了B样条基函数的p尺度细分方程中细分系数的计算公式及其性质,讨论了用p-nary细分生成非有理及有理B样条曲线的细分规则。采用该文的方法可方便而快速地在计算机上绘制有理B样条曲线。文章的结果可用于对一般p-nary曲线细分法收敛性及连续性的分析。     

Subdivision method to create furcating object with multibranches

A novel interim core scheme (ICS) is presented in this paper to construct a furcating object with multibranches. These M branches with arbitrary N-sided boundaries can be positioned freely but cannot be overlapped with each other. A furcating object can be built by blending these branches. The essence of the scheme is to construct a joint mesh that blends the initial control meshes of the M branches, and the smoothness of the resulting surfaces will only depend on the joint mesh and subdivision scheme applied. Some illustrative objects are given to verify the feasibility of ICS.     

A New Class of Guided C2 Subdivision Surfaces Combining Good Shape with Nested Refinement

Converting quadrilateral meshes to smooth manifolds, guided subdivision offers a way to combine the good highlight line distribution of recent G‐spline constructions with the refinability of subdivision surfaces. This avoids the complex refinement of G‐spline constructions and the poor shape of standard subdivision. Guided subdivision can then be used both to generate the surface and hierarchically compute functions on the surface. Specifically, we present a C² subdivision algorithm of polynomial degree bi‐6 and a curvature bounded algorithm of degree bi‐5. We prove that the common eigenstructure of this class of subdivision algorithms is determined by their guide and demonstrate that their eigenspectrum (speed of contraction) can be adjusted without harming the shape. For practical implementation, a finite number of subdivision steps can be completed by a high‐quality cap. Near irregular points this allows leveraging standard polynomial tools both for rendering of the surface and for approximately integrating functions on the surface.     

Non-Uniform Subdivision Surfaces with Sharp Features

Sharp features are important characteristics in surface modelling. However, it is still a significantly difficult task to create complex sharp features for Non-Uniform Rational B-Splines compatible subdivision surfaces. Current non-uniform subdivision methods produce sharp features generally by setting zero knot intervals, and these sharp features may have unpleasant visual effects. In this paper, we construct a non-uniform subdivision scheme to create complex sharp features by extending the eigen-polyhedron technique. The new scheme allows arbitrarily specifying sharp edges in the initial mesh and generates non-uniform cubic B-spline curves to represent the sharp features. Experimental results demonstrate that the present method can generate visually more pleasant sharp features than other existing approaches.     

Surface interpolation of meshes by geometric subdivision

Subdivision surfaces are generated by repeated approximation or interpolation from initial control meshes. In this paper, two new non-linear subdivision schemes, face based subdivision scheme and normal based subdivision scheme, are introduced for surface interpolation of triangular meshes. With a given coarse mesh more and more details will be added to the surface when the triangles have been split and refined. Because every intermediate mesh is a piecewise linear approximation to the final surface, the first type of subdivision scheme computes each new vertex as the solution to a least square fitting problem of selected old vertices and their neighboring triangles. Consequently, sharp features as well as smooth regions are generated automatically. For the second type of subdivision, the displacement for every new vertex is computed as a combination of normals at old vertices. By computing the vertex normals adaptively, the limit surface is G¹ smooth. The fairness of the interpolating surface can be improved further by using the neighboring faces. Because the new vertices by either of these two schemes depend on the local geometry, but not the vertex valences, the interpolating surface inherits the shape of the initial control mesh more fairly and naturally. Several examples are also presented to show the efficiency of the new algorithms.