Quadrature Prefilters for the Discrete Wavelet Transform. Bruce R. Johnson. James L. Kinsey. Abstract
|
|
- Shanon Beasley
- 5 years ago
- Views:
Transcription
1 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 to be reducible by initialization of input data. Prefilters based on Lagrange interpolants are derived here for biorthogonal compact support wavelet systems, providing exact subspace proection in cases of local polynomial smoothness. The resulting convergence acceleration in a non-polynomial test case is examined. Irregular sampling rates are also accommodated. * This work was supported by the Robert A. Welch Foundation and by Grant CHE from the National Science Foundation. * The authors are with the Department of Chemistry and the Rice Quantum Institute, MS #600, Rice University, Houston, TX EDICS Classification: SP I. Introduction The use of FIR filters based on compact support wavelet families allows for efficient digital multiresolution analysis of time series and other signals [1-4]. A general signal f ( t) may be expanded approximately in a finite Wavelet Series (WS) ( ) f t 0 φ k, f 0 φk ( t ) + ψ k, f ψ k ( t ), (1) k J 1 =0 where the synthesis scaling functions obey k
2 φ k ( t ) = 2 /2 φ( 2 /2 t k ), (2) and similar relations hold for the wavelets ψ k and the biorthogonal analysis functions φ k and ψ k. All of the expansion coefficients in subspaces J 1,J 2,...,0 can be generated from the proection integrals φ k J, f using the two-scale relations [5] φ k 1 t ( ) = ( t), ψ 1 k t k h k φ 2k+ k ( ) = ( t). (3) k g k ψ 2k+ k Thus it is the evaluation of the integrals φ k J, f from discrete input data which is addressed here. Their accurate calculation is important in order to exploit the computational efficiency of the Discrete Wavelet Transform (DWT) for signal analysis. For nonsingular functions f ( t), one may approximate the integrals on the finest octave by simple function samples, φ J k, f f k / 2 J ( ) = [ k]. (4) This underlies the standard Mallat Algorithm [2] [6], fast pyramidal filtering of an input data sequence. Shensa [7] and Rioul and Duhamel [8] clarified the point that, for general compact support wavelets, this procedure does not precisely produce the coefficients of the WS. The essential difficulty is that of ensuring accurate proection to V J, the subspace spanned by the φ k J, when limited to knowledge only of discrete data that is not strongly oversampled. The Shensa Algorithm [7] [8] focuses instead on a function ˆ ( t) = [ k] χ J ( t k ), (5) k which is defined in terms of the J level samples and a sampling function χ J designed to make ˆ in some manner close to f. The V J proections of this function are then given by φ J k, f ˆ J = [ l] dt φ J k t l ( ) χj t l ( ). (6) 2
3 J One approach takes χ J ( t l )=φ l ( t ), for which the right hand side of Eq. (6) reduces to [ k]; that is, ˆ ( t) is entirely in V J and its expansion coefficients are exactly given by the samples of f ( t) [7]. A different approach requires the samples of f ( t) and ˆ ( t) to coincide, using interpolating functions with the property χ J ( n) =δ n,0. Thus, for instance, the use of sinc functions has been explored [9] [10]. For this case, the integral in Eq. (6), which constitutes a prefiltering of the data, is more complicated. Other choices of prefiltering have also been examined, in some cases optimized to the specific data [11] [Chan, 1996 #29]. II. Lagrange Interpolation It is natural to adopt polynomial interpolation here since the scaling functions span loworder polynomials by construction. Such a procedure is directly connected to numerical quadrature estimation of the proection integrals [12-16]. The interpolants in this case are Lagrange polynomials, already familiar in wavelet-related contexts [17] [7] [Donoho, 1992 #30]. A unique level of generality may be obtained here since it is not necessary to assume uniform sample spacing for the derivation [16]. Thus, for some ordered set of distinct points t q in the neighborhood of a particular φ k J, ˆ r ( ) ( t) = f t q L r,q ( t ), (7) q=1 where the Lagrange polynomials are L r,q ( t ) = r q q t t q t q t q (8) If f ( t) is a polynomial of order less than r, then ˆ ( t) exactly equals f ( t); otherwise they are only equal at the sample points. A different set of points is generally to be used for the next value of k, though the sets may overlap; this reduces to a simple unit translation [cf., Eq. (5)] in the uniform-spacing case. The general-spacing analog of Eq. (6) is now 3
4 r ( ) dt φ J k, f ˆ J J = f t q φ k ( t ) L r,q t q=1 ( ) r ( ) = f t q q=1 r 1 p= 0 1 p! ( p) L r,q ( k /2 J ) dt t k /2 J ( ) p ( ) φ k J t r ( ) = ω q f t q, (9) q=1 where r 1 ω q = 2 J/2 Jp m p p! p= 0 ( p) L r,q ( k /2 J ), (10) m p = dt t p φ ( t ) = 1 2 p p 1 ( ) 1 p p m p µ p p, (11) p =0 µ p = 1 2 L 1 h k k p. (12) k=0 Equation (9) reflects the well-known connection between Lagrange interpolation of a function and numerical quadrature of its integrals. The quadrature weights ω q (related to those derived by Sweldens and Piessens [14]) are obtained by Taylor expansion of the L r,q around the position of φ k J. The scaling function moments m p can be iteratively calculated [Eq. (11)] using the discrete moments µ p as shown in [12]. Thus, the weights are broken into parts m p (calculated once) ( p) specific to the wavelet basis and L r,q specific to k and the associated points tq. If the latter have ( p) the same relative configuration for each k, as in the case of uniform sample spacing, the L r,q also only need to be calculated once. We will concentrate on the uniform case below. III. Simulation Results For testing purposes, the symmetric biorthogonal wavelet system with recursion lengths ( L, L) = (9, 7) described by Cohen, et al. [5] and Antonini, et al. [18] is used. Symmetry forces 4
5 the odd moments to be exactly zero here in contrast to the situation for orthogonal wavelet families [1], [12]. A damped and chirped sine wave, ( ) = e t 2 /8 sin t ( 1 + t ) f t [ ], (13) is chosen as a test function and illustrated in Fig. 1. There are several oscillations within the support, rendering inaccurate any few-point quadrature at the = 0 level. To avoid undue complexity in the quadrature computations, samples were chosen symmetrically distributed around the center of the support (2 J k) and with spacing 2 J. For odd r, these sets were of the form 2 J k, 2 J (k±1), 2 J (k±2), etc. For even r, they were of the form 2 J (k±1/2), 2 J (k±3/2), etc. Table 1 shows the values of the = k = 0 integral φ, f calculated by quadrature on different octaves J and iterating with the low-pass filter h k back to = 0. Natural sampling [cf. Eq. (4)] is equivalent to quadrature of order r = 1; a 7-point series of calculations is also shown. Both series ultimately lead to the same final result, but natural sampling proceeds at a very slow rate by comparison. The rate of convergence has been examined as a function of both quadrature order r and starting octave J. Figure 2 shows the relative error obtained for odd r with respect to the converged value For J = 0 (no recursion), the error is always large, as expected, except for the fortuitously accurate case r = 3; higher order r accomplishes nothing except to place new sample points outside of the range of importance of f(x) since the spacing is large on this scale. The picture changes dramatically with higher J, where strong dependence of the accuracy on r is found. One finds from detailed examination the following two patterns: for a given even value of r (not shown here), the asymptotic error decreases as 2 Jr with increasing J while, for a given odd value of r, the error decreases as 2 J(r +1). The former case is the expected asymptotic behavior [19], [14]. The accelerated convergence in the odd case can be understood from consideration of Eq. (9) in the case of monomials centered around k /2 J r ( ) p ω q t q k / 2 J = 2 J/2 Jp m p, 0 p r 1. (14) q =1 5
6 The symmetrical distribution of sample points and the vanishing of m r for odd r leads to automatic satisfaction of one higher equation (p = r ). (This situation is similar to an accelerated convergence found with a specific one-point quadrature for Daubechies scaling functions [12].) For r = 1(no prefiltering), the asymptotic convergence rate is 4 J rather than 2 J since the sample is taken at the center of the support. For higher r the convergence is very strong, drastically reducing the number of scales that must be traversed to meet a given convergence threshold. IV. Conclusion Where are these considerations important in practical applications? They are certainly critical in numerical analysis, solution of partial differential equations and some cases of spectral analysis where natural sampling would lead to the need for excessive refinement in pursuit of convergence. On the other hand, many signal and image processing applications do not seek such accuracy as a result of noise, inherent limits of precision, and/or lack of any known functions underlying discrete experimental data. In such cases one may use the standard DWT with natural sampling and regard the resulting multiresolution analysis as a sampling of the exact wavelet transform of the function ˆ J ( t)= [ k]φ k ( t ), as discussed above. This has the appearance of k reducing the matter to one of interpretation. In actuality, however, the behavior of ˆ output of its DWT are sensitive to the value o until the scale is so fine that all relevant frequencies are resolved (as in the first column of Table 1). The occasionally-encountered ( t) and the statement that single wavelets do not require use of prefilters is really a statement about the case for sufficiently large J. In many cases where the sampling rate has a hard upper limit, the proectionbased prefilter described above can be regarded as extrapolating the ordinary DWT algorithm to higher resolution. The order r 1of polynomials for which the Lagrange-polynomial-based prefilter gives exact subspace proection is not limited by the approximating power of the wavelets chosen. In this sense, r is a design parameter. To some extent, the same may be said of the sample positions used. Numerical experiments verify a certain level of robustness of the calculated WS coefficients 6
7 with respect to irregular sampling. This can be useful for specific applications, e.g., unevenly spaced antenna arrays, or for generic situations such as drop-outs in data sequences. As expected for any polynomial interpolation, there is reasonable stability as long as sample points do not approach each other too closely and the interpolated point is roughly within their range. On the other hand, numerical experiments including strong noise contamination show no greater stability than for natural sampling, indicating an avenue needing further exploration. Acknowledgment Many thanks are due to R. O. Wells, Jr., C. S. Burrus, J. E. Odegard, and other members of the Computational Mathematics Laboratory and the Digital Signal Processing group at Rice University for extended interactions. Useful conversations with W. Sweldens and G. Strang are also gratefully acknowledged. References [1] I. Daubechies, Ten Lectures on Wavelets, vol. 61. Philadelphia: SIAM Publications, [2] S. G. Mallat, A theory for multiresolution signal decomposition: The wavelet representation, IEEE Trans. Pattern Anal. Machine Intell., vol. 11, pp , [3] G. Strang and T. Nguyen, Wavelets and Filter Banks. Wellesley: Wellesley-Cambridge Press, [4] C. S. Burrus, R. A. Gopinath, and H. Guo, Introduction to Wavelets and Wavelet Transforms: Prentice-Hall, [5] A. Cohen, I. Daubechies, and J.-C. Feauveau, Biorthogonal Bases of Compactly Supported Wavelets, Comm. Pure Appl. Math., vol. 45, pp , [6] S. Mallat, Multifrequency channel decomposition of images and wavelet methods, IEEE Trans. Acoust. Speech, Signal Processing, vol. 37, pp , [7] M. J. Shensa, The Discrete Wavelet Transform: Wedding the À Trous and Mallat Algorithms, IEEE Trans. Signal Proc., vol. 40, pp , [8] O. Rioul and P. Duhamel, Fast Algorithms for Discrete and Continuous Wavelet Transforms, IEEE Trans. Inform. Theory, vol. 38, pp ,
8 [9] P. Abry and P. Flandrin, On the Initialization of the Discrete Wavelet Transform Algorithm, IEEE Sig. Proc. Lett., vol. 1, pp , [10] X.-P. Zhang, L.-S. Tian, and Y.-N. Peng, From the Wavelet Series to the Discrete Wavelet Transform - The Initialization, IEEE Trans. Signal Proc., vol. 44, pp , [11] X.-G. Xia, C.-C. J. Kuo, and Z. Zhang, Wavelet Coefficient Computation with Optimal Prefiltering, IEEE Trans. Signal Proc., vol. 42, pp , [12] R. A. Gopinath and C. S. Burris, On the moments of the scaling function φ, presented at Proceedings of the ISCAS-92, San Diego, [13] G. Beylkin, R. Coifman, and V. Rokhlin, Fast Wavelet Transforms and Numerical Algorithms I, Comm. Pure and Appl. Math, vol. 44, pp , [14] W. Sweldens and R. Piessens, Quadrature Formulae and Asymptotic Error Expansions for Wavelet Approximations of Smooth Functions, SIAM J. Numer. Anal., vol. 31, pp , [15] W. Sweldens and R. Piessens, Asymptotic error expansion of wavelet approximations of smooth functions II, Numerische Mathematik, vol. 68, pp. 377, [16] B. R. Johnson, J. P. Modisette, P. J. Nordlander, and J. L. Kinsey, Quadrature Integration for Orthogonal Wavelet Systems, J. Chem. Phys., vol. 110, pp , [17] G. Deslauriers and S. Dubuc, Symmetric Iterative Interpolation Processes, Constr. Approx., vol. 5, pp , [18] M. Antonini, M. Barlaud, P. Mathieu, and I. Daubechies, Image Coding Using Wavelet Transform, IEEE Trans. Image Proc., vol. 1, pp , [19] G. Strang, Wavelets and dilation equations: A brief introduction, SIAM Review, vol. 31, pp ,
9 Table Table 1. Integral φ 0 0, f for the symmetric (9,7) biorthogonal Wavelet System from Refs. [5] and [18] and the signal of Eq. (13), as calculated starting from various octaves J. J Natural Sampling Wavelet quadrature (r = 7) Figures t Figure 1. Damped and chirped sine wave. 9
10 J = 0 J = 1 J = 2 J = 3 J = Quadrature Order r Figure 2. Relative error φ [, f φ, f calc exact ] / φ, f exact for different quadrature orders and different finest scales used with recursion. The function f ( t) is as given in Fig
Introduction to Wavelets and Wavelet Transforms
Introduction to Wavelets and Wavelet Transforms A Primer C. Sidney Burrus, Ramesh A. Gopinath, and Haitao Guo with additional material and programs by Jan E. Odegard and Ivan W. Selesnick Electrical and
More informationQuadrature integration for orthogonal wavelet systems
Quadrature integration for orthogonal wavelet systems Bruce R. Johnson, a Jason P. Modisette, b Peter J. Nordlander, b and James L. Kinsey a a Department of Chemistry and Rice Quantum Institute Rice University
More informationA Higher-Density Discrete Wavelet Transform
A Higher-Density Discrete Wavelet Transform Ivan W. Selesnick Abstract In this paper, we describe a new set of dyadic wavelet frames with three generators, ψ i (t), i =,, 3. The construction is simple,
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 informationVARIOUS types of wavelet transform are available for
IEEE TRANSACTIONS ON SIGNAL PROCESSING A Higher-Density Discrete Wavelet Transform Ivan W. Selesnick, Member, IEEE Abstract This paper describes a new set of dyadic wavelet frames with two generators.
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 informationMoments conservation in adaptive Vlasov solver
12 Moments conservation in adaptive Vlasov solver M. Gutnic a,c, M. Haefele b,c and E. Sonnendrücker a,c a IRMA, Université Louis Pasteur, Strasbourg, France. b LSIIT, Université Louis Pasteur, Strasbourg,
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 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 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 informationThe Application of Legendre Multiwavelet Functions in Image Compression
Journal of Modern Applied Statistical Methods Volume 5 Issue 2 Article 3 --206 The Application of Legendre Multiwavelet Functions in Image Compression Elham Hashemizadeh Department of Mathematics, Karaj
More informationQUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS *
QUADRATURES INVOLVING POLYNOMIALS AND DAUBECHIES' WAVELETS * WEI-CHANG SHANN AND JANN-CHANG YAN Abstract. Scaling equations are used to derive formulae of quadratures involving polynomials and scaling/wavelet
More informationMultiresolution image processing
Multiresolution image processing Laplacian pyramids Some applications of Laplacian pyramids Discrete Wavelet Transform (DWT) Wavelet theory Wavelet image compression Bernd Girod: EE368 Digital Image Processing
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 informationSolution of Cartesian and Curvilinear Quantum. Equations via Multiwavelets on the Interval
Solution of Cartesian and Curvilinear Quantum Equations via Multiwavelets on the Interval Bruce R. Johnson, Jeffrey. Mackey, a and James. Kinsey Department of Chemistry and Rice Quantum Institute Rice
More informationLittlewood Paley Spline Wavelets
Proceedings of the 6th WSEAS International Conference on Wavelet Analysis & Multirate Systems, Bucharest, Romania, October 6-8, 26 5 Littlewood Paley Spline Wavelets E. SERRANO and C.E. D ATTELLIS Escuela
More informationLifting Parameterisation of the 9/7 Wavelet Filter Bank and its Application in Lossless Image Compression
Lifting Parameterisation of the 9/7 Wavelet Filter Bank and its Application in Lossless Image Compression TILO STRUTZ Deutsche Telekom AG, Hochschule für Telekommunikation Institute of Communications Engineering
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 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 informationLecture 16: Multiresolution Image Analysis
Lecture 16: Multiresolution Image Analysis Harvey Rhody Chester F. Carlson Center for Imaging Science Rochester Institute of Technology rhody@cis.rit.edu November 9, 2004 Abstract Multiresolution analysis
More informationOn the fast algorithm for multiplication of functions in the wavelet bases
Published in Proceedings of the International Conference Wavelets and Applications, Toulouse, 1992; Y. Meyer and S. Roques, edt., Editions Frontieres, 1993 On the fast algorithm for multiplication of functions
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing () Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids Subband coding
More informationMULTIRATE DIGITAL SIGNAL PROCESSING
MULTIRATE DIGITAL SIGNAL PROCESSING Signal processing can be enhanced by changing sampling rate: Up-sampling before D/A conversion in order to relax requirements of analog antialiasing filter. Cf. audio
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 informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
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 informationThe optimal filtering of a class of dynamic multiscale systems
Science in China Ser. F Information Sciences 2004 Vol.47 No.4 50 57 50 he optimal filtering of a class of dynamic multiscale systems PAN Quan, ZHANG Lei, CUI Peiling & ZHANG Hongcai Department of Automatic
More informationThe Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems
668 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 49, NO 10, OCTOBER 2002 The Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems Lei Zhang, Quan
More informationCOMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS
COMPLEX WAVELET TRANSFORM IN SIGNAL AND IMAGE ANALYSIS MUSOKO VICTOR, PROCHÁZKA ALEŠ Institute of Chemical Technology, Department of Computing and Control Engineering Technická 905, 66 8 Prague 6, Cech
More informationWavelet-Based Numerical Homogenization for Scaled Solutions of Linear Matrix Equations
International Journal of Discrete Mathematics 2017; 2(1: 10-16 http://www.sciencepublishinggroup.com/j/dmath doi: 10.11648/j.dmath.20170201.13 Wavelet-Based Numerical Homogenization for Scaled Solutions
More informationaxioms Construction of Multiwavelets on an Interval Axioms 2013, 2, ; doi: /axioms ISSN
Axioms 2013, 2, 122-141; doi:10.3390/axioms2020122 Article OPEN ACCESS axioms ISSN 2075-1680 www.mdpi.com/journal/axioms Construction of Multiwavelets on an Interval Ahmet Altürk 1 and Fritz Keinert 2,
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 informationTowards Global Design of Orthogonal Filter Banks and Wavelets
Towards Global Design of Orthogonal Filter Banks and Wavelets Jie Yan and Wu-Sheng Lu Department of Electrical and Computer Engineering University of Victoria Victoria, BC, Canada V8W 3P6 jyan@ece.uvic.ca,
More informationAnalysis of Redundant-Wavelet Multihypothesis for Motion Compensation
Analysis of Redundant-Wavelet Multihypothesis for Motion Compensation James E. Fowler Department of Electrical and Computer Engineering GeoResources Institute GRI Mississippi State University, Starville,
More informationWavelet transforms and compression of seismic data
Mathematical Geophysics Summer School, Stanford, August 1998 Wavelet transforms and compression of seismic data G. Beylkin 1 and A. Vassiliou 2 I Introduction Seismic data compression (as it exists today)
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 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 informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationWavelet Bases of the Interval: A New Approach
Int. Journal of Math. Analysis, Vol. 1, 2007, no. 21, 1019-1030 Wavelet Bases of the Interval: A New Approach Khaled Melkemi Department of Mathematics University of Biskra, Algeria kmelkemi@yahoo.fr Zouhir
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 informationIntroduction to Discrete-Time Wavelet Transform
Introduction to Discrete-Time Wavelet Transform Selin Aviyente Department of Electrical and Computer Engineering Michigan State University February 9, 2010 Definition of a Wavelet A wave is usually defined
More informationCOMPACTLY 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 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 informationHarmonic Wavelet Transform and Image Approximation
Harmonic Wavelet Transform and Image Approximation Zhihua Zhang and Naoki Saito Dept of Math, Univ of California, Davis, California, 95616, USA Email: zhangzh@mathucdavisedu saito@mathucdavisedu Abstract
More informationWavelet Transform. Figure 1: Non stationary signal f(t) = sin(100 t 2 ).
Wavelet Transform Andreas Wichert Department of Informatics INESC-ID / IST - University of Lisboa Portugal andreas.wichert@tecnico.ulisboa.pt September 3, 0 Short Term Fourier Transform Signals whose frequency
More informationDenoising and Compression Using Wavelets
Denoising and Compression Using Wavelets December 15,2016 Juan Pablo Madrigal Cianci Trevor Giannini Agenda 1 Introduction Mathematical Theory Theory MATLAB s Basic Commands De-Noising: Signals De-Noising:
More informationABSTRACT. Design of vibration inspired bi-orthogonal wavelets for signal analysis. Quan Phan
ABSTRACT Design of vibration inspired bi-orthogonal wavelets for signal analysis by Quan Phan In this thesis, a method to calculate scaling function coefficients for a new biorthogonal wavelet family derived
More informationToward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information
Toward a Realization of Marr s Theory of Primal Sketches via Autocorrelation Wavelets: Image Representation using Multiscale Edge Information Naoki Saito 1 Department of Mathematics University of California,
More informationA Wavelet-Based Technique for Identifying, Labeling, and Tracking of Ocean Eddies
MARCH 2002 LUO AND JAMESON 381 A Wavelet-Based Technique for Identifying, Labeling, and Tracking of Ocean Eddies JINGJIA LUO Department of Earth and Planetary Physics, Graduate School of Science, University
More informationRepresentation: Fractional Splines, Wavelets and related Basis Function Expansions. Felix Herrmann and Jonathan Kane, ERL-MIT
Representation: Fractional Splines, Wavelets and related Basis Function Expansions Felix Herrmann and Jonathan Kane, ERL-MIT Objective: Build a representation with a regularity that is consistent with
More informationWAVELET EXPANSIONS IN VOLUME INTEGRAL METHOD OF EDDY-CURRENT MODELING
WAVELET EXPANSIONS IN VOLUME INTEGRAL METHOD OF EDDY-CURRENT MODELING Bing Wang John P. Basart and John C. Moulder Center for NDE and Department of Electrical Engineering and Computer Engineering Iowa
More informationA Novel Fast Computing Method for Framelet Coefficients
American Journal of Applied Sciences 5 (11): 15-157, 008 ISSN 1546-939 008 Science Publications A Novel Fast Computing Method for Framelet Coefficients Hadeel N. Al-Taai Department of Electrical and Electronic
More informationA Comparative Study of Non-separable Wavelet and Tensor-product. Wavelet; Image Compression
Copyright c 007 Tech Science Press CMES, vol., no., pp.91-96, 007 A Comparative Study o Non-separable Wavelet and Tensor-product Wavelet in Image Compression Jun Zhang 1 Abstract: The most commonly used
More informationEFFICIENT IMAGE COMPRESSION ALGORITHMS USING EVOLVED WAVELETS
International Journal of ISSN 0974-2107 Systems and Technologies IJST Vol.4, No.2, pp 127-146 KLEF 2011 EFFICIENT IMAGE COMPRESSION ALGORITHMS USING EVOLVED WAVELETS 1.G.CHENCHU KRISHNAIAH* 2.T.JAYACHANDRAPRASAD
More informationJPEG2000 has recently been approved as an international
1080 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 12, NO. 9, SEPTEMBER 2003 Mathematical Properties of the JPEG2000 Wavelet Filters Michael Unser, Fellow, IEEE, and Thierry Blu, Member, IEEE Abstract The
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More informationIdentification and Classification of High Impedance Faults using Wavelet Multiresolution Analysis
92 NATIONAL POWER SYSTEMS CONFERENCE, NPSC 2002 Identification Classification of High Impedance Faults using Wavelet Multiresolution Analysis D. Cha N. K. Kishore A. K. Sinha Abstract: This paper presents
More information1 Introduction to Wavelet Analysis
Jim Lambers ENERGY 281 Spring Quarter 2007-08 Lecture 9 Notes 1 Introduction to Wavelet Analysis Wavelets were developed in the 80 s and 90 s as an alternative to Fourier analysis of signals. Some of the
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 informationWavelets in Image Compression
Wavelets in Image Compression M. Victor WICKERHAUSER Washington University in St. Louis, Missouri victor@math.wustl.edu http://www.math.wustl.edu/~victor THEORY AND APPLICATIONS OF WAVELETS A Workshop
More informationDesign of Orthonormal Wavelet Filter Banks Using the Remez Exchange Algorithm
Electronics and Communications in Japan, Part 3, Vol. 81, No. 6, 1998 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J80-A, No. 9, September 1997, pp. 1396 1402 Design of Orthonormal Wavelet
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 informationDesign and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform
Design and Implementation of Multistage Vector Quantization Algorithm of Image compression assistant by Multiwavelet Transform Assist Instructor/ BASHAR TALIB HAMEED DIYALA UNIVERSITY, COLLEGE OF SCIENCE
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 information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationREGULARITY AND CONSTRUCTION OF BOUNDARY MULTIWAVELETS
REGULARITY AND CONSTRUCTION OF BOUNDARY MULTIWAVELETS FRITZ KEINERT Abstract. The conventional way of constructing boundary functions for wavelets on a finite interval is to form linear combinations of
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationA CLASS OF M-DILATION SCALING FUNCTIONS WITH REGULARITY GROWING PROPORTIONALLY TO FILTER SUPPORT WIDTH
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 26, Number 2, December 998, Pages 350 3506 S 0002-9939(98)05070-9 A CLASS OF M-DILATION SCALING FUNCTIONS WITH REGULARITY GROWING PROPORTIONALLY
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 informationNinth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673
Ninth International Water Technology Conference, IWTC9 2005, Sharm El-Sheikh, Egypt 673 A NEW NUMERICAL APPROACH FOR THE SOLUTION OF CONTAMINANT TRANSPORT EQUATION Mohamed El-Gamel Department of Mathematical
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 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 informationWavelets and Signal Processing
Wavelets and Signal Processing John E. Gilbert Mathematics in Science Lecture April 30, 2002. Publicity Mathematics In Science* A LECTURE SERIES FOR UNDERGRADUATES Wavelets Professor John Gilbert Mathematics
More informationDigital Image Processing
Digital Image Processing Wavelets and Multiresolution Processing (Wavelet Transforms) Christophoros Nikou cnikou@cs.uoi.gr University of Ioannina - Department of Computer Science 2 Contents Image pyramids
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 informationOn the Hilbert Transform of Wavelets
On the Hilbert Transform of Wavelets Kunal Narayan Chaudhury and Michael Unser Abstract A wavelet is a localized function having a prescribed number of vanishing moments. In this correspondence, we provide
More informationA Tutorial on Wavelets and their Applications. Martin J. Mohlenkamp
A Tutorial on Wavelets and their Applications Martin J. Mohlenkamp University of Colorado at Boulder Department of Applied Mathematics mjm@colorado.edu This tutorial is designed for people with little
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More informationTime series denoising with wavelet transform
Paper Time series denoising with wavelet transform Bartosz Kozłowski Abstract This paper concerns the possibilities of applying wavelet analysis to discovering and reducing distortions occurring in time
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 informationFrom Fourier to Wavelets in 60 Slides
From Fourier to Wavelets in 60 Slides Bernhard G. Bodmann Math Department, UH September 20, 2008 B. G. Bodmann (UH Math) From Fourier to Wavelets in 60 Slides September 20, 2008 1 / 62 Outline 1 From Fourier
More informationModule 4 MULTI- RESOLUTION ANALYSIS. Version 2 ECE IIT, Kharagpur
Module 4 MULTI- RESOLUTION ANALYSIS Lesson Theory of Wavelets Instructional Objectives At the end of this lesson, the students should be able to:. Explain the space-frequency localization problem in sinusoidal
More informationTHE APPLICATION OF MULTIWAVELET FILTER BANKS TO IMAGE PROCESSING
THE APPLICATION OF MULTIWAVELET FILTER BANKS TO IMAGE PROCESSING V. Strela,, P. N. Heller, G. Strang, P. Topiwala, C. Heil Abstract Multiwavelets are a new addition to the body of wavelet theory. Realizable
More informationRobust Nonlinear Wavelet Transform based on Median-Interpolation
Robust Nonlinear Wavelet Transform based on Median-Interpolation David L. Donoho & Thomas P.Y. Yu Department of Statistics Stanford University Stanford, CA 943, USA donoho@stat.stanford.edu pyu@stat.stanford.edu
More informationDirect Design of Orthogonal Filter Banks and Wavelets
Direct Design of Orthogonal Filter Banks and Wavelets W.-S. Lu T. Hinamoto Dept. of Electrical & Computer Engineering Graduate School of Engineering University of Victoria Hiroshima University Victoria,
More informationWavelets in Scattering Calculations
Wavelets in Scattering Calculations W. P., Brian M. Kessler, Gerald L. Payne polyzou@uiowa.edu The University of Iowa Wavelets in Scattering Calculations p.1/43 What are Wavelets? Orthonormal basis functions.
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 informationMathematical Methods in Machine Learning
UMD, Spring 2016 Outline Lecture 2: Role of Directionality 1 Lecture 2: Role of Directionality Anisotropic Harmonic Analysis Harmonic analysis decomposes signals into simpler elements called analyzing
More informationPhysics 6303 Lecture 11 September 24, LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation
Physics 6303 Lecture September 24, 208 LAST TIME: Cylindrical coordinates, spherical coordinates, and Legendre s equation, l l l l l l. Consider problems that are no axisymmetric; i.e., the potential depends
More informationChapter 7 Wavelets and Multiresolution Processing. Subband coding Quadrature mirror filtering Pyramid image processing
Chapter 7 Wavelets and Multiresolution Processing Wavelet transform vs Fourier transform Basis functions are small waves called wavelet with different frequency and limited duration Multiresolution theory:
More informationHarmonic Analysis: from Fourier to Haar. María Cristina Pereyra Lesley A. Ward
Harmonic Analysis: from Fourier to Haar María Cristina Pereyra Lesley A. Ward Department of Mathematics and Statistics, MSC03 2150, 1 University of New Mexico, Albuquerque, NM 87131-0001, USA E-mail address:
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 informationPARAMETRIC OPTIMIZATION OF BIORTHOGONAL WAVELETS AND FILTERBANKS VIA PSEUDOFRAMES FOR SUBSPACES
PARAMETRIC OPTIMIZATION OF BIORTHOGONAL WAVELETS AND FILTERBANKS VIA PSEUDOFRAMES FOR SUBSPACES SHIDONG LI AND MICHAEL HOFFMAN Abstract. We present constructions of biorthogonal wavelets and associated
More informationLecture 7 Multiresolution Analysis
David Walnut Department of Mathematical Sciences George Mason University Fairfax, VA USA Chapman Lectures, Chapman University, Orange, CA Outline Definition of MRA in one dimension Finding the wavelet
More informationWAVELET RECONSTRUCTION OF NONLINEAR DYNAMICS
International Journal of Bifurcation and Chaos, Vol. 8, No. 11 (1998) 2191 2201 c World Scientific Publishing Company WAVELET RECONSTRUCTION OF NONLINEAR DYNAMICS DAVID ALLINGHAM, MATTHEW WEST and ALISTAIR
More informationApplications of filter coefficients and wavelets parametrized by moments
Radon Series Comp. Appl. Math 3, 9 24 c de Gruyter 27 Applications of filter coefficients and wavelets parametrized by moments Georg Regensburger Key words. Orthonormal wavelets, parametrized filter coefficients,
More informationResearch Statement. James Bremer Department of Mathematics, University of California, Davis
Research Statement James Bremer Department of Mathematics, University of California, Davis Email: bremer@math.ucdavis.edu Webpage: https.math.ucdavis.edu/ bremer I work in the field of numerical analysis,
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 informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationSome Aspects of Band-limited Extrapolations
October 3, 2009 I. Introduction Definition of Fourier transform: F [f ](ω) := ˆf (ω) := + f (t)e iωt dt, ω R (1) I. Introduction Definition of Fourier transform: F [f ](ω) := ˆf (ω) := + f (t)e iωt dt,
More information