1
|
Jha AK, Zhu Y, Arridge S, Wong DF, Rahmim A. Incorporating reflection boundary conditions in the Neumann series radiative transport equation: application to photon propagation and reconstruction in diffuse optical imaging. BIOMEDICAL OPTICS EXPRESS 2018; 9:1389-1407. [PMID: 29675291 PMCID: PMC5905895 DOI: 10.1364/boe.9.001389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/28/2017] [Revised: 11/28/2017] [Accepted: 12/26/2017] [Indexed: 05/11/2023]
Abstract
We propose a formalism to incorporate boundary conditions in a Neumann-series-based radiative transport equation. The formalism accurately models the reflection of photons at the tissue-external medium interface using Fresnel's equations. The formalism was used to develop a gradient descent-based image reconstruction technique. The proposed methods were implemented for 3D diffuse optical imaging. In computational studies, it was observed that the average root-mean-square error (RMSE) for the output images and the estimated absorption coefficients reduced by 38% and 84%, respectively, when the reflection boundary conditions were incorporated. These results demonstrate the importance of incorporating boundary conditions that model the reflection of photons at the tissue-external medium interface.
Collapse
Affiliation(s)
- Abhinav K. Jha
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD,
USA
- Department of Biomedical Engineering, Washington University in St. Louis, St. Louis, MO, USA
- Mallinckrodt Institute of Radiology, Washington University in St. Louis, St. Louis, MO, USA
| | - Yansong Zhu
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD,
USA
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD, USA
| | - Simon Arridge
- Department of Computer Science, University College London, Gower Street, London WC1E 6BT, UK
| | - Dean F. Wong
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD,
USA
- Department of Neuroscience, Johns Hopkins University, Baltimore, MD, USA
- Department of Psychiatry and Behavioral Sciences, Johns Hopkins University, Baltimore, MD, USA
- Department of Neurology, Johns Hopkins University, Baltimore, MD, USA
| | - Arman Rahmim
- Department of Radiology and Radiological Sciences, Johns Hopkins University, Baltimore, MD,
USA
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, MD, USA
| |
Collapse
|
2
|
Zhu Y, Jha AK, Dreyer JK, Le HND, Kang JU, Roland PE, Wong DF, Rahmim A. A three-step reconstruction method for fluorescence molecular tomography based on compressive sensing. PROCEEDINGS OF SPIE--THE INTERNATIONAL SOCIETY FOR OPTICAL ENGINEERING 2017; 10059. [PMID: 28596634 DOI: 10.1117/12.2252664] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Fluorescence molecular tomography (FMT) is a promising tool for real time in vivo quantification of neurotransmission (NT) as we pursue in our BRAIN initiative effort. However, the acquired image data are noisy and the reconstruction problem is ill-posed. Further, while spatial sparsity of the NT effects could be exploited, traditional compressive-sensing methods cannot be directly applied as the system matrix in FMT is highly coherent. To overcome these issues, we propose and assess a three-step reconstruction method. First, truncated singular value decomposition is applied on the data to reduce matrix coherence. The resultant image data are input to a homotopy-based reconstruction strategy that exploits sparsity via ℓ1 regularization. The reconstructed image is then input to a maximum-likelihood expectation maximization (MLEM) algorithm that retains the sparseness of the input estimate and improves upon the quantitation by accurate Poisson noise modeling. The proposed reconstruction method was evaluated in a three-dimensional simulated setup with fluorescent sources in a cuboidal scattering medium with optical properties simulating human brain cortex (reduced scattering coefficient: 9.2 cm-1, absorption coefficient: 0.1 cm-1) and tomographic measurements made using pixelated detectors. In different experiments, fluorescent sources of varying size and intensity were simulated. The proposed reconstruction method provided accurate estimates of the fluorescent source intensity, with a 20% lower root mean square error on average compared to the pure-homotopy method for all considered source intensities and sizes. Further, compared with conventional ℓ2 regularized algorithm, overall, the proposed method reconstructed substantially more accurate fluorescence distribution. The proposed method shows considerable promise and will be tested using more realistic simulations and experimental setups.
Collapse
Affiliation(s)
- Yansong Zhu
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, USA
| | - Abhinav K Jha
- Department of Radiology, Johns Hopkins University, Baltimore, USA
| | - Jakob K Dreyer
- Department of Neuroscience and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Hanh N D Le
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, USA
| | - Jin U Kang
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, USA
| | - Per E Roland
- Department of Neuroscience and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Dean F Wong
- Department of Radiology, Johns Hopkins University, Baltimore, USA
| | - Arman Rahmim
- Department of Electrical and Computer Engineering, Johns Hopkins University, Baltimore, USA.,Department of Radiology, Johns Hopkins University, Baltimore, USA
| |
Collapse
|