A New Semiblind Deconvolution Approach for Fourier-Based Image Restoration: An Application in Astronomy

The aim of this paper is to develop a new optimization algorithm for the restoration of an image starting from samples of its Fourier transform, when only partial information about the data frequencies is provided. The corresponding constrained optimization problem is approached with a cyclic block alternating scheme, in which projected gradient methods are used to find a regularized solution. Our algorithm is then applied to the imaging of high-energy radiation emitted during a solar flare through the analysis of the photon counts collected by the NASA Reuven Ramaty High Energy Solar Spectroscopic Imager satellite. Numerical experiments on simulated data show that, in both the presence and absence of statistical noise, the proposed approach provides some improvements in the reconstructions.