Issue |
A&A
Volume 688, August 2024
|
|
---|---|---|
Article Number | A100 | |
Number of page(s) | 20 | |
Section | Numerical methods and codes | |
DOI | https://doi.org/10.1051/0004-6361/202449325 | |
Published online | 08 August 2024 |
Swarm intelligence for full Stokes dynamic imaging reconstruction of interferometric data
1
Departament d’Astronomia i Astrofísica, Universitat de València,
C. Dr. Moliner 50,
46100 Burjassot,
València,
Spain
e-mail: alejandro.mus@uv.es
2
Observatori Astronòmic, Universitat de València,
Parc Científic, C. Catedràtico José Beltrán 2,
46980
Paterna,
València,
Spain
3
Departament de Matemàtiques i Informàtica, Universitat de les Illes Balears,
07071
Palma de Mallorca,
Spain
4
Max-Planck-Institut für Radioastronomie,
Auf dem Hügel 69,
53121
Bonn (Endenich),
Germany
e-mail: hmueller@mpifr-bonn.mpg.de
5
Jansky Fellow of National Radio Astronomy Observatory,
1011 Lopezville Rd,
Socorro,
NM
87801,
USA
Received:
23
January
2024
Accepted:
2
May
2024
Context. In very long baseline interferometry (VLBI), the combination of multiple antennas permits the synthesis of a virtual telescope with a larger diameter and consequently higher resolution than the individual antennas. However, due to the sparse nature of the array, recovering an image from the observed data is a challenging ill-posed inverse problem.
Aims. The VLBI community is interested in not only recovering an image in total intensity from interferometric data, but also in obtaining results in the polarimetric and the temporal domain. Only a few algorithms are able to work in all these domains simultaneously. In particular, the algorithms based on optimization that consider various penalty terms specific to static total intensity imaging, time-variability and polarimetry are restricted to grids in the domain of the objective function. In this work we present a novel algorithm, multiobjective particle swarm optimization (MO-PSO), that is able to recover the optimal weights without any space-gridding, and to obtain the marginal contribution of each of the playing terms.
Methods. To this end, we utilized multiobjective optimization together with particle swarm metaheuristics. We let the swarm of weights converge to the best position.
Results. We evaluate our algorithm with synthetic data sets that are representative for the main science targets and instrumental configuration of the Event Horizon Telescope Collaboration (EHTC) and its planned successors. We successfully recover the polarimetric, static, and time-dynamic signature of the ground truth movie' even with relative sparsity, and a set of realistic data corruptions.
Conclusions. We have built a novel, fast, hyperparameter space gridding-free algorithm that successfully recovers static and dynamic polarimetric reconstructions. Compared to regularized maximum likelihood (RML) methods, it avoids the need for parameter surveys, and it is not limited to the number of pixels, unlike recently proposed multiobjective imaging algorithms. Hence, this technique is a novel useful alternative tool to characterize full Stokes time-(in)dependent signatures in a VLBI data set robustly with a minimal set of user-based choices.
Key words: methods: numerical / techniques: high angular resolution / techniques: image processing / techniques: interferometric
© The Authors 2024
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.