Scalable Algorithm for Polarization Constrained Sparse Interferometric Stokes Imaging

摘要

Polarized emissions from various astrophysical sources provide an additional means to infer the properties of the source. This information can be extracted by imaging the Stokes parameters characterizing these emissions, from the incomplete Fourier sampling measurements provided by the radio interferometer. To solve the corresponding ill-posed inverse problem, we recently proposed the Polarized SARA method which exploits the polarization constraint along with imposing an average sparsity prior for each of the Stokes parameters, encompassed in a reweighting approach. Motivated by the large amounts of data produced by the next-generation radio interferometers, in the current work, we propose a scalable version of this method, combined with an accelerating strategy. We showcase the performance of this improved method on a real dataset acquired by the Very Large Array (VLA).