Image restoration for confocal microscopy: improving the limits of deconvolution, with application to the visualization of the mammalian hearing organ.

Deconvolution algorithms have proven very effective in conventional (wide-field) fluorescence microscopy. Their application to confocal microscopy is hampered, in biological experiments, by the presence of important levels of noise in the images and by the lack of a precise knowledge of the point spread function (PSF) of the system. We investigate the application of wavelet-based processing tools to deal with these problems, in particular wavelet denoising methods, which turn out to be very effective in application to three-dimensional confocal images. When used in combination with more classical deconvolution algorithms, these methods provide a robust and efficient restoration scheme allowing one to deal with difficult imaging conditions. To make our approach applicable in practical situations, we measured the PSF of a Biorad-MRC1024 confocal microscope under a large set of imaging conditions, including in situ acquisitions. As a specific biological application, we present several examples of restorations of three-dimensional confocal images acquired inside an intact preparation of the hearing organ. We also provide a quantitative assessment of the gain in quality achieved by wavelet-aided restorations over classical deconvolution schemes, based on a set of numerical experiments that we performed with test images.

[1]  S. Gibson,et al.  Experimental test of an analytical model of aberration in an oil-immersion objective lens used in three-dimensional light microscopy. , 1992, Journal of the Optical Society of America. A, Optics and image science.

[2]  D. Donoho Nonlinear Solution of Linear Inverse Problems by Wavelet–Vaguelette Decomposition , 1995 .

[3]  S. Hell,et al.  Aberrations in confocal fluorescence microscopy induced by mismatches in refractive index , 1993 .

[4]  S. Gibson,et al.  Experimental test of an analytical model of aberration in an oil-immersion objective lens used in three-dimensional light microscopy. , 1991, Journal of the Optical Society of America. A, Optics and image science.

[5]  Eero P. Simoncelli Modeling the joint statistics of images in the wavelet domain , 1999, Optics & Photonics.

[6]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[7]  Donald Geman,et al.  Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images , 1984 .

[8]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[9]  D. Donoho,et al.  Translation-Invariant De-Noising , 1995 .

[10]  K. Miller Least Squares Methods for Ill-Posed Problems with a Prescribed Bound , 1970 .

[11]  Jonathan Michael Pryce Statistical mechanics of image restoration , 1995 .

[12]  D. L. Donoho,et al.  Ideal spacial adaptation via wavelet shrinkage , 1994 .

[13]  P. Jansson,et al.  Resolution Enhancement of Spectra , 1970 .

[14]  M. Ulfendahl,et al.  Perilymphatic Fluid Compartments and Intercellular Spaces of the Inner Ear and the Organ of Corti , 2000, NeuroImage.

[15]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[16]  N O Petersen,et al.  Quantitation of membrane receptor distributions by image correlation spectroscopy: concept and application. , 1993, Biophysical journal.

[17]  Jean-Luc Starck,et al.  Filtering and deconvolution by the wavelet transform , 1994, Signal Process..

[18]  J. Pawley,et al.  Handbook of Biological Confocal Microscopy , 1990, Springer US.

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

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

[21]  Tomaso Poggio,et al.  Probabilistic Solution of Ill-Posed Problems in Computational Vision , 1987 .

[22]  S. Gull,et al.  Image reconstruction from incomplete and noisy data , 1978, Nature.

[23]  L. Shepp,et al.  Maximum Likelihood Reconstruction for Emission Tomography , 1983, IEEE Transactions on Medical Imaging.

[24]  I. Johnstone,et al.  Ideal spatial adaptation by wavelet shrinkage , 1994 .

[25]  B. R. Hunt,et al.  The Application of Constrained Least Squares Estimation to Image Restoration by Digital Computer , 1973, IEEE Transactions on Computers.

[26]  E. Candès,et al.  Recovering edges in ill-posed inverse problems: optimality of curvelet frames , 2002 .

[27]  M. Ulfendahl,et al.  Vital staining of the hearing organ: Visualization of cellular structure with confocal microscopy , 1998, Neuroscience.

[28]  Robert D. Nowak,et al.  A statistical multiscale framework for Poisson inverse problems , 2000, IEEE Trans. Inf. Theory.

[29]  A. N. Tikhonov,et al.  Solutions of ill-posed problems , 1977 .

[30]  N O Petersen,et al.  Scanning fluorescence correlation spectroscopy. I. Theory and simulation of aggregation measurements. , 1986, Biophysical journal.

[31]  I. Csiszár Why least squares and maximum entropy? An axiomatic approach to inference for linear inverse problems , 1991 .

[32]  Mats Ulfendahl,et al.  A temporal bone preparation for the study of cochlear micromechanics at the cellular level , 1989, Hearing Research.

[33]  D. M. Carlucci,et al.  Image restoration using the chiral Potts spin glass. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[34]  Sami K. Solanki,et al.  Noise reduction in astronomical spectra using wavelet packets , 1997 .

[35]  M Ulfendahl,et al.  An in vitro preparation to access cellular and neuronal components in the mouse inner ear , 2000, Journal of neurocytology.

[36]  Timothy J. Holmes Maximum-likelihood image restoration adapted for noncoherent optical imaging , 1988 .

[37]  I. Johnstone,et al.  Adapting to Unknown Smoothness via Wavelet Shrinkage , 1995 .

[38]  Mats Ulfendahl,et al.  Supporting Cells Contribute to Control of Hearing Sensitivity , 1999, The Journal of Neuroscience.

[39]  Mark A. Lukas,et al.  Methods for choosing the regularization parameter , 1992 .

[40]  Mario Bertero,et al.  Three‐dimensional image restoration and super‐resolution in fluorescence confocal microscopy , 1990 .

[41]  H. Nishimori,et al.  Statistical mechanics of image restoration and error-correcting codes. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.