Bin Han Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, Alberta, Canada
|
|
- Giles Clarke
- 5 years ago
- Views:
Transcription
1 Bivariate (Two-dimensional) Wavelets Bin Han Department of Mathematical and Statistical Sciences, University of Alberta, Edmonton, Alberta, Canada Article Outline Glossary I. Definitions II. Introduction III. Bivariate Refinable Functions and Their Properties IV. The Projection Method V. Bivariate Orthonormal and Biorthogonal Wavelets VI. Bivariate Riesz Wavelets VII. Pairs of Dual Wavelet Frames VIII. Future Directions. IX. Bibliography Glossary Dilation matrix A 2 2 matrix M is called a dilation matrix if all the entries of M are integers and all the eigenvalues of M are greater than one in modulus. Isotropic dilation matrix A dilation matrix M is said to be isotropic if M is similar to a diagonal matrix and all its eigenvalues have the same modulus. Wavelet system A wavelet system is a collection of square integrable functions that are generated from a finite set of functions (which are called wavelets) by using integer shifts and dilations. I. Definitions Throughout this article, R, C, Z denote the real line, the complex plane, and the set of all integers, respectively. For 1 p, L p (R 2 ) denotes the set of all Lebesgue measurable bivariate functions f such that f p L p (R 2 ) := R 2 f(x) p dx <. In particular, the space L 2 (R 2 ) of square integrable functions is a Hilbert space under the inner product f, g := f(x)g(x) dx, f, g L 2 (R 2 ), R 2 where g(x) denotes the complex conjugate of the complex number g(x). In applications such as image processing and computer graphics, the following are commonly used isotropic dilation matrices: M 2 = [ ], Q 2 = [ ] 2 1, M 3 =, di = d 0, (1) 0 d where d is an integer with d > 1. Using a dilation matrix M, a bivariate M-wavelet system is generated by integer shifts and dilates from a finite set {ψ 1,..., ψ L } of functions in L 2 (R 2 ). More precisely, the set of 1
2 2 all the basic wavelet building blocks of an M-wavelet system generated by {ψ 1,..., ψ L } is given by X M ({ψ 1,..., ψ L }) := {ψ l,m : j Z, k Z 2, l = 1,..., L}, (2) where ψ l,m (x) := detm j/2 ψ l (M j x k), x R 2. (3) One of the main goals in wavelet analysis is to find wavelet systems X M ({ψ 1,..., ψ L }) with some desirable properties such that any two-dimensional function or signal f L 2 (R 2 ) can be sparsely and efficiently represented under the M-wavelet system X M ({ψ 1,..., ψ L }): L f = h l (f)ψ l,m, (4) l=1 j Z k Z 2 where h l : L 2(R 2 ) C are linear functionals. There are many types of wavelet systems studied in the literature. In the following, let us outline some of the most important types of wavelet systems. Orthonormal wavelets. We say that {ψ 1,..., ψ L } generates an orthonormal M-wavelet basis in L 2 (R 2 ) if the system X M ({ψ 1,..., ψ L }) is an orthonormal basis of the Hilbert space L 2 (R 2 ). That is, the linear span of elements in X M ({ψ 1,..., ψ L }) is dense in L 2 (R 2 ) and ψ l,m,,m ψl j,k = δ l l δ j j δ k k, j, j Z, k, k Z 2, l, l = 1,..., L, (5) where δ denotes the Dirac sequence such that δ 0 = 1 and δ k = 0 for all k 0. For an orthonormal wavelet basis X M ({ψ 1,..., ψ L }), the linear functional h l in (4) is given by h l (f) = f, ψl,m and the representation in (4) becomes f = L l=1 j Z k Z 2 f, ψ l,m ψl,m, f L 2(R 2 ) (6) with the series converging in L 2 (R 2 ). Riesz wavelets. We say that {ψ 1,..., ψ L } generates a Riesz M-wavelet basis in L 2 (R 2 ) if the system X M ({ψ 1,..., ψ L }) is a Riesz basis of L 2 (R 2 ). That is, the linear span of elements in X M ({ψ 1,..., ψ L }) is dense in L 2 (R 2 ) and there exist two positive constants C 1 and C 2 such that L 2 C 1 c l 2 L c l ψ l,m L C 2 c l 2 l=1 j Z k Z 2 l=1 j Z k Z 2 l=1 j Z k Z 2 L 2 (R 2 ) for all finitely supported sequences {c l } j Z,k Z 2,l=1,...,L. Clearly, a Riesz M-wavelet generalizes an orthonormal M-wavelet by relaxing the orthogonality requirement in (5). For a Riesz basis X M ({ψ 1,..., ψ L }), it is well-known that there exists a dual Riesz basis { ψ l, : j Z, k Z 2, l = 1,..., L} of elements in L 2 (R 2 ) (this set is not necessarily generated by integer shifts and dilates from some finite set of functions) such that (5) still holds after replacing ψ l,m j,k by ψ l,j,k. For a Riesz wavelet basis, the linear functional in (4) becomes h l (f) = f, ψ l,. In fact, ψl, := F 1 (ψ l,m ), where F : L 2(R 2 ) L 2 (R 2 ) is defined to be F(f) := L l=1 j Z k Z 2 f, ψ l,m ψl,m, f L 2(R 2 ). (7)
3 3 Wavelet frames. A further generalization of a Riesz wavelet is a wavelet frame. We say that {ψ 1,..., ψ L } generates an M-wavelet frame in L 2 (R 2 ) if the system X M ({ψ 1,..., ψ L }) is a frame of L 2 (R 2 ). That is, there exist two positive constants C 1 and C 2 such that C 1 f 2 L 2 (R 2 ) L l=1 j Z k Z 2 f, ψ l,m 2 C 2 f 2 L 2 (R 2 ), f L 2(R 2 ). (8) It is not difficult to check that a Riesz M-wavelet is an M-wavelet frame. Then (8) guarantees that the frame operator F in (7) is a bounded and invertible linear operator. For an M-wavelet frame, the linear functional in (4) can be chosen to be h l (f) = f, F 1 (ψ l,m ) ; however, such functionals may not be unique. The fundamental difference between a Riesz wavelet and a wavelet frame lies in that for any given function f L 2 (R 2 ), the representation in (4) is unique under a Riesz wavelet while it may not be unique under a wavelet frame. The representation in (4) with the choice h l (f) = f, F 1 (ψ l,m ) is called the canonical representation of a wavelet frame. In other words, {F 1 (ψ l,m ) : j Z, k Z2, l = 1,..., L} is called the canonical dual frame of the given wavelet frame X M ({ψ 1,..., ψ L }). Biorthogonal wavelets. We say that ({ψ 1,..., ψ L }, { ψ 1,..., ψ L }) generates a pair of biorthogonal M-wavelet bases in L 2 (R 2 ) if each of X M ({ψ 1,..., ψ L }) and X M ({ ψ 1,..., ψ L }) is a Riesz basis of L 2 (R 2 ) and (5) still holds after replacing ψ l,m j,k by ψ l,m j,k. In other words, the dual Riesz basis of X M ({ψ 1,..., ψ L }) has the wavelet structure and is given by X M ({ ψ 1,..., ψ L }). For a biorthogonal wavelet, the wavelet representation in (4) becomes f = L l=1 f, j Z k Z 2 ψ l,m ψl,m, f L 2(R 2 ). (9) Obviously, {ψ 1,..., ψ L } generates an orthonormal M-wavelet basis in L 2 (R 2 ) if and only if ({ψ 1,..., ψ L }, {ψ 1,..., ψ L }) generates a pair of biorthogonal M-wavelet bases in L 2 (R 2 ). Dual wavelet frames. Similarly, we have the notion of a pair of dual wavelet frames. We say that ({ψ 1,..., ψ L }, { ψ 1,..., ψ L }) generates a pair of dual M-wavelet frames in L 2 (R 2 ) if each of X M ({ψ 1,..., ψ L }) and X M ({ ψ 1,..., ψ L }) is an M-wavelet frame in L 2 (R 2 ) and f, g = L l=1 f, j Z k Z 2 ψ l,m ψl,m, g, f, g L 2(R 2 ). It follows from the above identity that (9) still holds for a pair of dual M-wavelet frames. We say that {ψ 1,..., ψ L } generates a tight M-wavelet frame if ({ψ 1,..., ψ L }, {ψ 1,..., ψ L }) generates a pair of dual M-wavelet frames in L 2 (R 2 ). For a tight wavelet frame {ψ 1,..., ψ L }, the wavelet representation in (6) still holds. So, a tight wavelet frame is a generalization of an orthonormal wavelet basis. II. Introduction Bivariate (two-dimensional) wavelets are of interest in representing and processing two dimensional data such as images and surfaces. In this article we shall discuss some basic background and results on bivariate wavelets. We denote Π J the set of all bivariate polynomials of total degree at most J. For compactly supported functions ψ 1,..., ψ L, we say that {ψ 1,..., ψ L } has J vanishing moments if P, ψ l = 0 for all l = 1,..., L and all polynomials P Π J 1.
4 4 The advantages of wavelet representations largely lie in the following aspects: (1) There is a fast wavelet transform (FWT) for computing the wavelet coefficients h l (f) in the wavelet representation (4). (2) The wavelet representation has good time and frequency localization. Roughly speaking, the basic building blocks ψ l,m have good time localization and smoothness. (3) The wavelet representation is sparse. For a smooth function f, most wavelet coefficients are negligible. Generally, the wavelet coefficient h l (f) only depends on the information of f in a small neighborhood of the support of ψ l,m (more precisely, the support of ψl,m if h l l,m (f) = f, ψ ). If f in a small neighborhood of the support of ψ l,m is smooth or behaves like a polynomial to certain degree, then by the vanishing moments of the dual wavelet functions, h l (f) is often negligible. (4) For a variety of function spaces B such as Sobolev and Besov spaces, its norm f B, f B, is equivalent to certain weighted sequence norm of its wavelet coefficients {h l (f)} j Z,k Z 2,l=1,...,L. For more details on advantages and applications of wavelets, see [1, 4, 5, 6, 7, 10]. For a dilation matrix M = di 2, the easiest way for obtaining a bivariate wavelet is to use the tensor product method. In other words, all the functions in the generator set {ψ 1,..., ψ L } L 2 (R 2 ) take the form ψ l (x, y) = f l (x)g l (y), x, y R, where f l and g l are some univariate functions in L 2 (R). However, tensor product (also called separable) bivariate wavelets give preference to the horizontal and vertical directions which may not be desirable in applications such as image processing ([7, 10, 17, 52]). Also, for a nondiagonal dilation matrix, it is difficult or impossible to use the tensor product method to obtain bivariate wavelets. Therefore, nonseparable bivariate wavelets themselves are of importance and interest in both theory and application. In this article, we shall address several aspects of bivariate wavelets with a general two-dimensional dilation matrix. III. Bivariate Refinable Functions and Their Properties In order to have a fast wavelet transform to compute the wavelet coefficients in (4), the generators ψ 1,..., ψ L in a wavelet system are generally obtained from a refinable function φ via a multiresolution analysis ([4, 7, 10, 21]). Let M be a 2 2 dilation matrix. For a function φ in L 2 (R 2 ), we say that φ is M-refinable if it satisfies the following refinement equation φ(x) = detm k Z 2 a k φ(mx k), a.e. x R 2, (10) where a : Z 2 C is a finitely supported sequence on Z 2 satisfying k Z a 2 k = 1. Such a sequence a is often called a mask in wavelet analysis and computer graphics, or a low-pass filter in signal processing. Using the Fourier transform ˆφ of φ L 2 (R 2 ) L 1 (R 2 ) and the Fourier series â of a, which are defined to be ˆφ(ξ) := φ(x)e ix ξ dx and â(ξ) := a k e ik ξ, ξ R 2, (11) R 2 k Z 2 the refinement equation (10) can be equivalently rewritten as ˆφ(M T ξ) = â(ξ) ˆφ(ξ), ξ R 2, (12)
5 5 where M T denotes the transpose of the dilation matrix M. Since â(0) = 1 and â is a trigonometric polynomial, one can define a function ˆφ by ˆφ(ξ) := â((m T ) j ξ), ξ R 2, (13) j=1 with the series converging uniformly on any compact set of R 2. It is known ([3, 7]) that φ is a compactly supported tempered distribution and clearly φ satisfies the refinement equation in (12) with mask a. We call φ the standard refinable function associated with mask a and dilation M. Now the generators ψ 1,..., ψ L of a wavelet system are often obtained from the refinable function φ via ψ l (M T ξ) := âl (ξ) ˆφ(ξ), ξ R 2, l = 1,..., L, for some 2π-periodic trigonometric polynomials âl. Such sequences a l are called wavelet masks or high-pass filters. Except for very few special masks a such as masks for box spline refinable functions (see [8] for box splines), most refinable functions φ, obtained in (13) from mask a and dilation M, do not have explicit analytic expressions and a so-called cascade algorithm or subdivision scheme is used to approximate the refinable function φ. Let B be a Banach space of bivariate functions. Starting with a suitable initial function f B, we iteratively compute a sequence {Q n a,m f} n=0 of functions, where Q a,m f := detm k Z 2 a k f(m k). (14) If the sequence {Q n a,m f} n=0 of functions converges to some f B in the Banach space B, then we have Q a,m f = f. That is, as a fixed point of the cascade operator Q a,m, f is a solution to the refinement equation (10). A cascade algorithm plays an important role in the study of refinable functions in wavelet analysis and of subdivision surfaces in computer graphics. For more details on cascade algorithms and subdivision schemes, see [1, 3, 9, 24, 25, 33, 38, 55] and numerous references therein. One of the most important properties of a refinable function φ is its smoothness, which is measured by its L p smoothness exponent ν p (φ) and is defined to be ν p (φ) := sup{ν : φ W ν p (R 2 )}, 1 p, (15) where Wp ν (R 2 ) denotes the fractional Sobolev space of order ν > 0 in L p (R 2 ). The notion of sum rules of a mask is closely related to the vanishing moments of a wavelet system ([7, 46]). For a bivariate mask a, we say that a satisfies the sum rules of order J with the dilation matrix M if a j+mk P (j + Mk) = a Mk P (Mk), P Π J 1. (16) k Z 2 k Z 2 Or equivalently, µ 1 1 µ 2 2 â(2πγ) = 0 for all µ 1 + µ 2 < J and γ [(M T ) 1 Z 2 ]\Z 2, where 1 and 2 denote the partial derivatives in the first and second coordinate, respectively. Throughout the article, we denote sr(a, M) the highest order of sum rules satisfied by the mask a with the dilation matrix M. To investigate various properties of refinable functions, next we introduce a quantity ν p (a, M) in wavelet analysis (see [33]). For a bivariate mask a and a dilation matrix M,
6 6 we denote ν p (a, M) := log ρ(m) [ detm 1 1/p ρ(a, M, p)], 1 p, (17) where ρ(m) denotes the spectral radius of M and { } ρ(a, M, p) := max lim sup a n,(µ 1,µ 2 ) 1/n l p (Z n 2 ) : µ 1 + µ 2 = sr(a, M), µ 1, µ 2 N {0}, where the sequence a n,(µ 1,µ 2 ) is defined by its Fourier series as follows: for ξ = (ξ 1, ξ 2 ) R 2, a n,(µ 1,µ 2) (ξ) := (1 e iξ 1 ) µ 1 (1 e iξ 2 ) µ2â((m T ) n 1 ξ)â((m T ) n 2 ξ) â(m T ξ)â(ξ). It is known that ([33,38]) ν p (a, M) ν q (a, M) ν p (a, M) + (1/q 1/p) log ρ(m) detm for all 1 p q. Parts of the following result essentially appeared in various forms in many papers in the literature. For the study of refinable functions and cascade algorithms, see [3, 7, 9, 16, 19, 21, 25, 27, 31, 33, 47, 50, 56] and references therein. The following is from [33]. Theorem 1. Let M be a 2 2 isotropic dilation matrix and a be a finitely supported bivariate mask. Let φ denote the standard refinable function associated with mask a and dilation M. Then ν p (φ) ν p (a, M) for all 1 p. If the shifts of φ are stable, that is, {φ( k) : k Z 2 } is a Riesz system in L p (R 2 ), then ν p (φ) = ν p (a, M). Moreover, for every nonnegative integer J, the following statements are equivalent (1) For every compactly supported function f Wp J (R 2 ) (if p =, we require f C J (R 2 )) such that ˆf(0) = 1 and µ 1 1 µ 2 2 ˆf(2πk) = 0 for all µ 1 + µ 2 < J and k Z 2 \{0}, the cascade sequence {Q n a,m f} n=0 converges in the Sobolev space Wp J (R 2 ) (in fact, the limit function is φ). (2) ν p (a, M) > J. Symmetry of a refinable function is also another important property of a wavelet system in applications. For example, symmetry is desirable in order to handle the boundary of an image or to improve the visual quality of a surface ([1, 7, 9, 10, 34, 52]). Let M be a 2 2 dilation matrix and G be a finite set of 2 2 integer matrices. We say that G is a symmetry group with respect to M ([29, 31]) if G is a group under matrix multiplication and MEM 1 G for all E G. In dimension two, two commonly used symmetry groups are D 4 and D 6 : and D 6 := { ± D 4 := 1 0, ± 0 1 { ± [ ] 0 1, ± , ±, ± , ± , ± , ± 1 0 [ 0 ]} , ± 0 1 } For symmetric refinable functions, we have the following result ([34, Proposition 2.1]): Proposition 2. Let M be a 2 2 dilation matrix and G be a symmetry group with respect to M. Let a be a bivariate mask and φ be the standard refinable function associated with mask a and dilation M. Then a is G-symmetric with center c a : a E(k ca )+c a = a k for all k Z 2 and E G, if and only if, φ is G-symmetric with center c: φ(e( c) + c) = φ, E G with c := (M I 2 ) 1 c a. In the following, let us present some examples of bivariate refinable functions.
7 7 Example 3. Let M = 2I 2 and â(ξ 1, ξ 2 ) := cos 4 (ξ 1 /2) cos 4 (ξ 2 /2). Then a is a tensor product mask with sr(a, 2I 2 ) = 4 and a is D 4 -symmetric with center 0. The associated standard refinable function φ is the tensor product spline of order 4. The Catmull-Clark subdivision scheme for quadrilateral meshes in computer graphics is based on this mask a ([1]). Example 4. Let M = 2I 2 and â(ξ 1, ξ 2 ) := cos 2 (ξ 1 /2) cos 2 (ξ 2 /2) cos 2 (ξ 1 /2 + ξ 2 /2). Then sr(a, 2I 2 ) = 4 and a is D 6 -symmetric with center 0. The associated refinable function φ is the convolution of the three direction box spline with itself ([8]). The Loop subdivision scheme for triangular meshes in computer graphics is based on this mask a ([1]). Another important property of a refinable function is interpolation. We say that a bivariate function φ is interpolating if φ is continuous and φ(k) = δ k for all k Z 2. Theorem 5. Let M be a 2 2 dilation matrix and a be a finitely supported bivariate mask. Let φ denote the standard refinable function associated with mask a and dilation M. Then φ is an interpolating function if and only if ν (a, M) > 0 and a is an interpolatory mask with dilation M, that is, a 0 = detm 1 and a Mk = 0 for all k Z 2 \{0}. Example 6. Let M = 2I 2 and â(ξ 1, ξ 2 ) := cos(ξ 1 /2) cos(ξ 2 /2) cos(ξ 1 /2 + ξ 2 /2) [ cos(ξ 1 ) + 2 cos(ξ 2 ) + 2 cos(ξ 1 + ξ 2 ) cos(2ξ 1 + ξ 2 ) cos(ξ 1 + 2ξ 2 ) cos(ξ 1 ξ 2 )]/4. Then sr(a, 2I 2 ) = 4, a is D 6 -symmetric with center 0 and a is an interpolatory mask with dilation 2I 2. Since ν 2 (a, 2I 2 ) , so, ν (a, 2I 2 ) ν 2 (a, 2I 2 ) > 0. By Theorem 5, the associated standard refinable function φ is interpolating. The butterfly interpolatory subdivision scheme for triangular meshes is based on this mask a (see [25]). For more details on bivariate interpolatory masks and interpolating refinable functions, see [3, 9, 24, 25, 27, 28, 29, 33, 39, 40, 41, 60]. IV. The Projection Method In applications, one is interested in analyzing some optimal properties of multivariate wavelets. The projection method is useful for this purpose. Let r and s be two positive integers with r s. Let P be an r s real-valued matrix. For a compactly supported function φ and a finitely supported mask a in dimension s, we define the projected function P φ and projected mask P a in dimension r by (18) P φ(ξ) := ˆφ(P t T ξ) and P a(ξ) := â(p T ξ + 2πε j ), ξ R r, (19) where ˆφ and â are understood to be continuous and {ε 1,..., ε t } is a complete set of representatives of the distinct cosets of [P T R r ]/Z s. If P is an integer matrix, then P a(ξ) = â(p T ξ). Now we have the following result on projected refinable functions ([29, 30, 36, 37]). Theorem 7. Let N be an s s dilation matrix and M be an r r dilation matrix. Let P be an r s integer matrix such that P N = MP and P Z s = Z r. Let â be a 2πperiodic trigonometric polynomials in s-variables with â(0) = 1 and φ the standard N- refinable function associated with mask a. Then sr(a, N) sr(p a, M) and P φ(m T ξ) = P a(ξ) P φ(ξ). That is, P φ is M-refinable with mask P a. Moreover, for all 1 p, ν p (φ) ν p (P φ) and det M 1 1/p ρ(p a, M, p) det N 1 1/p ρ(a, N, p). j=1
8 8 If we further assume that ρ(m) = ρ(n), then ν p (a, N) ν p (P a, M). As pointed out in [36, 37], the projection method is closely related to box splines. For a given r s (direction) integer matrix Ξ of rank r with r s, the Fourier transform of its associated box spline M Ξ and its mask a Ξ are given by (see [8]) M Ξ (ξ) := k Ξ 1 e ik ξ ik ξ and â Ξ (ξ) := 1 + e ik ξ, ξ R r, (20) 2 k Ξ where k Ξ means that k is a column vector of Ξ and k goes through all the columns of Ξ once and only once. Let χ [0,1] s denote the characteristic function of the unit cube [0, 1] s. From (20), it is evident that the box spline M Ξ is just the projected function Ξχ [0,1] s, since Ξχ [0,1] s = M Ξ by Ξχ [0,1] s(ξ) = χ [0,1] s(ξ T ξ), ξ R r. Note that M Ξ is 2I r -refinable with the mask a Ξ, since M Ξ (2ξ) = â Ξ (ξ) M Ξ (ξ). As an application of the projection method in Theorem 7, we have ([27, Theorem 3.5]): Corollary 8. Let M = 2I s and a be an interpolatory mask with dilation M such that a is supported inside [ 3, 3] s. Then ν (a, 2I s ) 2 and therefore, φ C 2 (R s ), where φ is the standard refinable function associated with mask a and dilation 2I s. We use proof by contradiction. Suppose ν (a, 2I s ) > 2. Let P = [1, 0,..., 0] be a 1 s matrix. Then we must have sr(a, 2I s ) 3, which, combining with other assumptions on a, will force P a(ξ) = 1/2+9/16 cos(ξ) 1/16 cos(3ξ). Since ν (P a, 2) = 2, by Theorem 7, we must have ν (a, 2I s ) ν (P a, 2) = 2. A contradiction. So, ν (a, 2I s ) 2. In particular, the refinable function in the butterfly subdivision scheme in Example 6 is not C 2. The projection method can be used to construct interpolatory masks painlessly ([29]). Theorem 9. Let M be an r r dilation matrix. Then there is an r r integer matrix H such that MZ r = HZ r and H r = detm I r. Let P := detm 1 H. Then for any (tensor product) interpolatory mask a with dilation detm I r, the projected mask P a is an interpolatory mask with dilation M and sr(p a, M) sr(a, detm I r ). Example 10. For M = M 2 or M = Q 2, we can take H := M 2 in Theorem 9. For more details on the projection method, see [27, 29, 30, 34, 36, 37]. V. Bivariate Orthonormal and Biorthogonal Wavelets In this section, we shall discuss the analysis and construction of bivariate orthonormal and biorthogonal wavelets. For analysis of biorthogonal wavelets, the following result is well-known (see [7, 17, 18, 26, 28, 33, 48, 53, 56] and references therein): Theorem 11. Let M be a 2 2 dilation matrix and a, ã be two finitely supported bivariate masks. Let φ and φ be the standard M-refinable functions associated with masks a and ã, respectively. Then φ, φ L 2 (R 2 ) and satisfy the biorthogonality relation φ, φ( k) = δ k, k Z 2, (21) if and only if, ν 2 (a, M) > 0, ν 2 (ã, M) > 0 and (a, ã) is a pair of dual masks: â(ξ + 2πγ)ˆã(ξ + 2πγ) = 1, ξ R 2, (22) γ Γ M T
9 where Γ M T is a complete set of representatives of distinct cosets of [(M T ) 1 Z 2 ]/Z 2 with 0 Γ M T. Moreover, if (21) holds and there exist 2π-periodic trigonometric polynomials â 1,..., âm 1, ã 1,..., ã m 1 with m := detm such that where for ξ R 2, M [â, (ξ) := â 1,...,âm 1 ] with {γ 0,..., γ m 1 } := Γ M T ψ l (M T ξ) := âl (ξ) ˆφ(ξ) M [â, â 1,...,âm 1 ] (ξ)m [ˆã, ã 1,..., ã m 1 ] (ξ)t = I m, ξ R 2, (23) â(ξ + 2πγ 0 ) â 1 (ξ + 2πγ 0 ) â m 1 (ξ + 2πγ 0 ) â(ξ + 2πγ 1 ) â 1 (ξ + 2πγ 1 ) â m 1 (ξ + 2πγ 1 ) â(ξ + 2πγ m 1 ) â 1 (ξ + 2πγ m 1 ) â m 1 (ξ + 2πγ m 1 ) and γ 0 := 0. Define ψ 1,..., ψ m 1, ψ 1,..., ψ m 1 by and 9 (24) ψ l (M T ξ) := ã l (ξ) ˆ φ(ξ), l = 1,..., m 1. (25) Then ({ψ 1,..., ψ m 1 }, { ψ 1,..., ψ m 1 }) generates a pair of biorthogonal M-wavelet bases in L 2 (R 2 ). As a direct consequence of Theorem 11, for orthonormal wavelets, we have Corollary 12. Let M be a 2 2 dilation matrix and a be a finitely supported bivariate mask. Let φ be the standard refinable function associated with mask a and dilation M. Then φ has orthonormal shifts, that is, φ, φ( k) = δ k for all k Z 2, if and only if, ν 2 (a, M) > 0 and a is an orthogonal mask (that is, (22) is satisfied with ˆã being replaced by â). If in addition there exist 2π-periodic trigonometric polynomials â1,..., âm 1 with m := detm such that M [â, â 1,...,âm 1 ] (ξ)m [â,â1,...,âm 1 ](ξ) T = I m. Then {ψ 1,..., ψ m 1 }, defined in (25), generates an orthonormal M-wavelet basis in L 2 (R 2 ). The masks a and ã are called low-pass filters and the masks a 1,..., a m 1, ã 1,..., ã m 1 are called high-pass filters in the engineering literature. The equation in (23) is called the matrix extension problem in the wavelet literature (see [4, 7, 10, 48]). Given a pair of dual masks a and ã, though the existence of â1,..., âm 1, ã 1,..., ã m 1 (without symmetry) in the matrix extension problem in (23) is guaranteed by the Quillen-Suslin theorem, it is far from a trivial task to construct them (in particular, with symmetry) algorithmically. For the orthogonal case, the general theory of the matrix extension problem in (23) with ã = a and ã l = a l, l = 1,..., m 1, remains unanswered. However, when detm = 2, the matrix extension problem is trivial. In fact, letting γ Γ M T \{0} (that is, Γ M T = {0, γ}, one defines â1 (ξ) := e iη ξˆã(ξ + 2πγ) and ã 1 (ξ) := e iη ξ â(ξ + 2πγ), where η Z 2 satisfies η γ = 1/2. It is not an easy task to construct nonseparable orthogonal masks with desirable properties for a general dilation matrix. For example, for the dilation matrix Q 2 in (1), it is not known so far in the literature ([17]) whether there is a compactly supported C 1 Q 2-refinable function with orthonormal shifts. See [11, 12, 17, 27, 29, 32, 34, 52, 53, 56, 59] for more details on construction of orthogonal masks and orthonormal refinable functions. However, for any dilation matrix, separable orthogonal masks with arbitrarily high orders of sum rules can be easily obtained via the projection method. Let M be a 2 2 dilation matrix. Then M = Ediag(d 1, d 2 )F for some integer matrices E and F such that dete = detf = 1 and d 1, d 2 N. Let a be a tensor product orthogonal mask with the
10 10 diagonal dilation matrix diag(d 1, d 2 ) satisfying the sum rules of any preassigned order J. Then the projected mask Ea (that is, (Ea) k := a E 1 k, k Z 2 ) is an orthogonal mask with dilation M and sr(ea, M) J. See [29, Corollary 3.4] for more detail. Using the high-pass filters in the matrix extension problem, for a given mask a, dual masks ã (without any guaranteed order of sum rules of the dual masks) of a can be obtained by the lifting scheme in [64] and the stable completion method in [13]. Special families of dual masks with sum rules can be obtained by the convolution method in [26, Proposition 3.7] and [45], but the constructed dual masks generally have longer supports with respect to their orders of sum rules. Another method, which we shall cite here, for constructing all finitely supported dual masks with any preassigned orders of sum rules is the CBC (coset by coset) algorithm proposed in [14, 27, 28]. For µ = (µ 1, µ 2 ) and ν = (ν 1, ν 2 ), we say ν µ if ν 1 µ 1 and ν 2 µ 2. Also we denote µ = µ 1 + µ 2, µ! := µ 1!µ 2! and x µ := x µ 1 1 x µ 2 2 for x = (x 1, x 2 ). We denote Ω M a complete set of representatives of distinct cosets of Z 2 /[MZ 2 ] with 0 Ω M. The following result is from [27]. Theorem 13. (CBC Algorithm) Let M be a 2 2 dilation matrix and a be an interpolatory mask with dilation M. Let J be any preassigned positive integer. (1) Compute the quantities h a µ, µ < J, from the mask a by the recursive formula: h a µ := δ µ ( 1) µ ν µ! ν!(µ ν)! ha ν a k k µ ν. 0 ν<µ k Z 2 (2) For every nonzero coset γ Ω M \{0} and pre-selected subsets E γ Z 2, construct ã on the coset γ + MZ 2, with ã γ+mk = 0 for k Z 2 \E γ, satisfying the equation: k E γ ã γ+mk (γ + Mk) µ = detm 1 h a µ, µ < J. (3) Construct ã at the zero coset MZ 2 by ã Mj := detm 1 δ k a Mk Mj+γ ã γ+mk, j Z 2. k E γ γ Ω M \{0} Then mask ã is a dual mask of a and satisfies the sum rules of order J with dilation M. For a given interpolatory mask a and an arbitrary integer J, in fact any finitely supported dual mask ã of a, with ã satisfying the sum rules of order J, must come from the above CBC algorithm. The CBC algorithm also works well with symmetry. When a is symmetric, a dual mask ã with symmetry and sum rules can be easily obtained via the CBC algorithm by appropriately selecting the sets E γ. A general CBC algorithm is given in [28, Page 33 and Theorem 3.4] for constructing finitely supported dual masks ã with arbitrarily high orders of sum rules as long as the primal mask a possesses at least one finitely supported dual mask which may not have any order of sum rules. See [14, 27, 28] for more details on the CBC algorithm. The projection method can be used to obtain some optimal properties of bivariate orthonormal and biorthogonal wavelets ([27, 30, 36, 37]). Theorem 14. Let M = di 2 and (a, ã) be a pair of dual masks with dilation M. Let P = [1, 0]. Suppose that the mask a is projectable, that is, â(ξ, γπ/d) = 0 for all ξ R and γ = 1,..., d 1. Then (P a, P ã) is also a pair of dual masks with dilation d. If (21) holds
11 11 for φ and φ, where φ and φ are standard M-refinable functions with masks a and ã, then (21) still holds with φ and φ being replaced by P φ and P φ. Moreover, ν p (φ) ν p (P φ), ν p ( φ) ν p (P φ), ν p (a, di 2 ) ν p (P a, d) and ν p (ã, di 2 ) ν p (P ã, d) for all 1 p. Let us present a simple example to illustrate the CBC algorithm and the projection method on construction and analysis of bivariate biorthogonal wavelets. Example 15. Let M = 2I 2 and a be the interpolatory mask given in (18). Note that the mask a is projectable. Let P = [1, 0]. Then P a(ξ) = (2 cos(ξ)) cos 4 (ξ/2). By the CBC algorithm, we see that the shortest symmetric dual mask b of P a such that b satisfies the sum rules of order at least one is supported on [ 4, 4] and is uniquely given by ˆb(ξ) := cos 2 (ξ/2)[14 5 cos(ξ) 2 cos(2ξ) + cos(3ξ)]/8. By calculation, we have ν 2 (b, 2) Take J = 2 since sr(b, 2) = 2. By the CBC algorithm, any D 6 - symmetric dual mask ã of a, such that ã is supported inside [ 4, 4] 2 and sr(ã, 2I 2 ) = 2, must take the form t t 2 0 t t 3 t 3 t 3 t t 2 t 3 t 1 t 4 t 1 t 3 1 2t t 3 t t t 1 t 4 t 3 0 t 2 t 3 t t t 1 + 6t t 1 t 1 t 3 t 2 0 t 3 t t t 1 t 4 t t 2 t 3 t 1 t 4 t 1 t 3 1 2t t 3 t 3 t 3 t t t 2 0 t with t 3 := t 1 + 4t 2 and t 4 := 30 6t 1 4t 2, where ã (0,0) = 10 6t 1 + 6t 2. By Theorem 14 and P ã = b, we have ν 2 (ã, 2I 2 ) ν 2 (b, 2) Take t 1 = 17/4 and t 2 = 1/2, by calculation, we have ν 2 (ã, 2I 2 ) So, the conditions in Theorem 11 are satisfied and ã is near best in terms of the L 2 smoothness with respect to the support of the dual mask. VI. Bivariate Riesz Wavelets Wavelet-based numerical algorithms have been successfully used in numerical solutions to partial differential equations and integral equations ([5, 6]). The wavelets used in such wavelet-based methods are Riesz wavelets. For analysis on Riesz wavelets, see [16, 35, 42, 43, 44, 49, 58]. For applications based on Riesz wavelets, see [5, 6, 44, 51]. Theorem 16. Let M be a 2 2 dilation matrix with m := detm and a be a finitely supported bivariate mask. Let â1,..., âm 1 be some 2π-periodic trigonometric polynomials such that det M [â, â 1,...,âm 1 ] (ξ) 0 for all ξ R2 and â1 (0) = = âm 1 (0) = 0. Let φ be the standard refinable function associated with mask a and dilation M. If ν 2 (a, M) > 0 and ν 2 (ã, M) > 0, where ˆã is the (1, 1)-entry of the matrix [M [â, â 1,...,âm 1 ] (ξ)t ] 1, then {ψ 1,..., ψ m 1 }, which are defined in (25), generates a Riesz M-wavelet basis in L 2 (R 2 ). Example 17. Let M = 2I 2 and â(ξ 1, ξ 2 ) := cos 2 (ξ 1 /2) cos 2 (ξ 2 /2) cos 2 (ξ 1 /2 + ξ 2 /2) in Example 4. Then sr(a, 2I 2 ) = 4 and a is D 6 -symmetric with center 0. Define ([44, 61]) â 1 (ξ 1, ξ 2 ) := e i(ξ 1+ξ 2)â(ξ 1 + π, ξ 2 ), â 3 (ξ 1, ξ 2 ) := e iξ1â(ξ 1 + π, ξ 2 + π). â 2 (ξ 1, ξ 2 ) := e iξ2â(ξ 1, ξ 2 + π),
12 12 Then all the conditions in Theorem 16 are satisfied and {ψ 1, ψ 2, ψ 3 } generates a Riesz 2I 2 - wavelet basis in L 2 (R 2 ) ([44]). This Riesz wavelet derived from the Loop scheme has been used in [51] for mesh compression in computer graphics with impressive performance. VII. Pairs of Dual Wavelet Frames In this section, we mention a method for constructing bivariate dual wavelet frames. The following Oblique Extension Principle (OEP) has been proposed in [23] (and independently in [15]. Also see [22]) for constructing pairs of dual wavelet frames. Theorem 18. Let M be a 2 2 dilation matrix and a, ã be two finitely supported masks. Let φ, φ be the standard M-refinable functions with masks a and ã, respectively, such that φ, φ L 2 (R 2 ). If there exist 2π-periodic trigonometric polynomials Θ, â1,..., âl, ã 1,..., ã L such that Θ(0) = 1, â1 (0) = = âl (0) = ã 1 (0) = = ã L (0) = 0, and M [Θ(M T )â,â1,...,âl ] (ξ)m [ˆã, ã 1,..., ã L ] (ξ)t = diag(θ(ξ + γ 0 ),..., Θ(ξ + γ m 1 )), (26) where {γ 0,..., γ m 1 } = Γ M T with γ 0 = 0 in Theorem 11. Then ({ψ 1,..., ψ L }, { ψ 1,..., ψ L }), defined in (25), generates a pair of dual M-wavelet frames in L 2 (R 2 ). For dimension one, many interesting tight wavelet frames and pairs of dual wavelet frames in L 2 (R) have been constructed via the OEP method in the literature, for more details, see [15, 22, 23, 26, 32, 54, 62, 63] and references therein. The application of the OEP in high dimensions is much more difficult, mainly due to the matrix extension problem in (26). The projection method can also be used to obtain pairs of dual wavelet frames ([36, 37]). Theorem 19. Let M be an r r dilation matrix and N be an s s dilation matrix with r s. Let P be an r s integer matrix of rank r such that MP = P N and P T (Z r \[M T Z r ]) Z s \[N T Z s ]. Let ψ 1,..., ψ L, ψ 1,..., ψ L be compactly supported functions in L 2 (R s ) such that ν 2 (ψ l ) > 0 and ν 2 ( ψ l ) > 0 l = 1,..., L. If ({ψ 1,..., ψ L }, { ψ 1,..., ψ L }) generates a pair of dual N-wavelet frames in L 2 (R s ), then ({P ψ 1,..., P ψ L }, {P ψ 1,..., P ψ L }) generates a pair of dual M-wavelet frames in L 2 (R r ). Example 20. Let Ξ be an r s (direction) integer matrix such that Ξ T (Z r \[2Z r ]) Z s \[2Z 2 ]. Let M = 2I r and N = 2I s. Let {ψ 1,..., ψ 2s 1 } be the generators of the tensor product Haar orthonormal wavelet in dimension s, derived from the Haar orthonormal refinable function φ := χ [0,1] s. Then by Theorem 19, {Ξψ 1,..., Ξψ 2s 1 } generates a tight 2I r -wavelet frame in L 2 (R r ) and all the projected wavelet functions are derived from the refinable box spline function Ξφ = M Ξ. VIII. Future Directions. There are still many challenging problems on bivariate wavelets. In the following, we only mention a few here. (1) For any 2 2 dilation matrix M, e.g., M = Q 2 in (1), can one always construct a family of MRA compactly supported orthonormal M-wavelet bases with arbitrarily high smoothness (and with symmetry if detm > 2)? (2) The matrix extension problem in (23) for orthogonal masks. That is, for a given orthogonal mask a with dilation M, find finitely supported high-pass filters a 1,..., a m 1 (with symmetry if possible) such that (23) holds with ã = a and ã l = a l.
13 13 (3) The matrix extension problem in (23) for a given pair of dual masks with symmetry. That is, for a given pair of symmetric dual masks a and ã with dilation M, find finitely supported symmetric high-pass filters a 1,..., a m 1, ã 1,..., ã m 1 such that (23) is satisfied. (4) Directional bivariate wavelets. In order to handle edges of different orientations in images, directional wavelets are of interest in applications. See [20, 57] and many references on this topic. IX. Bibliography Books and Reviews 1. Subdivision for Modeling and Animation, course notes by T. DeRose, D. R. Forsey, L. Kobbelt, M. Lounsbery, J. Peters, P. Schröder, and D. Zorin, (1998) 2. C. Cabrelli, C. Heil and U. Molter, Self-similarity and multiwavelets in higher dimensions, Mem. Amer. Math. Soc. 170 (2004), no A. S. Cavaretta, W. Dahmen and C. A. Micchelli, Stationary subdivision, Mem. Amer. Math. Soc. 93 (1991), no C. K. Chui, An introduction to wavelets. Academic Press, Inc., Boston, MA, (1992). 5. A. Cohen, Numerical analysis of wavelet methods, North-Holland Publishing Co., Amsterdam, (2003). 6. W. Dahmen, Wavelet and multiscale methods for operator equations, Acta Numer. 6 (1997), I. Daubechies, Ten Lectures on Wavelets, CBMS-NSF Series, SIAM, Philadelphia, C. de Boor, K. H ollig and S. D. Riemenschneider, Box splines. Springer-Verlag, New York, (1993). 9. N. Dyn and D. Levin, Subdivision schemes in geometric modeling, Acta Numerica, 11 (2002), S. Mallat, A wavelet tour of signal processing, Academic Press, Inc., San Diego, CA, (1998). Primary Literature 11. A. Ayache, Some methods for constructing nonseparable, orthonormal, compactly supported wavelet bases, Appl. Comput. Harmon. Anal. 10 (2001), E. Belogay and Y. Wang, Arbitrarily smooth orthogonal nonseparable wavelets in R 2, SIAM J. Math. Anal. 30 (1999), J. M. Carnicer, W. Dahmen and J. M. Peña, Local decomposition of refinable spaces and wavelets, Appl. Comput. Harmon. Anal. 3 (1996), D. R. Chen, B. Han and S. D. Riemenschneider, Construction of multivariate biorthogonal wavelets with arbitrary vanishing moments, Adv. Comput. Math. 13 (2000), C. K. Chui, W. He and J. Stöckler, Compactly supported tight and sibling frames with maximum vanishing moments, Appl. Comput. Harmon. Anal. 13 (2002), A. Cohen and I. Daubechies, A new technique to estimate the regularity of refinable functions, Rev. Mat. Iberoamericana 12 (1996), A. Cohen and I. Daubechies, Nonseparable bidimensional wavelet bases, Rev. Mat. Iberoamericana 9 (1993), A. Cohen, I. Daubechies and J.C. Feauveau, Biorthogonal bases of compactly supported wavelets, Comm. Pure Appl. Math. 45 (1992), A. Cohen, K. Gröchenig and L. F. Villemoes, Regularity of multivariate refinable functions, Constr. Approx. 15 (1999),
14 A. Cohen and J. M. Schlenker, Compactly supported bidimensional wavelet bases with hexagonal symmetry, Constr. Approx. 9 (1993), I. Daubechies, Orthonormal bases of compactly supported wavelets, Comm. Pure Appl. Math. 41 (1988), I. Daubechies and B. Han, Pairs of dual wavelet frames from any two refinable functions, Constr. Approx. 20 (2004), I. Daubechies, B. Han, A. Ron, and Z. Shen, Framelets: MRA-based constructions of wavelet frames, Appl. Comput. Harmon. Anal. 14 (2003), S. Dahlke, K. Gröchenig and P. Maass, A new approach to interpolating scaling functions, Appl. Anal. 72 (1999), N. Dyn, J. A. Gregory and D. Levin, A butterfly subdivision scheme for surface interpolation with tension control, ACM Trans. Graphics, 9 (1990), B. Han, On dual wavelet tight frames, Appl. Comput. Harmon. Anal. 4 (1997), B. Han, Analysis and construction of optimal multivariate biorthogonal wavelets with compact support, SIAM Math. Anal. 31 (2000), B. Han, Approximation properties and construction of Hermite interpolants and biorthogonal multiwavelets, J. Approx. Theory 110 (2001), B. Han, Symmetry property and construction of wavelets with a general dilation matrix, Linear Algebra Appl. 353, (2002), B. Han, Projectable multivariate refinable functions and biorthogonal wavelets, Appl. Comput. Harmon. Anal. 13 (2002), B. Han, Computing the smoothness exponent of a symmetric multivariate refinable function, SIAM J. Matrix Anal. Appl. 24 (2003), B. Han, Compactly supported tight wavelet frames and orthonormal wavelets of exponential decay with a general dilation matrix, J. Comput. Appl. Math. 155 (2003), B. Han, Vector cascade algorithms and refinable function vectors in Sobolev spaces, J. Approx. Theory 124 (2003), B. Han, Symmetric multivariate orthogonal refinable functions, Appl. Comput. Harmon. Anal., 17 (2004), B. Han, On a conjecture about MRA Riesz wavelet bases, em Proc. Amer. Math. Soc., 134 (2006), B. Han, The projection method in wavelet analysis, Mod. Methods Math., G. Chen and M.J. Lai eds, Nashboro Press, Brentwood, TN, (2006), B. Han, Construction of wavelets and framelets by the projection method, International J. Math. Sci., 1 (2007), to appear. 38. B. Han and R. Q. Jia, Multivariate refinement equations and convergence of subdivision schemes, SIAM J. Math. Anal. 29 (1998), B. Han and R. Q. Jia, Optimal interpolatory subdivision schemes in multidimensional spaces, SIAM J. Numer. Anal. 36 (1999), B. Han and R. Q. Jia, Quincunx fundamental refinable functions and quincunx biorthogonal wavelets, Math. Comp. 71 (2002), B. Han and R. Q. Jia, Optimal C 2 two-dimensional interpolatory ternary subdivision schemes with two-ring stencils, Math. Comp. 75 (2006), B. Han and R. Q. Jia, Characterization of Riesz bases of wavelets generated from multiresolution analysis, Appl. Comput. Harmon. Anal., to appear.
15 B. Han and Z. Shen, Wavelets with short support, SIAM J. Math. Anal. 38 (2006), B. Han and Z. Shen, Wavelets from the Loop scheme, J. Fourier Anal. Appl. 11 (2005), H. Ji, S. D. Riemenschneider and Z. Shen, Multivariate compactly supported fundamental refinable functions, duals, and biorthogonal wavelets, Stud. Appl. Math. 102 (1999), R. Q. Jia, Approximation properties of multivariate wavelets, Comp. Math. 67 (1998), R. Q. Jia, Characterization of smoothness of multivariate refinable functions in Sobolev spaces, Trans. Amer. Math. Soc. 351 (1999), R. Q. Jia and C. A. Micchelli, Using the refinement equation for the construction of pre-wavelets II: Power of two, Curves and Surfaces (P.J. Laurent, A. Le Méhauté and L. L. Schumaker eds.), Academic Press, New York, (1991), R. Q. Jia, J. Z. Wang, and D. X. Zhou, Compactly supported wavelet bases for Sobolev spaces, Appl. Comput. Harmon. Anal. 15 (2003), Q. T. Jiang, On the regularity of matrix refinable functions, SIAM J. Math. Anal. 29 (1998), A. Khodakovsky, P. Schröder, and W. Sweldens, Progressive geometry compression, Proceedings of SIGGRAPH J. Kovačević and M. Vetterli, Nonseparable multidimensional perfect reconstruction filter banks and wavelet bases for R n, IEEE Trans. Inform. Theory 38 (1992), M. J. Lai, Construction of multivariate compactly supported orthonormal wavelets, Adv. Comput. Math. 25 (2006), M. J. Lai and J. Stöckler, Construction of multivariate compactly supported tight wavelet frames, Appl. Comput. Harmon. Anal. 21 (2006), W. Lawton, S. L. Lee and Z. Shen, Convergence of multidimensional cascade algorithm, Numer. Math. 78 (1998), W. Lawton, S. L. Lee and Z. Shen, Stability and orthonormality of multivariate refinable functions, SIAM J. Math. Anal. 28 (1997), E. Le Pennec and S. Mallat, Sparse geometric image representations with bandelets, IEEE Trans. Image Process. 14 (2005), R. Lorentz and P. Oswald, Criteria for hierarchical bases in Sobolev spaces, Appl. Comput. Harmon. Anal. 8 (2000), P. Maass, Families of orthogonal two-dimensional wavelets, SIAM J. Math. Anal. 27 (1996), S. D. Riemenschneider and Z. Shen, Multidimensional interpolatory subdivision schemes, SIAM J. Numer. Anal. 34 (1997), S. D. Riemenschneider and Z. W. Shen, Wavelets and pre-wavelets in low dimensions, J. Approx. Theory, 71 (1992), A. Ron and Z. Shen, Affine systems in L 2 (R d ): the analysis of the analysis operator, J. Funct. Anal. 148 (1997), A. Ron and Z. Shen, Affine systems in L 2 (R d ) II: dual systems, J. Fourier Anal. Appl. 3 (1997), W. Sweldens, The lifting scheme: a custom-design construction of biorthogonal wavelets, Appl. Comput. Harmon. Anal. 3 (1996),
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY
COMPACTLY SUPPORTED ORTHONORMAL COMPLEX WAVELETS WITH DILATION 4 AND SYMMETRY BIN HAN AND HUI JI Abstract. In this paper, we provide a family of compactly supported orthonormal complex wavelets with dilation
More informationBAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS
BAND-LIMITED REFINABLE FUNCTIONS FOR WAVELETS AND FRAMELETS WEIQIANG CHEN AND SAY SONG GOH DEPARTMENT OF MATHEMATICS NATIONAL UNIVERSITY OF SINGAPORE 10 KENT RIDGE CRESCENT, SINGAPORE 119260 REPUBLIC OF
More informationSMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN
SMALL SUPPORT SPLINE RIESZ WAVELETS IN LOW DIMENSIONS BIN HAN, QUN MO, AND ZUOWEI SHEN Abstract. In [B. Han and Z. Shen, SIAM J. Math. Anal., 38 (006), 530 556], a family of univariate short support Riesz
More informationConstruction of Multivariate Compactly Supported Orthonormal Wavelets
Construction of Multivariate Compactly Supported Orthonormal Wavelets Ming-Jun Lai Department of Mathematics The University of Georgia Athens, GA 30602 April 30, 2004 Dedicated to Professor Charles A.
More informationEXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS
EXAMPLES OF REFINABLE COMPONENTWISE POLYNOMIALS NING BI, BIN HAN, AND ZUOWEI SHEN Abstract. This short note presents four examples of compactly supported symmetric refinable componentwise polynomial functions:
More informationConstruction of Biorthogonal Wavelets from Pseudo-splines
Construction of Biorthogonal Wavelets from Pseudo-splines Bin Dong a, Zuowei Shen b, a Department of Mathematics, National University of Singapore, Science Drive 2, Singapore, 117543. b Department of Mathematics,
More informationAvailable at ISSN: Vol. 2, Issue 2 (December 2007) pp (Previously Vol. 2, No.
Available at http://pvamu.edu.edu/pages/398.asp ISSN: 193-9466 Vol., Issue (December 007) pp. 136 143 (Previously Vol., No. ) Applications and Applied Mathematics (AAM): An International Journal A New
More informationMRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces
Chapter 6 MRA Frame Wavelets with Certain Regularities Associated with the Refinable Generators of Shift Invariant Spaces Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan University,
More informationBiorthogonal Spline Type Wavelets
PERGAMON Computers and Mathematics with Applications 0 (00 1 0 www.elsevier.com/locate/camwa Biorthogonal Spline Type Wavelets Tian-Xiao He Department of Mathematics and Computer Science Illinois Wesleyan
More informationHow smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron
How smooth is the smoothest function in a given refinable space? Albert Cohen, Ingrid Daubechies, Amos Ron A closed subspace V of L 2 := L 2 (IR d ) is called PSI (principal shift-invariant) if it is the
More informationSOME SMOOTH COMPACTLY SUPPORTED TIGHT WAVELET FRAMES WITH VANISHING MOMENTS
SOME SMOOTH COMPACTLY SUPPORTED TIGHT WAVELET FRAMES WITH VANISHING MOMENTS A. SAN ANTOLÍN AND R. A. ZALIK Abstract. Let A R d d, d 1 be a dilation matrix with integer entries and det A = 2. We construct
More informationBin Han Department of Mathematical Sciences University of Alberta Edmonton, Canada T6G 2G1
SYMMETRIC ORTHONORMAL SCALING FUNCTIONS AND WAVELETS WITH DILATION FACTOR d = Bin Han Department of Mathematical Sciences University of Alberta Edmonton, Canada T6G 2G1 email: bhan@math.ualberta.ca Abstract.
More informationWAVELETS WITH SHORT SUPPORT
WAVELETS WITH SHORT SUPPORT BIN HAN AND ZUOWEI SHEN Abstract. This paper is to construct Riesz wavelets with short support. Riesz wavelets with short support are of interests in both theory and application.
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationNonstationary Subdivision Schemes and Totally Positive Refinable Functions
Nonstationary Subdivision Schemes and Totally Positive Refinable Functions Laura Gori and Francesca Pitolli December, 2007 Abstract In this paper we construct a class of totally positive refinable functions,
More informationApplied and Computational Harmonic Analysis 11, (2001) doi: /acha , available online at
Applied and Computational Harmonic Analysis 11 305 31 (001 doi:10.1006/acha.001.0355 available online at http://www.idealibrary.com on LETTER TO THE EDITOR Construction of Multivariate Tight Frames via
More informationLINEAR INDEPENDENCE OF PSEUDO-SPLINES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 134, Number 9, September 006, Pages 685 694 S 000-9939(06)08316-X Article electronically published on March 3, 006 LINEAR INDEPENDENCE OF PSEUDO-SPLINES
More informationWavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing
Wavelet Bi-frames with Uniform Symmetry for Curve Multiresolution Processing Qingtang Jiang Abstract This paper is about the construction of univariate wavelet bi-frames with each framelet being symmetric.
More informationConstruction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities
Illinois Wesleyan University From the SelectedWorks of Tian-Xiao He 007 Construction of Biorthogonal B-spline Type Wavelet Sequences with Certain Regularities Tian-Xiao He, Illinois Wesleyan University
More information446 SCIENCE IN CHINA (Series F) Vol. 46 introduced in refs. [6, ]. Based on this inequality, we add normalization condition, symmetric conditions and
Vol. 46 No. 6 SCIENCE IN CHINA (Series F) December 003 Construction for a class of smooth wavelet tight frames PENG Lizhong (Λ Π) & WANG Haihui (Ξ ) LMAM, School of Mathematical Sciences, Peking University,
More informationPairs of Dual Wavelet Frames From Any Two Refinable Functions
Pairs of Dual Wavelet Frames From Any Two Refinable Functions Ingrid Daubechies Bin Han Abstract Starting from any two compactly supported refinable functions in L 2 (R) with dilation factor d, we show
More informationUNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES. Tight compactly supported wavelet frames of arbitrarily high smoothness
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES Tight compactly supported wavelet frames of arbitrarily high smoothness Karlheinz Gröchenig Amos Ron Department of Mathematics U-9 University
More informationThe Method of Virtual Components in the Multivariate Setting
The Method of Virtual Components in the Multivariate Setting Ming-Jun Lai and Alexander Petukhov May 21, 2007 Abstract We describe the so-called method of virtual components for tight wavelet framelets
More informationProperties of Dual Pseudo-Splines
Properties of Dual Pseudo-Splines Bin Dong, Nira Dyn, Kai Hormann IfI Technical Report Series IfI-09-03 Impressum Publisher: Institut für Informatik, Technische Universität Clausthal Julius-Albert Str.
More informationA NOTE ON MATRIX REFINEMENT EQUATIONS. Abstract. Renement equations involving matrix masks are receiving a lot of attention these days.
A NOTE ON MATRI REFINEMENT EQUATIONS THOMAS A. HOGAN y Abstract. Renement equations involving matrix masks are receiving a lot of attention these days. They can play a central role in the study of renable
More informationApproximation Properties and Construction of Hermite Interpolants and Biorthogonal Multiwavelets. Bin Han
Approximation Properties and Construction of Hermite Interpolants and Biorthogonal Multiwavelets Bin Han Abstract. Multiwavelets are generated from refinable function vectors by using multiresolution analysis.
More informationc 2000 Society for Industrial and Applied Mathematics
SIAM J MATH ANAL Vol 32, No 2, pp 420 434 c 2000 Society for Industrial and Applied Mathematics DISTRIBUTIONAL SOLUTIONS OF NONHOMOGENEOUS DISCRETE AND CONTINUOUS REFINEMENT EQUATIONS RONG-QING JIA, QINGTANG
More informationNonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames
Nonlinear Approximation Schemes Associated With Nonseparable Wavelet Bi-frames Martin Ehler September 26, 2008 Abstract In the present paper, we study nonlinear approximation properties of multivariate
More informationDISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS
DISCRETE CDF 9/7 WAVELET TRANSFORM FOR FINITE-LENGTH SIGNALS D. Černá, V. Finěk Department of Mathematics and Didactics of Mathematics, Technical University in Liberec Abstract Wavelets and a discrete
More informationConstruction of wavelets. Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam
Construction of wavelets Rob Stevenson Korteweg-de Vries Institute for Mathematics University of Amsterdam Contents Stability of biorthogonal wavelets. Examples on IR, (0, 1), and (0, 1) n. General domains
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationINFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS. Youngwoo Choi and Jaewon Jung
Korean J. Math. (0) No. pp. 7 6 http://dx.doi.org/0.68/kjm.0...7 INFINITUDE OF MINIMALLY SUPPORTED TOTALLY INTERPOLATING BIORTHOGONAL MULTIWAVELET SYSTEMS WITH LOW APPROXIMATION ORDERS Youngwoo Choi and
More information2 W. LAWTON, S. L. LEE AND ZUOWEI SHEN is called the fundamental condition, and a sequence which satises the fundamental condition will be called a fu
CONVERGENCE OF MULTIDIMENSIONAL CASCADE ALGORITHM W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. Necessary and sucient conditions on the spectrum of the restricted transition operators are given for the
More informationBand-limited Wavelets and Framelets in Low Dimensions
Band-limited Wavelets and Framelets in Low Dimensions Likun Hou a, Hui Ji a, a Department of Mathematics, National University of Singapore, 10 Lower Kent Ridge Road, Singapore 119076 Abstract In this paper,
More informationConstruction of Multivariate Compactly Supported Tight Wavelet Frames
Construction of Multivariate Compactly Supported Tight Wavelet Frames Ming-Jun Lai and Joachim Stöckler April 5, 2006 Abstract Two simple constructive methods are presented to compute compactly supported
More informationWAVELETS WITH COMPOSITE DILATIONS
ELECTRONIC RESEARCH ANNOUNCEMENTS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Pages 000 000 (Xxxx XX, XXXX S 1079-6762(XX0000-0 WAVELETS WITH COMPOSITE DILATIONS KANGHUI GUO, DEMETRIO LABATE, WANG-Q
More informationMultiresolution analysis by infinitely differentiable compactly supported functions. N. Dyn A. Ron. September 1992 ABSTRACT
Multiresolution analysis by infinitely differentiable compactly supported functions N. Dyn A. Ron School of of Mathematical Sciences Tel-Aviv University Tel-Aviv, Israel Computer Sciences Department University
More informationA short introduction to frames, Gabor systems, and wavelet systems
Downloaded from orbit.dtu.dk on: Mar 04, 2018 A short introduction to frames, Gabor systems, and wavelet systems Christensen, Ole Published in: Azerbaijan Journal of Mathematics Publication date: 2014
More informationDecomposition of Riesz frames and wavelets into a finite union of linearly independent sets
Decomposition of Riesz frames and wavelets into a finite union of linearly independent sets Ole Christensen, Alexander M. Lindner Abstract We characterize Riesz frames and prove that every Riesz frame
More information. Introduction Let M be an s s dilation matrix, that is, an s s integer matrix whose eigenvalues lie outside the closed unit disk. The compactly suppo
Multivariate Compactly Supported Fundamental Renable Functions, Duals and Biorthogonal Wavelets y Hui Ji, Sherman D. Riemenschneider and Zuowei Shen Abstract: In areas of geometric modeling and wavelets,
More informationSmoothness of Nonlinear and Non-Separable Subdivision Schemes
Book Title Book Editors IOS Press, 2003 1 Smoothness of Nonlinear and Non-Separable Subdivision Schemes Basarab Matei a,1, Sylvain Meignen b and Anastasia Zakharova b a LAGA Laboratory, Paris XIII University,
More informationAN ALGORITHM FOR MATRIX EXTENSION AND WAVELET CONSTRUCTION
MATHEMATICS OF COMPUTATION Volume 65, Number 14 April 1996, Pages 73 737 AN ALGORITHM FOR MATRIX EXTENSION AND WAVELET CONSTRUCTION W LAWTON, S L LEE AND ZUOWEI SHEN Abstract This paper gives a practical
More informationSpectral properties of the transition operator associated to a multivariate re nement equation q
Linear Algebra and its Applications 292 (1999) 155±178 Spectral properties of the transition operator associated to a multivariate re nement equation q Rong-Qing Jia *, Shurong Zhang 1 Department of Mathematical
More informationIsotropic Multiresolution Analysis: Theory and Applications
Isotropic Multiresolution Analysis: Theory and Applications Saurabh Jain Department of Mathematics University of Houston March 17th 2009 Banff International Research Station Workshop on Frames from first
More informationShift-Invariant Spaces and Linear Operator Equations. Rong-Qing Jia Department of Mathematics University of Alberta Edmonton, Canada T6G 2G1.
Shift-Invariant Spaces and Linear Operator Equations Rong-Qing Jia Department of Mathematics University of Alberta Edmonton, Canada T6G 2G1 Abstract In this paper we investigate the structure of finitely
More informationApproximation by Multiple Refinable Functions
Approximation by Multiple Refinable Functions Rong-Qing Jia, S. D. Riemenschneider, and Ding-Xuan Zhou Department of Mathematical Sciences University of Alberta Edmonton, Canada T6G 2G1 Abstract We consider
More informationMultiscale Frame-based Kernels for Image Registration
Multiscale Frame-based Kernels for Image Registration Ming Zhen, Tan National University of Singapore 22 July, 16 Ming Zhen, Tan (National University of Singapore) Multiscale Frame-based Kernels for Image
More informationDORIN ERVIN DUTKAY AND PALLE JORGENSEN. (Communicated by )
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 OVERSAMPLING GENERATES SUPER-WAVELETS arxiv:math/0511399v1 [math.fa] 16 Nov 2005 DORIN ERVIN DUTKAY
More informationNonseparable multivariate wavelets. Ghan Shyam Bhatt. A dissertation submitted to the graduate faculty
Nonseparable multivariate wavelets by Ghan Shyam Bhatt A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Major: Applied
More informationSparse Multidimensional Representation using Shearlets
Sparse Multidimensional Representation using Shearlets Demetrio Labate a, Wang-Q Lim b, Gitta Kutyniok c and Guido Weiss b, a Department of Mathematics, North Carolina State University, Campus Box 8205,
More informationUNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES. On the construction of multivariate (pre)wavelets
UNIVERSITY OF WISCONSIN-MADISON CENTER FOR THE MATHEMATICAL SCIENCES On the construction of multivariate (pre)wavelets Carl de Boor 1, Ronald A. DeVore 2, and Amos Ron 1 Technical Summary Report #92-09
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. MATH. ANAL. Vol. 30, No. 3, pp. 678 697 c 999 Society for Industrial and Applied Mathematics ARBITRARILY SMOOTH ORTHOGONAL NONSEPARABLE WAVELETS IN R EUGENE BELOGAY AND YANG WANG Abstract. For
More informationOn Dual Wavelet Tight Frames
On Dual Wavelet Tight Frames Bin Han Department of Mathematical Sciences University of Alberta Edmonton, AB T6G 2G, Canada Email: bhan@vega.math.ualberta.ca Abstract. A characterization of multivariate
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationApproximately dual frames in Hilbert spaces and applications to Gabor frames
Approximately dual frames in Hilbert spaces and applications to Gabor frames Ole Christensen and Richard S. Laugesen October 22, 200 Abstract Approximately dual frames are studied in the Hilbert space
More informationarxiv: v1 [cs.oh] 3 Oct 2014
M. Prisheltsev (Voronezh) mikhail.prisheltsev@gmail.com ADAPTIVE TWO-DIMENSIONAL WAVELET TRANSFORMATION BASED ON THE HAAR SYSTEM 1 Introduction arxiv:1410.0705v1 [cs.oh] Oct 2014 The purpose is to study
More informationAffine and Quasi-Affine Frames on Positive Half Line
Journal of Mathematical Extension Vol. 10, No. 3, (2016), 47-61 ISSN: 1735-8299 URL: http://www.ijmex.com Affine and Quasi-Affine Frames on Positive Half Line Abdullah Zakir Husain Delhi College-Delhi
More informationA Riesz basis of wavelets and its dual with quintic deficient splines
Note di Matematica 25, n. 1, 2005/2006, 55 62. A Riesz basis of wavelets and its dual with quintic deficient splines F. Bastin Department of Mathematics B37, University of Liège, B-4000 Liège, Belgium
More informationORTHONORMAL SAMPLING FUNCTIONS
ORTHONORMAL SAMPLING FUNCTIONS N. KAIBLINGER AND W. R. MADYCH Abstract. We investigate functions φ(x) whose translates {φ(x k)}, where k runs through the integer lattice Z, provide a system of orthonormal
More informationConstruction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling
Construction of Orthonormal Quasi-Shearlets based on quincunx dilation subsampling Rujie Yin Department of Mathematics Duke University USA Email: rujie.yin@duke.edu arxiv:1602.04882v1 [math.fa] 16 Feb
More informationFunctions: A Fourier Approach. Universitat Rostock. Germany. Dedicated to Prof. L. Berg on the occasion of his 65th birthday.
Approximation Properties of Multi{Scaling Functions: A Fourier Approach Gerlind Plona Fachbereich Mathemati Universitat Rostoc 1851 Rostoc Germany Dedicated to Prof. L. Berg on the occasion of his 65th
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationDigital Affine Shear Filter Banks with 2-Layer Structure
Digital Affine Shear Filter Banks with -Layer Structure Zhihua Che and Xiaosheng Zhuang Department of Mathematics, City University of Hong Kong, Tat Chee Avenue, Kowloon Tong, Hong Kong Email: zhihuache-c@my.cityu.edu.hk,
More informationComposite Dilation Wavelets with High Degrees
Illinois Wesleyan University From the SelectedWorks of Summer August 18, 2014, Illinois Wesleyan University Available at: https://works.bepress.com/tian_xiao_he/48/ Dept. Math, Illinois Wesleyan University
More informationCompactly Supported Tight Frames Associated with Refinable Functions 1. Communicated by Guido L. Weiss Received July 27, 1999
Applied and Computational Harmonic Analysis 8, 93 319 (000) doi:10.1006/acha.000.0301, available online at http://www.idealibrary.com on LETTER TO THE EDITOR Compactly Supported Tight Frames Associated
More information2 W. LAWTON, S. L. LEE AND ZUOWEI SHEN Denote by L 1 (R d ) and L 2 (R d ) the spaces of Lebesque integrable and modulus square integrable functions d
STABILITY AND ORTHONORMALITY OF MULTIVARIATE REFINABLE FUNCTIONS W. LAWTON, S. L. LEE AND ZUOWEI SHEN Abstract. This paper characterizes the stability and orthonormality of the shifts of a multidimensional
More informationWavelets and Image Compression. Bradley J. Lucier
Wavelets and Image Compression Bradley J. Lucier Abstract. In this paper we present certain results about the compression of images using wavelets. We concentrate on the simplest case of the Haar decomposition
More informationJournal of Mathematical Analysis and Applications. Properties of oblique dual frames in shift-invariant systems
J. Math. Anal. Appl. 356 (2009) 346 354 Contents lists available at ScienceDirect Journal of Mathematical Analysis and Applications www.elsevier.com/locate/jmaa Properties of oblique dual frames in shift-invariant
More informationC -Algebra B H (I) Consisting of Bessel Sequences in a Hilbert Space
Journal of Mathematical Research with Applications Mar., 2015, Vol. 35, No. 2, pp. 191 199 DOI:10.3770/j.issn:2095-2651.2015.02.009 Http://jmre.dlut.edu.cn C -Algebra B H (I) Consisting of Bessel Sequences
More informationWavelets and Image Compression Augusta State University April, 27, Joe Lakey. Department of Mathematical Sciences. New Mexico State University
Wavelets and Image Compression Augusta State University April, 27, 6 Joe Lakey Department of Mathematical Sciences New Mexico State University 1 Signals and Images Goal Reduce image complexity with little
More informationQuadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
Quadrature Prefilters for the Discrete Wavelet Transform Bruce R. Johnson James L. Kinsey Abstract Discrepancies between the Discrete Wavelet Transform and the coefficients of the Wavelet Series are known
More informationThe Construction of Smooth Parseval Frames of Shearlets
Math. Model. Nat. Phenom. Vol., No., 01 The Construction of Smooth Parseval Frames of Shearlets K. Guo b and D. Labate a1 a Department of Mathematics, University of Houston, Houston, Texas 7704 USA b Department
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationWavelets: Theory and Applications. Somdatt Sharma
Wavelets: Theory and Applications Somdatt Sharma Department of Mathematics, Central University of Jammu, Jammu and Kashmir, India Email:somdattjammu@gmail.com Contents I 1 Representation of Functions 2
More informationON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES
ON SAMPLING RELATED PROPERTIES OF B-SPLINE RIESZ SEQUENCES SHIDONG LI, ZHENGQING TONG AND DUNYAN YAN Abstract. For B-splineRiesz sequencesubspacesx span{β k ( n) : n Z}, thereis an exact sampling formula
More informationWe have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma and (3.55) with j = 1. We can write any f W 1 as
88 CHAPTER 3. WAVELETS AND APPLICATIONS We have to prove now that (3.38) defines an orthonormal wavelet. It belongs to W 0 by Lemma 3..7 and (3.55) with j =. We can write any f W as (3.58) f(ξ) = p(2ξ)ν(2ξ)
More informationConstruction of scaling partitions of unity
Construction of scaling partitions of unity arxiv:1710.08290v1 [math.fa] 23 Oct 2017 Ole Christensen and Say Song Goh Abstract Partitions of unity in R d formed by (matrix) scales of a fixed function appear
More informationNonlinear Stationary Subdivision
Nonlinear Stationary Subdivision Michael S. Floater SINTEF P. O. Box 4 Blindern, 034 Oslo, Norway E-mail: michael.floater@math.sintef.no Charles A. Micchelli IBM Corporation T.J. Watson Research Center
More informationFourier-like Transforms
L 2 (R) Solutions of Dilation Equations and Fourier-like Transforms David Malone December 6, 2000 Abstract We state a novel construction of the Fourier transform on L 2 (R) based on translation and dilation
More informationParametrizing orthonormal wavelets by moments
Parametrizing orthonormal wavelets by moments 2 1.5 1 0.5 0 0.5 1 1.5 0 0.5 1 1.5 2 2.5 3 1.4 1.2 1 0.8 0.6 0.4 0.2 0 0.2 0.4 0 1 2 3 4 5 Georg Regensburger Johann Radon Institute for Computational and
More informationINVARIANCE OF A SHIFT-INVARIANT SPACE
INVARIANCE OF A SHIFT-INVARIANT SPACE AKRAM ALDROUBI, CARLOS CABRELLI, CHRISTOPHER HEIL, KERI KORNELSON, AND URSULA MOLTER Abstract. A shift-invariant space is a space of functions that is invariant under
More informationRing-like structures of frequency domains of wavelets
Ring-like structures of frequency domains of wavelets Zhihua Zhang and Naoki aito Dept. of Math., Univ. of California, Davis, California, 95616, UA. E-mail: zzh@ucdavis.edu saito@math.ucdavis.edu Abstract.
More informationMoment Computation in Shift Invariant Spaces. Abstract. An algorithm is given for the computation of moments of f 2 S, where S is either
Moment Computation in Shift Invariant Spaces David A. Eubanks Patrick J.Van Fleet y Jianzhong Wang ẓ Abstract An algorithm is given for the computation of moments of f 2 S, where S is either a principal
More informationDAVID FERRONE. s k s k 2j = δ 0j. s k = 1
FINITE BIORTHOGONAL TRANSFORMS AND MULTIRESOLUTION ANALYSES ON INTERVALS DAVID FERRONE 1. Introduction Wavelet theory over the entire real line is well understood and elegantly presented in various textboos
More informationON SPECTRAL CANTOR MEASURES. 1. Introduction
ON SPECTRAL CANTOR MEASURES IZABELLA LABA AND YANG WANG Abstract. A probability measure in R d is called a spectral measure if it has an orthonormal basis consisting of exponentials. In this paper we study
More informationApplications of Polyspline Wavelets to Astronomical Image Analysis
VIRTUAL OBSERVATORY: Plate Content Digitization, Archive Mining & Image Sequence Processing edited by M. Tsvetkov, V. Golev, F. Murtagh, and R. Molina, Heron Press, Sofia, 25 Applications of Polyspline
More information290 J.M. Carnicer, J.M. Pe~na basis (u 1 ; : : : ; u n ) consisting of minimally supported elements, yet also has a basis (v 1 ; : : : ; v n ) which f
Numer. Math. 67: 289{301 (1994) Numerische Mathematik c Springer-Verlag 1994 Electronic Edition Least supported bases and local linear independence J.M. Carnicer, J.M. Pe~na? Departamento de Matematica
More informationOctober 7, :8 WSPC/WS-IJWMIP paper. Polynomial functions are renable
International Journal of Wavelets, Multiresolution and Information Processing c World Scientic Publishing Company Polynomial functions are renable Henning Thielemann Institut für Informatik Martin-Luther-Universität
More informationConstruction of orthogonal multiscaling functions and multiwavelets with higher approximation order based on the matrix extension algorithm
J. Math. Kyoto Univ. (JMKYAZ 6- (6, 75 9 Construction of orthogonal multiscaling functions and multiwavelets with higher approximation order based on the matrix extension algorithm By Shouzhi Yang and
More informationSemi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform
NTMSCI 6, No., 175-183 018) 175 New Trends in Mathematical Sciences http://dx.doi.org/10.085/ntmsci.018.83 Semi-orthogonal wavelet frames on positive half-line using the Walsh Fourier transform Abdullah
More informationWAVELET EXPANSIONS OF DISTRIBUTIONS
WAVELET EXPANSIONS OF DISTRIBUTIONS JASSON VINDAS Abstract. These are lecture notes of a talk at the School of Mathematics of the National University of Costa Rica. The aim is to present a wavelet expansion
More informationTwo-Dimensional Orthogonal Filter Banks with Directional Vanishing Moments
Two-imensional Orthogonal Filter Banks with irectional Vanishing Moments Jianping Zhou and Minh N. o epartment of Electrical and Computer Engineering University of Illinois at Urbana-Champaign, Urbana,
More informationWavelets and modular inequalities in variable L p spaces
Wavelets and modular inequalities in variable L p spaces Mitsuo Izuki July 14, 2007 Abstract The aim of this paper is to characterize variable L p spaces L p( ) (R n ) using wavelets with proper smoothness
More informationarxiv: v2 [math.fa] 27 Sep 2016
Decomposition of Integral Self-Affine Multi-Tiles Xiaoye Fu and Jean-Pierre Gabardo arxiv:43.335v2 [math.fa] 27 Sep 26 Abstract. In this paper we propose a method to decompose an integral self-affine Z
More informationSchemes. Philipp Keding AT15, San Antonio, TX. Quarklet Frames in Adaptive Numerical. Schemes. Philipp Keding. Philipps-University Marburg
AT15, San Antonio, TX 5-23-2016 Joint work with Stephan Dahlke and Thorsten Raasch Let H be a Hilbert space with its dual H. We are interested in the solution u H of the operator equation Lu = f, with
More informationLectures notes. Rheology and Fluid Dynamics
ÉC O L E P O L Y T E C H N IQ U E FÉ DÉR A L E D E L A U S A N N E Christophe Ancey Laboratoire hydraulique environnementale (LHE) École Polytechnique Fédérale de Lausanne Écublens CH-05 Lausanne Lectures
More informationFrame Diagonalization of Matrices
Frame Diagonalization of Matrices Fumiko Futamura Mathematics and Computer Science Department Southwestern University 00 E University Ave Georgetown, Texas 78626 U.S.A. Phone: + (52) 863-98 Fax: + (52)
More informationA real algebra perspective on multivariate tight wavelet frames
A real algebra perspective on multivariate tight wavelet frames arxiv:1202.3596v4 [math.fa] 9 Jul 2012 Maria Charina, Mihai Putinar, Claus Scheiderer and Joachim Stöckler June 15, 2018 Abstract Recent
More informationBoundary functions for wavelets and their properties
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 009 Boundary functions for wavelets and their properties Ahmet Alturk Iowa State University Follow this and additional
More information