Regular Article: Computing Fourier Transforms and Convolutions on the 2-Sphere
Advances in Applied Mathematics
Journal of Multivariate Analysis
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,
Handbook of Mathematical Functions, With Formulas, Graphs, and Mathematical Tables,
Consistent and efficient density estimation
ICCSA'03 Proceedings of the 2003 international conference on Computational science and its applications: PartI
Convolution on the n-sphere with application to PDF modeling
IEEE Transactions on Signal Processing
Geometric structures arising from kernel density estimation on Riemannian manifolds
Journal of Multivariate Analysis
Hi-index | 0.00 |
This paper is on density estimation on the 2-sphere, S2, using the orthogonal series estimator corresponding to spherical harmonics. In the standard approach of truncating the Fourier series of the empirical density, the Fourier transform is replaced with a version of the discrete fast spherical Fourier transform, as developed by Driscoll and Healy. The fast transform only applies to quantitative data on a regular grid. We will apply a kernel operator to the empirical density, to produce a function whose values at the vertices of such a grid will be the basis for the density estimation. The proposed estimation procedure also contains a deconvolution step, in order to reduce the bias introduced by the initial kernel operator. The main issue is to find necessary conditions on the involved discretization and the bandwidth of the kernel operator, to preserve the rate of convergence that can be achieved by the usual computationally intensive Fourier transform. Density estimation is considered in L2(S2) and more generally in Sobolev spaces Hυ(S2), any υ ≥ 0, with the regularity assumption that the probability density to be estimated belongs to Hs(S2) for some s υ. The proposed technique to estimate the Fourier transform of an unknown density keeps computing cost down to order O(n), where n denotes the sample size.