Space-Varying Iterative Restoration of Diffuse Optical Tomograms Reconstructed by the Photon Average Trajectories Method

The possibility of improving the spatial resolution of diffuse optical tomograms reconstructed by the photon average trajectories (PAT) method is substantiated. The PAT method recently presented by us is based on a concept of an average statistical trajectory for transfer of light energy, the photon average trajectory (PAT). The inverse problem of diffuse optical tomography is reduced to a solution of an integral equation with integration along a conditional PAT. As a result, the conventional algorithms of projection computed tomography can be used for fast reconstruction of diffuse optical images. The shortcoming of the PAT method is that it reconstructs the images blurred due to averaging over spatial distributions of photons which form the signal measured by the receiver. To improve the resolution, we apply a spatially variant blur model based on an interpolation of the spatially invariant point spread functions simulated for the different small subregions of the image domain. Two iterative algorithms for solving a system of linear algebraic equations, the conjugate gradient algorithm for least squares problem and the modified residual norm steepest descent algorithm, are used for deblurring. It is shown that a gain in spatial resolution can be obtained.

[1]  Vladimir V. Lyubimov Physical foundations of the strongly scattering media laser tomography , 1996, Other Conferences.

[2]  Gene H. Golub,et al.  Matrix computations (3rd ed.) , 1996 .

[3]  Dianne P. O'Leary,et al.  Restoring Images Degraded by Spatially Variant Blur , 1998, SIAM J. Sci. Comput..

[4]  M. Schweiger,et al.  The finite element method for the propagation of light in scattering media: boundary and source conditions. , 1995, Medical physics.

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

[6]  B. Pogue,et al.  Multiwavelength three-dimensional near-infrared tomography of the breast: initial simulation, phantom, and clinical results. , 2003, Applied optics.

[7]  A. Papoulis Systems and transforms with applications in optics , 1981 .

[8]  Dianne P. O'Leary,et al.  The Use of the L-Curve in the Regularization of Discrete Ill-Posed Problems , 1993, SIAM J. Sci. Comput..

[9]  M. Hanke Conjugate gradient type methods for ill-posed problems , 1995 .

[10]  C. Vogel Computational Methods for Inverse Problems , 1987 .

[11]  J. Nagy,et al.  KRONECKER PRODUCT AND SVD APPROXIMATIONS IN IMAGE RESTORATION , 1998 .

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

[13]  Alexander D. Klose,et al.  Gradient-based iterative image reconstruction scheme for time-resolved optical tomography , 1999, IEEE Transactions on Medical Imaging.

[14]  Thomas S. Huang,et al.  Inverse filtering for linear shift-variant imaging systems , 1972 .

[15]  E. Miller,et al.  Tomographic optical breast imaging guided by three-dimensional mammography. , 2003, Applied optics.

[16]  Britton Chance,et al.  Monte Carlo simulations of photon migration path distributions in multiple scattering media , 1993, Photonics West - Lasers and Applications in Science and Engineering.

[17]  B. Chance,et al.  Spectroscopy and Imaging with Diffusing Light , 1995 .

[18]  Mario Bertero,et al.  Introduction to Inverse Problems in Imaging , 1998 .

[19]  C. W. Groetsch,et al.  The theory of Tikhonov regularization for Fredholm equations of the first kind , 1984 .

[20]  Misha Elena Kilmer,et al.  Choosing Regularization Parameters in Iterative Methods for Ill-Posed Problems , 2000, SIAM J. Matrix Anal. Appl..

[21]  Wolfgang Osten,et al.  Introduction to Inverse Problems in Imaging , 1999 .

[22]  Bobby R. Hunt,et al.  Approximate shift-invariance by warping shift-variant systems , 1994, Optics & Photonics.

[23]  Linda Kaufman,et al.  Maximum likelihood, least squares, and penalized least squares for PET , 1993, IEEE Trans. Medical Imaging.

[24]  James G. Nagy,et al.  Iterative Methods for Image Deblurring: A Matlab Object-Oriented Approach , 2004, Numerical Algorithms.

[25]  Alexander B. Konovalov,et al.  Application of transform algorithms to high-resolution image reconstruction in optical diffusion tomography of strongly scattering media , 2003, J. Electronic Imaging.

[26]  J. Navarro-Pedreño Numerical Methods for Least Squares Problems , 1996 .

[27]  Alexander B. Konovalov,et al.  Influence of fast reconstruction algorithms on spatial resolution of optical diffuse tomography by photon average trajectories method , 2002, Saratov Fall Meeting.

[28]  J. Nagy,et al.  Steepest Descent, CG, and Iterative Regularization of Ill-Posed Problems , 2003 .

[29]  Lothar Reichel,et al.  Non-negativity and iterative methods for ill-posed problems , 2004 .

[30]  Eric C. Frey,et al.  Comparison between ML-EM and WLS-CG algorithms for SPECT image reconstruction , 1991 .

[31]  M. Glas,et al.  Principles of Computerized Tomographic Imaging , 2000 .

[32]  S. Arridge Optical tomography in medical imaging , 1999 .

[33]  Alexander A. Sawchuk,et al.  Space-variant image restoration by coordinate transformations* , 1974 .

[34]  Alexander B. Konovalov,et al.  Application of integral transform algorithms to high-resolution reconstruction of tissue inhomogeneities in medical diffuse optical tomography , 2002, SPIE/COS Photonics Asia.

[35]  D. Sandwell BIHARMONIC SPLINE INTERPOLATION OF GEOS-3 AND SEASAT ALTIMETER DATA , 1987 .

[36]  Olga V. Kravtsenyuk,et al.  Calculation of shadows induced by macroinhomogeneities located inside a strongly scattering object using integration over the average photon path , 1998, European Conference on Biomedical Optics.

[37]  H. Joel Trussell,et al.  Identification and restoration of spatially variant motion blurs in sequential images , 1992, IEEE Trans. Image Process..

[38]  Robert J. Hanisch,et al.  The restoration of HST images and spectra - II , 2015 .

[39]  V. V. Lyubimov,et al.  Application of the photon average trajectories method to real-time reconstruction of tissue inhomogeneities in diffuse optical tomography of strongly scattering media. , 2002, Physics in medicine and biology.

[40]  Ming Jiang,et al.  Blind deblurring of spiral CT images , 2003, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[41]  Simon R. Arridge,et al.  Investigation of alternative data types for time-resolved optical tomography , 1998 .

[42]  J. Nagy,et al.  Enforcing nonnegativity in image reconstruction algorithms , 2000, SPIE Optics + Photonics.

[43]  James G. Nagy,et al.  Fast iterative image restoration with a spatially varying PSF , 1997, Optics & Photonics.

[44]  E. V. Appleton,et al.  Optics and spectroscopy: An english translation of the russian , 1960 .

[45]  Alexander B. Konovalov,et al.  Application of photon average trajectories method for separate mapping of absorbing and scattering macroinhomogeneities using time-domain measurements technique , 2001, European Conference on Biomedical Optics.

[46]  D. A. Fish,et al.  Scanning singular-value-decomposition method for restoration of images with space-variant blur , 1996 .