Stochastic Gauss-Newton Method for Estimating Absorption and Scattering in Optical Tomography with the Monte Carlo Method for Light Transport

Overview

Journal Biomed Opt Express

Specialty Radiology

Date 2024 Sep 30

PMID 39347007

Authors

Jonna Kangasniemi

Meghdoot Mozumder

Aki Pulkkinen

Tanja Tarvainen

Affiliations

Soon will be listed here.

Abstract

Image reconstruction in optical tomography in the so-called transport regime, where the diffusion approximation is not valid, requires modeling of light transport using the radiative transfer equation. In this work, we approach this problem by utilizing the Monte Carlo method for light transport. In this work, we propose a methodology for absolute imaging of absorption and scattering in this regime utilizing a Monte Carlo method for light transport. The image reconstruction problem is formulated as a minimization problem that is solved using a stochastic Gauss-Newton method. In the construction of the Jacobian matrix for scattering, a perturbation approximation for Monte Carlo is utilized. The approach is evaluated with numerical simulations using an adaptive approach where the number of photon packets is adjusted during the iterations, and with different fixed numbers of photon packets. The simulations show that the Monte Carlo method for light transport can be utilized in the absolute imaging problem of optical tomography and that the absorption and scattering parameters can be estimated simultaneously with good accuracy.

Citing Articles

Diffuse Optical Spectroscopy: Technology and Applications: introduction to the feature issue.

OSullivan T, Dehghani H, Re R Biomed Opt Express. 2024; 15(11):6516-6520.

PMID: 39553860 PMC: 11563346. DOI: 10.1364/BOE.542635.

References

Gibson A, Hebden J, Arridge S . Recent advances in diffuse optical imaging. Phys Med Biol. 2005; 50(4):R1-43. DOI: 10.1088/0031-9155/50/4/r01. View

Yao R, Intes X, Fang Q . Direct approach to compute Jacobians for diffuse optical tomography using perturbation Monte Carlo-based photon "replay". Biomed Opt Express. 2018; 9(10):4588-4603. PMC: 6179418. DOI: 10.1364/BOE.9.004588. View

Sassaroli A, Martelli F . Equivalence of four Monte Carlo methods for photon migration in turbid media. J Opt Soc Am A Opt Image Sci Vis. 2012; 29(10):2110-7. DOI: 10.1364/JOSAA.29.002110. View

Montejo L, Jia J, Kim H, Netz U, Blaschke S, Muller G . Computer-aided diagnosis of rheumatoid arthritis with optical tomography, Part 2: image classification. J Biomed Opt. 2013; 18(7):076002. PMC: 3710916. DOI: 10.1117/1.JBO.18.7.076002. View

Hirvi P, Kuutela T, Fang Q, Hannukainen A, Hyvonen N, Nissila I . Effects of atlas-based anatomy on modelled light transport in the neonatal head. Phys Med Biol. 2023; 68(13). PMC: 10460200. DOI: 10.1088/1361-6560/acd48c. View

Fang Q, Boas D . Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units. Opt Express. 2009; 17(22):20178-90. PMC: 2863034. DOI: 10.1364/OE.17.020178. View

Hanninen N, Pulkkinen A, Arridge S, Tarvainen T . Adaptive stochastic Gauss-Newton method with optical Monte Carlo for quantitative photoacoustic tomography. J Biomed Opt. 2022; 27(8). PMC: 8993421. DOI: 10.1117/1.JBO.27.8.083013. View

Darne C, Lu Y, Sevick-Muraca E . Small animal fluorescence and bioluminescence tomography: a review of approaches, algorithms and technology update. Phys Med Biol. 2013; 59(1):R1-64. DOI: 10.1088/0031-9155/59/1/R1. View

Amendola C, Maffeis G, Farina A, Spinelli L, Torricelli A, Pifferi A . Application limits of the scaling relations for Monte Carlo simulations in diffuse optics. Part 1: theory. Opt Express. 2024; 32(1):125-150. DOI: 10.1364/OE.507646. View

10.

Grosenick D, Rinneberg H, Cubeddu R, Taroni P . Review of optical breast imaging and spectroscopy. J Biomed Opt. 2016; 21(9):091311. DOI: 10.1117/1.JBO.21.9.091311. View

11.

Leino A, Lunttila T, Mozumder M, Pulkkinen A, Tarvainen T . Perturbation Monte Carlo Method for Quantitative Photoacoustic Tomography. IEEE Trans Med Imaging. 2020; 39(10):2985-2995. DOI: 10.1109/TMI.2020.2983129. View

12.

Durduran T, Choe R, Baker W, Yodh A . Diffuse Optics for Tissue Monitoring and Tomography. Rep Prog Phys. 2015; 73(7). PMC: 4482362. DOI: 10.1088/0034-4885/73/7/076701. View

13.

Klose A, Hielscher A . Iterative reconstruction scheme for optical tomography based on the equation of radiative transfer. Med Phys. 1999; 26(8):1698-707. DOI: 10.1118/1.598661. View

14.

Macdonald C, Arridge S, Powell S . Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models. J Biomed Opt. 2020; 25(8). PMC: 7426481. DOI: 10.1117/1.JBO.25.8.085002. View

15.

Hoshi Y, Yamada Y . Overview of diffuse optical tomography and its clinical applications. J Biomed Opt. 2016; 21(9):091312. DOI: 10.1117/1.JBO.21.9.091312. View

16.

Schweiger M, Arridge S, Nissila I . Gauss-Newton method for image reconstruction in diffuse optical tomography. Phys Med Biol. 2005; 50(10):2365-86. DOI: 10.1088/0031-9155/50/10/013. View

17.

Wang L, Jacques S, Zheng L . MCML--Monte Carlo modeling of light transport in multi-layered tissues. Comput Methods Programs Biomed. 1995; 47(2):131-46. DOI: 10.1016/0169-2607(95)01640-f. View

18.

Montejo L, Jia J, Kim H, Netz U, Blaschke S, Muller G . Computer-aided diagnosis of rheumatoid arthritis with optical tomography, Part 1: feature extraction. J Biomed Opt. 2013; 18(7):076001. PMC: 3710917. DOI: 10.1117/1.JBO.18.7.076001. View

19.

Kumar Y, Vasu R . Reconstruction of optical properties of low-scattering tissue using derivative estimated through perturbation Monte-Carlo method. J Biomed Opt. 2004; 9(5):1002-12. DOI: 10.1117/1.1778733. View

20.

Zhu C, Liu Q . Review of Monte Carlo modeling of light transport in tissues. J Biomed Opt. 2013; 18(5):50902. DOI: 10.1117/1.JBO.18.5.050902. View