Multilevel regularization of wavelet based fitting of scattered data ? some experiments

• Published in: Numerical algorithms 2005-07, Vol.39 (1-3), p.81-96
• Main Authors: Casta o, Daniel, Kunoth, Angela
• Format: Article
• Language: English
• Subjects: Algorithms; Geophysics; Least squares; Regularization; Sobolev space
• ISSN: 1017-1398; 1572-9265
• DOI: 10.1007/s11075-004-3622-0

In [6], an adaptive method to approximate unorganized clouds of points by smooth surfaces based on wavelets has been described. The general fitting algorithm operates on a coarse-to-fine basis. It selects on each refinement level in a first step a reduced number of wavelets which are appropriate to represent the features of the data set. In a second step, the fitting surface is constructed as the linear combination of the wavelets which minimizes the distance to the data in a least squares sense. This is followed by a thresholding procedure on the wavelet coefficients to discard those which are too small to contribute much to the surface representation. In this paper, we firstly generalize this strategy to a classically regularized least squares functional by adding a Sobolev norm, taking advantage of the capability of wavelets to characterize Sobolev spaces of even fractional order. After recalling the usual cross-validation technique to determine the involved smoothing parameters, some examples of fitting severely irregularly distributed data, synthetically produced and of geophysical origin, are presented. In order to reduce computational costs, we then introduce a multilevel generalized cross-validation technique which goes beyond the Sobolev formulation and exploits the hierarchical setting based on wavelets. We illustrate the performance of the new strategy on some geophysical data.