Optical tomography with the discontinuous Galerkin formulation of the radiative transfer equation in frequency domain

Abstract Optical tomography is an inverse method of probing semi-transparent media with the help of light sources. The reconstruction of the optical properties usually employs finite volumes or continuous finite elements formulations of light transport as a forward model for the predictions. In a previous study, we have introduced a generalization of the inversion approach with finite elements formulations by using an integral form of the objective function. The novelty is that the surfaces of the detectors are taken into account in the reconstruction and compatibility is obtained for all finite element formulations. This present paper illustrates this new approach by developing a Discontinuous Galerkin formulation as a forward model for an optical tomography application in the frequency domain framework. Numerical tests are performed to gauge the accuracy of the method in recovering optical properties distribution with a gradient-based algorithm where the adjoint method is used to fastly compute the objective function gradient. It is seen that the reconstruction is accurate and can be affected by noise on the measurements as expected. Filtering of the gradient at each iteration of the reconstruction is used to cope with the ill-posed nature of the inverse problem and to improves the quality and accuracy of the reconstruction.

[1]  Todd E. Peterson,et al.  A note on the convergence of the discontinuous Galerkin method for a scalar hyperbolic equation , 1991 .

[2]  André Charette,et al.  New developments in frequency domain optical tomography. Part II: Application with a L-BFGS associated to an inexact line search , 2011 .

[3]  Dimitris Gorpas,et al.  A three-dimensional finite elements approach for the coupled radiative transfer equation and diffusion approximation modeling in fluorescence imaging , 2010 .

[4]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[5]  R. Alcouffe,et al.  Comparison of finite-difference transport and diffusion calculations for photon migration in homogeneous and heterogeneous tissues. , 1998, Physics in medicine and biology.

[6]  X. Cui,et al.  Discontinuous finite element solution of 2-D radiative transfer with and without axisymmetry , 2005 .

[7]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[8]  Suhas V. Patankar,et al.  Improved Treatment of Scattering Using the Discrete Ordinates Method , 1994 .

[9]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[10]  José Luis Morales,et al.  A numerical study of limited memory BFGS methods , 2002, Appl. Math. Lett..

[11]  D. Pelletier,et al.  Development of an adaptive Discontinuous-Galerkin finite element method for advection–reaction equations , 2007 .

[12]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[13]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[14]  Ben Q. Li.  Discontinuous Finite Elements in Fluid Dynamics and Heat Transfer , 2005 .

[15]  Hongkai Zhao,et al.  A Fast-Forward Solver of Radiative Transfer Equation , 2009 .

[16]  S Arridge,et al.  A gradient-based optimisation scheme foroptical tomography. , 1998, Optics express.

[17]  Simon R. Arridge,et al.  Variable order spherical harmonic expansion scheme for the radiative transport equation using finite elements , 2011, J. Comput. Phys..

[18]  Li-Ming Ruan,et al.  Finite element method for radiative heat transfer in absorbing and anisotropic scattering media , 2005 .

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

[20]  A. Crosbie,et al.  Radiative transfer in a two-dimensional rectangular medium exposed to diffuse radiation , 1984 .

[21]  Guillaume Bal,et al.  Frequency Domain Optical Tomography Based on the Equation of Radiative Transfer , 2006, SIAM J. Sci. Comput..

[22]  D. Marceau,et al.  A least square finite element formulation of the collimated irradiation in frequency domain for optical tomography applications , 2010 .

[23]  Alberto Cardona,et al.  Solution of non‐linear thermal transient problems by a reduction method , 1986 .

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

[25]  André Charette,et al.  New developments in frequency domain optical tomography. Part I: Forward model and gradient computation , 2011 .

[26]  S. Arridge,et al.  Optical imaging in medicine: II. Modelling and reconstruction , 1997, Physics in medicine and biology.

[27]  Li-Ming Ruan,et al.  On the discrete ordinates method for radiative heat transfer in anisotropically scattering media , 2002 .

[28]  Simon R. Arridge,et al.  Gauss–Newton reconstruction method for optical tomography using the finite element solution of the radiative transfer equation , 2008 .

[29]  Least-squares finite element formulations for one-dimensional radiative transfer , 2005 .

[30]  François Bay,et al.  Induction heating processes optimization a general optimal control approach , 2003 .

[31]  O. Balima,et al.  Comparison of light transport models based on finite element and the discrete ordinates methods in view of optical tomography applications , 2010, J. Comput. Appl. Math..

[32]  Tae-Kuk Kim,et al.  Radiative transfer in two-dimensional anisotropic scattering media with collimated incidence , 1989 .

[33]  Joan Boulanger,et al.  An overview on recent radiation transport algorithm development for optical tomography imaging , 2008 .

[34]  Claes Johnson.  Numerical solution of partial differential equations by the finite element method , 1988 .

[35]  W. Fiveland.  Discrete-Ordinates Solutions of the Radiative Transport Equation for Rectangular Enclosures , 1984 .

[36]  Thomas J. R. Hughes,et al.  A comparison of discontinuous and continuous Galerkin methods bases on error estimates, conservation, robustness and efficiency , 2000 .

[37]  Y. Rouizi,et al.  Numerical model reduction of 2D steady incompressible laminar flows: Application on the flow over a backward-facing step , 2009, J. Comput. Phys..