Sparse inpainting techniques are gaining in popularity as a tool for cosmological data analysis, in particular for handling data which present masked regions and missing observations. We investigate here the relationship between sparse inpainting techniques using the spherical harmonic basis as a dictionary and the isotropy properties of cosmological maps, as for instance those arising from cosmic microwave background (CMB) experiments. In particular, we investigate the possibility that inpainted maps may exhibit anisotropies in the behaviour of higher-order angular polyspectra. We provide analytic computations and simulations of inpainted maps for a Gaussian isotropic model of CMB data, suggesting that the resulting angular trispectrum may exhibit small but non-negligible deviations from isotropy.

Convex regularization techniques are now widespread tools for solving inverse problems in a variety of different frameworks. In some cases, the functions to be reconstructed are naturally viewed as realizations from random processes; an important question is thus whether such regularization techniques preserve the properties of the underlying probability measures. We focus here on a case which has produced a very lively debate in the cosmological literature, namely Gaussian and isotropic spherical random fields, and we prove that Gaussianity and isotropy are not conserved in general under convex regularization over a Fourier dictionary, such as the orthonormal system of spherical harmonics.

This paper provides quantitative Central Limit Theorems for nonlinear transforms of spherical random fields, in the high frequency limit. The sequences of fields that we consider are represented as smoothed averages of spherical Gaussian eigenfunctions and can be viewed as random coefficients from continuous wavelets/needlets; as such, they are of immediate interest for spherical data analysis. In particular, we focus on so-called needlets polyspectra, which are popular tools for nonGaussianity analysis in the astrophysical community, and on the area of excursion sets. Our results are based on Stein-Maliavin approximations for nonlinear transforms of Gaussian fields, and on an explicit derivation on the high-frequency limit of their variances, which may have some independent interest.