Localized bases for kernel spaces on the unit sphere

From MaRDI portal



Abstract: Approximation/interpolation from spaces of positive definite or conditionally positive definite kernels is an increasingly popular tool for the analysis and synthesis of scattered data, and is central to many meshless methods. For a set of N scattered sites, the standard basis for such a space utilizes N emph{globally} supported kernels; computing with it is prohibitively expensive for large N. Easily computable, well-localized bases, with "small-footprint" basis elements - i.e., elements using only a small number of kernels -- have been unavailable. Working on sphere, with focus on the restricted surface spline kernels (e.g. the thin-plate splines restricted to the sphere), we construct easily computable, spatially well-localized, small-footprint, robust bases for the associated kernel spaces. Our theory predicts that each element of the local basis is constructed by using a combination of only mathcalO((logN)2) kernels, which makes the construction computationally cheap. We prove that the new basis is Lp stable and satisfies polynomial decay estimates that are stationary with respect to the density of the data sites, and we present a quasi-interpolation scheme that provides optimal Lp approximation orders. Although our focus is on mathbbS2, much of the theory applies to other manifolds - mathbbSd, the rotation group, and so on. Finally, we construct algorithms to implement these schemes and use them to conduct numerical experiments, which validate our theory for interpolation problems on mathbbS2 involving over one hundred fifty thousand data sites.


For important special cases of radial basis functions, such as thin-plate splines and other surface splines of the Duchon type, and with arguments restricted to the unit sphere, the authors derive stable ``small footprint basis functions. Small footprint means that they are constructed with few coefficients, so their \(\ell_0\)-``norm is small. These function space bases are quite easily computable and stable with respect to \(L_p\) norms/spaces. They are local -- not usually compactly supported which is normally not possible here, but quickly decaying -- and admit useful approximations by quasi-interpolation using them. The approach is working via Lagrange functions and as mentioned for the two-dimensional unit-sphere (in parts generalizable as remarked by the authors). Numerical examples are included, too, as is a discussion how to use preconditioning of the Gram (interpolation) matrix in this context.




Cited in
(38)








This page was built for publication: Localized bases for kernel spaces on the unit sphere

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2870616)