This paper is concerned with the construction of biorthogonal multiresolution analyses on [0; 1] such that the corresponding wavelets realize any desired order of moment conditions throughout the interval. Our starting point is the family of biorthogonal pairs consisting of cardinal B-splines and compactly supported dual generators on IR developed by Cohen, Daubechies and Feauveau. In contrast to previous investigations we preserve the full degree of polynomial reproduction also for the dual multiresolution and prove in general that the corresponding modifications of dual generators near the end points of the interval still permit the biorthogonalization of the resulting bases. The subsequent construction of compactly supported biorthogonal wavelets is based on the concept of stable completions. As a first step we derive an initial decomposition of the spline spaces where the complement spaces between two successive levels are spanned by compactly supported splines which form uniformly...

This paper is concerned with the construction of biorthogonal wavelet bases defined on a union of parametric images of the unit d-cube. These bases are to satisfy certain requirements imposed by applications to a class of operator equations acting on such domains. This covers also elliptic boundary value problems although this study is primarily motivated by our previous analysis of wavelet methods for pseudo-differential equations with special emphasis on boundary integral equations. In this case it is natural to model the boundary surface as a union of parametric images of the unit cube. It will be shown how to construct wavelet bases on the surface which are composed of wavelet bases defined on each surface patch. Here the relevant properties are the validity of norm equivalences in certain ranges of Sobolev scales as well as appropriate moment conditions. Key Words: Biorthogonal wavelets, norm equivalences, boundary element methods, composite multiresolution, multiscale methods fo...

The potential of wavelets as a discretization tool for the numerical treatment of operator equations hinges on the validity of norm equivalences for Besov or Sobolev spaces in terms of weighted sequence norms of wavelet expansion coefficients and on certain cancellation properties. These features are crucial for the construction of optimal preconditioners, for matrix compression based on sparse representations of functions and operators as well as for the design and analysis of adaptive solvers. However, for realistic domain geometries the relevant properties of wavelet bases could so far only be realized to a limited extent. This paper is concerned with concepts that aim at expanding the applicability of wavelet schemes in this sense. The central issue is to construct wavelet bases with the desired properties on manifolds which can be represented as the disjoint union of smooth parametric images of the standard cube. The approach considered here is conceptually different though from o...

Recently, a framework for multiscale stochastic modeling was introduced based on coarse-to-fine scale-recursive dynamics defined on trees. This model class has some attractive characteristics which lead to extremely efficient, statistically optimal signal and image processing algorithms. In this paper, we show that this model class is also quite rich. In particular, we describe how 1-D Markov processes and 2-D Markov random fields (MRF’s) can be represented within this framework. The recursive structure of 1-D Markov processes makes them simple to analyze, and generally leads to computationally efficient algorithms for statistical inference. On the other hand, 2-D MRF’s are well known to be very difficult to analyze due to their noncausal structure, and thus their use typically leads to computationally intensive algorithms for smoothing and parameter identification. In contrast, our multiscale representations are based on scale-recursive models and thus lead naturally to scale-recursive algorithms, which can be substantially more efficient computationally than those associated with MRF models. In 1-D, the multiscale representation is a generalization of the midpoint deflection construction of Brownian motion. The representation of 2-D MRF’s is based on a further generalization to a “midline ” deflection construction. The exact representations of 2-D MRF’s are used to motivate a class of multiscale approximate MRF models based on one-dimensional wavelet transforms. We demonstrate the use of these latter models in the context of texture representation and, in particular, we show how they can be used as approximations for or alternatives to well-known MRF texture models.

. The Wavelet Element Method (WEM) combines biorthogonal wavelet systems with the philosophy of Spectral Element Methods in order to obtain a biorthogonal wavelet system on fairly general bounded domains in some IR n . The domain of interest is split into subdomains which are mapped to a simple reference domain, here n--dimensional cubes. Thus, one has to construct appropriate biorthogonal wavelets on the reference domain such that mapping them to each subdomain and matching along the interfaces leads to a wavelet system on the domain. In this paper we use adapted biorthogonal wavelet systems on the interval in such a way, that tensor products of these functions can be used for the construction of wavelet bases on the reference domain. We describe the matching procedure in any dimension n in order to impose continuity and prove that it leads to a construction of a biorthogonal wavelet system on the domain. These wavelet systems characterize Sobolev spaces measuring both piecewise and...

We describe an algorithm to estimate a discrete signal from its noisy observation, using a library of orthonormal bases (consisting of various wavelets, wavelet packets, and local trigonometric bases) and the information-theoretic criterion called minimum description length (MDL). The key to effective random noise suppression is that the signal component in the data may be represented efficiently by one or more of the bases in the library, whereas the noise component cannot be represented efficiently by any basis in the library. The MDL criterion gives the best compromise between the fidelity of the estimation result to the data (noise suppression) and the efficiency of the representation of the estimated signal (signal compression): it selects the "best" basis and the "best" number of terms to be retained out of various bases in the library in an objective manner. Because of the use of the MDL criterion, our algorithm is free from any parameter setting or subjective judgments. This ...

Abstract. This note is a very basic introduction to wavelets. It starts with an orthogonal basis of piecewise constant functions, constructed by dilation and translation. The “wavelet transform ” maps each f(x) to its coefficients with respect to this basis. The mathematics is simple and the transform is fast (faster than the Fast Fourier Transform, which we briefly explain), but approximation by piecewise constants is poor. To improve this first wavelet, we are led to dilation equations and their unusual solutions. Higher-order wavelets are constructed, and it is surprisingly quick to compute with them — always indirectly and recursively. We comment informally on the contest between these transforms in signal processing, especially for video and image compression (including high-definition television). So far the Fourier Transform — or its 8 by 8 windowed version, the Discrete Cosine Transform — is often chosen. But wavelets are already competitive, and they are ahead for fingerprints. We present a sample of this developing theory. 1. The Haar wavelet To explain wavelets we start with an example. It has every property we hope for, except one. If that one defect is accepted, the construction is simple and the computations are fast. By trying to remove the defect, we are led to dilation equations and recursively defined functions and a small world of fascinating new problems — many still unsolved. A sensible person would stop after the first wavelet, but fortunately mathematics goes on. The basic example is easier to draw than to describe: Figure 1. Scaling function φ(x), wavelet W(x), and the next level of detail.

We present a multiresolution construction for efficiently computing, compressing and applying large powers of operators that have high powers with low numerical rank. This allows the fast computation of functions of the operator, notably the associated Green’s function, in compressed form, and their fast application. Classes of operators satisfying these conditions include diffusion-like operators, in any dimension, on manifolds, graphs, and in non-homogeneous media. In this case our construction can be viewed as a far-reaching generalization of Fast Multipole Methods, achieved through a different point of view, and of the non-standard wavelet representation of Calderón-Zygmund and pseudodifferential operators, achieved through a different multiresolution analysis adapted to the operator. We show how the dyadic powers of an operator can be used to induce a multiresolution analysis, as in classical Littlewood-Paley and wavelet theory, and we show how to construct, with fast and stable algorithms, scaling function and wavelet bases associated to this multiresolution analysis, and the corresponding downsampling operators, and use them to compress the corresponding powers of the operator. This allows to extend multiscale signal processing to general spaces (such as manifolds and graphs) in a very natural way, with corresponding fast algorithms.

CT 06520 0

We show how importance-driven refinement and a wavelet basis can be combined to provide an efficient solution to the global illumination problem with glossy and diffuse reflections. Importance is used to focus the computation on the interactions having the greatest impact on the visible solution. Wavelets are used to provide an efficient representation of radiance, importance, and the transport operator. We discuss a number of choices that must be made when constructing a finite element algorithm for glossy global illumination. Our algorithm is based on the standard wavelet decomposition of the transport operator and makes use of a four-dimensional wavelet representation for spatially- and angularly-varying radiance distributions. We use a final gathering step to improve the visual quality of the solution. Features of our implementation include support for curved surfaces as well as texture-mapped anisotropic emission and reflection functions. 1 Introduction Radiosity algorithms assum...