Diffuse optical techniques have been extensively employed during the last years to retrieve the optical properties of tissue with and without inclusions. The usual approach is to use the diffusion approximation to the radiative transfer equation. However, if low- or non-diffusive regions are inside the studied volume, the diffusion approximation does not hold and the radiative transfer equation needs to be solved, which is computationally much more demanding. In this contribution, the problem of determining the optical properties of the whole volume of a turbid host medium containing both diffusive and non-diffusive inhomegeneities is examined. The situation reproduces clinical cases in which tumors and cysts can be present inside a studied tissue. To achieve this, an extended Kalman filter with compensation by Bayesian error modeling approach was adopted. Applying this technique, the diffusion approximation is used for calculations, reducing computation time, and discrepancies in the non-diffusive regions are compensated by the radiative transfer equation, thus keeping accuracy over the whole volume. The proposal is validated by phantom experiments showing very good results.
AlacamBYaziciBIntesX, et al.Extended Kalman filtering for the modeling and analysis of ICG pharmacokinetics in cancerous tumors using NIR optical methods. IEEE Trans Biomed Eng2006; 53: 1861–1871.
2.
AlacamBYaziciBIntesX, et al.Pharmacokinetic-rate images of indocyanine green for breast tumors using near-infrared optical methods. Phys Med Biol2008; 53: 837–859.
3.
AncoraDQiuLZacharakisG, et al.Noninvasive optical estimation of CSF thickness for brain-atrophy monitoring. Biomed Opt Express2018; 9: 4094–4112.
4.
AncoraDZacharopoulosARipollJ, et al.Fluorescence diffusion in the presence of optically clear tissues in a mouse head model. IEEE Trans Med Imaging2017; 36: 1086–1093.
5.
AncoraDZacharopoulosARipollJ, et al.Light propagation through weakly scattering media: a study of Monte Carlo vs. diffusion theory with application to neuroimaging. Proc SPIE2015; 9538: 95380G.
6.
ArridgeSR. Optical tomography in medical imaging. Inverse Probl1999; 15: R41.
7.
BaezGRPomaricoJAElicabeGE. An improved extended Kalman filter for diffuse optical tomography. Biomed Phys Eng Express2017; 3: 015013.
8.
BalGRenK. Generalized diffusion model in optical tomography with clear layers. J Opt Soc Am A2003; 20: 2355–2364.
9.
BolognaniSTubianaLZigliottoM. Extended Kalman filter tuning in sensorless PMSM drives. IEEE Trans Ind Appl2003; 39: 1741–1747.
10.
ChenXSunFYangD, et al.Hybrid simplified spherical harmonics with diffusion equation for light propagation in tissues. Phys Med Biol2015; 60: 6305–6322.
11.
ChenXZhangQYangD, et al.Hybrid radiosity-SP3 equation based bioluminescence tomography reconstruction for turbid medium with low- and non-scattering regions. J Appl Phys2014; 115: 024702.
12.
CoelhoPJCrouseillesNPereiraP, et al.Multi-scale methods for the solution of the radiative transfer equation. J Quant Spectrosc Radiat Transf2016; 172: 36–49.
13.
ContiniDMartelliFZaccantiG. Photon migration through a turbid slab described by a model based on diffusion approximation. I. Theory. Appl Opt1997; 36: 4587–4599.
14.
CubedduRPifferiATaroniP, et al.A solid tissue phantom for photon migration studies. Phys Med Biol1997; 42: 1971–1979.
15.
DiamondSGHuppertTJKolehmainenV, et al.Physiological system identification with the Kalman filter in diffuse optical tomography. Med Image Comput Comput Assist Interv2005; 8: 649–656.
16.
DehghaniHEamesMEYalavarthyPK, et al.Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction. Commun Numer Methods Eng2008; 25: 711–732.
17.
Di RoccoHOIriarteDIPomaricoJA, et al.Light transmittance in turbid media slabs used to model biological tissues; scaling laws. Optik Int J Light Electron Opt2008; 121: 435–441.
18.
FangQBoasDA. Monte Carlo simulation of photon migration in 3d turbid media accelerated by graphics processing units. Opt Express2009; 17: 20178–20190.
19.
FirbankMArridgeSRSchweigerM, et al.An investigation of light transport through scattering bodies with non-scattering regions. Phys Med Biol1996; 41: 767–783.
20.
GaoHZhaoH. A fast-forward solver of radiative transfer equation. Transp Theory Stat Phys2009; 38: 149–192.
21.
GibsonAPHebdenJCRileyJ, et al.Linear and nonlinear reconstruction for optical tomography of phantoms with nonscattering regions. Appl Opt2005; 44: 3925.
22.
GorpasDAndersson-EngelsS. Evaluation of a radiative transfer equation and diffusion approximation hybrid forward solver for fluorescence molecular imaging. J Biomed Opt2012; 17: 126010.
23.
HänninenNPulkkinenATarvainenT. Image reconstruction with reliability assessment in quantitative photoacoustic tomography. J Imaging2018; 4: 148.
24.
Heijblom M, Piras D, Maartens E, et al. Appearance of breast cysts in planar geometry photoacoustic mammography using 1064-nm excitation. Jour Biomed opt 2013 18(12): 126009/0-11.
25.
HoshiyaMSaitoE. Structural identification by extended Kalman filter. J Eng Mech1985; 110: 1757–1770.
26.
HumpherysJReddPWestJ. A fresh look at the Kalman filter. SIAM Rev2012; 54: 801–823.
27.
IshimaruA. Wave propagation and scattering in random media, New York: Academic Press, 1978.
28.
JacquesSL. Optical properties of biological tissues: a review. Phys Med Biol2013; 58: R37–R61.
29.
KaipioJSomersaloE. Statistical and computational inverse problems, Dordrecht: Springer, 2005.
30.
KaipioJSomersaloE. Statistical inverse problems: discretization, model reduction and inverse crimes. J Comput Appl Math2007; 198: 493–504.
31.
KalmanRE. A new approach to linear filtering and prediction problems. J Fluids Eng1960; 82: 35–45.
32.
KienleA. Light diffusion through a turbid parallelepiped. J Opt Soc Am A2005; 22: 1883–1888.
33.
KienleAPattersonMS. Improved solutions of the steady-state and the time-resolved diffusion equations for reflectance from a semi-infinite turbid medium. J Opt Soc Am A1997; 14: 246–254.
34.
LiemertAKienleA. Exact and efficient solution of the radiative transport equation for the semi-infinite medium. Sci Rep2013; 3: 1–7.
35.
Kolehmainen V, Arridge SR, Kaipio JP, et al. Approximation errors and model reduction in optical tomography. In: Annual international conference of the IEEE engineering in medicine and biology—proceedings, New York, NY, 30 August–3 September 2006, pp. 2659–2662. Piscataway, NJ: IEEE Engineering in Medicine and Biology Society.
36.
LehtikangasOTarvainenT. Hybrid forward-peaked-scattering-diffusion approximations for light propagation in turbid media with low-scattering regions. J Quant Spectrosc Radiat Transf2013; 116: 132–144.
37.
LiemertAKienleA. Analytical solution of the radiative transfer equation for infinite-space fluence. Phys Rev A2011; 83: 015804.
38.
NissilINoponenTKotilahtiK, et al.Instrumentation and calibration methods for the multichannel measurement of phase and amplitude in optical tomography. Rev Sci Instrum2005; 76: 044302.
39.
PardiniPPomaricoJAIriarteDI. A novelty technique for the fabrication of biomedical optics phantoms with cyst-mimicking inclusions. J Near Infrared Spectrosc2017; 0: 0967033516686045.
40.
RipollJNieto-VesperinasMArridgeSR, et al.Boundary conditions for light propagation in diffusive media with nonscattering regions. J Opt Soc Am A2000; 17: 1671–1681.
41.
RogerMCaliotCCrouseillesN, et al.A hybrid transport-diffusion model for radiative transfer in absorbing and scattering media. J Comput Phys2014; 275: 346–362.
42.
SabatiniAM. Quaternion-based extended Kalman filter for determining orientation by inertial and magnetic sensing. IEEE Rev Biomed Eng2006; 53: 1346–1356.
43.
SchweigerMArridgeSR. The Toast++ software suite for forward and inverse modeling in optical tomography. J Biomed Opt2014; 19: 040801.
44.
SchweigerMArridgeSRNissilaI. Gauss–Newton method for image reconstruction in diffuse optical tomography. Phys Med Biol2005; 50: 2365.
45.
SimonD. Optimal state rstimation: Kalman, H Infinity, and nonlinear approaches, Hoboken: John Wiley & Sons, 2006.
46.
TarvainenTKolehmainenVArridgeSR, et al.Image reconstruction in diffuse optical tomography using the coupled radiative transport-diffusion model. J Quant Spectrosc Radiat Transf2011; 112: 2600–2608.
47.
TarvainenTKolehmainenVPulkkinenA, et al.An approximation error approach for compensating for modelling errors between the radiative transfer equation and the diffusion approximation in diffuse optical tomography. Inverse Probl2010; 26: 15005–15018.
48.
TarvainenTVauhkonenMKolehmainenV, et al.Coupled radiative transfer equation and diffusion approximation model for photon migration in turbid medium with low-scattering and non-scattering regions. Phys Med Biol2005; 50: 4913–4930.
49.
SerraMVWCarboneNDi RoccoHO, et al.Diffuse light transmission profiles obtained using CW: A comparative analysis with time resolved experiments. Optik Int J Light Electron Opt2014; 125: 3507–3513.
50.
XuYZhangQJiangH. Optical image reconstruction of non-scattering and low scattering heterogeneities in turbid media based on the diffusion approximation model. J Opt A Pure Appl Opt2003; 6: 29–35.
51.
YangDChenXCaoX, et al.Performance investigation of SP3 and diffusion approximation for three-dimensional whole-body optical imaging of small animals. Med Biol Eng Comput2015; 53: 805–814.
52.
YangDChenXPengZ, et al.Light transport in turbid media with non-scattering, low-scattering and high absorption heterogeneities based on hybrid simplified spherical harmonics with radiosity model. Biomed Optics Exp2013; 4: 2209–2223.
53.
YangDChenXRenS, et al.Influence investigation of a void region on modeling light propagation in a heterogeneous medium. Appl Opt2013; 52: 400–408.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
