Spatial-spectral concentration on the ball

Post by Jason McEwen

In cosmology, observations are made over the celestial sphere, giving rise to observations defined on the two-dimensional sphere, i.e. a spherical surface.  If depth information is also available (i.e. redshift), then observations are defined on the three-dimensional ball, i.e. on a spherical volume.  For example, observations of the cosmic microwave background (CMB) are made on the sphere, while observations of the galaxy distribution that traces the large-scale structure (LSS) are made on the ball.

Often, however, observations cannot be made over the full sky.  For example, we must look through our galaxy, which contaminates observations. Foreground contamination can sometimes be modelled and reduced, however regions of significant contamination must be removed altogether. In addition, telescopes often simply cannot see the entire sky.

Dealing with partial-sky coverage can be difficult.  Wavelets are a powerful method to do this due to their dual spatial and spectral localisation properties.  Alternatively, one can build a basis concentrated in the observed region.  This is a well-studied problem in signal processing and is known as the Slepian spatial-spectral concentration problem.  Although this problem has been solved in the Euclidean setting, and also on the sphere, it has not been solved on the ball.  We recently submitted a paper solving the Slepian spatial-spectral concentration problem on the ball.

The abstract of our submission is reproduced below and you can find the full paper on the arXiv.

“We formulate and solve the Slepian spatial-spectral concentration problem on the three-dimensional ball. Both the standard Fourier-Bessel and also the Fourier-Laguerre spectral domains are considered since the latter exhibits a number of practical advantages (spectral decoupling and exact computation). The Slepian spatial and spectral concentration problems are formulated as eigenvalue problems, the eigenfunctions of which form an orthogonal family of concentrated functions. Equivalence between the spatial and spectral problems is shown. The spherical Shannon number on the ball is derived, which acts as the analog of the space-bandwidth product in the Euclidean setting, giving an estimate of the number of concentrated eigenfunctions and thus the dimension of the space of functions that can be concentrated in both the spatial and spectral domains simultaneously. Various symmetries of the spatial region are considered that reduce considerably the computational burden of recovering eigenfunctions, either by decoupling the problem into smaller subproblems or by affording analytic calculations. The family of concentrated eigenfunctions forms a Slepian basis that can be used be represent concentrated signals efficiently. We illustrate our results with numerical examples and show that the Slepian basis indeeds permits a sparse representation of concentrated signals.”

In additional to considering the standard Fourier-Bessel basis on the ball, we also consider the Fourier-Laguerre basis, which exhibits a number of practical advantages.  The first few Slepian functions concentrated within an example region are shown in the following plots for each basis on the ball.


Fourier-Bessel spatially concentrated Slepian functions


Fourier-Laguerre spatially concentrated Slepian functions



Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

%d bloggers like this: