Free access
Astron. Astrophys. Suppl. Ser.
Volume 140, Number 2, December I 1999
Page(s) 235 - 246
DOI: 10.1051/aas:1999420

Astron. Astrophys. Suppl. Ser. 140, 235-246

Comparison between ISRA and RLA algorithms. Use of a Wiener Filter based stopping criterion

H. Lantéri - R. Soummer - C. Aime

UMR 6525, "Astrophysique", Université de Nice, Sophia Antipolis, Faculté des Sciences, Parc Valrose, 06108 Nice Cedex 2, France

Received January 29; accepted July 23, 1999


This paper consists of two parts. In the first one, we make a comparative analysis of the algorithms of Richardson-Lucy (RLA) and the Image Space Reconstruction Algorithm (ISRA), for image deconvolution of astronomical images. These iterative algorithms keep the reconstructed image non-negative while maximizing the likelihood for a Poisson Process (RLA) or a Gaussian process (ISRA). Their comparison is made easier when these algorithms are rewritten as descent algorithms; the additive forms evidence the role of the variances of the noise and allow a better understanding of these algorithms. A numerical illustration is performed from simulated images. In practice, the results obtained by the two algorithms appeared to be very similar, independently of the statistics of the noise. In the second part of the paper, we propose a new objective-stopping technique that makes use of a comparison of the results of these algorithms with that of the Wiener filter. The comparison is made in the Fourier plane, computing the Euclidean distances between modules of the spatial frequencies components of the images. We propose then to use the results of ISRA and RLA at the iteration number corresponding to the minimum of that distance. The technique is checked in a numerical simulation, for which the optimal iteration numbers can be easily determined.

A good agreement is obtained between best ISRA and RLA results and the Wiener filter, in particular for rather noisy images. However ISRA and RLA produce very high frequency components, outside the cut-off frequency of the instrument; limiting the iteration number alone cannot allow a perfect agreement with the Wiener approach, stressing the need of an explicit regularization for astronomical applications.

Key words: methods: numerical -- techniques: image processing

Copyright The European Southern Observatory (ESO)

What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access.

An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.

  • If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
  • You can define your own OpenURL resolver with your EDPS Account.
    In this case your choice will be given priority over that of your library.
  • You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see You should disable this module if you wish to use the OpenURL server that you or your library have defined.