Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)

Jasleen Birdi, Audrey Repetti, Yves Wiaux

Research output: Contribution to journalArticlepeer-review

54 Downloads (Pure)


We develop a novel algorithm for sparse Stokes parameters imaging in radio interferometry under the polarization constraint. The latter is a physical non-linear relation between the Stokes parameters, imposing the polarization intensity as a lower bound on the total intensity. To solve the joint inverse Stokes imaging problem including this bound, we leverage epigraphical projection techniques in convex optimization and design a primal-dual method offering a highly flexible and parallelizable structure. In addition, we propose to regularize each Stokes parameter map through an average sparsity prior in the context of a reweighted analysis approach (SARA). The resulting approach is dubbed Polarized SARA. We demonstrate on simulated observations of M87 with the Event Horizon Telescope that imposing the polarization constraint leads to superior image quality. For the considered datasets, the results also indicate better performance of the average sparsity prior in comparison with the widely used Cotton-Schwab CLEAN algorithm and other total variation based priors for polarimetric imaging.
Original languageEnglish
Pages (from-to)4442-4463
Number of pages22
JournalMonthly Notices of the Royal Astronomical Society
Issue number4
Early online date4 Jul 2018
Publication statusPublished - 21 Aug 2018


Dive into the research topics of 'Sparse interferometric Stokes imaging under polarization constraint (Polarized SARA)'. Together they form a unique fingerprint.

Cite this