IViS: Interferometric Visibility-domain inversion Software - A GPU-accelerated Python framework for joint deconvolution with ASKAP
astro-ph.IM, astro-ph.GA
Submitted: 2026-09-15
Updated: 2026-09-15
Comments: 20 pages, 21 figures; Accepted for publication in A&A
Code: https://github.com/antoinemarchal/fBms
License: http://creativecommons.org/licenses/by-nc-sa/4.0/
The gist: Wide-field spectral-line imaging with modern radio interferometers remains challenging when the emission is extended, multiscale, and distributed over many overlapping pointings.
Terminology
Abstract
Wide-field spectral-line imaging with modern radio interferometers remains challenging when the emission is extended, multiscale, and distributed over many overlapping pointings. Accurate reconstruction requires joint treatment of the calibrated visibilities, control of image-domain regularity, and recovery of missing short spacings from single-dish data. We present IViS (Interferometric Visibility-domain Inversion Software), a GPU-accelerated Python framework for visibility-domain joint deconvolution, and assess its performance for wide-field HI imaging with the Australian Square Kilometre Array Pathfinder (ASKAP). The current model, Classic3D, reconstructs a non-parametric sky cube directly from calibrated visibilities by minimizing a regularized least-squares criterion. The forward model relies on a non-uniform fast Fourier transform (NUFFT), supports joint deconvolution across mosaic pointings and positivity constraints, and can incorporate single-dish data through a fusion term. We validate the method using point-source, noise-only, and multiscale diffuse-emission simulations, and apply it to ASKAP observations toward the Large Magellanic Cloud (LMC). The simulations show that IViS recovers point-source flux accurately without regularization and that regularization sets the trade-off between noise suppression and effective resolution. Joint deconvolution can recover more large-scale power than linear mosaicking of independently deconvolved pointings. With single-dish information, the reconstructed power spectra closely reproduce the input sky statistics. Applied to ASKAP data, IViS yields 10h mosaics with an effective resolution of 22'' and reduced residual side-lobe structure. Compared with ASKAPSoft multi-scale CLEAN for a single observing block, Classic3D exhibits fewer residual side-lobes and recovers substantially more power at low spatial frequencies.
Sources
- Aliasing error of the exp$(\beta \sqrt{1-z^2})$ kernel in the nonuniform fast Fourier transform
- Intermittent process analysis with scattering moments
- Multi-Scale CLEAN deconvolution of radio synthesis images
- PyTorch: An Imperative Style, High-Performance Deep Learning Library
- cuFINUFFT: a load-balanced GPU library for general-purpose nonuniform FFTs
Related papers
- A signal dedispersion algorithm for imaging-based transient searches
- AVICA: A fully automated CASA pipeline for large volume VLBI data calibration
- Spectral Map Making with SPHEREx
- Long-Integration Magnetar Burst Observatory (LIMBO): Instrument Summary and Early FRB Rate Constraints
- Towards independent event horizon imaging of the supermassive black holes in M87 and the Milky Way
- A PINK update: Improvements to the CELEBI fast radio burst data reduction and analysis pipeline