Fast directional correlation on the sphere with steerable filters

Abstract

A fast algorithm is developed for the directional correlation of scalar band-limited signals and band-limited steerable filters on the sphere. The asymptotic complexity associated to it through simple quadrature is of order O(L^5), where 2L stands for the square-root of the number of sampling points on the sphere, also setting a band limit L for the signals and filters considered. The filter steerability allows to compute the directional correlation uniquely in terms of direct and inverse scalar spherical harmonics transforms, which drive the overall asymptotic complexity. The separation of variables technique for the scalar spherical harmonics transform produces an O(L^3) algorithm independently of the pixelization. On equi-angular pixelizations, a sampling theorem introduced by Driscoll and Healy implies the exactness of the algorithm. The equi-angular and HEALPix implementations are compared in terms of memory requirements, computation times, and numerical stability. The computation times for the scalar transform, and hence for the directional correlation, of maps of several megapixels on the sphere (L~10^3) are reduced from years to tens of seconds in both implementations on a single standard computer. These generic results for the scale-space signal processing on the sphere are specifically developed in the perspective of the wavelet analysis of the cosmic microwave background (CMB) temperature (T) and polarization (E and B) maps of the WMAP and Planck experiments. As an illustration, we consider the computation of the wavelet coefficients of a simulated temperature map of several megapixels with the second Gaussian derivative wavelet.

Version accepted in APJ. 14 pages, 2 figures, Revtex4 (emulateapj). Changes include (a) a presentation of the algorithm as directly built on blocks of standard spherical harmonics transforms, (b) a comparison between the HEALPix and equi-angular implementations

References in corpus (9)

• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters
• HEALPix -- a Framework for High Resolution Discretization, and Fast Analysis of Data Distributed on the Sphere
• Detection of non-Gaussianity in the WMAP 1-year data using spherical wavelets
• First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Interpretation of the TT and TE Angular Power Spectrum Peaks
• Gauss--Legendre Sky Pixelization (GLESP) for CMB maps
• A high-significance detection of non-Gaussianity in the WMAP 1-year data using directional spherical wavelets
• Wavelets and WMAP non-Gaussianity
• Correspondence principle between spherical and Euclidean wavelets
• Full-sky correlation functions for CMB experiments with asymmetric window functions

Cited by in corpus (23)

• A novel sampling theorem on the sphere
• Exact reconstruction with directional wavelets on the sphere
• A Panorama on Multiscale Geometric Representations, Intertwining Spatial, Directional and Frequency Selectivity
• Cosmological applications of a wavelet analysis on the sphere
• Exact Wavelets on the Ball
• Alignment and signed-intensity anomalies in WMAP data
• S2LET: A code to perform fast wavelet analysis on the sphere
• Fast spin +-2 spherical harmonics transforms and application in cosmology
• Localisation of directional scale-discretised wavelets on the sphere
• Probing dark energy with steerable wavelets through correlation of WMAP and NVSS local morphological measures
• Complex data processing: fast wavelet analysis on the sphere
• Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression
• On the computation of directional scale-discretized wavelet transforms on the sphere
• Non-Gaussianity analysis on local morphological measures of WMAP data
• Directional spin wavelets on the sphere
• A pseudospectral matrix method for time-dependent tensor fields on a spherical shell
• A novel sampling theorem on the rotation group
• Wavelet reconstruction of E and B modes for CMB polarisation and cosmic shear analyses
• Fast directional spatially localized spherical harmonic transform
• Parallel Spherical Harmonic Transforms on heterogeneous architectures (GPUs/multi-core CPUs)
• Compressed convolution
• Cosmic Microwave Background anisotropies: the power spectrum and beyond
• A phase-sensitive method for filtering on the sphere