Wavelets for approximate Fourier transform and data compression

Author

Guo, Haitao

Date

1997

Advisor

Burrus, C. Sidney

Degree

Doctor of Philosophy

Abstract

This dissertation has two parts. In the first part, we develop a wavelet-based fast approximate Fourier transform algorithm. The second part is devoted to the developments of several wavelet-based data compression techniques for image and seismic data.
We propose an algorithm that uses the discrete wavelet transform (DWT) as a tool to compute the discrete Fourier transform (DFT). The classical Cooley-Tukey FFT is shown to be a special case of the proposed algorithm when the wavelets in use are trivial. The main advantage of our algorithm is that the good time and frequency localization of wavelets can be exploited to approximate the Fourier transform for many classes of signals, resulting in much less computation. Thus the new algorithm provides an efficient complexity versus accuracy tradeoff. When approximations are allowed, under certain sparsity conditions, the algorithm can achieve linear complexity, i.e. O(N). The proposed algorithm also has built-in noise reduction capability.
For waveform and image compression, we propose a novel scheme using the recently developed Burrows-Wheeler transform (BWT). We show that the discrete wavelet transform (DWT) should be used before the Burrows-Wheeler transform to improve the compression performance for many natural signals and images. We demonstrate that the simple concatenation of the DWT and BWT coding performs comparably as the embedded zerotree wavelet (EZW) compression for images. Various techniques that significantly improve the performance of our compression scheme are also discussed.
The phase information is crucial for seismic data processing. However, traditional compression schemes do not pay special attention to preserving the phase of the seismic data, resulting in the loss of critical information. We propose a lossy compression method that preserves the phase as much as possible. The method is based on the self-adjusting wavelet transform that adapts to the locations of the significant signal components. The elegant method of embedded zerotree wavelet compression is modified and incorporated into our compression scheme. Our method can be applied to both one dimensional seismic signals and two dimensional seismic images.