500 likes | 644 Views
Probing the Prostate with Diffusive Photons (2): Procession of Photon Diffuse Propagation in Cylinder-Applicator Geometry. Daqing ( Daching ) Piao, PhD Assistant Professor of Bioengineering School of Electrical and Computer Engineering Oklahoma State University, Stillwater, OK 74078-5032.
E N D
Probing the Prostate with Diffusive Photons (2): Procession of Photon Diffuse Propagation in Cylinder-Applicator Geometry Daqing (Daching) Piao, PhD Assistant Professor of Bioengineering School of Electrical and Computer Engineering Oklahoma State University, Stillwater, OK 74078-5032
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
Diffuse Optical Imaging (Tomgoraphy) Source Detector Multiple source-detector pairs----tomography reconstruction
Image reconstruction Measurement ? Forward problem Medium property ? ? Inverse problem
Homogenous ? Source: Cameron Musgrove
Tissue scattering dominates in NIR spectrum Photon propagation becomes essentially diffused after ~5mm for most soft tissue Source: Brian. Pogue
Radiative Transport Equation The equation of radiative transfer simply says that as a beam of radiation travels, it loses energy to absorption, gains energy by emission, and redistributes energy by scattering. Extinction Divergence L: Radiance Source Scattering
Photon Diffusion Equation Fluence rate Current density Photon diffusion equation is a second order partial differential equation describing the time behavior of photon fluence rate distribution in a low-absorption high-scattering medium.
Treatment of the Boundary Condition Dirichiliet Boundary Condition Where vanishes on the boundary Neumann Boundary Condition Where the normal gradient of vanishes on the boundary Robin Boundary Condition The linear combination of the type-I and type-II conditions
Non-planar Interface • Literatures • S.R. Arridge, M. Cope, D.T. Delpy, “The theoretical basis for the determination of optical pathlengths in tissue: temporal and frequency analysis,” Phys Med Biol., 37(7):1531-60 (1992). • B.W. Pogue, and M.S. Patterson, “Frequency-domain optical absorption spectroscopy of finite tissue volumes using diffusion theory,” Phys Med Biol., 39(7):1157-80 (1994). • A. Sassaroli, F. Martelli, G. Zaccanti, and Y. Yamada, "Performance of fitting procedures in curved geometry for retrieval of the optical properties of tissue from time-resolved measurements," Appl. Opt. 40, 185-197 (2001).
Non-planar Interface Probe
Motivation Joseph M. Lasker, etc Journal of Biomedical Optics 12(5), 052001 September/October 2007 Probe
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
Boltzmann (Radiative) Transport Equation Photon Diffusion Equation Time-resolved Steady-state Point source Homogenous medium Isotropic detector The equation for Green’s function
Boltzmann (Radiative) Transport Equation Photon Diffusion Equation Time-resolved Steady-state Point source Spherical coordinates Isotropic detector
Point source Spherical coordinates Isotropic detector
Cylindrical-coordinate solution to the steady-state photon diffusion in an infinite medium unknown
Cylindrical-coordinate solution to the steady-state photon diffusion in an infinite medium The equation for the radial Green’s function : where or . Following the approach for solving Poisson equation introduced by Jaskson, John David in Classical Electrodynamics Chapter3 , we find: where and are the modified Bessel functions, and and indicate the smaller and larger radial coordinates of the source and the detector.
Cylindrical-coordinate solution to the steady-state photon diffusion in an infinite medium where
Boltzmann (Radiative) Transport Equation Photon Diffusion Equation Time-resolved Steady-state Point source Spherical coordinates Isotropic detector Cylindrical coordinates
Numerical Validation to the cylindrical-coordinate solution of the steady-state photon diffusion in an infinite medium Spherical Coordinates ? Cylindrical Coordinates
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
When there is a boundary ? Image source Point source Isotropic detector
Steady-state photon diffusion in a medium with planner boundary Geometry introduced by Fantini et al. where
Steady-state photon diffusion in a medium with convex boundary Convex boundary: a radius of the equivalent isotropic source the detector is at real boundary The photon intensity associated with the real source is: known unknown
Steady-state photon diffusion in a medium with convex boundary where is to be determined by the boundary condition. Besides, the factor is introduced to simplify the determination of .
Steady-state photon diffusion in a medium with convex boundary At the extrapolated boundary, and Evaluating the same order of m
Steady-state photon diffusion in a medium with convex boundary For detector point which is located at the real boundary, and Contribution due to the real source Contribution due to the image source
Steady-state photon diffusion in a medium with concave boundary Concave boundary: At the extrapolated boundary,
Steady-state photon diffusion in a medium with concave boundary For detector point , and
Comparison between the solution under convex boundary and concave boundary Convex boundary: Concave boundary: physical boundary extrapolated boundary Position of the equivalent real source
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
Longitudinal direction … …………………….Part 1 Journal of Optical Society of America, A, Vol. 27, No. 3, pp. 648-662 (2010)
Outline • Motivation --------- Improve the model for trans-lumenalimaging • Steady-state photon diffusion in an infinite medium --------- Solution in cylindrical-coordinateand its numerical validation • Steady-state photon diffusion with a cylindrical boundary------ Analytic solutionfor trans-lumenal (convex & concave) geometry • Qualitative numerical evaluation of the analytic results • Quantitative examinations---Finite Element, Experiment, Monte Carlo
Results … …………………….Part 2 Journal of Optical Society of America, A, in preparation
Zhang A, Piao D, Bunting CF, Pogue BW, “Photon diffusion in a homogeneous medium bounded externally or internally by an infinitely long circular cylindrical applicator ---- Part I: steady-state theory,” Journal of Optical Society of America, A, accepted.
Summary • The Green’s function for steady-state photon diffusion in infinite medium is expanded in cylindrical coordinates. • The solutions of steady-state photon diffusion is derived for a circular concave or convex boundary. • The analytic solutions are qualitatively evaluated. • The analytic solutions are quantitatively examined
Acknowledgement • Anqi (Andrew) Zhang • Guan (Gary) Xu • Dr. Gang (Gary) Yao • Dr. Brian W. Pogue • Dr. Charles F. Bunting Oklahoma Center for the Advancement of Science and Technology (OCAST)