For the representation of spin-$s$ band-limited functions on the sphere, we propose a sampling scheme with optimal number of samples equal to the number of degrees of freedom of the function in harmonic space. In comparison to the existing sampling designs, which require ${\sim}2L^2$ samples for the representation of spin-$s$ functions band-limited at $L$, the proposed scheme requires $N_o=L^2-s^2$ samples for the accurate computation of the spin-$s$ spherical harmonic transform~($s$-SHT). For the proposed sampling scheme, we also develop a method to compute the $s$-SHT. We place the samples in our design scheme such that the matrices involved in the computation of $s$-SHT are well-conditioned. We also present a multi-pass $s$-SHT to improve the accuracy of the transform. We also show the proposed sampling design exhibits superior geometrical properties compared to existing equiangular and Gauss-Legendre sampling schemes, and enables accurate computation of the $s$-SHT corroborated through numerical experiments.

In this paper, we develop an analytical formulation for the Slepian spatial-spectral concentration problem on the sphere for a limited colatitude-longitude spatial region on the sphere, defined as the Cartesian product of a range of positive colatitudes and longitudes. The solution of the Slepian problem is a set of functions that are optimally concentrated and orthogonal within a spatial or spectral region. These properties make them useful for applications where measurements are taken within a spatially limited region of the sphere and/or a signal is only to be analyzed within a region of the sphere. To support localized spectral/spatial analysis, and estimation and sparse representation of localized data in these applications, we exploit the expansion of spherical harmonics in the complex exponential basis to develop an analytical formulation for the Slepian concentration problem for a limited colatitude-longitude spatial region. We also extend the analytical formulation for spatial regions that are comprised of a union of rotated limited colatitude-longitude subregions. By exploiting various symmetries of the proposed formulation, we design a computationally efficient algorithm for the implementation of the proposed analytical formulation. Such a reduction in computation time is demonstrated through numerical experiments. We present illustrations of our results with the help of numerical examples and show that the representation of a spatially concentrated signal is indeed sparse in the Slepian basis.

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 this paper, we present an optimal filter for the enhancement or estimation of signals on the 2-sphere corrupted by noise, when both the signal and noise are realizations of anisotropic processes on the 2-sphere. The estimation of such a signal in the spatial or spectral domain separately can be shown to be inadequate. Therefore, we develop an optimal filter in the joint spatio-spectral domain by using a framework recently presented in the literature --- the spatially localized spherical harmonic transform --- enabling such processing. Filtering of a signal in the spatio-spectral domain facilitates taking into account anisotropic properties of both the signal and noise processes. The proposed spatio-spectral filtering is optimal under the mean-square error criterion. The capability of the proposed filtering framework is demonstrated with by an example to estimate a signal corrupted by an anisotropic noise process.

We propose a transform for signals defined on the sphere that reveals their localized directional content in the spatio-spectral domain when used in conjunction with an asymmetric window function. We call this transform the directional spatially localized spherical harmonic transform (directional SLSHT) which extends the SLSHT from the literature whose usefulness is limited to symmetric windows. We present an inversion relation to synthesize the original signal from its directional-SLSHT distribution for an arbitrary window function. As an example of an asymmetric window, the most concentrated band-limited eigenfunction in an elliptical region on the sphere is proposed for directional spatio-spectral analysis and its effectiveness is illustrated on the synthetic and Mars topographic data-sets. Finally, since such typical data-sets on the sphere are of considerable size and the directional SLSHT is intrinsically computationally demanding depending on the band-limits of the signal and window, a fast algorithm for the efficient computation of the transform is developed. The floating point precision numerical accuracy of the fast algorithm is demonstrated and a full numerical complexity analysis is presented.