1
|
Jia M, Jiang J, Ma W, Li C, Wang S, Zhao H, Gao F. Accelerating nonlinear reconstruction in laminar optical tomography by use of recursive SVD inversion. BIOMEDICAL OPTICS EXPRESS 2017; 8:4275-4293. [PMID: 28966864 PMCID: PMC5611940 DOI: 10.1364/boe.8.004275] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/27/2017] [Revised: 08/21/2017] [Accepted: 08/28/2017] [Indexed: 06/07/2023]
Abstract
Image reconstruction in the most model-based biophotonic imaging modalities essentially poses an ill-posed nonlinear inverse problem, which has been effectively tackled in the diffusion-approximation-satisfied scenarios such as diffuse optical tomography. Nevertheless, a nonlinear implementation in high-resolution laminar optical tomography (LOT) is normally computationally-costly due to its strong dependency on a dense source-detector configuration and a physically-rigorous photon-transport model. To circumvent the adversity, we herein propose a practical nonlinear LOT approach to the absorption reconstruction. The scheme takes advantage of the numerical stability of the singular value decomposition (SVD) for the ill-posed linear inversion, and is accelerated by adopting an explicitly recursive strategy for the time-consuming repeated SVD inversion, which is based on a scaled expression of the sensitivity matrix. Experiments demonstrate that the proposed methodology can perform as well as the traditional nonlinear one, while the computation time of the former is merely 26.27% of the later on average.
Collapse
Affiliation(s)
- Mengyu Jia
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, 300072, China
| | - Jingying Jiang
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, 300072, China
| | - Wenjuan Ma
- Tianjin Medical University Cancer Institute & Hospital Tianjin, 300072, China
| | - Chenxi Li
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, 300072, China
| | - Shuang Wang
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, 300072, China
| | - Huijuan Zhao
- College of Precision Instrument and Optoelectronics Engineering, Tianjin University, Tianjin, 300072, China
| | - Feng Gao
- Tianjin Key Laboratory of Biomedical Detecting Techniques and Instruments, Tianjin, 300072, China
| |
Collapse
|
2
|
Ancora D, Zacharopoulos A, Ripoll J, Zacharakis G. Fluorescence Diffusion in the Presence of Optically Clear Tissues in a Mouse Head Model. IEEE TRANSACTIONS ON MEDICAL IMAGING 2017; 36:1086-1093. [PMID: 28055860 DOI: 10.1109/tmi.2016.2646518] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Diffuse Optical Tomography commonly neglects or assumes as insignificant the presence of optically clear regions in biological tissues, estimating their contribution as a small perturbation to light transport. The inaccuracy introduced by this practice is examined in detail in the context of a complete, based on realistic geometry, virtual fluorescence Diffuse Optical Tomography experiment where a mouse head is imaged in the presence of cerebral spinal fluid. Despite the small thickness of such layer, we point out that an error is introduced when neglecting it from the model with possibly reduction in the accuracy of the reconstruction and localization of the fluorescence distribution within the brain. The results obtained in the extensive study presented here suggest that fluorescence diffuse neuroimaging studies can be improved in terms of quantitative and qualitative reconstruction by accurately taking into account optically transparent regions especially in the cases where the reconstruction is aided by the prior knowledge of the structural geometry of the specimen. Thus, this has only recently become an affordable choice, thanks to novel computation paradigms that allow to run Monte Carlo photon propagation on a simple graphic card, hence speeding up the process a thousand folds compared to CPU-based solutions.
Collapse
|
3
|
Detectability of absorption and reduced scattering coefficients in frequency-domain measurements using a realistic head phantom. SENSORS 2012; 13:152-64. [PMID: 23262479 PMCID: PMC3574670 DOI: 10.3390/s130100152] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/07/2012] [Revised: 12/17/2012] [Accepted: 12/18/2012] [Indexed: 11/17/2022]
Abstract
Detection limits of the changes in absorption and reduced scattering coefficients were investigated using a frequency-domain near-infrared system in a realistic head phantom. The results were quantified in terms of the maximum detectable depth for different activation volumes in the range of 0.8-20 microliters. The non-linear relation between the maximum detectable depth and the magnitude of changes in the absorption coefficient conform well with the Born approximation to the diffusion equation. The minimal detectable changes in the reduced scattering coefficient measured in terms of the phase signal were found to be approximately twice as large as that of the absorption coefficient using the AC signal for the same volume and at the same depth. The phase delay, which can be used to quantify the fast neuronal optical response in the human brain, showed a linear dependence on the reciprocal of the reduced scattering coefficient, as predicted by the Rytov approximation.
Collapse
|
4
|
Bedside optical imaging of occipital resting-state functional connectivity in neonates. Neuroimage 2011; 59:2529-38. [PMID: 21925609 DOI: 10.1016/j.neuroimage.2011.08.094] [Citation(s) in RCA: 76] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2011] [Revised: 08/07/2011] [Accepted: 08/30/2011] [Indexed: 11/20/2022] Open
Abstract
Resting-state networks derived from temporal correlations of spontaneous hemodynamic fluctuations have been extensively used to elucidate the functional organization of the brain in adults and infants. We have previously developed functional connectivity diffuse optical tomography methods in adults, and we now apply these techniques to study functional connectivity in newborn infants at the bedside. We present functional connectivity maps in the occipital cortices obtained from healthy term-born infants and premature infants, including one infant with an occipital stroke. Our results suggest that functional connectivity diffuse optical tomography has potential as a valuable clinical tool for the early detection of functional deficits and for providing prognostic information on future development.
Collapse
|
5
|
Depth sensitivity analysis of functional near-infrared spectroscopy measurement using three-dimensional Monte Carlo modelling-based magnetic resonance imaging. Lasers Med Sci 2010; 25:431-8. [DOI: 10.1007/s10103-010-0754-4] [Citation(s) in RCA: 39] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2009] [Accepted: 01/04/2010] [Indexed: 10/19/2022]
|
6
|
Heiskala J, Hiltunen P, Nissilä I. Significance of background optical properties, time-resolved information and optode arrangement in diffuse optical imaging of term neonates. Phys Med Biol 2009; 54:535-54. [PMID: 19124950 DOI: 10.1088/0031-9155/54/3/005] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
The significance of accurate knowledge of background optical properties and time-resolved information in reconstructing images of hemodynamic changes in the neonatal brain from diffuse optical imaging data was studied using Monte Carlo (MC) simulation. A segmented anatomical magnetic resonance (MR) image and literature-derived optical properties for each tissue type were used to create a voxel-based anatomical model. Small absorbing perturbations were introduced into the anatomical model to simulate localized hemodynamic responses related to brain activation. Perturbation MC (pMC) was used as the primary method of image reconstruction. For comparison, reconstructions were also performed using the finite element method (FEM) to solve the diffusion approximation (DA) to the radiative transfer equation (RTE). The effect of optode layout was investigated using three different grids. Of the factors studied, the density of the optode grid was found to have the greatest effect on image quality. The use of time-resolved information significantly improved the spatial accuracy with all optode grids. Adequate knowledge and modeling of the optical properties of the background was found to significantly improve the spatial accuracy of the reconstructed images and make the recovery of contrast of absorption changes more consistent over simplified modeling. Localization accuracy of small perturbations was found to be 2-3 mm with accurate a priori knowledge of the background optical properties, when a grid with high optode density (>1 optode cm(-2)) was used.
Collapse
Affiliation(s)
- J Heiskala
- BioMag Laboratory, HUSLAB, Helsinki University Central Hospital, PO Box 340, FI-00029 HUS, Finland.
| | | | | |
Collapse
|
7
|
Patwardhan SV, Culver JP. Quantitative diffuse optical tomography for small animals using an ultrafast gated image intensifier. JOURNAL OF BIOMEDICAL OPTICS 2008; 13:011009. [PMID: 18315358 DOI: 10.1117/1.2830656] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The quantitative accuracy of fluorescence and bioluminescence imaging of small animals can be improved by knowledge of the in situ optical properties of each animal. Obtaining in situ optical property maps is challenging, however, due to short propagation distances, requirements for high dynamic range, and the need for dense spatial, temporal, and spectral sampling. Using an ultrafast gated image intensifier and a pulsed laser source, we have developed a small animal diffuse optical tomography system with multiple synthetic modulation frequencies up to >1 GHz. We show that amplitude and phase measurements with useful contrast-to-noise ratios can be obtained for modulation frequencies over the range of approximately 250 to 1250 MHz. Experiments with tissue simulating phantoms demonstrate the feasibility of reconstructing the absorption and scattering optical properties in a small animal imaging system.
Collapse
Affiliation(s)
- Sachin V Patwardhan
- Washington University School of Medicine, Department of Radiology, Mallinckrodt Institute of Radiology, St. Louis, Missouri 63110, USA
| | | |
Collapse
|
8
|
Lasker JM, Fong CJ, Ginat DT, Dwyer E, Hielscher AH. Dynamic optical imaging of vascular and metabolic reactivity in rheumatoid joints. JOURNAL OF BIOMEDICAL OPTICS 2007; 12:052001. [PMID: 17994887 DOI: 10.1117/1.2798757] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/20/2023]
Abstract
Dynamic optical imaging is increasingly applied to clinically relevant areas such as brain and cancer imaging. In this approach, some external stimulus is applied and changes in relevant physiological parameters (e.g., oxy- or deoxyhemoglobin concentrations) are determined. The advantage of this approach is that the prestimulus state can be used as a reference or baseline against which the changes can be calibrated. Here we present the first application of this method to the problem of characterizing joint diseases, especially effects of rheumatoid arthritis (RA) in the proximal interphalangeal finger joints. Using a dual-wavelength tomographic imaging system together with previously implemented model-based iterative image reconstruction schemes, we have performed initial dynamic imaging case studies on a limited number of healthy volunteers and patients diagnosed with RA. Focusing on three cases studies, we illustrated our major finds. These studies support our hypothesis that differences in the vascular reactivity exist between affected and unaffected joints.
Collapse
Affiliation(s)
- Joseph M Lasker
- Columbia University, Department of Biomedical Engineering, 500 West 120th Street, ET351 Mudd Building, MC8904, New York, New York 10027, USA
| | | | | | | | | |
Collapse
|
9
|
Gibson AP, Austin T, Everdell NL, Schweiger M, Arridge SR, Meek JH, Wyatt JS, Delpy DT, Hebden JC. Three-dimensional whole-head optical tomography of passive motor evoked responses in the neonate. Neuroimage 2005; 30:521-8. [PMID: 16246586 DOI: 10.1016/j.neuroimage.2005.08.059] [Citation(s) in RCA: 103] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2005] [Revised: 07/08/2005] [Accepted: 08/30/2005] [Indexed: 11/20/2022] Open
Abstract
Optical tomography has been used to reconstruct three-dimensional images of the entire neonatal head during motor evoked responses. Data were successfully acquired during passive movement of each arm on four out of six infants examined, from which eight sets of bilateral images of hemodynamic parameters were reconstructed. Six out of the eight images showed the largest change in total hemoglobin in the region of the contralateral motor cortex. The mean distance between the peak response in the image and the estimated position of the contralateral motor cortex was 10.8 mm. These results suggest that optical tomography may provide an appropriate technique for non-invasive cot-side imaging of brain function.
Collapse
Affiliation(s)
- A P Gibson
- Department of Medical Physics and Bioengineering, University College London, Malet Place Engineering Building, London WC1E 6BT, UK.
| | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Diamond SG, Huppert TJ, Kolehmainen V, Franceschini MA, Kaipio JP, Arridge SR, Boas DA. Dynamic physiological modeling for functional diffuse optical tomography. Neuroimage 2005; 30:88-101. [PMID: 16242967 PMCID: PMC2670202 DOI: 10.1016/j.neuroimage.2005.09.016] [Citation(s) in RCA: 82] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2005] [Accepted: 09/14/2005] [Indexed: 11/20/2022] Open
Abstract
Diffuse optical tomography (DOT) is a noninvasive imaging technology that is sensitive to local concentration changes in oxy- and deoxyhemoglobin. When applied to functional neuroimaging, DOT measures hemodynamics in the scalp and brain that reflect competing metabolic demands and cardiovascular dynamics. The diffuse nature of near-infrared photon migration in tissue and the multitude of physiological systems that affect hemodynamics motivate the use of anatomical and physiological models to improve estimates of the functional hemodynamic response. In this paper, we present a linear state-space model for DOT analysis that models the physiological fluctuations present in the data with either static or dynamic estimation. We demonstrate the approach by using auxiliary measurements of blood pressure variability and heart rate variability as inputs to model the background physiology in DOT data. We evaluate the improvements accorded by modeling this physiology on ten human subjects with simulated functional hemodynamic responses added to the baseline physiology. Adding physiological modeling with a static estimator significantly improved estimates of the simulated functional response, and further significant improvements were achieved with a dynamic Kalman filter estimator (paired t tests, n=10, P<0.05). These results suggest that physiological modeling can improve DOT analysis. The further improvement with the Kalman filter encourages continued research into dynamic linear modeling of the physiology present in DOT. Cardiovascular dynamics also affect the blood-oxygen-dependent (BOLD) signal in functional magnetic resonance imaging (fMRI). This state-space approach to DOT analysis could be extended to BOLD fMRI analysis, multimodal studies and real-time analysis.
Collapse
Affiliation(s)
- Solomon Gilbert Diamond
- Massachusetts General Hospital, Athinoula A. Martinos Center for Biomedical Imaging, Charlestown, MA 02129, USA.
| | | | | | | | | | | | | |
Collapse
|