# Distributed and parallel sparse convex optimization for radio interferometry with PURIFY

@article{Pratley2019DistributedAP, title={Distributed and parallel sparse convex optimization for radio interferometry with PURIFY}, author={Luke Pratley and Jason D. McEwen and Mayeul d'Avezac and Xiaohao Cai and David P{\'e}rez-Su{\'a}rez and Ilektra A. Christidi and Roland Guichard}, journal={arXiv: Instrumentation and Methods for Astrophysics}, year={2019} }

Next generation radio interferometric telescopes are entering an era of big data with extremely large data sets. While these telescopes can observe the sky in higher sensitivity and resolution than before, computational challenges in image reconstruction need to be overcome to realize the potential of forthcoming telescopes. New methods in sparse image reconstruction and convex optimisation techniques (cf. compressive sensing) have shown to produce higher fidelity reconstructions of simulations… Expand

#### Supplemental Code

Github Repo

Sparse OPTimisation using state-of-the-art convex optimisation algorithms.

#### Paper Mentions

#### 5 Citations

Radio astronomy image reconstruction in the Big Data era

- Computer Science
- 2019

It is shown that PURIFY’s methods provide higher dynamic range than traditional image reconstruction methods, providing a more accurate and detailed sky model for real observations, and sets the stage for state-of-the-art image Reconstruction methods to be distributed and applied to next generation interferometric telescopes. Expand

Inference in optical interferometry: a quick review of modeling and imaging techniques

- Engineering, Computer Science
- Astronomical Telescopes + Instrumentation
- 2020

These topics are covered in an attempt to predict what the future holds for inference in optical interferometry in the next two decades. Expand

Offline and online reconstruction for radio interferometric imaging

- Physics, Computer Science
- 2020 XXXIIIrd General Assembly and Scientific Symposium of the International Union of Radio Science
- 2020

This article reviews methods proposed recently to resolve the ill-posed inverse problem of imaging the raw visibilities acquired by RI telescopes in the big-data scenario and focuses on the recently proposed online reconstruction method. Expand

Sparse Image Reconstruction for the SPIDER Optical Interferometric Telescope

- Computer Science, Physics
- PHOTOPTICS
- 2021

Algorithms from radio interferometric imaging and sparse image reconstruction are described and their application to the SPIDER concept telescope is demonstrated through simulated observation and reconstruction of the optical sky. Expand

w-Stacking w-projection hybrid algorithm for wide-field interferometric imaging: implementation details and improvements

- Physics, Computer Science
- Publications of the Astronomical Society of Australia
- 2020

It is shown that the use of conjugate symmetry can increase the w-stacking efficiency, decrease the time required to construction, and apply w-projection kernels for large data sets. Expand

#### References

SHOWING 1-10 OF 68 REFERENCES

Scalable splitting algorithms for big-data interferometric imaging in the SKA era

- Physics
- 2016

In the context of next generation radio telescopes, like the Square Kilometre Array, the efficient processing of large-scale datasets is extremely important. Convex optimisation tasks under the… Expand

Robust sparse image reconstruction of radio interferometric observations with PURIFY

- Physics
- 2018

Next-generation radio interferometers, such as the Square Kilometre Array, will revolutionize our understanding of the Universe through their unprecedented sensitivity and resolution. However, to… Expand

Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging

- Computer Science, Physics
- 2015

The theory behind interferometric imaging, sparse representations and convex optimization, and their application with numerical tests with SASIR, an implementation of the FISTA, a Forward-Backward splitting algorithm hosted in a LOFAR imager are exposed. Expand

Distributed multi-frequency image reconstruction for radio-interferometry

- Computer Science, Physics
- 2016 24th European Signal Processing Conference (EUSIPCO)
- 2016

Numerical results show that the IUWT-based version of MUFFIN can be successfully implemented at large scale with performances comparable to union of bases, and opens the possibility of comparing these efficient dictionaries, and others, in spatio-spectral reconstruction. Expand

Online radio interferometric imaging: assimilating and discarding visibilities on arrival

- Computer Science, Physics
- Monthly Notices of the Royal Astronomical Society
- 2019

It is anticipated online imaging techniques, as proposed here, will be critical in scaling RI imaging to the emerging big-data era of radio astronomy. Expand

MORESANE: MOdel REconstruction by Synthesis-ANalysis Estimators - A sparse deconvolution algorithm for radio interferometric imaging

- Computer Science, Physics
- 2014

This paper presents a new radio deconvolution algorithm named MORESANE and its application to fully realistic simulated data of MeerKAT, one of the SKA precursors and indicates that MORESAN provides competitive results in terms of both total flux/surface brightness conservation and fidelity of the reconstructed model. Expand

PURIFY: a new approach to radio-interferometric imaging

- Physics
- 2014

In a recent article series, the authors have promoted convex optimization algorithms for radio-interferometric imaging in the framework of compressed sensing, which leverages sparsity regularization… Expand

Cygnus A super-resolved via convex optimization from VLA data

- Physics
- 2018

We leverage the Sparsity Averaging Reweighted Analysis (SARA) approach for interferometric imaging, that is based on convex optimisation, for the super-resolution of Cyg A from observations at the… Expand

Simulating full‐sky interferometric observations

- Physics
- 2008

Aperture array interferometers, such as that proposed for the Square Kilometre Array, will see the entire sky, hence the standard approach to simulating visibilities will not be applicable since it… Expand

An improved method for polarimetric image restoration in interferometry

- Physics
- 2016

Interferometric radio astronomy data require the effects of limited coverage in the Fourier plane to be accounted for via a deconvolution process. For the last 40 years this process, known as… Expand