MiRA: image reconstruction for optical interferometry
MiRA is an algorithm for image reconstruction from data provided by optical interferometers. I wrote this algorithm to cope with problems specific to this type of data (sparseness of the measurements, missing of Fourier phase information, etc.) and not solved by existing methods, e.g., developed for radio-astronomy. MiRA won the 2008' Interferometric Imaging Beauty Contest [6] organized by IAU to compare the image synthesis algorithms designed for optical interferometry. An overview of the existing image synthesis algorithms in radio and optical interferometry is provided in [3]; more specifically, MiRA algorithm has been briefly described in [1, 5] and MiRA has been compared to Wisard method in [4].
In a nutshell, MiRA proceeds by direct minimization of a penalized likelihood. This penalty is the sum of two terms: a likelihood term (typically a χ2) which enforces agreement of the model with the data, plus a regularization term to account for priors. The priors are required to lever the many degeneracies due to the sparseness of the spatial frequency sampling. MiRA implements many different regularizations (quadratic or edge-preserving smoothness, total variation, maximum entropy, etc.) and let the user defines his own priors. The likelihood penalty is modular and designed to account for available data of any kind (complex visibilities, powerspectra and/or closure phase). One of the strength of MiRA is that it is purely based on an inverse problem approach and can therefore cope with incomplete data set; for instance, MiRA can build an image without any Fourier phase information.
MiRA is a research tool, it is however free software usable by others (though you need some basic knowledge in optical interferometry). For instance, MiRA has been used to teach image reconstruction in interferometry to the students of the 2008' VLTI Summer School held in Keszthely [3]. MiRA has been used to obtained the results published in a number of scientific journals [2, 7, 8] or conferences.
For now, input data must be in OI-FITS format, extension to other formats (e.g. UV-tables) is on the way. The algorithm is currently modified to achieve true 3-D reconstruction of multi-wavelength images and to account for other interferometers such as ALMA and next generation VLTI instruments (MATIS).
Installation and usage
You can get the last version of MiRA at GitHub
See INSTALL.md for prerequisites and instructions to install MiRA.
See USAGE.md for a short introduction about using MiRA.