Fast Sky to Sky Interpolation for Radio Interferometric Imaging

Nicolas Monnier (1) , François Orieux (1) , Nicolas Gac (1) , Cyril Tasse (2) , Erwan Raffin (3) , David Guibert (3)

1 L2S - Laboratoire des signaux et systèmes
2 GEPI - Galaxies, Etoiles, Physique, Instrumentation
3 Atos/Bull Center for Excellence in Parallel Programming

Reconstruction of radio interferometric images requires the processing of data in Fourier space that dot not have regular coordinates, preventing direct use of the Fast Fourier Transform. The most common solution is to rely on interpolation algorithms, called gridding, that are computationally expensive. In this paper, we propose an algorithmic reinterpretation, named sky to sky method, to reduce the computation cost of the gridding operation and its adjoint, the degridding, when used successively. We analyze the impact of interpolation step size regarding the computation cost and the reconstruction error. We also illustrate this optimization on a full reconstruction with gradient descent and CLEAN algorithm. Finally, we obtain acceleration factors between 1.2 and 16.4 without additional approximation.

Mots clés: Interpolation; inverse problem; gridding; Computation acceleration

Domaines: Traitement du signal et de l'image [eess.SP]

Nicolas Monnier, François Orieux, Nicolas Gac, Cyril Tasse, Erwan Raffin, et al.. Fast Sky to Sky Interpolation for Radio Interferometric Imaging. IEEE International Conference on Image Processing (ICIP), Oct 2022, Bordeaux, France. ⟨hal-03725824v1⟩