Draft:Slepian function

Slepian functions are a class of spatio-spectrally concentrated functions (that is, space- or time-concentrated while spectrally bandlimited, or spectral-band-concentrated while space- or time-limited) that form an orthogonal basis for bandlimited or spacelimited spaces. They are widely used as basis functions for approximation and in linear inverse problems, and as apodization tapers or window functions in quadratic problems of spectral density estimation.

Slepian functions exist in discrete and continuous varieties, in one, two, and three dimensions, in Cartesian and spherical geometry, on surfaces and in volumes, and in scalar, vector, and tensor forms.

Without reference to any of these particularities, let $f$ be a square-integrable function of physical space, and let ${\mathcal {H}}$ represent Fourier transformation, such that $F={\mathcal {H}}f$ and ${\mathcal {H}}^{-1}F=f$. Let the operators ${\mathcal {R}}$ and ${\mathcal {L}}$ project onto the space of spacelimited functions, ${\mathcal {S}}_{R}$, and the space of bandlimited functions, ${\mathcal {S}}_{L}$, respectively, whereby $R$ is an arbitrary nontrivial subregion of all of physical space, and $L$ an arbitrary nontrivial subregion of spectral space. Thus, the operator ${\mathcal {R}}$ acts to spacelimit, and the operator ${\mathcal {H}}^{-1}{\mathcal {L}}{\mathcal {H}}$ acts to bandlimit the function $f$.

Slepian's quadratic spectral concentration problem aims to maximize the concentration of spectral power to a target region $L$, for a function that is spatially limited to a target region $R$. Conversely, Slepian's spatial concentration problem maximizes the spatial concentration to $R$ of a function bandlimited to $L$. Using $\langle \cdot ,\cdot \rangle$ for the inner product both in the space and the spectral domains, both problems are stated equivalently using Rayleigh quotients in the form

$\lambda ={\frac {\langle {\mathcal {R}}{\mathcal {H}}^{-1}{\mathcal {L}}\,F,{\mathcal {R}}{\mathcal {H}}^{-1}{\mathcal {L}}\,F\,\rangle }{\langle {\mathcal {H}}^{-1}{\mathcal {L}}\,F,{\mathcal {H}}^{-1}{\mathcal {L}}\,F\,\rangle }}={\frac {\langle {\mathcal {L}}{\mathcal {H}}{\mathcal {R}}f,{\mathcal {L}}{\mathcal {H}}{\mathcal {R}}f\rangle }{\langle {\mathcal {H}}{\mathcal {R}}f,{\mathcal {H}}{\mathcal {R}}f\rangle }}={\mbox{maximum}}.$

The equivalent spectral-domain and spatial-domain eigenvalue equations are

$({\mathcal {L}}{\mathcal {H}}{\mathcal {R}}{\mathcal {H}}^{-1\!}{\mathcal {L}})({\mathcal {L}}\,F\,)=\lambda ({\mathcal {L}}\,F\,)$ and $({\mathcal {R}}{\mathcal {H}}^{-1\!}{\mathcal {L}}{\mathcal {H}}{\mathcal {R}})({\mathcal {R}}f)=\lambda ({\mathcal {R}}f),$

given that ${\mathcal {H}}$ and ${\mathcal {H}}^{-1}$ are each others' adjoints, and that ${\mathcal {R}}$ and ${\mathcal {L}}$ are self-adjoint and idempotent.

The Slepian functions are solutions to either of these types of equations with positive-definite kernels, that is, they are bandlimited functions $G={\mathcal {L}}F$, concentrated to the spatial domain within $R$, or spacelimited functions of the form $h={\mathcal {R}}f$, concentrated to the spectral domain within $L$.

Scalar Slepian functions in one dimension

Let $g(t)$ and its Fourier transform $G(\omega )$ be strictly bandlimited in angular frequency between $[-W,W]$. Attempting to concentrate $g(t)$ in the time domain, to be contained within the time interval $[-T,T]$, amounts to maximizing

$\lambda ={\frac {\int _{-T}^{T}g^{2}(t)\,dt}{\int _{-\infty }^{\infty }g^{2}(t)\,dt}}={\text{maximum}},$

which is equivalent to solving either, in the frequency domain, the convolutional integral eigenvalue (Fredholm) equation

$\int _{-W}^{W}D_{T}(\omega ,\omega ')G(\omega ')\,d\omega '=\lambda G(\omega ),\qquad D_{T}(\omega ,\omega ')={\frac {\sin T(\omega -\omega ')}{\pi (\omega -\omega ')}},\qquad |\omega |\leq W$,

or the time- or space-domain version

$\int _{-T}^{T}D_{W}(t,t')g(t')\,dt'=\lambda g(t),\qquad D_{W}(t,t')={\frac {\sin W(t-t')}{\pi (t-t')}}=(2\pi )^{-1}\int _{-W}^{W}e^{i\omega (t-t')}d\omega ,\qquad t\in \mathbb {R}$.

(a) Slepian functions in the time domain. (b) Slepian functions in the frequency domain. Shown is the square of the absolute value of the Fourier transform of the Slepian functions shown in (a). (c) Concentration factors associated with the successive Slepian functions shown in (a). (d) Cumulative energy by summation the square of the Slepian functions shown in (a).

Either of these can be transformed and rescaled to the dimensionless

$\int _{-1}^{1}D(x,x')g(x')\,dx'=\lambda g(x),\qquad D(x,x')={\frac {\sin TW(x-x')}{\pi (x-x')}}$.

The trace of the positive definite kernel is the sum of the infinite number of real and positive eigenvalues,

$N=\sum _{\alpha =1}^{\infty }\lambda _{\alpha }=\int _{-1}^{1}D(x,x')\,dx={\frac {2TW}{\pi }},$

that is, the area of the concentration domain in time-frequency space (a time-bandwidth product).

One-dimensional scalar Slepian functions are the workhorse of the Thomson multitaper method of spectral density estimation.

Scalar Slepian functions in two Cartesian dimensions

We use $g(\mathbf {x} )$ and its Fourier transform $G(\mathbf {k} )$ to denote a function that is strictly bandlimited to ${\mathcal {K}}$, an arbitrary subregion of the spectral space of spatial wave vectors. Seeking to concentrate $g(\mathbf {x} )$ into a finite spatial region $R\in \mathbb {R} ^{2}$, of area $A$, we must find the unknown functions for which

$\lambda ={\frac {\int _{R}g^{2}(\mathbf {x} )\,d\mathbf {x} }{\int _{-\infty }^{\infty }g^{2}(\mathbf {x} )\,d\mathbf {x} }}={\mbox{maximum}}.$

Maximizing this Rayleigh quotient requires solving the Fredholm integral equation

$\int _{\mathcal {K}}D_{R}(\mathbf {k} ,\mathbf {k} ')\,G(\mathbf {k} ')\,d\mathbf {k} '=\lambda G(\mathbf {k} ),\qquad D_{R}(\mathbf {k} ,\mathbf {k} ')=(2\pi )^{-2}\int _{R}e^{i(\mathbf {k} '-\mathbf {k} )\cdot \mathbf {x} }\,d\mathbf {x} ,\qquad \mathbf {k} \in {\mathcal {K}}.$

The corresponding problem in the spatial domain is

$\int _{R}\!D_{\mathcal {K}}(\mathbf {x} ,\mathbf {x} ')\,g(\mathbf {x} ')\,d\mathbf {x} '=\lambda g(\mathbf {x} ),\qquad D_{\mathcal {K}}(\mathbf {x} ,\mathbf {x} ')=(2\pi )^{-2}\int _{\mathcal {K}}e^{i\mathbf {k} \cdot (\mathbf {x} -\mathbf {x} ')}\,d\mathbf {k} ,\qquad \mathbf {x} \in \mathbb {R} ^{2}.$

Concentration to the disk-shaped spectral band ${\mathcal {K}}=\{\mathbf {k} :\|\mathbf {k} \|\leq K\}$ allows us to rewrite the spatial kernel as

Slepian functions concentrated to a cat-like spatial (top row; rank $\alpha$ and concentration eigenvalue $\lambda$) and a duck-like spectral domain (bottom row; shown is the square of the absolute value of the Fourier transform of the functions shown in the top row).

$D_{\mathcal {K}}(\mathbf {x} ,\mathbf {x} ')={\frac {KJ_{1}(K\|\mathbf {x} -\mathbf {x} '\|)}{2\pi \|\mathbf {x} -\mathbf {x} '\|}},$

with $J_{1}$ a Bessel function of the first kind, from which we may derive that

$N=\sum _{\alpha =1}^{\infty }\lambda _{\alpha }=\int _{R}D(\mathbf {x} ,\mathbf {x} ')\,d\mathbf {x} =K^{2}{\frac {A}{4\pi }},$

in other words, again the area of the concentration domain in space-frequency space (a space-bandwidth product).

Scalar Slepian functions on the surface of a sphere

We denote $g(\mathbf {\hat {r}} )$ a function on the unit sphere $\Omega$ and its spherical harmonic transform coefficient $g_{lm}$ at the degree $l$ and order $m$, respectively, and we consider bandlimitation to spherical harmonic degree $L$, that is, $g\in {\mathcal {S}}_{L}$. Maximizing the quadratic energy ratio within the spatial subdomain $R\subset \Omega$ via

$\lambda ={\frac {\int _{R}g^{2}(\mathbf {\hat {r}} )\,d\Omega }{\int _{\Omega }g^{2}(\mathbf {\hat {r}} )\,d\Omega }}={\mbox{maximum}}$

amounts in the spectral domain to solving the algebraic eigenvalue equation

$\sum _{l'=0}^{L}\sum _{m'=-l'}^{l'}D_{lm,l'm'}g_{l'm'}=\lambda g_{lm},\qquad D_{lm,l'm'}=\int _{R}Y_{lm}(\mathbf {\hat {r}} )Y_{l'm'}(\mathbf {\hat {r}} )\,d\Omega$,

with $Y_{lm}$ the spherical harmonic at degree $l$ and order $m$. The equivalent spatial-domain equation, $\int _{R}D(\mathbf {\hat {r}} ,\mathbf {\hat {r}} ')g(\mathbf {\hat {r}} )\,d\Omega =\lambda g(\mathbf {\hat {r}} ),\qquad D(\mathbf {\hat {r}} ,\mathbf {\hat {r}} ')=\sum _{l=0}^{L}\sum _{m=-l}^{l}Y_{lm}(\mathbf {\hat {r}} )Y_{lm}(\mathbf {\hat {r}} ')=\sum _{l=0}^{L}\left({\frac {2l+1}{4\pi }}\right)P_{l}(\mathbf {\hat {r}} \cdot \mathbf {\hat {r}} '),$

Spherical Slepian functions of spherical-harmonic bandwidth 18, and of spherical-harmonic order 0 (that is, only made of zonal spherical harmonics), either very well (top row) or very poorly (bottom row) concentrated, as indicated by the concentration ratio $\lambda$ to the North-polar cap of opening angle 40$^{\circ }$.

is a homogeneous Fredholm integral equation of the second kind, with a finite-rank, symmetric, separable kernel.

The last equality is a consequence of the spherical harmonic addition theorem which involves $P_{l}$, the Legendre polynomial. The trace of this kernel is given by

$N=\sum _{\alpha =1}^{(L+1)^{2}}\lambda _{\alpha }=\int _{R}D(\mathbf {\hat {r}} ,\mathbf {\hat {r}} )\,d\Omega =\sum _{l=0}^{L}\sum _{m=-l}^{l}D_{lm,lm}=(L+1)^{2}{\frac {A}{4\pi }}$,

that is, once again a space-bandwidth product, of the dimension of ${\mathcal {S}}_{L}$ and the fractional area of $R$ on the unit sphere $\Omega$, namely $A/(4\pi )$.

References

I. Daubechies. Ten Lectures on Wavelets. SIAM, 1992, ISBN 0-89871-274-2

V. Michel. Spherical Slepian Functions, in Lectures on Constructive Approximation. Birkhäuser, 2012, doi:10.1007/978-0-8176-8403-7_8

V. Michel, A. Plattner, and K. Seibert. A unified approach to scalar, vector, and tensor Slepian functions on the sphere and their construction by a commuting operator. Analysis and Applications, 2022, doi:10.1142/S0219530521500317

C. T. Mullis and L. L. Scharf. Quadratic estimators of the power spectrum, in Advances in Spectrum Analysis and Array Processing, Vol. 1, chap. 1, pp. 1–57, ed. S. Haykin. Prentice-Hall, 1991, ISBN 978-0130074447

W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery. Numerical Recipes. The Art of Scientific Computing (Third Edition). Cambridge, 2007, ISBN 978-0-521-88068-8

F. J. Simons. Slepian functions and their use in signal estimation and spectral analysis. Handbook of Geomathematics, 2010, doi:10.1007/978-3-642-01546-5_30

F. J. Simons, M. A. Wieczorek and F. A. Dahlen. Spatiospectral concentration on a sphere. SIAM Review, 2006, doi:10.1137/S0036144504445765.

F. J. Simons and D. V. Wang. Spatiospectral concentration in the Cartesian plane. Int. J. Geomath, 2011, doi:10.1007/s13137-011-0016-z.

F. J. Simons and A. Plattner. Scalar and vector Slepian functions, spherical signal estimation and spectral analysis. Handbook of Geomathematics, 2015, doi:10.1007/978-3-642-54551-1_30

D. J. Thomson. Spectrum estimation and harmonic analysis. Proceedings of the IEEE, 1982, doi:10.1109/PROC.1982.12433