BIPP: An efficient HPC implementation of the Bluebild algorithm for radio astronomy

Abstract

The Bluebild algorithm is a new technique for image synthesis in radio astronomy which decomposes the sky into distinct energy levels using functional principal component analysis. These levels can be linearly combined to construct a least-squares estimate of the radio sky, i.e. minimizing the residuals between measured and predicted visibilities. This approach is particularly useful for deconvolution-free imaging or for scientific applications that need to filter specific energy levels. We present an HPC implementation of the Bluebild algorithm for radio-interferometric imaging: Bluebild Imaging++ (BIPP). The library features interfaces to C++, C and Python and is designed with seamless GPU acceleration in mind. We evaluate the accuracy and performance of BIPP on simulated observations of the upcoming Square Kilometer Array Observatory and real data from the Low-Frequency Array (LOFAR) telescope. We find that BIPP offers accurate wide-field imaging and has competitive execution time with respect to the interferometric imaging libraries CASA and WSClean for images with $\le 10^6$ pixels. Furthermore, due to the energy level decomposition, images produced with BIPP can reveal information about faint and diffuse structures before any cleaning iterations. BIPP does not perform any regularization, but we suggest methods to integrate the output of BIPP with CLEAN. The source code of BIPP is publicly released.

21 pages, 13 figures

References in corpus (14)

• CASA, the Common Astronomy Software Applications for Radio Astronomy
• WSClean: an implementation of a fast, generic wide-field imager for radio astronomy
• Multi-Scale CLEAN deconvolution of radio synthesis images
• The non-coplanar baselines effect in radio interferometry: The W-Projection algorithm
• An optimized algorithm for multi-scale wideband deconvolution of radio astronomical images
• Revisiting the radio interferometer measurement equation. I. A full-sky Jones formalism
• Faceting for direction-dependent spectral deconvolution
• Advances in Calibration and Imaging Techniques in Radio Interferometry
• Image Domain Gridding: a fast method for convolutional resampling of visibilities
• Deep Learning-based Imaging in Radio Interferometry
• Propagation Effects in Quiet Sun Observations at Meter Wavelengths
• A Generic and Efficient E-field Parallel Imaging Correlator for Next-Generation Radio Telescopes
• Scalable precision wide-field imaging in radio interferometry: I. uSARA validated on ASKAP data
• HVOX: Scalable Interferometric Synthesis and Analysis of Spherical Sky Maps