Application of the photon average trajectories method to real-time reconstruction of tissue inhomogeneities in diffuse optical tomography of strongly scattering media.

The possibility of application of the photon average trajectories (PAT) method to real-time reconstruction of tissue inhomogeneities in diffuse optical tomography of strongly scattering media has been substantiated. By this method, the inverse problem is reduced to solution of the integral equation with integration along a conditional PAT. Such an approach allows the standard fast algebraic algorithms commonly used in projection computed tomography to be applied to diffuse optical image reconstruction. To demonstrate the capabilities of the PAT method, a numerical experiment on cross-sectional reconstruction of cylindrical strongly scattering objects with absorbing inhomogeneities has been done. Relative shadows caused by inhomogeneities are simulated via numerical solution of the non-stationary diffusion equation. To solve the inverse problem, the QR-factorization least-squares algorithm and the multiplicative algebraic reconstruction technique are used. The results are compared with those obtained by a well-known software package for temporal optical absorption and scattering tomography based on multiple solution of the diffusion equation. It is shown that the PAT method allows reconstruction of the optical structure of objects with comparable accuracy while saving reconstruction time considerably.

[1]  Simon R. Arridge,et al.  A method for three‐dimensional time‐resolved optical tomography , 2000 .

[2]  I. Good,et al.  The Maximum Entropy Formalism. , 1979 .

[3]  Y. Censor.  Finite series-expansion reconstruction methods , 1983, Proceedings of the IEEE.

[4]  Arridge,et al.  THE USE OF MULTIPLE DATA-TYPES IN TIME-RESOLVED OPTICAL-ABSORPTION AND SCATTERING TOMOGRAPHY (TOAST) , 1993 .

[5]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

[6]  S Arridge,et al.  Optical tomography using the SCIRun problem solving environment: Preliminary results for three-dimensional geometries and parallel processing. , 1999, Optics express.

[7]  Gabor T. Herman,et al.  Image reconstruction from projections : the fundamentals of computerized tomography , 1980 .

[8]  J. Melissen,et al.  Tomographic image reconstruction from optical projections in light-diffusing media. , 1997, Applied optics.

[9]  S R Arridge,et al.  Comparison of two- and three-dimensional reconstruction methods in optical tomography. , 1998, Applied optics.

[10]  S. Helgason.  The Radon Transform , 1980 .

[11]  C. Gardiner.  Handbook of Stochastic Methods , 1983 .

[12]  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.

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

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

[15]  Simon R. Arridge,et al.  Use of multiple data types in time-resolved optical absorption and scattering tomography , 1993, Optics & Photonics.

[16]  G. Herman,et al.  Algebraic reconstruction techniques (ART) for three-dimensional electron microscopy and x-ray photography. , 1970, Journal of theoretical biology.

[17]  M. Schweiger,et al.  A general framework for iterative reconstruction algorithms in optical tomography using a finite element method , 1997 .

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

[19]  R. Lewitt.  Reconstruction algorithms: Transform methods , 1983, Proceedings of the IEEE.

[20]  Simon R. Arridge,et al.  Iterative reconstruction of near-infrared absorption images , 1992, Optics & Photonics.

[21]  Matthias Orlt,et al.  Reconstruction algorithm for near-infrared imaging in turbid media by means of time-domain data , 1997 .

[22]  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.

[23]  R. J. Gaudette,et al.  A comparison study of linear reconstruction techniques for diffuse optical tomographic imaging of absorption coefficient. , 2000, Physics in medicine and biology.

[24]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .