Influence of source-detector separation on diffuse correlation spectroscopy measurements of cerebral blood flow with a multilayered analytical model

Study of Time-Resolved Dynamics in Turbid Medium Using a Single-Cavity Dual-Comb Laser

In measuring cerebral blood flow (CBF) noninvasively using optical techniques, diffusing-wave spectroscopy is often combined with near-infrared spectroscopy to obtain a reliable blood flow index. Measuring the blood flow index at a determined depth remains the ultimate goal. In this study, we present a simple approach using dual-comb lasers where we simultaneously measure the absorption coefficient (μa), the reduced scattering coefficient (μs '), and dynamic properties. This system can also effectively differentiate dynamics from various depths, which is crucial for analyzing multilayer dynamics. For CBF measurements, this capability is particularly valuable as it helps mitigate the influence of the scalp and skull, thereby enhancing the specificity of deep tissue.

Modified Beer-Lambert algorithm to measure pulsatile blood flow, critical closing pressure, and intracranial hypertension

We introduce a frequency-domain modified Beer-Lambert algorithm for diffuse correlation spectroscopy to non-invasively measure flow pulsatility and thus critical closing pressure (CrCP). Using the same optical measurements, CrCP was obtained with the new algorithm and with traditional nonlinear diffusion fitting. Results were compared to invasive determination of intracranial pressure (ICP) in piglets (n = 18). The new algorithm better predicted ICP elevations; the area under curve (AUC) from logistic regression analysis was 0.85 for ICP ≥ 20 mmHg. The corresponding AUC for traditional analysis was 0.60. Improved diagnostic performance likely results from better filtering of extra-cerebral tissue contamination and measurement noise.

A Device-on-Chip Solution for Real-Time Diffuse Correlation Spectroscopy Using FPGA

Diffuse correlation spectroscopy (DCS) is a non-invasive technology for the evaluation of blood perfusion in deep tissue. However, it requires high computational resources for data analysis, which poses challenges in its implementation for real-time applications. To address the unmet need, we developed a novel device-on-chip solution that fully integrates all the necessary computational components needed for DCS. It takes the output of a photon detector and determines the blood flow index (BFI). It is implemented on a field-programmable gate array (FPGA) chip including a multi-tau correlator for the calculation of the temporal light intensity autocorrelation function and a DCS analyzer to perform the curve fitting operation that derives the BFI at a rate of 6000 BFIs/s. The FPGA DCS system was evaluated against a lab-standard DCS system for both phantom and cuff ischemia studies. The results indicate that the autocorrelation of the light correlation and BFI from both the FPGA DCS and the reference DCS matched well. Furthermore, the FPGA DCS system was able to achieve a measurement rate of 50 Hz and resolve pulsatile blood flow. This can significantly lower the cost and footprint of the computational components of DCS and pave the way for portable, real-time DCS systems.

Comparing the performance potential of speckle contrast optical spectroscopy and diffuse correlation spectroscopy for cerebral blood flow monitoring using Monte Carlo simulations in realistic head geometries

Significance

The non-invasive measurement of cerebral blood flow based on diffuse optical techniques has seen increased interest as a research tool for cerebral perfusion monitoring in critical care and functional brain imaging. Diffuse correlation spectroscopy (DCS) and speckle contrast optical spectroscopy (SCOS) are two such techniques that measure complementary aspects of the fluctuating intensity signal, with DCS quantifying the temporal fluctuations of the signal and SCOS quantifying the spatial blurring of a speckle pattern. With the increasing interest in the use of these techniques, a thorough comparison would inform new adopters of the benefits of each technique.

Aim

We systematically evaluate the performance of DCS and SCOS for the measurement of cerebral blood flow.

Approach

Monte Carlo simulations of dynamic light scattering in an MRI-derived head model were performed. For both DCS and SCOS, estimates of sensitivity to cerebral blood flow changes, coefficient of variation of the measured blood flow, and the contrast-to-noise ratio of the measurement to the cerebral perfusion signal were calculated. By varying complementary aspects of data collection between the two methods, we investigated the performance benefits of different measurement strategies, including altering the number of modes per optical detector, the integration time/fitting time of the speckle measurement, and the laser source delivery strategy.

Results

Through comparison across these metrics with simulated detectors having realistic noise properties, we determine several guiding principles for the optimization of these techniques and report the performance comparison between the two over a range of measurement properties and tissue geometries. We find that SCOS outperforms DCS in terms of contrast-to-noise ratio for the cerebral blood flow signal in the ideal case simulated here but note that SCOS requires careful experimental calibrations to ensure accurate measurements of cerebral blood flow.

Conclusion

We provide design principles by which to evaluate the development of DCS and SCOS systems for their use in the measurement of cerebral blood flow.

Quantification of blood flow index in diffuse correlation spectroscopy using a robust deep learning method

Significance

Diffuse correlation spectroscopy (DCS) is a powerful, noninvasive optical technique for measuring blood flow. Traditionally the blood flow index (BFi) is derived through nonlinear least-square fitting the measured intensity autocorrelation function (ACF). However, the fitting process is computationally intensive, susceptible to measurement noise, and easily influenced by optical properties (absorption coefficient μ a and reduced scattering coefficient μ s ' ) and scalp and skull thicknesses.

Aim

We aim to develop a data-driven method that enables rapid and robust analysis of multiple-scattered light's temporal ACFs. Moreover, the proposed method can be applied to a range of source-detector distances instead of being limited to a specific source-detector distance.

Approach

We present a deep learning architecture with one-dimensional convolution neural networks, called DCS neural network (DCS-NET), for BFi and coherent factor (β ) estimation. This DCS-NET was performed using simulated DCS data based on a three-layer brain model. We quantified the impact from physiologically relevant optical property variations, layer thicknesses, realistic noise levels, and multiple source-detector distances (5, 10, 15, 20, 25, and 30 mm) on BFi and β estimations among DCS-NET, semi-infinite, and three-layer fitting models.

Results

DCS-NET shows a much faster analysis speed, around 17,000-fold and 32-fold faster than the traditional three-layer and semi-infinite models, respectively. It offers higher intrinsic sensitivity to deep tissues compared with fitting methods. DCS-NET shows excellent anti-noise features and is less sensitive to variations of μ a and μ s ' at a source-detector separation of 30 mm. Also, we have demonstrated that relative BFi (rBFi) can be extracted by DCS-NET with a much lower error of 8.35%. By contrast, the semi-infinite and three-layer fitting models result in significant errors in rBFi of 43.76% and 19.66%, respectively.

Conclusions

DCS-NET can robustly quantify blood flow measurements at considerable source-detector distances, corresponding to much deeper biological tissues. It has excellent potential for hardware implementation, promising continuous real-time blood flow measurements.

Pulsatile microvascular cerebral blood flow waveforms change with intracranial compliance and age

Significance

Diffuse correlation spectroscopy (DCS) is an optical method to measure relative changes in cerebral blood flow (rCBF) in the microvasculature. Each heartbeat generates a pulsatile signal with distinct morphological features that we hypothesized to be related to intracranial compliance (ICC).

Aim

We aim to study how three features of the pulsatile rCBF waveforms: the augmentation index (AIx), the pulsatility index, and the area under the curve, change with respect to ICC. We describe ICC as a combination of vascular compliance and extravascular compliance.

Approach

Since patients with Chiari malformations (CM) (n = 30 ) have been shown to have altered extravascular compliance, we compare the morphology of rCBF waveforms in CM patients with age-matched healthy control (n = 30 ).

Results

AIx measured in the supine position was significantly less in patients with CM compared to healthy controls (p < 0.05 ). Since physiologic aging also leads to changes in vessel stiffness and intravascular compliance, we evaluate how the rCBF waveform changes with respect to age and find that the AIx feature was strongly correlated with age (R healthy subjects = - 0.63 , R preoperative   CM   patient = - 0.70 , and R postoperative   CM   patients = - 0.62 , p < 0.01 ).

Conclusions

These results suggest that the AIx measured in the cerebral microvasculature using DCS may be correlated to changes in ICC.

Comparison of diffuse correlation spectroscopy analytical models for measuring cerebral blood flow in adults

Significance

Although multilayer analytical models have been proposed to enhance brain sensitivity of diffuse correlation spectroscopy (DCS) measurements of cerebral blood flow, the traditional homogeneous model remains dominant in clinical applications. Rigorous in vivo comparison of these analytical models is lacking.

Aim

We compare the performance of different analytical models to estimate a cerebral blood flow index (CBFi) with DCS in adults.

Approach

Resting-state data were obtained on a cohort of 20 adult patients with subarachnoid hemorrhage. Data at 1 and 2.5 cm source-detector separations were analyzed with the homogenous, two-layer, and three-layer models to estimate scalp blood flow index and CBFi. The performance of each model was quantified via fitting convergence, fit stability, brain-to-scalp flow ratio (BSR), and correlation with transcranial Doppler ultrasound (TCD) measurements of cerebral blood flow velocity in the middle cerebral artery (MCA).

Results

The homogeneous model has the highest pass rate (100%), lowest coefficient of variation (CV) at rest (median [IQR] at 1 Hz of 0.18 [0.13, 0.22]), and most significant correlation with MCA blood flow velocities (R s = 0.59 , p = 0.010 ) compared with both the two- and three-layer models. The multilayer model pass rate was significantly correlated with extracerebral layer thicknesses. Discarding datasets with non-physiological BSRs increased the correlation between DCS measured CBFi and TCD measured MCA velocities for all models.

Conclusions

We found that the homogeneous model has the highest pass rate, lowest CV at rest, and most significant correlation with MCA blood flow velocities. Results from the multilayer models should be taken with caution because they suffer from lower pass rates and higher coefficients of variation at rest and can converge to non-physiological values for CBFi. Future work is needed to validate these models in vivo, and novel approaches are merited to improve the performance of the multimodel models.

Deep-learning-based separation of shallow and deep layer blood flow rates in diffuse correlation spectroscopy

Diffuse correlation spectroscopy faces challenges concerning the contamination of cutaneous and deep tissue blood flow. We propose a long short-term memory network to directly quantify the flow rates of shallow and deep-layer tissues. By exploiting the different contributions of shallow and deep-layer flow rates to auto-correlation functions, we accurately predict the shallow and deep-layer flow rates (RMSE = 0.047 and 0.034 ml/min/100 g of simulated tissue, R2 = 0.99 and 0.99, respectively) in a two-layer flow phantom experiment. This approach is useful in evaluating the blood flow responses of active muscles, where both cutaneous and deep-muscle blood flow increase with exercise.

scatterBrains: an open database of human head models and companion optode locations for realistic Monte Carlo photon simulations

Significance

Monte Carlo (MC) simulations are currently the gold standard in the near-infrared and diffuse correlation spectroscopy (NIRS/DCS) communities for generating light transport paths through tissue. However, realistic and diverse models that capture complex tissue layers are not easily available to all; moreover, manually placing optodes on such models can be tedious and time consuming. Such limitations may hinder the adoption of representative models for basic simulations and the use of these models for large-scale simulations, e.g., for training machine learning algorithms.

Aim

We aim to provide the NIRS/DCS communities with an open-source, user-friendly database of morphologically and optically realistic head models, as well as a succinct software pipeline to prepare these models for mesh-based Monte Carlo simulations of light transport.

Approach

Sixteen anatomical models were created from segmented T1-weighted magnetic resonance imaging (MRI) head scans and converted to tetrahedral mesh volumes. Approximately 800 companion scalp surface locations were distributed on each model, comprising full head coverage. A pipeline was created to place custom source and optical detectors at each location, and guidance is provided on how to use these parameters to set up MC simulations.

Results

The models, head surface locations, and all associated code are freely available under the scatterBrains project on Github.

Conclusions

The NIRS/DCS community benefits from having shared resources for conducting MC simulations on realistic head geometries. We hope this will make MRI-based head models and virtual optode placement easily accessible to all. Contributions to the database are welcome and encouraged.

Influence of oversimplifying the head anatomy on cerebral blood flow measurements with diffuse correlation spectroscopy

Significance

Diffuse correlation spectroscopy (DCS) is an emerging optical modality for non-invasive assessment of an index of regional cerebral blood flow. By the nature of this noninvasive measurement, light must pass through extracerebral layers (i.e., skull, scalp, and cerebral spinal fluid) before detection at the tissue surface. To minimize the contribution of these extracerebral layers to the measured signal, an analytical model has been developed that treats the head as a series of three parallel and infinitely extending slabs (mimicking scalp, skull, and brain). The three-layer model has been shown to provide a significant improvement in cerebral blood flow estimation over the typically used model that treats the head as a bulk homogenous medium. However, the three-layer model is still a gross oversimplification of the head geometry that ignores head curvature, the presence of cerebrospinal fluid (CSF), and heterogeneity in layer thickness.

Aim

Determine the influence of oversimplifying the head geometry on cerebral blood flow estimated with the three-layer model.

Approach

Data were simulated with Monte Carlo in a four-layer slab medium and a three-layer sphere medium to isolate the influence of CSF and curvature, respectively. Additionally, simulations were performed on magnetic resonance imaging (MRI) head templates spanning a wide-range of ages. Simulated data were fit to both the homogenous and three-layer model for CBF. Finally, to mitigate the errors in potential CBF estimation due to the difficulty in defining layer thickness, we investigated an approach to identify an equivalent, "optimized" thickness via a pressure modulation.

Results

Both head curvature and failing to account for CSF lead to significant errors in the estimation of CBF. However, the effect of curvature and CSF on relative changes in CBF is minimal. Further, we found that CBF was underestimated in all MRI-templates, although the magnitude of these underestimations was highly influenced by small variations in the source and detector optode positioning. The optimized thickness obtained from pressure modulation did not improve estimation accuracy of CBF, although it did significantly improve the estimation accuracy of relative changes in CBF.

Conclusions

In sum, these findings suggest that the three-layer model holds promise for improving estimation of relative changes in cerebral blood flow; however, estimations of absolute cerebral blood flow with the approach should be viewed with caution given that it is difficult to account for appreciable sources of error, such as curvature and CSF.

High-density diffuse correlation tomography with enhanced depth localization and minimal surface artefacts

A spatially weighted filter applied to both the measurement and the Jacobian is proposed for high-density diffuse correlation tomography (DCT) to remove unwanted extracerebral interferences and artefacts along with better depth localization in the reconstructed blood flow images. High-density DCT is implemented by appropriate modification of recently introduced Multi-speckle Diffuse Correlation Spectroscopy (M-DCS) system. Additionally, we have used autocorrelation measurements at multiple delay-times in an iterative manner to improve the reconstruction results. The proposed scheme has been validated by simulations, phantom experiments and in-vivo human experiments.