Wavelets and multirate filter banks : theory, structure, design, and applications

Thesis (Ph. D.)--Massachusetts Institute of Technology, Dept. of Civil and Environmental Engineering, 2004.Includes bibliographical references (p. 219-230) and index.Wavelets and filter banks have revolutionized signal processing with their ability to process data at multiple temporal and spatial resolutions. Fundamentally, continuous-time wavelets are governed by discrete-time filter banks with properties such as perfect reconstruction, linear phase and regularity. In this thesis, we study multi-channel filter bank factorization and parameterization strategies, which facilitate designs with specified properties that are enforced by the actual factorization structure. For M-channel filter banks (M =/> 2), we develop a complete factorization, M-channel lifting factorization, using simple ladder-like structures as predictions between channels to provide robust and efficient implementation; perfect reconstruction is structurally enforced, even under finite precision arithmetic and quantization of lifting coefficients. With lifting, optimal low-complexity integer wavelet transforms can thus be designed using a simple and fast algorithm that incorporates prescribed limits on hardware operations for power-constrained environments. As filter bank regularity is important for a variety of reasons, an aspect of particular interest is the structural imposition of regularity onto factorizations based on the dyadic form uvt. We derive the corresponding structural conditions for regularity, for which M-channel lifting factorization provides an essential parameterization. As a result, we are able to design filter banks that are exactly regular and amenable to fast implementations with perfect reconstruction, regardless of the choice of free parameters and possible finite precision effects. Further constraining u = v ensures regular orthogonal filter banks,(cont.) whereas a special dyadic form is developed that guarantees linear phase. We achieve superior coding gains within 0.1% of the optimum, and benchmarks conducted on image compression applications show clear improvements in perceptual and objective performance. We also consider the problem of completing an M-channel filter bank, given only its scaling filter. M-channel lifting factorization can efficiently complete such biorthogonal filter banks. On the other hand, an improved scheme for completing paraunitary filter banks is made possible by a novel order-one factorization which allows greater design flexibility, resulting in improved frequency selectivity and energy compaction over existing state of the art methods. In a dual setting, the technique can be applied to transmultiplexer design to achieve higher-rate data transmissions.by Ying-Jui Chen.Ph.D

Multi-scale edge-guided image gap restoration

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University London.The focus of this research work is the estimation of gaps (missing blocks) in digital images. To progress the research two main issues were identified as (1) the appropriate domains for image gap restoration and (2) the methodologies for gap interpolation. Multi-scale transforms provide an appropriate framework for gap restoration. The main advantages are transformations into a set of frequency and scales and the ability to progressively reduce the size of the gap to one sample wide at the transform apex. Two types of multi-scale transform were considered for comparative evaluation; 2-dimensional (2D) discrete cosines (DCT) pyramid and 2D discrete wavelets (DWT). For image gap estimation, a family of conventional weighted interpolators and directional edge-guided interpolators are developed and evaluated. Two types of edges were considered; ‘local’ edges or textures and ‘global’ edges such as the boundaries between objects or within/across patterns in the image. For local edge, or texture, modelling a number of methods were explored which aim to reconstruct a set of gradients across the restored gap as those computed from the known neighbourhood. These differential gradients are estimated along the geometrical vertical, horizontal and cross directions for each pixel of the gap. The edge-guided interpolators aim to operate on distinct regions confined within edge lines. For global edge-guided interpolation, two main methods explored are Sobel and Canny detectors. The latter provides improved edge detection. The combination and integration of different multi-scale domains, local edge interpolators, global edge-guided interpolators and iterative estimation of edges provided a variety of configurations that were comparatively explored and evaluated. For evaluation a set of images commonly used in the literature work were employed together with simulated regular and random image gaps at a variety of loss rate. The performance measures used are the peak signal to noise ratio (PSNR) and structure similarity index (SSIM). The results obtained are better than the state of the art reported in the literature

Finite element modeling of transient ultrasonic waves in linear viscoelastic media

Linear viscoelasticity offers a minimal framework within which to construct a causal model for wave propagation in absorptive media. Viscoelastic media are often described as media with memory, that is, the present state of stress is dependent on the present strain and the complete time history of strain weighted by time convolution with an appropriate time-dependent stress relaxation modulus. An axisymmetric, displacement based finite element method for modeling pulsed ultrasonic waves in linear, homogeneous viscoelastic media is developed that does not require storage of the complete time history of displacement at every node. This is accomplished by modeling the stress relaxation moduli as discrete or continuous spectra of decaying exponentials. The viscoelastic finite element method serves as a test bed for studying three inverse methods for recovering time dependent longitudinal moduli from pulsed ultrasonic waves transmitted through a slab of viscoelastic material with properties known a priori. Specifically, two existing inverse methods called propagator methods, denoted here as the two-slab method and slab-substitution method, are modeled and compared to show relative advantages and disadvantages of both. Both methods require attenuation and wave speed as a function of frequency derived from transmitted wave data for inversion and recovery of modulus data. Several different variables such as measurement location and source radius are varied to discern those variables that have greatest influence on accuracy of reconstructed moduli. It is found that an increase in source aperture radius causes the greatest improvement in modulus accuracy. Another novel inverse method known as wave splitting is applied to numerical data generated by the finite element test bed. Wave splitting requires a time-dependent transmission kernel for recovery of a viscoelastic modulus rather than frequency-dependent attenuation and wave speed. It is shown that in principle wave splitting can recover the material modulus with a data derived from a simulated ultrasonic experiment, but it is not as robust as the other two frequency-domain inverse methods studied. Its main drawback is that transmission kernel data required for inversion must be known for the same thickness of viscoelastic slab implying that pulses with relatively high center frequencies must be propagated through slabs whose thickness is only appropriate for low frequency measurement. Material attenuation quickly reduces transmitted waves at high frequencies to unacceptable low levels when propagated through thick slabs appropriate for pulses centered at lower frequencies. In general, the finite element method has been utilized as an effective tool for comparing alternative inverse methods

Nonlinear optics

Nonlinear light-matter interactions have been drawing attention of physicists since the 1960's. Quantum mechanics played a significant role in their description and helped to derive important formulas showing the dependence on the intensity of the electromagnetic field. High intensity light is able to generate second and third harmonics which translates to generation of electromagnetic field with multiples of the original frequency. In comparison with the linear behaviour of light, the nonlinear interactions are smaller in scale. This makes perturbation methods well suited for obtaining solutions to equations in nonlinear optics. In particular, the method of multiple scales is deployed in paper 3, where it is used to solve nonlinear dispersive wave equations. The key difference in our multiple scale solution is the linearity of the amplitude equation and a complex valued frequency of the mode. Despite the potential ill-posedness of the amplitude equation, the multiple scale solution remained a valid approximation of the solution to the original model. The results showed great potential of this method and its promising wider applications. Other methods use pseudo-spectral methods which require an orthogonal set of eigenfunctions (modes) used to create a substitute for the usual Fourier transform. This mode transform is only useful if it succeeds to represent target functions well. Papers 1 and 2 deal with investigating such modes called resonant and leaky modes and their ability to construct a mode transform. The modes in the first paper are the eigenvalues for a quantum mechanical system where an external radiation field is used to excite an electron trapped in an electrical potential. The findings show that the resonant mode expansion converges inside the potential independently of its depth. Equivalently, leaky modes are obtained in paper 2 which are in close relation to resonant modes. Here, the modes emerge from a system where a channel is introduced with transparent boundaries for simulation of one-directional optical beam propagation. Artificial index material is introduced outside the channel which gives rise to leaky modes associated with such artificial structure. The study is showing that leaky modes are well suited for function representation and thus solving the nonlinear version of this problem. In addition, the transparent boundary method turns out to be useful for spectral propagators such as the unidirectional pulse propagation equation in contrast to a perfectly matched layer

Proceedings of the Workshop on Applications of Distributed System Theory to the Control of Large Space Structures

Two general themes in the control of large space structures are addressed: control theory for distributed parameter systems and distributed control for systems requiring spatially-distributed multipoint sensing and actuation. Topics include modeling and control, stabilization, and estimation and identification

Bibliography of Lewis Research Center technical publications announced in 1987

This compilation of abstracts describes and indexes the technical reporting that resulted from the scientific and engineering work performed and managed by the Lewis Research Center in 1987. All the publications were announced in the 1987 issues of STAR (Scientific and Technical Aerospace Reports) and/or IAA (International Aerospace Abstracts). Included are research reports, journal articles, conference presentations, patents and patent applications, and theses

Laboratory directed research and development. FY 1995 progress report

This document presents an overview of Laboratory Directed Research and Development Programs at Los Alamos. The nine technical disciplines in which research is described include materials, engineering and base technologies, plasma, fluids, and particle beams, chemistry, mathematics and computational science, atmic and molecular physics, geoscience, space science, and astrophysics, nuclear and particle physics, and biosciences. Brief descriptions are provided in the above programs