Application of Regularized Richardson-Lucy Algorithm for Deconvolution of Confocal Microscopy Images

While confocal microscopes have considerably smaller contribution of out-of-focus light than widefield microscopes, the confocal images can still be enhanced by deconvolution if the optical and data acquisition effects are accounted for.Several deconvolution algorithms have been proposed for 3D microscopy. In this work we analyze the Richardson-Lucy iterative algorithm that is derived for Poisson noise and combined with total variation (TV) regularization. The influence of TV regularization on deconvolution process is determined by one parameter. However, the choice of regularization parameters is often unknown while it has considerable effect on the result of deconvolution process.The aims of this work were: to find good estimates of regularization parameter from the input; to develop an open source software package that would allow testing different deconvolution algorithms and that would be easy to use in practice. For that, we derived a formula to estimate this regularization parameter automatically from the images as the algorithm progresses. To assess the effectiveness of this algorithm, synthetic images were composed on the basis of confocal images of rat cardiomyocytes. From the analysis of deconvolved results, we have determined under which conditions our estimation of TV regularization parameter gives good results. The estimated TV regularization parameter can be monitored during deconvolution process and used as a stopping criterion. As a result, we propose a practical method to deconvolve confocal microscope images that uses estimated regularization parameter depending on the input image.We applied the deconvolution algorithm to study mitochondrial organization in rat cardiomyocytes. An open source software for deconvolving 3D images is available in http://sysbio.ioc.ee/software/.

[1]  Marko Vendelin,et al.  ADP Compartmentation Analysis Reveals Coupling between Pyruvate Kinase and ATPases in Heart Muscle , 2010, Biophysical journal.

[2]  G. B. Avinash,et al.  Data‐driven, simultaneous blur and image restoration in 3‐D fluorescence microscopy , 1996 .

[3]  Zhiping Lin,et al.  Noise suppression of point spread functions and its influence on deconvolution of three‐dimensional fluorescence microscopy image sets , 2005, Journal of microscopy.

[4]  Thomas J. Asaki,et al.  A Variational Approach to Reconstructing Images Corrupted by Poisson Noise , 2007, Journal of Mathematical Imaging and Vision.

[5]  Eric Jones,et al.  SciPy: Open Source Scientific Tools for Python , 2001 .

[6]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[7]  Mats Ulfendahl,et al.  Image-adaptive deconvolution for three-dimensional deep biological imaging. , 2003, Biophysical journal.

[8]  F S Fay,et al.  Superresolution three-dimensional images of fluorescence in cells with minimal light exposure. , 1995, Science.

[9]  Pearu Peterson,et al.  F2PY: a tool for connecting Fortran and Python programs , 2009, Int. J. Comput. Sci. Eng..

[10]  Josiane Zerubia,et al.  Richardson–Lucy algorithm with total variation regularization for 3D confocal microscope deconvolution , 2006, Microscopy research and technique.

[11]  Guido van Rossum,et al.  Python Programming Language , 2007, USENIX Annual Technical Conference.

[12]  S. Lockett,et al.  Calibration and standardization of the emission light path of confocal microscopes , 2006, Journal of microscopy.

[13]  Mark B. Cannell,et al.  Image Enhancement by Deconvolution , 2006 .

[14]  David M. Beazley,et al.  An Extensible Compiler for Creating Scriptable Scientific Software , 2002, International Conference on Computational Science.

[15]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[16]  J. Zerubia,et al.  3D Microscopy Deconvolution using Richardson-Lucy Algorithm with Total Variation Regularization , 2004 .

[17]  L. J. Thomas,et al.  Regularized linear method for reconstruction of three-dimensional microscopic objects from optical sections. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[18]  S. Riad,et al.  Filtering capabilities and convergence of the Van-Cittert deconvolution technique , 1991 .

[19]  J. Conchello,et al.  Parametric blind deconvolution: a robust method for the simultaneous estimation of image and blur. , 1999, Journal of the Optical Society of America. A, Optics, image science, and vision.

[20]  J. Lichtman,et al.  Optical sectioning microscopy , 2005, Nature Methods.

[21]  Chrysanthe Preza,et al.  Depth-variant maximum-likelihood restoration for three-dimensional fluorescence microscopy. , 2004, Journal of the Optical Society of America. A, Optics, image science, and vision.

[22]  S. Inoué,et al.  Foundations of Confocal Scanned Imaging in Light Microscopy , 2006 .

[23]  David M. Beazley,et al.  Automated scientific software scripting with SWIG , 2003, Future Gener. Comput. Syst..

[24]  Steven G. Johnson,et al.  The Design and Implementation of FFTW3 , 2005, Proceedings of the IEEE.

[25]  Michael I. Miller,et al.  Image reconstruction for 3D light microscopy with a regularized linear method incorporating a smoothness prior , 1993, Electronic Imaging.

[26]  D. Rawlins,et al.  Three-dimensional fluorescence microscopy. , 1991, Progress in biophysics and molecular biology.

[27]  G. M. P. VAN KEMPEN,et al.  A quantitative comparison of image restoration methods for confocal microscopy , 1997 .

[28]  D. Rawlins,et al.  The point‐spread function of a confocal microscope: its measurement and use in deconvolution of 3‐D data , 1991 .

[29]  F Ianzini,et al.  An open‐source deconvolution software package for 3‐D quantitative fluorescence microscopy imaging , 2009, Journal of microscopy.

[30]  Richard G. Baraniuk,et al.  ForWaRD: Fourier-wavelet regularized deconvolution for ill-conditioned systems , 2004, IEEE Transactions on Signal Processing.

[31]  L. Lucy.  An iterative technique for the rectification of observed distributions , 1974 .

[32]  P. Shaw.  Comparison of Widefield/Deconvolution and Confocal Microscopy for Three-Dimensional Imaging , 2006 .

[33]  J Boutet de Monvel,et al.  Image restoration for confocal microscopy: improving the limits of deconvolution, with application to the visualization of the mammalian hearing organ. , 2001, Biophysical journal.

[34]  T J Holmes,et al.  Blind deconvolution of quantum-limited incoherent imagery: maximum-likelihood approach. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[35]  Marko Vendelin,et al.  Anisotropic diffusion of fluorescently labeled ATP in rat cardiomyocytes determined by raster image correlation spectroscopy , 2008, American journal of physiology. Cell physiology.

[36]  Van Kempen,et al.  The influence of the regularization parameter and the first estimate on the performance of tikhonov regularized non-linear image restoration algorithms , 2000 .

[37]  J A Conchello,et al.  Superresolution and convergence properties of the expectation-maximization algorithm for maximum-likelihood deconvolution of incoherent images. , 1998, Journal of the Optical Society of America. A, Optics, image science, and vision.

[38]  D. Agard.  Optical sectioning microscopy: cellular architecture in three dimensions. , 1984, Annual review of biophysics and bioengineering.

[39]  Fang Li,et al.  Selection of regularization parameter in total variation image restoration. , 2009, Journal of the Optical Society of America. A, Optics, image science, and vision.

[40]  Marko Vendelin,et al.  Three-dimensional mitochondrial arrangement in ventricular myocytes: from chaos to order. , 2006, American journal of physiology. Cell physiology.

[41]  J. Zerubia,et al.  Parametric Blind Deconvolution for Confocal Laser Scanning Microscopy , 2007, 2007 29th Annual International Conference of the IEEE Engineering in Medicine and Biology Society.