Matrix-valued 4-point Spline and 3-point Non-spline Interpolatory Curve Subdivision Schemes
|
|
- Carmel Carroll
- 5 years ago
- Views:
Transcription
1 Matrix-valued 4-point Spline and -point Non-spline Interpolatory Curve Subdivision Schemes Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis, MO 6, U.S.A. Abstract The objective of this paper is to study and construct matrix-valued templates for interpolatory curve subdivision. Since our investigation of this problem was motivated by the need of such subdivision stencils as boundary templates for interpolatory surface subdivision, we provide both spline and non-spline templates that are necessarily symmetric, due to the lack of direction-orientation in carrying out surface subdivision in general. For example, the minimum-supported Hermite interpolatory C cubic spline curve subdivision scheme, with the skew-symmetric basis function for interpolating first derivatives, does not meet the symmetry specification. Non-spline C interpolatory templates constructed in this paper are particularly important, due to their smaller support needed to minimize undesirable surface oscillations, when adopted as boundary templates for interpolatory C surface subdivision. The curve subdivision templates introduced in this paper are adopted as boundary stencils for interpolatory surface subdivision with matrix-valued templates.. Introduction It is well-known that subdivision schemes provide the most efficient way for generating curves and surfaces, see e. g. [5,. To construct a smooth curve in the plane or in the -dimensional space, the subdivision process is carried out by taking weighted averages of the vertices iteratively, using the weights as specified by some subdivision templates. The initial vertices, called control points, are vertices of a given (control) polygon, and the number of vertices is increased after each iterative step, thereby generating a sequence of finer and finer polygons that eventually converge to a desirable limiting curve, called subdivision curve. For instance, when a dyadic subdivision scheme is employed, a finer polygon with twice as many vertices is generated from the previous polygon for each iterative step. This is accomplished by creating one new vertex (called odd numbered vertex) between every two consecutive old vertices (called even numbered vertices) of the coarser polygon (i.e. the polygon before the iterative step is carried out), while the positions of the old vertices may or may not be updated, again according to Research supported by ARO Grant #W9NF This author is also with the Department of Statistics, Stanford University, Stanford, CA 945. Research supported by UM Research Board /5 and UMSL Research Award /6.
2 the subdivision templates. If the even numbered vertices are not to be updated, then the vertices of each coarser polygon are among the vertices of the finer polygon. This subdivision process is then called an interpolatory subdivision scheme. Otherwise, it is called an approximation subdivision scheme. The subdivision templates of a dyadic curve subdivision scheme are constructed by following some local averaging rule, usually derived from certain refinement equation φ(x) = k Z p k φ(x k), x IR, (.) where φ(x) is called a refinable function, with dilation factor and corresponding (subdivision) refinement mask {p k }, which is necessarily a finite sequence. For any control polygon with vertices vk, an application of the refinement equation (.), immediately yields the local averaging rule: v m+ j = k v m k p j k, m =,,, (.) where, for each m =,,, the point-set {vk m} k denotes the set of vertices obtained after taking m iterations according to this local averaging rule. To facilitate implementational convenience, the local averaging rule can also be described in terms of two subdivision templates, by considering the index j in (.) to be even and odd integers, respectively. For example, the local averaging rule for the subdivision mask {,, 9,, 9,, }, introduced independently in [4 and [5, can be re-formulated as: v m+ i+ = 6 ( vm i + 9vi m + 9vi+ m vi+), m v m+ i = vi m, m ; and therefore may be described by the two subdivision templates shown in the middle and right of the three diagrams in Fig.. Here, solid circles represent the old vertices (i.e. before the iterative subdivision step is carried out) and the hollow circles represent the new vertices obtained after each iterative subdivision step. Observe that in this example the first subdivision template (shown in the middle) is a single solid circle. This means that the old vertices are not moved, and the subdivision scheme is interpolatory. On the other hand, the second template (shown on the right) consists of four solid circles and one hollow circle (for the new vertex). Hence, this subdivision scheme is called a 4-point interpolatory scheme Figure : Templates of 4-point interpolatory subdivision As illustrated by the above example, a (scalar) subdivision scheme is interpolatory, if and only if the subdivision mask {p k } satisfies the condition: p =, p k =, k,
3 which, under certain mild condition, is equivalent to the interpolatory property φ(k) = δ k, k Z of the scaling function φ(x). It is well-known, as proved in [5, that the above 4-point interpolatory scheme is a C scheme, meaning that the (limiting) subdivision curve is C in terms of the parameter on the real line. On the other hand, to construct C interpolatory schemes, it is necessary to admit larger templates (that is, more solid circles in the diagram on the right of Fig. ), unless a larger integer dilation factor is considered. In this regard, the interested reader is referred to [ for some 6-point C interpolatory scheme, and [6 for some 4-point C ternary interpolatory scheme. However, the subdivision curves so obtained are not spline curves; that is, the corresponding scaling functions are not (piecewise polynomial) spline functions. More recently, in the construction of multi-wavelets, the refinement equation (.) is extended to vector refinement, which naturally leads to matrix-valued subdivision masks. For example, in [, we studied matrix-valued interpolatory subdivision schemes for surface design and illustrated that matrix-valued C interpolatory schemes with small subdivision templates can be constructed. In fact, the construction in [ includes matrixvalued compactly supported spline scaling functions that assure the subdivision surfaces (at regular vertices) to be spline surfaces. The main objective of this paper is to apply the method introduced in [ to construct matrix-valued 4-point C spline interpolatory schemes and -point C (non-spline) interpolatory schemes for curve design. In addition to being interpolatory, the curve subdivision templates constructed in this paper are symmetric, and therefore can be readily adopted as boundary subdivision templates for matrix-valued interpolatory surface subdivision. This paper is organized as follows. In Section, we briefly recall the definition and some basic properties of matrix-valued interpolatory schemes introduced in [. Construction of some 4-point C quadratic-spline and C cubic-spline interpolatory schemes, both with symmetric subdivision templates, is demonstrated in Section. A 4-point C cubic-spline interpolatory scheme, again with symmetric subdivision templates, is introduced in Section 4. In Section 5, a -point C non-spline symmetric interpolatory scheme is presented. Finally, these curve subdivision schemes are adopted as boundary schemes for matrix-valued interpolatory surface subdivision in Section 6.. Matrix-valued interpolatory schemes for curve design Matrix-valued subdivision schemes for curve generation are associated with some vector refinement equation Φ(x) = k Z P k Φ(x k), x IR, (.) with an r-dimensional vector-valued refinable function Φ = [φ,, φ r T (also called refinable function vector), and refinement (or two-scale) sequence of r-dimensional square matrices {P k }, called a refinement (or subdivision) mask. The refinement equation (.)
4 can also be written as Φ(ω) = P ( ω) Φ( ω), where P (ω) = k P ke ikω, is called the (twoscale) symbol of Φ. Here and throughout, Φ = [ φ,, φ r T, where φ l denotes the Fourier transform of φ l. In the following, Φ is always assumed to satisfy the property of partition of unity, meaning that there exists a nonzero constant r-dimensional vector y such that y Φ(x k) = const, x IR. (.) k Z When the vector refinement equation (.) is applied to curve subdivision, the local averaging rule (.) is extended to the matrix setting: v m+ k = j v m j P k j, m =,,, (.) where v m k =: [v m k, s m k,, s m k,r (.4) are row-vectors with r components of vectors vk m, sm k,i, i =,, r. For the vector y = [,,,, as chosen in our earlier work [,,, we may use the first components vk m to denote the vertices of the -D subdivision polygons generated after the m-th iterative step, with initial vertices vk being the control points of the curve subdivision. The other components s k,,, s k,r of v k, can be used to control the geometric shape of the limiting curve. Then, it can be shown by applying the property of partition of unity, as in [, that the vertices vk m provide an accurate discrete approximation of the target subdivision curve, formulated by the series representation: F (x) = vkφ (x k) + ( ) s k,φ (x k) + + s k,r φ r (x k). k k For the matrix-valued surface subdivision, the other components s k,,, s k,r of v k are called the shape control pararemeters in [,,. It was demonstrated in [ that the shape control parameters can change the shape of the limiting surface significantly. The local averaging rule for a matrix-valued scheme can also be described in terms of subdivision templates. For example, the two templates for the local averaging rule of any subdivision mask {P k } k=,, are shown in Fig.. P P P P P P P Figure : Templates of matrix-valued 4-point subdivision scheme Following [ in the study of surface subdivision schemes, we have the following definition for matrix-valued interpolatory curve subdivision schemes. Definition A matrix-valued (dyadic) subdivision scheme with matrix-valued mask {P k } is said to be interpolatory, if for any given control vectors vk, the first components vk m of vm k, m =,,, in (.) remain the same as the corresponding first components of coarser vectors v m k ; i.e., vk m = vm k, m =,,. 4
5 The following result on characterization of interpolatory subdivision masks for the matrix-valued setting can be easily established. Proposition A (dyadic) subdivision scheme with matrix-valued mask {P k } is interpolatory if and only if P j satisfies P =..., P k =..., k, k Z; (.5) that is, the first column of P must be the unit vector [,,, T and the first columns of the other P k (with even indices) must be the zero vector. Let Φ = [φ, φ,, φ r T be a refinable function vector with subdivision mask {P k }. Then under certain mild conditions (see, for example, [ and [4), the subdivision scheme associated with {P k } is interpolatory if and only if Φ satisfies: φ (k) = δ(k), φ (k) =,, φ r (k) =, k Z. (.6). Matrix-valued 4-point C quadratic and cubic spline interpolatory scheme This section is devoted to the study of cardinal spline interpolation by compactly supported, symmetric, and refinable piecewise-polynomial functions, where cardinal interpolation means interpolation of the data δ k, k Z. The statement of this problem itself demands an explanation, since the problem of spline interpolation is well, and perhaps, overly studied in the published literature in both mathematics and engineering (see particularly [9,, for quadratic spline interpolation). However, the specification of the problem, as governed by the adjectives, compactly supported, symmetric, and refinable, does impose some challenge, and this particular specification must be met in constructing interpolatory curve subdivision templates, that can be adopted as boundary stencils for interpolatory surface subdivision, (the topic of discussion in Section 6). In particular, in application to connecting individually designed subdivision surface patches in the -space, symmetry of the (curve) subdivision templates would compensate for the lack of orientation of the surface subdivision schemes. For example, the minimum-supported Hermite interpolatory C cubic spline curve subdivision scheme, with the skew-symmetric basis function for interpolating first derivatives, does not meet the symmetry specification; and the issue of refinability of basis functions were not considered in the quadratic spline interpolation papers [9,. For spline-based subdivision, we only need to study quadratic and cubic spline interpolation. Let B (x) be the normalized C quadratic cardinal B-spline with support [, 5
6 . In the following, we consider some spline function ϕ with support [-, and symmetric with respect to x =. (In this regard, observe that the spline function B (x+) with support [, lacks symmetry). Let ϕ (x) := B (x + ) + B (x + ). (.) Note that ϕ is not refinable, in the sense that ϕ ( x ) is not a (finite) linear combination of ϕ (x k), k Z. Yet we may still write ˆϕ (ω) = 8 ( + eiω + e iω + e i ω + e i ω ) ˆϕ (ω), (.) which follows directly from the definition of ϕ and the refinement of B (x) : ˆB (ω) = 8 ( + e iω ) ˆB (ω). Next, we introduce ϕ (x) = ϕ (x) + (ϕ (x + ) + ϕ (x )), namely: ϕ (x) := B (x + ) + B (x + ) + (B (x + ) + B (x)). (.) The Bézier coefficients for ϕ and ϕ are shown in Figs. and 4 respectively. In the following proposition, we will see that the vector-valued function Φ a := [ϕ, ϕ T is refinable. / /.5.5 Figure : Support and Bézier coefficients for ϕ / / 5/6 5/6 / / Figure 4: Support and Bézier coefficients for ϕ Proposition Let ϕ and ϕ be the compactly supported C quadratic spline functions defined by (.) and (.) respectively. Then Φ a = [ϕ, ϕ T is refinable with refinement mask {P k } given by P = [ , P = P = [ 4 4 [, P = P = 6, P k =, k >. (.4) We remark that the matrices P k in (.4) do not satisfy the conditions in Proposition, but they will be later modified to satisfy these conditions. Proof of Proposition. From ϕ (x) = ϕ (x) + (ϕ (x + ) + ϕ (x )), or equivalently, ˆϕ (ω) = ˆϕ (ω) + ω (ei + e i ω ) ˆϕ (ω), 6
7 and in view of (.), we have On the other hand, since ˆϕ (ω) = 8 ( + eiω + e iω ) ˆϕ (ω) + 4 ( ˆϕ (ω) ˆϕ (ω)) = 8 ( + eiω + e iω ) ˆϕ (ω) + 8 ˆϕ (ω). (.5) ˆϕ (ω) = ˆϕ (ω) + (eiω + e iω ) ˆϕ (ω) = ( + (eiω + e iω )){ 8 ( + eiω + e iω ) ˆϕ (ω) + 8 ˆϕ (ω)} = 4 (5 + eiω + e iω + e iω + e iω ) ˆϕ (ω) + 8 ( + eiω + e iω ) ˆϕ (ω),(.6) It follows, by combining (.5) and (.6), that Φ a is refinable with refinement mask given by (.4). The refinement mask of Φ a yields the subdivision templates for the local averaging rule, as shown in Fig. 5 with G = P, H = P, L = P. L G L H H Figure 5: Local averaging rule Since P, P do not satisfy (.5), this scheme is not interpolatory. However, since ϕ is supported on [-,, we can modify ϕ by subtracting from ϕ (x) a suitable finite linear combination of ϕ (x k). Observe that ϕ () =, ϕ () = ϕ ( ) =. So, if we let ϕ (x) := ϕ (x) ϕ (x) (ϕ (x ) + ϕ (x + )), then ϕ satisfies (.6). Therefore, [ϕ, ϕ (x) T will yield a matrix-valued interpolatory scheme, provided that it is refinable. However, although it is indeed refinable, we will consider instead, Φ a := [ϕ + ϕ 8, ϕ 8 T, in the following, since Φ a satisfies the condition of generalized partition of unity (.) with y = [,. Proposition Let ϕ and ϕ be the compactly supported C quadratic spline functions defined by (.) and (.) resp. Set ϕ (x) = 5 8 ϕ (x) + 8 ϕ (x) 8 (ϕ (x ) + ϕ (x + )); ϕ (x) = 8 (ϕ (x) ϕ (x)) 8 (ϕ (x ) + ϕ (x + )). Then Φ a = [ ϕ, ϕ T is refinable with refinement mask { P k } given by P = [ 6 8, P = 6 8 P = [ 9, (.7) 6 7 P = P = [ 4, P = 6 4 P = [, Pk =, k >. 6 7
8 Proof. With z = e iω, let U(z) = [ (z + z ) (z + z ) 8. Then Φa (ω) = U(z) Φ a (ω). Thus, the two-scale symbol P (ω)=u(z )P (ω)u(z) of Φ a is [ z 6z + 9(z 4 + z ) (z 6 + ) 8z (z 4 + z ) 4(z 5 + z) + (z 6 + ) (z 4 + z ) (z 6 + ) 8z + 7(z 4 + z ) 4(z 5 + z) + (z 6. + ) Therefore, Φ a is refinable with refinable mask given by (.7). Z W Z Y X X Y Figure 6: Matrix-valued 4-point interpolatory scheme Clearly, the mask of Φ a immediately yields an interpolatory scheme with templates shown in Fig. 6 with W = P, X = P, Y = P, Z = P. We call such a scheme a matrix-valued 4-point interpolatory scheme. Figure 7: Support and Bézier coefficients for f We can also construct the matrix-valued C cubic 4-point interpolatory scheme by the same method. To be more precise, let f (x) be the C cubic spline supported on [-, with its Bézier coefficients shown in Fig. 7. Note that f is not refinable. So, let f (x) = f ( x), and by applying the C -smoothing formula, we have the Bézier representations f (x) and f ( x ), as shown in Fig. 8 and Fig. 9, respectively. By evaluating the Bézier coefficients of f ( x) and f (x k), f (x k), we can show that Φ b (x) := [f (x), f (x) T is refinable and its refinement mask can be easily computed. Furthermore, we can show that Φ b (x) := [ f (x), f (x) T yields a matrix-valued 4-point interpolatory scheme, where { f (x) = f (x) + f (x) (f (x ) + f (x + )), f (x) = (f (x) f (x)) (f (.8) (x ) + f (x + )). More precisely, we have the following result. /4 / /4 /4 / /4 Figure 8: Support and Bézier coefficients for f (x)(= f ( x )) 8
9 5/6 7/ /4 5/8 / /8 /4 5/ /6 4 Figure 9: Bézier coefficients for f ( x )(= f ( x 4 )) on [, 4 Proposition 4 Suppose f is the compactly supported C cubic spline functions with Bézier coefficients shown in Fig. 7. Let f (x) = f ( x), and f, f be the functions defined by (.8). Then both Φ b = [f, f T and Φ b = [ f, f T are refinable, with refinement masks, (G k ) k and ( G k ), given respectively by [ [ G =, G = G =, and 4 4 [ G = G = 8 G = [ G = G = 6 8 7, G = G =, G = G = 6 [ [ [ 9 7, G = G = 6, [, G k =, k >,, Gk =, k >. Again, the mask of Φ b yields a matrix-valued 4-point interpolatory scheme with templates shown in Fig. 6, where W = G, X = G, Y = G, Z = G. 4. Matrix-valued 4-point C cubic spline interpolatory scheme In this section we consider C cubic spline functions and show that the above approach can be applied to construct a matrix-valued 4-point C cubic-spline interpolatory scheme. Let B 4 (x) be the normalized cardinal C cubic B-spline with support [, 4, and φ (x) = B 4(x + ). Then φ is supported in [-, and symmetric with respect to x =. However, φ is not refinable. More precisely, from the definition of φ and the refinement of B 4 (x) : ˆB 4 (ω) = 6 ( + e iω ) 4 ˆB4 (ω), we have ˆφ (ω) = 6 (6 + eiω + e iω ) ˆφ (ω) + 4 (ei ω + e i ω ) ˆφ (ω). (4.) Next, we define φ (x) = 7φ 9 (x) + 4(φ 9 (x + ) + φ (x )), namely: φ (x) := 7 6 B 4(x + ) + (B 4(x + ) + B 4 (x + )). (4.) The Bézier coefficients for φ and φ are shown in Figs. and resp., and the vectorvalued function Φ c := [φ, φ T is refinable. More precisely, we have the following result. 9
10 /4 / / /4.5.5 Figure : Support and Bézier coefficients for φ /9 /9 4/9 /6 5/6 5/6 /6 4/9 /9 / Figure : Support and Bézier coefficients for φ Proposition 5 Let φ = B 4(x + ) and φ be defined by (4.). Then Φ c = [φ, φ T is refinable and its associated refinement mask {H k } is given by [ [ 9 [ H = 8 8, H = H = 8, H = H =, H k =, k >.(4.) Proof. From φ (x) = 7φ 9 (x) + 4(φ 9 (x + ) + φ (x )), we have ˆφ (ω) = 7 9 ˆφ (ω) (ei ω + e i ω ) ˆφ (ω). This, together with (4.), leads to ˆφ (ω) = 6 (6 + eiω + e iω ) ˆφ (ω) + 4 (9 4 ˆφ (ω) 7 4 ˆφ (ω)) = 6 ( + eiω + e iω ) ˆφ (ω) ˆφ (ω). (4.4) On the other hand, we have ˆφ (ω) = 9 (7 + 4eiω + 4e iω ) ˆφ (ω) = 9 (7 + 4eiω + 4e iω )){ 6 ( + eiω + e iω ) ˆφ (ω) ˆφ (ω)} = 44 ( + eiω + e iω + 4e iω + 4e iω ) ˆφ (ω) + 6 (7 + 4eiω + 4e iω ) ˆφ (ω).(4.5) From (4.4) and (4.5), we see that Φ c is refinable and its associated refinement mask is given by (4.). The refinement mask of Φ c yields the templates for the local averaging rule as shown in Fig. 5 with G = H, H = H, L = H. In this case the vector y for the partition unity property of Φ c is [, 9. In practice, if such a matrix-valued C cubic approximation subdivision scheme is used for curve design, we would recommend the choice of Φ d := [φ +9φ, 9φ +9φ T, since the vector y for the generalized partition of unity property of Φ d is [,, and Φ d is refinable with refinement mask, denoted as (H (d) ) k, given by H (d) = [ 9 H (d) = H (d) = [, H (d) = H (d) = [ 8 8, H (d) k =, k >., k
11 The scheme derived from Φ c is not interpolatory, since H, H, H do not satisfy (.5). However, we can modify φ as in the above section such that the resulting refinable function vector will yields an interpolatory scheme. Observe that φ () =, φ () = φ ( ) = 9. Then φ (x) := 4 (φ (x) φ (x)) (φ (x ) + φ (x + )) (4.6) satisfies (.6). We can show that [φ, φ T is refinable, and therefore, its associated mask yields a matrix-valued interpolatory scheme. However, the vector y for the property of generalized partition of unity (.) is [,. Hence, in the following we consider, instead, Φ c := [ φ, φ T so that the vector y for generalized partition of unity (.) of Φ c is [,, where φ = φ + φ, namely, φ (x) := 4 φ (x) + 4 φ (x) (φ (x ) + φ (x + )). (4.7) Proposition 6 Suppose φ = B 4(x + ) and φ is defined by (4.). Let φ and φ be the C cubic spline functions defined by (4.7) and (4.6). Then Φ c = [ φ, φ T is refinable and the refinement mask is given by H = [ 48, H = 48 H = [ 5, (4.8) 48 H = H = [ 6, H = 48 6 H = [, Hk =, k >. 48 Proof. Let [ U (z) = (z + 4 z ) 4 (z + 4 z ) 4, where z = e iω. Then Φc (ω) = U (z) Φ c (ω). Thus, the two-scale symbol H(ω)=U (z )H(ω)U (z) of Φ c is [ z 48z + 5(z 4 + z ) (z 6 + ) z (z 4 + z ) 6(z 5 + z) + (z 6 + ) 96 (z 4 + z ) (z 6 + ) z + (z 4 + z ) 6(z 5 + z) + (z 6. + ) Therefore Φ c is refinable with refinement mask given by (4.8), as desired; and Φ c yields a matrix-valued 4-point interpolatory scheme with templates shown in Fig. 6, where W = H, X = H, Y = H, Z = H. 5. Matrix-valued -point C non-spline interpolatory scheme In the previous two sections, we have constructed matrix-valued spline interpolatory schemes; namely, schemes with spline functions as refinable functions. In this section, we will construct -point non-spline C interpolatory schemes with templates shown in Fig.
12 . More precisely, we consider the refinable function vector Φ e (x) := [φ e (x), φ e (x) T that satisfies the refinement relation: with Q = [ t 4 (t + t ) Φ e (x) = Q k Φ e (x k), k= [, Q = Q = 8 8 [ t, Q = Q = t, t, t IR. The mask {Q k } k=,,,, satisfies the so-called sum rule of order 4 (refer to e. g., [7), which implies that the integer shifts φ e (x k), φ e (x k), k Z reproduce polynomials of degrees up to. We can determine t, t by the Sobolev smoothness estimate formula in [8 for the refinable function vectors so that the resulting φ e, φ e are C. There are many choices of t, t for C φ e, φ e. For example, if we choose t = 64, t = 7 64, the resulting φ e, φ e are in the Sobolev space W.48 (IR). For such t, t, the corresponding Q k, k, are Q = W, Q = Q = X, Q = Q = Z, where [ [ [ W =, X =, Z = (5.) Thus, we have a matrix-valued -point C interpolatory scheme shown in Fig.. Z W Z X X Figure : Matrix-valued -point interpolatory scheme 6. Application to surface design In this section, we will apply the matrix-valued C interpolatory schemes for curve subdivision constructed in Sections -5 as boundary schemes for interpolatory surface design. However, for simplicity, we will only consider the triangular mesh and the interior boundary vertices which are not corner vertices. Furthermore, we will only illustrate our approach by considering the non-spline -ring C interpolatory surface subdivision schemes constructed in [ and [ for the interior vertices. The subdivision templates for regular vertices are shown in Fig., where the matrix-valued weights are given by [ [ [ [ 45 A = , D = 5 45, B = , C = (6.)
13 D D C D A D B B D D C Figure : Templates of -ring C interpolatory scheme for regular old ( even ) vertices (left) and new ( odd ) vertices (right) An (interior) boundary vertex on a triangular mesh is called a regular boundary vertex if its valence (number of its neighbors) is 4; otherwise, it is called an extraordinary boundary vertex. We will use the templates of the averaging rule in Fig. 6 for both regular and extraordinary boundary vertices. Using the same templates for both regular and extraordinary boundary vertices has the advantage in connecting two individually designed subdivision surface patches, as pointed out in [5. Here, we show that the eigenvalues λ, λ, (in the order of non-increasing magnitudes) of the subdivision matrices near boundary extraordinary vertices satisfy certain basic properties. In this regard, we need the following lemma. Lemma Let M n (s, t) be a -banded n n symmetric matrix defined by s M n (s, t) = t t s t t t s n n, (6.) where s, t IR. Then the eigenvalues of M n (s, t) are given by s t cos jπ n+, j n. The eigenvalues of a matrix of the form (6.) have been discussed in [. For completeness, we will provide the proof of Lemma in the Appendix.
14 u w 4 v u w v u n y 4 z y u w n v n f e z d d y 4 n f x c e d n 4 x (b ) (b 4 ) c d z n x e n n f (b n ) a 4 a x (b c ) f n n a n a c e y z w v u n y d n x n (b n ) a ( a ) a n n x x (b ) d y u Figure 4: -D representation of initial vectors and refined vectors after one step of matrixvalued subdivision iteration near the extraordinary boundary vertex x Let x be an initial vector associated with an extraordinary boundary vertex with valence n, and let x j, y j, u j, z k, v k, w k, j n, k n, be the initial vectors associated with vertices in the neighborhood of this extraordinary vertex. Also, let a, a j, b j, d j, c k, e k, f k, j n, k n, be the corresponding vectors after carrying out one iteration of the subdivision process with the averaging rules given in Figs. 6 and. See Fig. 4 for the -D representation of these initial vectors and refined vectors. Consider the subdivision matrix S, defined by [ [a, a, a n, b, b n, d, d n, [a, a,, a n, [c, c,, c n, [b, b,, b n, [d, d,, d n, [e, e,, e n, [f, f,, f n [ = [x, x, x n, y, y n, u, u n, [x, x,, x n, [z, z,, z n, [y, y,, y n, [u, u,, u n, [v, v,, v n, [w, w,, w n S. Then S is an upper-triangular block matrix with nonzero diagonal blocks S j,j, j 4
15 given by S, = S, = W X X Z Z Y Y Z X Y W X Z Y X W X Y Z X, S, = Y Z X Y Y C D C..., S, = C (n ) (n ) B C C B C D... D C C B (n ) (n ) (n ) (n )., Here, S, has (nonzero) eigenvalues and 5, each with multiplicity n, and 8 5 S, has (nonzero) eigenvalue 45 with multiplicity n. For S 5,, let us consider 5S,, and introduce the notation [ [ 9 64 B = 5B =, C = 5C = Then by interchanging the second and third rows, followed by interchanging the second and third columns of the matrix 5S,, we have the matrix S, = B C. C B C C C B Next, by moving the 5th row to become the third row, followed by moving the 5th 5
16 column to become the third column, of the matrix S,, we have the matrix S, = 64 9 C C B C C C B In general, this procedure is repeated by moving the (j + )th row to become the (j +)th row, followed by moving the (j +)th column to become the (j +)th column, for j n, to arrive at the block matrix: [ Mn (9, 64) M n ( 47, 7) M n (69, 5) where M n (s, t) is the matrix defined by (6.). Thus, S, is similar to [ Mn (, ) 8 8 M n ( 47, 7 ) M 5 5 n ( 69, 5 ). 5 5 Therefore, it follows from Lemma, that S, has eigenvalues 8 4 cos jπ n, cos jπ, j n, (6.) n among which = + cos π is the largest in magnitude. 8 4 n 8 4 n If we apply the matrix-valued 4-point C spline interpolatory scheme for the regular boundary vertices, namely, X = H, Y = H, W = H, Z = H with H k, k =,,, given in (4.8), then the nonzero eigenvalues of S, are cos (n )π,, 4, 8, 8. Hence, the eigenvalues λ k of S satisfy λ = > λ λ. Also, if we use the matrix-valued -point C non-spline interpolatory scheme for the regular boundary vertices, namely, X = X, Y =, W = W, Z = Z with W, X, Z given in (5.), then the nonzero eigenvalues of S, are,, 4, 8, 8, 6. 6,
17 Figure 5: Initial open mesh (left) and limiting surface (right) Thus, the eigenvalues λ k of S also satisfy λ = > λ λ. On the right of Fig. 5, we show the limiting surface by applying this scheme with control parameters s k = v k to an initial (open) mesh shown on the left of Fig. 5. Observe that when n > 7, the first two largest eigenvalues + cos π and n 8 π cos of S 4 n, listed in (6.) are bigger than, an eigenvalue of S. So, it is likely that the resulting surface is not C at boundary extraordinary vertices with large valences. To force the largest eigenvalue of S, to be, one can use a modified template for the new vertices next to an extraordinary boundary vertex, as suggested in [5. More precisely, for an extraordinary boundary vertex x with valence n, we may apply the template in Fig. 6 to a, a,, a n in Fig. 4, where B (n), B (n) are matrices depending on n, and C is given in (6.). C B (n) x B (n) C Figure 6: Template for new ( odd ) vertices next to an extraordinary vertex x In this case the subdivision matrix S is still an upper-triangular block matrix with exactly the same nonzero matrices S j,j except S,, which is now given by B (n) C C B (n) C S, = C C B (n) If we choose B (n) = γ n B, B (n) = ( γ n )B, 7 (n ) (n )
18 where γ n is a real number and B is given by (6.), then the nonzero eigenvalues of S, become 8 γ n 4 cos jπ n, 69 5 γ n cos jπ n, j n. Therefore, if γ n is so chosen that γ 8 n + cos π =, namely, 4 n γ n = 4 cos π n, (6.4) then is an eigenvalue of S, and it is the largest eigenvalue in magnitude of S,. Hence, with γ n as given in (6.4), the eigenvalues of S are,,, λ, λ 4, with > λ λ 4 for either the 4-point C spline scheme or -point C non-spline scheme to be used as the boundary scheme. Appendix Proof of Lemma. For t =, the result is trivial, since M n (s, t) becomes a diagonal matrix. In the following, we consider the case t. Let p n (λ) =det(m n (s, t) λi n ) be the polynomial in λ with degree n. Clearly, p (λ) = s λ, and with the notation p (λ) :=, we have p n = (s t)p n t p n, n. To determine the zeros of p n (λ), write s λ = t(z + z ), where z = e iθ for some θ IR\{}. Then we have Therefore, Hence, we have p n p n = (tz + tz )p n t p n. p n tz p n = tzp n t p n = tz(p n tz p n ) = (tz) (p n tz p n ) = = (tz) n (p tz p ) = (tz) n (s λ tz ) = (tz) n (tz + tz tz ) = (tz) n. = tz p n + (tz) n = tz (tz p n + (tz) n ) + (tz) n = (tz ) p n + t n z n + t n z n = = (tz ) n p + t n z n (n ) + + t n z n + t n z n = t n (z n + z n (n ) + + z n + z n ) = t n zn+ z (n+) n sin(n + )θ = t. z z sin θ Clearly, when θ = θ j = jπ sin(n+)θ, j n, we have =. Thus, λ n+ sin θ j = s t cos θ j are the zeros of p n ; and hence, they are the eigenvalues of M n (s, t). In addition, for t, the n values λ j, j n are distinct. Thus, λ j, j n, are the n eigenvalues of M n (s, t). 8
19 References [ C.K. Chui and Q.T. Jiang, Matrix-valued symmetric templates for interpolatory surface subdivisions I. Regular vertices, Appl. Comput. Harmonic Anal. 9 (5), 9. [ C.K. Chui and Q.T. Jiang, Matrix-valued subdivision schemes for generating surfaces with extraordinary vertices, Comput. Aided Geom. Design (6), [ C.K. Chui and Q.T. Jiang, From extension of Loop s approximation scheme to interpolatory subdivisions, Comput. Aided Geom. Design 5 (8), [4 S. Dubuc, Interpolation through an iterative scheme, J. Math. Anal. Appl. 4 (986), [5 N. Dyn, J.A. Gregory and D. Levin, A 4-point interpolatory subdivision scheme for curve designs, Comput. Aided Geom. Design 4 (987), [6 M.F. Hassan, I. P. Ivrissimitzis, N. A. Dodgson, and M. A. Sabin, An interpolating 4-point C ternary stationary subdivision scheme, Comput. Aided Geom. Design 9 (), 8. [7 R.Q. Jia, and Q.T. Jiang, Approximation power of refinable vectors of functions, in: Wavelet Analysis and Applications, AMS/IP Stud. Adv. Math., Vol. 5, Amer. Math. Soc., Providence, RI,, pp [8 R.Q. Jia and Q.T. Jiang, Spectral analysis of transition operators and its applications to smoothness analysis of wavelets, SIAM J. Matrix Anal. Appl. 4 (), 7 9. [9 S.S. Rana, Quadratic spline interpolation, J. Approx. Theory 57 (989), 5. [ U. Reif and J. Peters, Structural analysis of subdivision surfaces a summary, in Topics in Multivariate Approximation and Interpolation, K. Jetter et al. (Eds.), 5. [ L.L. Schumaker, On shape preserving quadratic spline interpolation, SIAM J. Numer. Anal. (98), [ J.E. Schweitzer, Analysis and application of subdivision surfaces, Ph.D. Thesis, University of Washington, Dept. of Computer Science and Engineering, Seattle, 996. [ A. Weissman, A 6-point C interpolatory subdivision scheme for curve design, M. Sc. Thesis, Tel-Aviv University, 989. [4 D.X. Zhou, Multiple refinable Hermite interpolants, J. Approx. Theory (),
20 [5 D. Zorin and P. Schröder, A. DeRose, L. Kobbelt, A. Levin, and W. Sweldens, Subdivision for Modeling and Animation, SIGGRAPH Course Notes.
Refinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions
Refinable bivariate quartic and quintic C 2 -splines for quadrilateral subdivisions Charles K. Chui, Qingtang Jiang Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis,
More informationTangents and curvatures of matrix-valued subdivision curves and their applications to curve design
Tangents and curvatures of matrix-valued subdivision curves and their applications to curve design Qingtang Jiang, James J. Smith Department of Mathematics and Computer Science University of Missouri St.
More informationFrom Extension of Loop s Approximation Scheme to Interpolatory Subdivisions
From Extension of Loop s Approximation Scheme to Interpolatory Subdivisions Charles K. Chui, Qingtang Jiang 2 Department of Mathematics and Computer Science University of Missouri St. Louis St. Louis,
More informationDepartment of Mathematics and Computer Science University of Missouri St. Louis St. Louis, MO 63121
Triangular 7 and Quadrilateral 5 Subdivision Schemes: Regular Case Charles K. Chui 1, Qingtang Jiang 2 3, Rokhaya N. Ndao Department of Mathematics and Computer Science University of Missouri St. Louis
More informationAn interpolating 4-point C 2 ternary stationary subdivision scheme
Computer Aided Geometric Design 9 (2002) 8 www.elsevier.com/locate/comaid An interpolating 4-point C 2 ternary stationary subdivision scheme M.F Hassan a,, I.P. Ivrissimitzis a, N.A. Dodgson a,m.a.sabin
More informationHighly Symmetric Bi-frames for Triangle Surface Multiresolution Processing
Highly Symmetric Bi-frames for Triangle Surface Multiresolution Processing Qingtang Jiang and Dale K. Pounds Abstract In this paper we investigate the construction of dyadic affine (wavelet) bi-frames
More informationQuad/triangle Subdivision, Nonhomogeneous Refinement Equation and Polynomial Reproduction
Quad/triangle Subdivision, Nonhomogeneous Refinement Equation and Polynomial Reproduction Qingtang Jiang a,, Baobin Li b, a Department of Mathematics and Computer Science, University of Missouri St. Louis,
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn Michael S. Floater Kai Hormann Abstract. We present a new four-point subdivision scheme that generates C 2 curves.
More informationSubdivision Curves and Surfaces
Subdivision Surfaces or How to Generate a Smooth Mesh?? Subdivision Curves and Surfaces Subdivision given polyline(2d)/mesh(3d) recursively modify & add vertices to achieve smooth curve/surface Each iteration
More informationMatrix-valued Subdivision Schemes for Generating Surfaces with Extraordinary Vertices
Matrix-valued Subdivision Schemes for Generating Surfaces with Extraordinary Vertices Charles K. Chui, Qingtang Jiang Department of Mathematics & Computer Science University of Missouri St. Louis St. Louis,
More informationu 0+u 2 new boundary vertex
Combined Subdivision Schemes for the design of surfaces satisfying boundary conditions Adi Levin School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel. Email:fadilev@math.tau.ac.ilg
More informationNormals of subdivision surfaces and their control polyhedra
Computer Aided Geometric Design 24 (27 112 116 www.elsevier.com/locate/cagd Normals of subdivision surfaces and their control polyhedra I. Ginkel a,j.peters b,,g.umlauf a a University of Kaiserslautern,
More informationFourier transform of Bernstein Bézier polynomials
Illinois Wesleyan University From the SelectedWorks of Tian-Xiao He 27 Fourier transform of Bernstein Bézier polynomials Tian-Xiao He, Illinois Wesleyan University Charles K. Chui Qingtang Jiang Available
More informationTernary Butterfly Subdivision
Ternary Butterfly Subdivision Ruotian Ling a,b Xiaonan Luo b Zhongxian Chen b,c a Department of Computer Science, The University of Hong Kong b Computer Application Institute, Sun Yat-sen University c
More informationSubdivision on Arbitrary Meshes: Algorithms and Theory
Subdivision on Arbitrary Meshes: Algorithms and Theory Denis Zorin New York University 719 Broadway, 12th floor, New York, USA E-mail: dzorin@mrl.nyu.edu Subdivision surfaces have become a standard geometric
More informationA C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions
A C 2 Four-Point Subdivision Scheme with Fourth Order Accuracy and its Extensions Nira Dyn School of Mathematical Sciences Tel Aviv University Michael S. Floater Department of Informatics University of
More informationSubdivision of Curves and Surfaces: An Overview
Subdivision of Curves and Surfaces: An Overview Ben Herbst, Karin M Hunter, Emile Rossouw Applied Mathematics, Department of Mathematical Sciences, University of Stellenbosch, Private Bag X1, Matieland,
More informationFourier Transform of Bernstein-Bézier Polynomials
Fourier Transform of Bernstein-Bézier Polynomials Charles. K. Chui, Tian-Xiao He 2, and Qingtang Jiang 3 Department of Mathematics and Computer Science University of Missouri-St. Louis, St. Louis, MO 632,
More informationSubdivision based Interpolation with Shape Control
Subdivision based Interpolation with Shape Control Fengtao Fan University of Kentucky Deparment of Computer Science Lexington, KY 40506, USA ffan2@uky.edu Fuhua (Frank) Cheng University of Kentucky Deparment
More informationEdge Detections Using Box Spline Tight Frames
Edge Detections Using Box Spline Tight Frames Ming-Jun Lai 1) Abstract. We present a piece of numerical evidence that box spline tight frames are very useful for edge detection of images. Comparsion with
More informationA NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME
J. Appl. Math. & Computing Vol. 18(2005), No. 1-2, pp. 553-562 A NEW PROOF OF THE SMOOTHNESS OF 4-POINT DESLAURIERS-DUBUC SCHEME YOUCHUN TANG, KWAN PYO KO* AND BYUNG-GOOK LEE Abstract. It is well-known
More informationInterpolatory 3-Subdivision
EUROGRAPHICS 2000 / M. Gross and F.R.A. Hopgood (Guest Editors) Volume 19 (2000), Number 3 Interpolatory 3-Subdivision U. Labsik G. Greiner Computer Graphics Group University of Erlangen-Nuremberg Am Weichselgarten
More informationResearch Article A Family of Even-Point Ternary Approximating Schemes
International Scholarly Research Network ISRN Applied Mathematics Volume, Article ID 97, pages doi:.5//97 Research Article A Family of Even-Point Ternary Approximating Schemes Abdul Ghaffar and Ghulam
More informationSubdivision overview
Subdivision overview CS4620 Lecture 16 2018 Steve Marschner 1 Introduction: corner cutting Piecewise linear curve too jagged for you? Lop off the corners! results in a curve with twice as many corners
More informationOn Smooth Bicubic Surfaces from Quad Meshes
On Smooth Bicubic Surfaces from Quad Meshes Jianhua Fan and Jörg Peters Dept CISE, University of Florida Abstract. Determining the least m such that one m m bi-cubic macropatch per quadrilateral offers
More informationNormals of subdivision surfaces and their control polyhedra
Normals of subdivision surfaces and their control polyhedra I. Ginkel, a, J. Peters b, and G. Umlauf a, a University of Kaiserslautern, Germany b University of Florida, Gainesville, FL, USA Abstract For
More informationQuadratic and cubic b-splines by generalizing higher-order voronoi diagrams
Quadratic and cubic b-splines by generalizing higher-order voronoi diagrams Yuanxin Liu and Jack Snoeyink Joshua Levine April 18, 2007 Computer Science and Engineering, The Ohio State University 1 / 24
More informationA new 8-node quadrilateral spline finite element
Journal of Computational and Applied Mathematics 195 (2006) 54 65 www.elsevier.com/locate/cam A new 8-node quadrilateral spline finite element Chong-Jun Li, Ren-Hong Wang Institute of Mathematical Sciences,
More informationUsing Semi-Regular 4 8 Meshes for Subdivision Surfaces
Using Semi-Regular 8 Meshes for Subdivision Surfaces Luiz Velho IMPA Instituto de Matemática Pura e Aplicada Abstract. Semi-regular 8 meshes are refinable triangulated quadrangulations. They provide a
More informationA Study on Subdivision Scheme-Draft. Kwan Pyo Ko Dongseo University Busan, South Korea
A Study on Subdivision Scheme-Draft Kwan Pyo Ko Dongseo University Busan, South Korea April 30, 007 Contents Introduction 9 Subdivision of Univariate Data 3. Definitions and Basic Results.............................
More informationSubdivision Curves and Surfaces: An Introduction
Subdivision Curves and Surfaces: An Introduction Corner Cutting De Casteljau s and de Boor s algorithms all use corner-cutting procedures. Corner cutting can be local or non-local. A cut is local if it
More informationSubdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches
Subdivision Depth Computation for Extra-Ordinary Catmull-Clark Subdivision Surface Patches Fuhua Frank Cheng,GangChen, and Jun-Hai Yong University of Kentucky, Lexington, KY, USA Tsinghua University, Beijing,
More informationSubdivision Surfaces
Subdivision Surfaces 1 Geometric Modeling Sometimes need more than polygon meshes Smooth surfaces Traditional geometric modeling used NURBS Non uniform rational B-Spline Demo 2 Problems with NURBS A single
More informationC 1 Quintic Spline Interpolation Over Tetrahedral Partitions
C 1 Quintic Spline Interpolation Over Tetrahedral Partitions Gerard Awanou and Ming-Jun Lai Abstract. We discuss the implementation of a C 1 quintic superspline method for interpolating scattered data
More informationUNIVERSITY OF CALGARY. Subdivision Surfaces. Advanced Geometric Modeling Faramarz Samavati
Subdivision Surfaces Surfaces Having arbitrary Topologies Tensor Product Surfaces Non Tensor Surfaces We can t find u-curves and v-curves in general surfaces General Subdivision Coarse mesh Subdivision
More informationTechnical Report. Removing polar rendering artifacts in subdivision surfaces. Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin.
Technical Report UCAM-CL-TR-689 ISSN 1476-2986 Number 689 Computer Laboratory Removing polar rendering artifacts in subdivision surfaces Ursula H. Augsdörfer, Neil A. Dodgson, Malcolm A. Sabin June 2007
More informationGeneralizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that g
Generalizing the C 4 Four-directional Box Spline to Surfaces of Arbitrary Topology Luiz Velho Abstract. In this paper we introduce a new scheme that generalizes the four-directional box spline of class
More informationParameterization of triangular meshes
Parameterization of triangular meshes Michael S. Floater November 10, 2009 Triangular meshes are often used to represent surfaces, at least initially, one reason being that meshes are relatively easy to
More informationNon-uniform interpolatory curve subdivision with edge parameters built-upon compactly supported cardinal splines
Non-uniform interpolatory curve subdivision with edge parameters built-upon compactly supported cardinal splines Carolina Beccari a, Giulio Casciola a, Lucia Romani b a Department of Mathematics, University
More informationSmooth Patching of Refined Triangulations
Smooth Patching of Refined Triangulations Jörg Peters July, 200 Abstract This paper presents a simple algorithm for associating a smooth, low degree polynomial surface with triangulations whose extraordinary
More informationREGULAR GRAPHS OF GIVEN GIRTH. Contents
REGULAR GRAPHS OF GIVEN GIRTH BROOKE ULLERY Contents 1. Introduction This paper gives an introduction to the area of graph theory dealing with properties of regular graphs of given girth. A large portion
More informationCurve Corner Cutting
Subdivision ision Techniqueses Spring 2010 1 Curve Corner Cutting Take two points on different edges of a polygon and join them with a line segment. Then, use this line segment to replace all vertices
More informationEvaluation of Loop Subdivision Surfaces
Evaluation of Loop Subdivision Surfaces Jos Stam Alias wavefront, Inc. 8 Third Ave, 8th Floor, Seattle, WA 980, U.S.A. jstam@aw.sgi.com Abstract This paper describes a technique to evaluate Loop subdivision
More informationSubdivision. Outline. Key Questions. Subdivision Surfaces. Advanced Computer Graphics (Spring 2013) Video: Geri s Game (outside link)
Advanced Computer Graphics (Spring 03) CS 83, Lecture 7: Subdivision Ravi Ramamoorthi http://inst.eecs.berkeley.edu/~cs83/sp3 Slides courtesy of Szymon Rusinkiewicz, James O Brien with material from Denis
More informationSubdivision Surfaces
Subdivision Surfaces CS 4620 Lecture 31 Cornell CS4620 Fall 2015 1 Administration A5 due on Friday Dreamworks visiting Thu/Fri Rest of class Surfaces, Animation, Rendering w/ prior instructor Steve Marschner
More informationINF3320 Computer Graphics and Discrete Geometry
INF3320 Computer Graphics and Discrete Geometry More smooth Curves and Surfaces Christopher Dyken, Michael Floater and Martin Reimers 10.11.2010 Page 1 More smooth Curves and Surfaces Akenine-Möller, Haines
More informationStationary and Non-Stationary Univariate Subdivision Schemes
Punjab University Journal of Mathematics (ISSN 6-56) Vol. 5(3)(8) pp. 5- Stationary and Non-Stationary Univariate Subdivision Schemes Muhammad Asghar Department of Mathematics, The Islamia University of
More informationUniversity of Twente. Faculty of Mathematical Sciences. Convexity preservation of the four-point interpolatory subdivision scheme
Faculty of Mathematical Sciences University of Twente University for Technical and Social Sciences P.O. Box 17 7500 AE Enschede The Netherlands Phone: +31-53-4893400 Fax: +31-53-4893114 Email: memo@math.utwente.nl
More informationG 2 Bezier Crust on Quad Subdivision Surfaces
Pacific Graphics (2013) B. Levy, X. Tong, and K. Yin (Editors) Short Papers G 2 Bezier Crust on Quad Subdivision Surfaces paper 1348 Figure 1: Two examples of Bezier crust applied on Catmull-Clark subdivision
More informationNon-uniform interpolatory curve subdivision with edge parameters built upon compactly supported fundamental splines
BIT manuscript No. (will be inserted by the editor) Non-uniform interpolatory curve subdivision with edge parameters built upon compactly supported fundamental splines Carolina Vittoria Beccari Giulio
More informationApproximating Catmull-Clark Subdivision Surfaces with Bicubic Patches
Approximating Catmull-Clark Subdivision Surfaces with Bicubic Patches Charles Loop Microsoft Research Scott Schaefer Texas A&M University April 24, 2007 Technical Report MSR-TR-2007-44 Microsoft Research
More informationG 2 Interpolation for Polar Surfaces
1 G 2 Interpolation for Polar Surfaces Jianzhong Wang 1, Fuhua Cheng 2,3 1 University of Kentucky, jwangf@uky.edu 2 University of Kentucky, cheng@cs.uky.edu 3 National Tsinhua University ABSTRACT In this
More informationOn the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations
On the graphical display of Powell-Sabin splines: a comparison of three piecewise linear approximations Hendrik Speleers Paul Dierckx Stefan Vandewalle Report TW515, January 008 Ò Katholieke Universiteit
More information08 - Designing Approximating Curves
08 - Designing Approximating Curves Acknowledgement: Olga Sorkine-Hornung, Alexander Sorkine-Hornung, Ilya Baran Last time Interpolating curves Monomials Lagrange Hermite Different control types Polynomials
More informationPolar Embedded Catmull-Clark Subdivision Surface
Polar Embedded Catmull-Clark Subdivision Surface Anonymous submission Abstract In this paper, a new subdivision scheme with Polar embedded Catmull-Clark mesh structure is presented. In this new subdivision
More informationNon-Uniform Recursive Doo-Sabin Surfaces (NURDSes)
Non-Uniform Recursive Doo-Sabin Surfaces Zhangjin Huang 1 Guoping Wang 2 1 University of Science and Technology of China 2 Peking University, China SIAM Conference on Geometric and Physical Modeling Doo-Sabin
More informationCurve Subdivision in SE(2)
Curve Subdivision in SE(2) Jan Hakenberg, ETH Zürich 2018-07-26 Figure: A point in the special Euclidean group SE(2) consists of a position in the plane and a heading. The figure shows two rounds of cubic
More informationRobustness of Boolean operations on subdivision-surface models
Robustness of Boolean operations on subdivision-surface models Di Jiang 1, Neil Stewart 2 1 Université de Montréal, Dép t. Informatique CP6128, Succ. CentreVille, Montréal, H3C 3J7, Qc, Canada jiangdi@umontreal.ca
More informationComputer Graphics Curves and Surfaces. Matthias Teschner
Computer Graphics Curves and Surfaces Matthias Teschner Outline Introduction Polynomial curves Bézier curves Matrix notation Curve subdivision Differential curve properties Piecewise polynomial curves
More informationElement Quality Metrics for Higher-Order Bernstein Bézier Elements
Element Quality Metrics for Higher-Order Bernstein Bézier Elements Luke Engvall and John A. Evans Abstract In this note, we review the interpolation theory for curvilinear finite elements originally derived
More informationOn-Line Geometric Modeling Notes REFINEMENT
On-Line Geometric Modeling Notes REFINEMENT Kenneth I Joy Visualization and Graphics Research Group Department of Computer Science University of California, Davis Overview Bézier curves, B-spline curves
More informationCOMPUTER AIDED ENGINEERING DESIGN (BFF2612)
COMPUTER AIDED ENGINEERING DESIGN (BFF2612) BASIC MATHEMATICAL CONCEPTS IN CAED by Dr. Mohd Nizar Mhd Razali Faculty of Manufacturing Engineering mnizar@ump.edu.my COORDINATE SYSTEM y+ y+ z+ z+ x+ RIGHT
More informationTriangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property
Triangle Mesh Subdivision with Bounded Curvature and the Convex Hull Property Charles Loop cloop@microsoft.com February 1, 2001 Technical Report MSR-TR-2001-24 The masks for Loop s triangle subdivision
More informationA Subdivision-Regularization Framework for Preventing Over Fitting of Data by a Model
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 178-190 Applications and Applied Mathematics: An International Journal (AAM) A Subdivision-Regularization
More informationAs a consequence of the operation, there are new incidences between edges and triangles that did not exist in K; see Figure II.9.
II.4 Surface Simplification 37 II.4 Surface Simplification In applications it is often necessary to simplify the data or its representation. One reason is measurement noise, which we would like to eliminate,
More informationExample: Loop Scheme. Example: Loop Scheme. What makes a good scheme? recursive application leads to a smooth surface.
Example: Loop Scheme What makes a good scheme? recursive application leads to a smooth surface 200, Denis Zorin Example: Loop Scheme Refinement rule 200, Denis Zorin Example: Loop Scheme Two geometric
More informationBernstein-Bezier Splines on the Unit Sphere. Victoria Baramidze. Department of Mathematics. Western Illinois University
Bernstein-Bezier Splines on the Unit Sphere Victoria Baramidze Department of Mathematics Western Illinois University ABSTRACT I will introduce scattered data fitting problems on the sphere and discuss
More informationComputation of interpolatory splines via triadic subdivision
Computation of interpolatory splines via triadic subdivision Valery A. Zheludev and Amir Z. Averbuch Abstract We present an algorithm for computation of interpolatory splines of arbitrary order at triadic
More informationTo appear in Computer-Aided Design Revised June 18, J-splines
To appear in Computer-Aided Design Revised June 18, 2008 J-splines Jarek Rossignac School of Interactive Computing, College of Computing, Georgia Institute of Technology, Atlanta, GA http://www.gvu.gatech.edu/~jarek
More informationGeneralized barycentric coordinates
Generalized barycentric coordinates Michael S. Floater August 20, 2012 In this lecture, we review the definitions and properties of barycentric coordinates on triangles, and study generalizations to convex,
More informationGraph Adjacency Matrix Automata Joshua Abbott, Phyllis Z. Chinn, Tyler Evans, Allen J. Stewart Humboldt State University, Arcata, California
Graph Adjacency Matrix Automata Joshua Abbott, Phyllis Z. Chinn, Tyler Evans, Allen J. Stewart Humboldt State University, Arcata, California Abstract We define a graph adjacency matrix automaton (GAMA)
More informationGauss curvature. curvature lines, needle plot. vertices of progenitor solid are unit size. Construction Principle. edge. edge. vertex. Patch.
The Platonic Spheroids Jorg Peters and Leif Kobbelt Abstract. We present a gallery of simple curvature continuous surfaces that possess the topological structure of the Platonic solids. These spherelike
More information3. Lifting Scheme of Wavelet Transform
3. Lifting Scheme of Wavelet Transform 3. Introduction The Wim Sweldens 76 developed the lifting scheme for the construction of biorthogonal wavelets. The main feature of the lifting scheme is that all
More informationDeficient Quartic Spline Interpolation
International Journal of Computational Science and Mathematics. ISSN 0974-3189 Volume 3, Number 2 (2011), pp. 227-236 International Research Publication House http://www.irphouse.com Deficient Quartic
More informationJoe Warren, Scott Schaefer Rice University
Joe Warren, Scott Schaefer Rice University Polygons are a ubiquitous modeling primitive in computer graphics. Their popularity is such that special purpose graphics hardware designed to render polygons
More informationGeneralized Lane Riesenfeld algorithms
Generalized Lane Riesenfeld algorithms Thomas J. Cashman a, Kai Hormann a, Ulrich Reif b a Università della Svizzera italiana, Via Giuseppe Buffi 13, CH-6904 Lugano b TU Darmstadt, Schlossgartenstrasse
More informationModified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation
Modified Catmull-Clark Methods for Modelling, Reparameterization and Grid Generation Karl-Heinz Brakhage RWTH Aachen, 55 Aachen, Deutschland, Email: brakhage@igpm.rwth-aachen.de Abstract In this paper
More informationRemoving Polar Rendering Artifacts in Subdivision Surfaces
This is an electronic version of an article published in Journal of Graphics, GPU, and Game Tools, Volume 14, Issue 2 pp. 61-76, DOI: 10.1080/2151237X.2009.10129278. The Journal of Graphics, GPU, and Game
More informationAn introduction to interpolation and splines
An introduction to interpolation and splines Kenneth H. Carpenter, EECE KSU November 22, 1999 revised November 20, 2001, April 24, 2002, April 14, 2004 1 Introduction Suppose one wishes to draw a curve
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationInterpolation and Splines
Interpolation and Splines Anna Gryboś October 23, 27 1 Problem setting Many of physical phenomenona are described by the functions that we don t know exactly. Often we can calculate or measure the values
More informationApplications of Wavelets and Framelets
Applications of Wavelets and Framelets Bin Han Department of Mathematical and Statistical Sciences University of Alberta, Edmonton, Canada Present at 2017 International Undergraduate Summer Enrichment
More informationAlmost Curvature Continuous Fitting of B-Spline Surfaces
Journal for Geometry and Graphics Volume 2 (1998), No. 1, 33 43 Almost Curvature Continuous Fitting of B-Spline Surfaces Márta Szilvási-Nagy Department of Geometry, Mathematical Institute, Technical University
More informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationCS 450 Numerical Analysis. Chapter 7: Interpolation
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationConvergence of C 2 Deficient Quartic Spline Interpolation
Advances in Computational Sciences and Technology ISSN 0973-6107 Volume 10, Number 4 (2017) pp. 519-527 Research India Publications http://www.ripublication.com Convergence of C 2 Deficient Quartic Spline
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J. MATRIX ANAL. APPL. Vol. 26, No. 2, pp. 390 399 c 2004 Society for Industrial and Applied Mathematics HERMITIAN MATRICES, EIGENVALUE MULTIPLICITIES, AND EIGENVECTOR COMPONENTS CHARLES R. JOHNSON
More informationRemark. Jacobs University Visualization and Computer Graphics Lab : ESM4A - Numerical Methods 331
Remark Reconsidering the motivating example, we observe that the derivatives are typically not given by the problem specification. However, they can be estimated in a pre-processing step. A good estimate
More informationconsisting of compact sets. A spline subdivision scheme generates from such
Spline Subdivision Schemes for Compact Sets with Metric Averages Nira Dyn and Elza Farkhi Abstract. To dene spline subdivision schemes for general compact sets, we use the representation of spline subdivision
More informationHarmonic Spline Series Representation of Scaling Functions
Harmonic Spline Series Representation of Scaling Functions Thierry Blu and Michael Unser Biomedical Imaging Group, STI/BIO-E, BM 4.34 Swiss Federal Institute of Technology, Lausanne CH-5 Lausanne-EPFL,
More information1. Introduction. 2. Parametrization of General CCSSs. 3. One-Piece through Interpolation. 4. One-Piece through Boolean Operations
Subdivision Surface based One-Piece Representation Shuhua Lai Department of Computer Science, University of Kentucky Outline. Introduction. Parametrization of General CCSSs 3. One-Piece through Interpolation
More informationSubdivision Scheme Tuning Around Extraordinary Vertices
Subdivision Scheme Tuning Around Extraordinary Vertices Loïc Barthe Leif Kobbelt Computer Graphics Group, RWTH Aachen Ahornstrasse 55, 52074 Aachen, Germany Abstract In this paper we extend the standard
More informationCurvature of Approximating Curve Subdivision Schemes
Curvature of Approximating Curve Subdivision Schemes Kȩstutis Karčiauskas 1 and Jörg Peters 2 1 Vilnius University, 2 University of Florida Abstract. The promise of modeling by subdivision is to have simple
More informationHierarchical Grid Conversion
Hierarchical Grid Conversion Ali Mahdavi-Amiri, Erika Harrison, Faramarz Samavati Abstract Hierarchical grids appear in various applications in computer graphics such as subdivision and multiresolution
More informationAdaptive and Smooth Surface Construction by Triangular A-Patches
Adaptive and Smooth Surface Construction by Triangular A-Patches Guoliang Xu Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences, Beijing, China Abstract
More information4 Parametrization of closed curves and surfaces
4 Parametrization of closed curves and surfaces Parametrically deformable models give rise to the question of obtaining parametrical descriptions of given pixel or voxel based object contours or surfaces,
More informationFinite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 24
Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 24 So in today s class, we will look at quadrilateral elements; and we will
More informationTwo Algorithms for Approximation in Highly Complicated Planar Domains
Two Algorithms for Approximation in Highly Complicated Planar Domains Nira Dyn and Roman Kazinnik School of Mathematical Sciences, Tel-Aviv University, Tel-Aviv 69978, Israel, {niradyn,romank}@post.tau.ac.il
More informationQuasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters
Quasi-Quartic Trigonometric Bézier Curves and Surfaces with Shape Parameters Reenu Sharma Assistant Professor, Department of Mathematics, Mata Gujri Mahila Mahavidyalaya, Jabalpur, Madhya Pradesh, India
More informationOn the Dimension of the Bivariate Spline Space S 1 3( )
On the Dimension of the Bivariate Spline Space S 1 3( ) Gašper Jaklič Institute of Mathematics, Physics and Mechanics University of Ljubljana Jadranska 19, 1000 Ljubljana, Slovenia Gasper.Jaklic@fmf.uni-lj.si
More information